<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=utf-8">

</head>

<body>

<div dir="ltr">

<div class="gmail_quote">

<div dir="ltr">Junchao :<br>

</div>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<div>

<div dir="ltr">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?

<div>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?</div>

</div>

</div>

</blockquote>

<div><br>

</div>

<div>If you suspect something is wrong, produce an example that validate your guess. We then work from there and improve PETSc. </div>

<div>Hong </div>

</div>

</div>

</body>

</html>