[petsc-dev] SNESMonitorVI fix: maint or master?

Dmitry Karpeyev karpeev at mcs.anl.gov
Wed Oct 7 17:48:46 CDT 2015


Now that we allow NaNs bubble up through the solver,
this can trip mysterious-looking errors in SNESVI:
PETSC_ERR_PLIB, "Can never get here"
is thrown from SNESMonitorVI(), because a NaN in
the residual can defeat all of the seemingly-exhaustive
if-then-else branches counting the number of active constraints.

I think the fix should be to count the number of NaNs separately
and report them alongside the legitimate active bounds to give
the user as much useful information as possible.  Since this entails
a substantial difference to the output format of -snes_vi_monitor,
should the fix go to maint or master?

Dmitry.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20151007/a9f24bad/attachment.html>


More information about the petsc-dev mailing list