[petsc-dev] NaN constant?

[Adrian Croucher]

hi

The reason I'm using NaNs is that John O'Sullivan asked on here some
time ago about how to make a SNES abort when something goes wrong in
the function evaluation (in our case, probably the solution going out
of range of thermodynamic functions). The answer he got was you could
put a NaN or Inf anywhere in the function vector to make the SNES return
SNES_DIVERGED_FNORM_NAN.

Cheers, Adrian
 
-- 
Dr Adrian Croucher
Department of Engineering Science
University of Auckland
New Zealand
tel 64-9-373-7599 ext 84611