[petsc-users] Automatically re-solving after MUMPS error

[Matthew Knepley]

On Tue, Sep 15, 2015 at 4:30 AM, David Knezevic <david.knezevic at akselos.com>
wrote:

> In some cases, I get MUMPS error -9, i.e.:
> [2]PETSC ERROR: Error reported by MUMPS in numerical factorization phase:
> INFO(1)=-9, INFO(2)=98927
>
> This is easily fixed by re-running the executable with -mat_mumps_icntl_14
> on the commandline.
>
> However, I would like to update my code in order to do this automatically,
> i.e. detect the -9 error and re-run with the appropriate option. Is there a
> recommended way to do this? It seems to me that I could do this with a
> PETSc error handler (e.g. PetscPushErrorHandler) in order to call a
> function that sets the appropriate option and solves again, is that right?
> Are there any examples that illustrate this type of thing?
>

I would not use the error handler. I would just check the ierr return code
from the solver. I think you need the
INFO output, for which you can use MatMumpsGetInfo().

  Thanks,

    Matt


> Thanks,
> David
>
>  <http://www.akselos.com>
>
>


-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150915/f7e7d0c9/attachment.html>