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.
