Dear all, 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?
Thanks a lot, BehZad -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120604/2cf0debe/attachment.html>
