<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Jun 6, 2016 at 4:16 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
David,<br>
<br>
Actually all the routines for accessing the MUMPS information already exist. I have added this to the manual page for MATSOLVERMUMPS<br>
<br>
Notes: When a MUMPS factorization fails inside a KSP solve, for example with a KSP_DIVERGED_PCSETUP_FAILED, one can find the MUMPS information about the failure by calling<br>
$ KSPGetPC(ksp,&pc);<br>
$ PCFactorGetMatrix(pc,&mat);<br>
$ MatMumpsGetInfo(mat,....);<br>
$ MatMumpsGetInfog(mat,....); etc.<br>
<br>
.seealso: PCFactorSetMatSolverPackage(), MatSolverPackage, MatMumpsSetICntl(), MatMumpsGetIcntl(), MatMumpsSetCntl(), MatMumpsGetCntl(), MatMumpsGetInfo(), MatMumpsGetInfog(), MatMumpsGetRinfo(), MatMumpsGetRinfog(), KSPGetPC(), PCGetFactor(), PCFactorGetMatrix()<br>
<br>
Please let us know if this does not work<br></blockquote><div><br></div><div><br></div><div>Thanks! That looks good. I've changed my code to use this, and I'll let you know if I have any issues with it.</div><div><br></div><div>Thanks,</div><div>David</div><div><br></div><div><br></div><div><br></div><div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="HOEnZb"><div class="h5">> On Jun 5, 2016, at 12:45 PM, David Knezevic <<a href="mailto:david.knezevic@akselos.com">david.knezevic@akselos.com</a>> wrote:<br>
><br>
> On Sun, Jun 5, 2016 at 1:42 PM, Barry Smith <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> wrote:<br>
><br>
> > On Jun 5, 2016, at 12:12 PM, David Knezevic <<a href="mailto:david.knezevic@akselos.com">david.knezevic@akselos.com</a>> wrote:<br>
> ><br>
> > In PETSc-3.6, when I hit an error using the MUMPS solver, I would get an error like this:<br>
> ><br>
> > [0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------<br>
> > [0]PETSC ERROR: Error in external library<br>
> > [0]PETSC ERROR: Error reported by MUMPS in numerical factorization phase: INFO(1)=-9, INFO(2)=351513<br>
> ><br>
> > In PETSc-3.7, the solver now reports "KSP_DIVERGED_PCSETUP_FAILED" instead. The new behavior is definitely better, so thanks for the update! However, I'd still be curious to be able to see what error codes MUMPS returned (which was previously provided in INFO(1) and INFO(2)), since that can be helpful for diagnosing the reason for the solver failure. Is that info still available somehow?<br>
><br>
> David,<br>
><br>
> There are a couple of ways to obtain them both clunky, we should also provide a simple function call interface such as MatMumpsGetInfo(Mat,int *infog,int *info);<br>
><br>
> Anyways right now you can run with -info and it prints the information (along with too much other information) or you can run with -ksp_error_if_not_converged and then it reverts back to the old behavior of stopping immediately and printing the error message and info values.<br>
><br>
><br>
> OK, got it, thanks!<br>
><br>
> David<br>
<br>
</div></div></blockquote></div><br></div></div>