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
