Error in subnormal computation

Bernhard Fischer rep.nop@aon.at
Sat Oct 14 19:29:00 GMT 2006


On Sun, Oct 01, 2006 at 04:21:51PM +0200, Dominique Dhumieres wrote:
>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
>

Can't reproduce this.
gcc version 4.2.0 20061012 (experimental) gives

 -4.2038954E-45

What version are you using?



More information about the Fortran mailing list