http://gcc.gnu.org/bugzilla/show_bug.cgi?id=61002
Dominique d'Humieres <dominiq at lps dot ens.fr> changed: What |Removed |Added ---------------------------------------------------------------------------- Status|UNCONFIRMED |WAITING Last reconfirmed| |2014-05-01 Ever confirmed|0 |1 --- Comment #3 from Dominique d'Humieres <dominiq at lps dot ens.fr> --- Another way to get Infinity or NaN can be achieved with the following code (which traps at run time when compiled with -ffpe-trap=invalid) real :: x, y, z x = 0.0 y = 1.0/x z = log(-y) print *, y, z end If there is no objection I'll close this PR as INVALID in a couple of days.