On Wed, Dec 14, 2011 at 13:50, Mark F. Adams <mark.adams at columbia.edu>wrote:

> What I am doing is classic inexact Newton (or it looks like these folks
> are doing Picard).


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.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: 
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20111214/4cef83c9/attachment.html>

Reply via email to