How to determine that 2.0 == (float) 2 in a front-end.

[Toon Moene]

Robert Dewar wrote:

> I wrote:

> > x**6.0 can be rendered as x**6 (which is y=x*x, z=y*y, z=z*y).
> 
> Note that this is liikely less accurate than calling a well written
> pow routine.
> 
> Note that for example x*x*x*x is generally more accurate than x**2**2
> 
> I don't necesarily argue for the more accurate form (generally Fortran
> operates in the "who cares that much about last bit accuracy anyway")
> but it is good to be aware of these things.

Yes, that's why it will only be allowed under unsafe_math_optimizations
(once I get around figgerin' out how to do it, of course).

-- 
Toon Moene - mailto:toon@moene.indiv.nluug.nl - phoneto: +31 346 214290
Saturnushof 14, 3738 XG  Maartensdijk, The Netherlands
Maintainer, GNU Fortran 77: http://gcc.gnu.org/onlinedocs/g77_news.html
Join GNU Fortran 95: http://g95.sourceforge.net/ (under construction)