------- Comment #5 from pinskia at gcc dot gnu dot org 2006-11-01 04:30 ------- (In reply to comment #3) > This is now fixed on mainline provided the user specifies -ffast-math. > There are some complications where imagpart(z*~z) can be non-zero, if > imagpart(z) is non-finite, such as an Inf or a NaN. It's unclear from > the fortran-95 standard whether gfortran is allowed to optimize this > even without -ffast-math.
Well Fortran 95 does not have non-finite except with a TR and the IEEE fp types, I forgot the number and the name -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=23452