<div class="gmail_quote">On Wed, Dec 14, 2011 at 13:50, Mark F. Adams <span dir="ltr"><<a href="mailto:mark.adams@columbia.edu">mark.adams@columbia.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
What I am doing is classic inexact Newton (or it looks like these folks are doing Picard).</blockquote></div><br><div>It appears that they are doing Picard in one-shot mode instead of in defect correction mode. If you write it as a defect correction, then the initial guess will be zero and the standard tolerance is the correct one. I maintain that defect correction is better in every way for solvers, and it naturally gives you a nonlinear residual so that you can use a reliable nonlinear solver tolerance.</div>