<div class="gmail_extra">On Mon, Dec 3, 2012 at 9:56 AM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div id=":29r">   Ok, clearly we still don't understand Picard and what is called Picard in PETSc is probably all wrong.<br>
<br>
   Start with F(x) = 0  Introduce any "reasonable" matrix function G(x) and write the equivalent nonlinear system G(x) x = -F(x) + G(x)x<br>
non-defect correction Picard iteration is solve G(x^n) x^{n+1} = -F(x^n) + G(x^n) x^n write this in defect correction style as<br>
solve G(x^n)(x^{n+1} - x^n) = -F(x^n).  When G(x^n) is the Jacobian of F() this is Newton's method. Correct?<br>
<br>
   So a user can do (defect correction) "Picard" for their problem simply by providing any G() they want and the original F() using the usual PETSc interface, correct?<br>
<br>
   Given this introduction it is wrong to say there are many methods _between_ Newton and Picard. It is really that (defect correction) Picard is a spectrum of methods (any reasonable G()) and Newton is a special case. Correct?<br>

<br>
    The way "Picard" comes into PETSc is when the user does not have F(x) instead they have A(x) x = f(x). They want to supply A(x) and f(x) and A(x) may not be particularly near the Jacobian of A(x)x - f(x) (since f() can have all kinds of stuff in it).<br>

<br>
Thus Picard in PETSc means the user provides A(x) and f(x) and Newton means the user provides F(x) and any "Jacobian" they want.</div></blockquote></div><br></div><div class="gmail_extra">Fine, but I want to encourage people to "do Picard" by writing F(x) and A(x) because then they can incrementally add terms to A without breaking the structure. If they write to the "Picard" (A(x), f(x)) interface, they are constraining themselves and probably won't try using other methods. Regardless, we still want to support the equivalent of -snes_mf_operator (Newton) despite them using that silly "Picard" interface.</div>