8587 – (double complex zero)**(double prec number) = NAN instead of zero (original) (raw)

There are a couple situations (produced in the small attached test source file) when a double complex variable set to zero, and raised to the power of another variable, is not equaling zero - it's coming out as NAN instead. This all works fine in Sun F77 if that's useful for comparison... The case where a double complex variable equaling zero is raised to a real (as opposed to double precision) constant works here but is producing slightly off values in our program when the double complex variable does not equal zero.

Release: GNU Fortran 0.5.26 20000731 (Red Hat Linux 7.2 2.96-112.7.2) and also GNU Fortran 0.5.25 19991024

Environment: GNU Fortran 0.5.26 20000731 on RedHat i386/7.3 on athlon GNU Fortran 0.5.25 19991024 on SunOS 5.8 on Sun-Blade-1000

How-To-Repeat: compile and run the attached dprec.f

Comment 1 aganse 2002-11-14 17:46:01 UTC

Fix: ?

Comment 2 Toon Moene 2002-11-18 14:14:54 UTC

Responsible-Changed-From-To: unassigned->toon Responsible-Changed-Why: Fortran Maintainer.

Comment 3 Toon Moene 2002-11-18 14:14:54 UTC

State-Changed-From-To: open->analyzed State-Changed-Why: Also fails on Debian 3.0+ (testing) using 3.2 prerelease on powerpc-unknown-linux-gnu

Comment 4 Toon Moene 2002-11-19 15:25:10 UTC

State-Changed-From-To: analyzed->closed State-Changed-Why: Fixed for GCC/g77-3.3