<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body>
<div dir="ltr">
<div dir="ltr">OK, I found -ksp_error_if_not_converged will trigger PETSc to fail in this case.</div>
<div dir="ltr"><br clear="all">
<div>
<div dir="ltr" class="gmail_signature">
<div dir="ltr">--Junchao Zhang</div>
</div>
</div>
<br>
</div>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr">On Tue, Oct 9, 2018 at 3:38 PM Junchao Zhang <<a href="mailto:jczhang@mcs.anl.gov">jczhang@mcs.anl.gov</a>> wrote:<br>
</div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div>I met a case where MUMPS returned an out-of-memory code but PETSc continued to run. When PETSc calls MUMPS, it checks if (A->erroriffailure). I added -mat_error_if_failure, but it did not work since it was overwritten by MatSetErrorIfFailure(pc->pmat,pc->erroriffailure)</div>
<div>Does it suggest we should add a new option -pc_factor_error_if_failure and check it in PCSetFromOptions_Factor()?</div>
<div><br>
</div>
<div>--Junchao Zhang</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr">On Fri, Oct 5, 2018 at 8:12 PM Zhang, Hong <<a href="mailto:hzhang@mcs.anl.gov" target="_blank">hzhang@mcs.anl.gov</a>> wrote:<br>
</div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div dir="ltr">
<div class="gmail_quote">
<div dir="ltr">Mike:<br>
</div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div>Hello PETSc team:<br>
</div>
<div><br>
</div>
<div>I am trying to solve a PDE problem with high-order finite elements. The matrix is getting denser and my experience is that MUMPS just outperforms iterative solvers.</div>
<div><br>
</div>
<div>For certain problems, MUMPS just fail in the middle for no clear reason. I just wander if there is any suggestion to improve the robustness of MUMPS? Or in general, any suggestion for interative solver with very high-order finite elements?</div>
</div>
</blockquote>
<div> </div>
<div>What error message do you get when MUMPS fails? Out of memory, zero pivoting, or something?</div>
<div> Hong</div>
</div>
</div>
</div>
</blockquote>
</div>
</blockquote>
</div>
</body>
</html>