[Bug libfortran/33469] Default formats for real input are not precise enough

[dominiq at lps dot ens dot fr]

------- 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