<div class="gmail_quote">On Mon, Apr 23, 2012 at 08:11, Dominik Szczerba <span dir="ltr">&lt;<a href="mailto:dominik@itis.ethz.ch">dominik@itis.ethz.ch</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="im">What is &#39;func&#39; on the SNESSetPicard manual page. It only says<br>
&quot;function evaluation routine&quot;. What function? Do you mean Ax-b?<br></div></blockquote><div><br></div><div>Just b(x)</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="im">
<br>
&gt; All SNESSetPicard does is evaluate both A and b in &quot;residual evaluation&quot;,<br>
&gt; compute F(x) = A(x)x - b(x), and internally cache the matrix A(x). Using it<br>
&gt; is completely optional. If you are willing to write residual evaluation for<br>
&gt; the whole F(x), then you will benefit from having less expensive residual<br>
&gt; evaluations, making line searches and matrix-free Newton (-snes_mf_operator)<br>
&gt; affordable.<br>
<br>
</div>I understood from the documentation sec. 5.1 and 5.6 that<br>
-snes_mf_operator is not well suited for unstructured problems, rather<br>
for structured problems with known stencil. Is it not so?</blockquote></div><br><div>You must be thinking of coloring.</div>