This is the mail archive of the
gcc-bugs@gcc.gnu.org
mailing list for the GCC project.
fortran/8587: (double complex zero)**(double prec number) = NAN instead of zero
- From: aganse at apl dot washington dot edu
- To: gcc-gnats at gcc dot gnu dot org
- Date: 15 Nov 2002 01:41:22 -0000
- Subject: fortran/8587: (double complex zero)**(double prec number) = NAN instead of zero
- Reply-to: aganse at apl dot washington dot edu
>Number: 8587
>Category: fortran
>Synopsis: (double complex zero)**(double prec number) = NAN instead of zero
>Confidential: no
>Severity: serious
>Priority: medium
>Responsible: unassigned
>State: open
>Class: sw-bug
>Submitter-Id: net
>Arrival-Date: Thu Nov 14 17:46:01 PST 2002
>Closed-Date:
>Last-Modified:
>Originator: Andy Ganse, Applied Physics Laboratory, Seattle
>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
>Organization:
>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
>Description:
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.
>How-To-Repeat:
compile and run the attached dprec.f
>Fix:
?
>Release-Note:
>Audit-Trail:
>Unformatted:
----gnatsweb-attachment----
Content-Type: application/octet-stream; name="dprec.f"
Content-Transfer-Encoding: base64
Content-Disposition: attachment; filename="dprec.f"
CiAgICAgIGRvdWJsZSBjb21wbGV4IGRjMSwgZGMyLCBkYzMsIGRjNCwgZGM1LCBkYzYKICAgICAg
ZG91YmxlIHByZWNpc2lvbiBkMSwgZDIKICAgICAgcmVhbCByMQogICAgICBkYzE9ZGNtcGx4KDAu
ZDAsMC5kMCkKICAgICAgZDE9MgogICAgICBkMj0yLmQwCiAgICAgIHIxPTIKICAgICAgZGMyPWRj
MSoqZDEKICAgICAgZGMzPWRjMSoqZDIKICAgICAgZGM0PWRjMSoqcjEKICAgICAgZGM1PWRjMSoq
MgogICAgICBkYzY9ZGMxKioyLmQwCiAgICAgIHByaW50ICosIGRjMSwgZGMyLCBkYzMsIGRjNCwg
ZGM1LCBkYzYKCmMgICAgIGlzIHByb2R1Y2luZzoKYyAgICAgKDAuLDAuKSAoTkFOLE5BTikgKE5B
TixOQU4pIChOQU4sTkFOKSAoMC4sMC4pIChOQU4sTkFOKQpjICAgICBpbnN0ZWFkIG9mOgpjICAg
ICAoMC4sMC4pICgwLiwwLikgKDAuLDAuKSAoMC4sMC4pICgwLiwwLikgKDAuLDAuKQoKICAgICAg
ZW5kIAo=