------- Comment #6 from jv244 at cam dot ac dot uk 2006-01-13 18:15 ------- another example: REAL*8 :: a,b read(6,*) a b=a**(3.D0/2.D0) write(6,*) b END gets computed by ifort as b=a*sqrt(a) but this is also turned into pow by gfortran at -O3 -ffast-math -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=25620