From 9407d5e44b5f395b623e63fe75eec84d2ad47361 Mon Sep 17 00:00:00 2001 From: toon Date: Tue, 19 Nov 2002 23:17:36 +0000 Subject: 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 --- libf2c/ChangeLog | 5 +++++ libf2c/libF77/pow_zz.c | 6 ++++++ 2 files changed, 11 insertions(+) (limited to 'libf2c') diff --git a/libf2c/ChangeLog b/libf2c/ChangeLog index e06717f03bd..7e10295040d 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 522b639bc7e..d5cfbf3d931 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); -- cgit v1.2.1