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