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>
