------- Comment #4 from tkoenig at gcc dot gnu dot org 2006-06-01 19:13 ------- (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.
I'll ask con comp.lang.fortran :-) -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=23452
