[petsc-users] sources of floating point randomness in JFNK in serial

[Mark Lohry]

I'm running multiple iterations of newtonls with an MFFD/JFNK nonlinear
solver where I give it the sparsity. PC asm, KSP gmres, with
SNESSetLagJacobian -2 (compute once and then frozen jacobian).

I'm seeing slight (<1%) but nonzero differences in residuals from run to
run. I'm wondering where randomness might enter here -- does the jacobian
coloring use a random seed?
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20230503/ea100732/attachment.html>