On Aug 15, 2013, at 8:40 PM, Karl Rupp <[email protected]> wrote: > 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
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 Barry > > > Best regards, > Karli
