[petsc-dev] How to know MatSolve() was successful?
Zhang, Junchao
jczhang at mcs.anl.gov
Fri Oct 12 11:33:22 CDT 2018
I learned users should call KSPGetConvergedReason to check if a KSPSolve was successful. But if users directly call MatSolve() or MatMatSolve etc, how can they know it was successful?
I see in MatSolve_MUMPS, it first checks factor errors. If there is, it sets the solution to infinity; But it does not do that in MatMatSolve. Also, in MatSolve_MUMPS, if the factorization was fine but the solve failed, the code just aborts. Is it wrong in a KSP context?
--Junchao Zhang
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20181012/7891daad/attachment.html>
More information about the petsc-dev
mailing list