>>/ Thank you Barry for your suggestions. />/> //The error I get is now KSP_DIVERGED_INDEFINITE_PC. The matrix that I try to invert is actually symmetric and positive definite. I was not aware that this can lead to a indefinite preconditioner. / > Absolutely. Many preconditioners do not retain this feature even in exact > precision and with numerical effects it can even appear unexpected. > By default BoomAMG doesn't retain this.
Is there then a possibility to tell the BlockJacobi preconditioner to be positive definit as there is for the BoomerAMG preconditioner via -pc_hypre_boomeramg_relax_type_all symmetric-SOR/Jacobi? Bests Daniel // -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20121219/40490170/attachment.html>
