[petsc-dev] -mat_error_if_fpe?

[Dmitry Karpeyev]

Currently -pc_type lu will not detect a zero pivot (it will happily produce
an Inf or a NaN) and the flag it checks whether to see whether
to do it is erroriffpe.   KSP_DIVERGED_NANORINF is too generic.  What if
the user wants to verify there is a zero pivot?

Also, irrespective of the use case, if we have a XXXSetYYY(), isn't our
policy to have a command-line way of doing it?

On Tue, Jun 2, 2015 at 2:21 PM Barry Smith <bsmith at mcs.anl.gov> wrote:

>
> > On Jun 2, 2015, at 3:02 PM, Dmitry Karpeyev <karpeev at mcs.anl.gov> wrote:
> >
> > Should we have a way to turn on/off MatErrorIfFPE() on the command line?
>
>   Sure, one does it with -ksp_error_if_not_converged or
> -snes_error_if_not_converged  It only adds to user confusion, but no user
> benefit, to have a -mat_error_if_fpe thing.
>
>   What user case do think there is for a -mat_error_if_fpe ?
>
>   Barry
>
>
>
> > It is now turned off by default to allow Inf/NaNs to propagate and
> currently the only way to toggle it is programmatic.
>
>
> >
> > Dmitry.
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20150602/530f7453/attachment.html>