<div class="gmail_quote">On Mon, May 28, 2012 at 1:48 AM, behzad baghapour <span dir="ltr"><<a href="mailto:behzad.baghapour@gmail.com" target="_blank">behzad.baghapour@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
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.</blockquote></div><br><div>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.</div>