[petsc-users] Initial estimation on SNES and KSP

Jed Brown jed at jedbrown.org
Thu Sep 4 12:28:30 CDT 2014


Miguel Angel Salazar de Troya <salazardetroya at gmail.com> writes:

> Dear all
>
> SNES uses internally a KSP to solve the linear system of equations right?
> Now the case that we had a linear system of equations that we are solving
> with SNES, how could we set the initial estimation for the KSP? If we just
> included the option -ksp_initial_guess_nonzero, the KSP will grab the
> vector X we passed to the SNES? 

You definitely don't want this for Newton-type methods.  If you have a
nonzero guess for the solution of the linear system, you should have
evaluated the Jacobian at that point.  A zero initial guess is optimal
for Newton-type methods.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 818 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140904/f756154c/attachment.pgp>


More information about the petsc-users mailing list