aboutsummaryrefslogtreecommitdiff
path: root/libgfortran
diff options
context:
space:
mode:
authorSteven G. Kargl <kargls@comcast.net>2005-10-11 23:35:27 +0000
committerSteven G. Kargl <kargl@gcc.gnu.org>2005-10-11 23:35:27 +0000
commita2694f689cbb6c4780c2aba88c8462c8059cc182 (patch)
tree3f30663653c6b162b724455b8cf3773c387ba347 /libgfortran
parent55792875db8332861380d621fc5a2877c5cc2617 (diff)
downloadgcc-a2694f689cbb6c4780c2aba88c8462c8059cc182.zip
gcc-a2694f689cbb6c4780c2aba88c8462c8059cc182.tar.gz
gcc-a2694f689cbb6c4780c2aba88c8462c8059cc182.tar.bz2
re PR libfortran/24313 (complex sqrt function does not return principal value)
PR libgfortran/24313 * c99_functions.c (csqrtf, csqrt): Fix choice of branch cut. Note csqrt{f} were imported from glibc, and this bug is still present there. glibc PR is 1146. From-SVN: r105274
Diffstat (limited to 'libgfortran')
-rw-r--r--libgfortran/ChangeLog7
-rw-r--r--libgfortran/intrinsics/c99_functions.c4
2 files changed, 9 insertions, 2 deletions
diff --git a/libgfortran/ChangeLog b/libgfortran/ChangeLog
index d5df1d3..cca3bbe 100644
--- a/libgfortran/ChangeLog
+++ b/libgfortran/ChangeLog
@@ -1,3 +1,10 @@
+2005-10-11 Steven G. Kargl <kargls@comcast.net>
+
+ PR libgfortran/24313
+ * c99_functions.c (csqrtf, csqrt): Fix choice of branch cut. Note
+ csqrt{f} were imported from glibc, and this bug is still present
+ there. glibc PR is 1146.
+
2005-10-07 Janne Blomqvist <jblomqvi@cc.hut.fi>
PR fortran/16339
diff --git a/libgfortran/intrinsics/c99_functions.c b/libgfortran/intrinsics/c99_functions.c
index 2417e69..906a85a 100644
--- a/libgfortran/intrinsics/c99_functions.c
+++ b/libgfortran/intrinsics/c99_functions.c
@@ -669,7 +669,7 @@ csqrtf (float complex z)
r = sqrtf (0.5 * fabsf (im));
- COMPLEX_ASSIGN (v, copysignf (r, im), r);
+ COMPLEX_ASSIGN (v, r, copysignf (r, im));
}
else
{
@@ -722,7 +722,7 @@ csqrt (double complex z)
r = sqrt (0.5 * fabs (im));
- COMPLEX_ASSIGN (v, copysign (r, im), r);
+ COMPLEX_ASSIGN (v, r, copysign (r, im));
}
else
{