[petsc-users] How to use jacobi or Gauss-seilder iteration to solve problem

[Karl Rupp]

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