Hey,
>>
KSPSetType(ksp,?);
Jacobi and Gauss-Seidel are available as preconditioners. To use them as
standalone solvers, make the linear solver run the preconditioner only via
KSPPREONLY from here:
http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/KSP/KSPType.html
and then set the preconditioner (e.g. PCJACOBI) via PCSetType():
http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/PC/PCType.html
This won't work for jacobi (it will only do one iteration) and will only
work for sor if you use PCSetSORIterations() instead of -ksp_max_its
Hmm, should have tried it first... :-/ Thanks for the correction.
Best regards,
Karli
