Hi, > > 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 Best regards, Karli
