Error in subnormal computation

Dominique Dhumieres dominiq@lps.ens.fr
Sun Oct 1 14:21:00 GMT 2006


The following code

real x
x = -3.0*(nearest(tiny(x),1.0)-tiny(x))
print *, x
end

gives

   0.000000    

instead of something like (xlf)

 -0.4203895393E-44

TIA

Dominique



More information about the Fortran mailing list