Miguel Angel Salazar de Troya <[email protected]> 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.
pgpz7SreS0X9D.pgp
Description: PGP signature
