[Bug fortran/40070] Some math expressions containing exponents fail on a Windows 64 build

["kargl at gcc dot gnu dot org" <gcc-bugzilla at gcc dot gnu dot org>]

------- Comment #8 from kargl at gcc dot gnu dot org  2009-05-08 15:42 -------
(In reply to comment #0)

>       Subroutine myfunc(a,b,c,d,e,res1,res2) 
>       Real a,b,c,d,e
>       Parameter (g = 273.15)
>       res1 = a*(b/a)**((c-d-g)/(e-d))
>       res2 = (b/a)**((c-d-g)/(e-d))
>       res2 = res2 * a
>       Return
>       End

Can you compile the above subroutine with -fdump-tree-original
and post the *.original file?

>     a=1513.25;
>     b=238.149994;
>     c=-16.500040;
>     d=3.2;
>     e=0;
>     f=9.806650;
>     g=20.964399;
>     h=7614.320313;
>     i=-0.1065;

Have you read Goldberg's paper?

What does 

   program a
   print *, precision(1.e0), digits(1.e0)
   end program a

produce?


-- 


http://gcc.gnu.org/bugzilla/show_bug.cgi?id=40070