[petsc-users] Gauss-Seidel

[Jed Brown]

On Mon, Apr 23, 2012 at 20:27, tribur <burckhardt at itis.ethz.ch> wrote:

> I would like to do 10 iterations of a point-by-point Gauss-Seidel
> procedure for the purpose of testing.
>
> Is it correct if I type "mpiexec -n 1 ... -ksp_type preonly -pc_type sor
> -pc_sor_its 10"?
>
> How can I define the initial guess nonzero?
>

With a nonzero initial guess, you are doing Richardson (not preonly). You
can set the nonzero initial guess as usual (KSPSetInitialGuessNonzero()).
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120423/2a115d14/attachment.htm>