On Mon, May 28, 2012 at 1:48 AM, behzad baghapour < behzad.baghapour at gmail.com> wrote:
> Besides, Is the SNES context call residual evaluation then Jacobian or > vise versa? I just concern that the limiting procedure affects both > residual and Jacobian in the same iteration. The Jacobian is, by definition, the derivative of the residual. If you are doing Newton's method, then you are solving with the Jacobian, therefore the Jacobian also "sees" the limiter. -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120528/99b89922/attachment.html>
