"Toon Moene" <toon@xxxxxxxxxxxxxxxxxxxx> wrote in message news:422F54B8.4090904@xxxxxxxxxxxxxxxxxxxxxxx > Lionel B wrote: > > > Is there a Fortran 77 intrinsic equivalent to the libc isnan() function for detecting IEEE fp NaNs? Otherwise any other > > reliable means of detection for my platform (gcc 3.3.3 cygwin on Win 2K)? Testing X.EQ.X doesn't appear to work > > correctly... > > Not surprising - you should test X.NE.X :-) Of course - that was a typo ... but it still doesn't work for me ... > Oh, and don't use -ffast-math, because then you're telling the compiler > you're not interested in exact IEEE 754 conformant floating point math. Yes, that was it. Many thanks, -- Lionel B