[petsc-dev] Domain error in SNESComputeJacobian()?

[Dmitry Karpeyev]

We don't check for domainerror in SNESComputeJacobian(), but it may very
well occur there. It would cost a reduction, of course, but the alternative
is to set a Jacobian entry to Inf and wait for the nearest MatMult(),
VecNorm() combination to propagate it.

MOOSE right now works around that by setting an Inf into the residual  -
 a side effect.

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