Dear all,<br><br>I want to compute Jacobian matrix via Matrix-Free manner in SNES procedure. I should consider the time-marching term as M*dQ/dt with is not the part of finite-differencing and should be added into diagonal blocks of the matrix. So, how should be right to calculate the Jacobian using the prepared Petsc routines?<br>
<br>Thanks a lot,<br>BehZad<br>