------- Comment #11 from dominiq at lps dot ens dot fr 2007-10-03 19:34 -------
BTW I have forgotten to explain why I have to use an auxiliary variable 'z': if
I usenearest(0.0_8,1.0_8); I get
default_format_1_db.f90:70.29:
if (abs(x-y)>nearest(0.0_8,1.0_8)) print *, x, y, x-y
1
Error: Result of NEAREST underflows its kind at (1)
default_format_1_db.f90:84.29:
if (abs(x-y)>nearest(0.0_8,1.0_8)) print *, x, y, x-y
1
Error: Result of NEAREST underflows its kind at (1)
default_format_1_db.f90:93.25:
use test_default_format
1
Fatal Error: Can't open module file 'test_default_format.mod' for reading at
(1): No such file or directory
Producing a hard error for an underflow is indeed an incredibely clever idea!-(
Reminds me of a compiler which crunched NaNs and Infs with delight, but aborted
on underflow, very handy for waves simulations!
--
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=33469
