<div dir="ltr">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.<div><br></div><div>MOOSE right now works around that by setting an Inf into the residual  -</div><div> a side effect.  <br><div><br></div><div>Dmitry.</div></div></div>