From 9407d5e44b5f395b623e63fe75eec84d2ad47361 Mon Sep 17 00:00:00 2001 From: toon Date: Tue, 19 Nov 2002 23:17:36 +0000 Subject: [PATCH] 2002-11-19 Toon Moene PR fortran/8587 * libF77/pow_zz.c: Handle (0.0, 0.0) ** power. 2002-11-19 Toon Moene PR fortran/8587 * news.texi: Show PR fortran/8587 fixed. git-svn-id: svn+ssh://gcc.gnu.org/svn/gcc/trunk@59280 138bc75d-0d04-0410-961f-82ee72b054a4 --- gcc/f/ChangeLog | 5 +++++ gcc/f/news.texi | 4 +++- libf2c/ChangeLog | 5 +++++ libf2c/libF77/pow_zz.c | 6 ++++++ 4 files changed, 19 insertions(+), 1 deletion(-) diff --git a/gcc/f/ChangeLog b/gcc/f/ChangeLog index 4960f2f..2142b4c 100644 --- a/gcc/f/ChangeLog +++ b/gcc/f/ChangeLog @@ -1,3 +1,8 @@ +2002-11-19 Toon Moene + + PR fortran/8587 + * news.texi: Show PR fortran/8587 fixed. + 2002-11-19 Jason Thorpe * g77spec.c (lang_specific_spec_functions): New. diff --git a/gcc/f/news.texi b/gcc/f/news.texi index 24cec0a..2d0dc0f 100644 --- a/gcc/f/news.texi +++ b/gcc/f/news.texi @@ -10,7 +10,7 @@ @c in the standalone derivations of this file (e.g. NEWS). @set copyrights-news 1995,1996,1997,1998,1999,2000,2001,2002 -@set last-update-news 2002-11-02 +@set last-update-news 2002-11-19 @include root.texi @@ -181,6 +181,8 @@ gave wrong results) (libf2c) DATE_AND_TIME milliseconds field inactive on Windows @item 7388 Incorrect output with 0-based array of characters +@item 8587 +Double complex zero ** double precision number -> NaN instead of zero @end table @item Richard Henderson (@email{rth@@redhat.com}) analysed and improved the handling diff --git a/libf2c/ChangeLog b/libf2c/ChangeLog index e06717f..7e10295 100644 --- a/libf2c/ChangeLog +++ b/libf2c/ChangeLog @@ -1,3 +1,8 @@ +2002-11-19 Toon Moene + + PR fortran/8587 + * libF77/pow_zz.c: Handle (0.0, 0.0) ** power. + 2002-09-23 Zack Weinberg * libF77/Version.c: Rename junk to __LIBF77_VERSION__. Add diff --git a/libf2c/libF77/pow_zz.c b/libf2c/libF77/pow_zz.c index 522b639..d5cfbf3 100644 --- a/libf2c/libF77/pow_zz.c +++ b/libf2c/libF77/pow_zz.c @@ -8,6 +8,12 @@ pow_zz (doublecomplex * r, doublecomplex * a, doublecomplex * b) { double logr, logi, x, y; + if (a->r == 0.0 && a->i == 0.0) + { + /* Algorithm below doesn't cope. */ + r->r = r->i = 0.0; + return; + } logr = log (f__cabs (a->r, a->i)); logi = atan2 (a->i, a->r); -- 2.7.4