<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Sep 4, 2014 at 12:28 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="">Miguel Angel Salazar de Troya <<a href="mailto:salazardetroya@gmail.com">salazardetroya@gmail.com</a>> writes:<br>
<br>
> Dear all<br>
><br>
> SNES uses internally a KSP to solve the linear system of equations right?<br>
> Now the case that we had a linear system of equations that we are solving<br>
> with SNES, how could we set the initial estimation for the KSP? If we just<br>
> included the option -ksp_initial_guess_nonzero, the KSP will grab the<br>
> vector X we passed to the SNES?<br>
<br>
</div>You definitely don't want this for Newton-type methods.  If you have a<br>
nonzero guess for the solution of the linear system, you should have<br>
evaluated the Jacobian at that point.  A zero initial guess is optimal<br>
for Newton-type methods.<br>
</blockquote></div><br>Notice that Newton is solving for the correction, not the solution itself.</div><div class="gmail_extra"><br></div><div class="gmail_extra">  Matt<br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>