aboutsummaryrefslogtreecommitdiff
path: root/gcc
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 /gcc
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 'gcc')
-rw-r--r--gcc/testsuite/ChangeLog5
-rw-r--r--gcc/testsuite/gfortran.dg/csqrt_2.f19
2 files changed, 24 insertions, 0 deletions
diff --git a/gcc/testsuite/ChangeLog b/gcc/testsuite/ChangeLog
index 39ea01a..d21780a 100644
--- a/gcc/testsuite/ChangeLog
+++ b/gcc/testsuite/ChangeLog
@@ -1,3 +1,8 @@
+2005-10-11 Steven G. Kargl <kargls@comcast.net>
+
+ PR libgfortran/24313
+ gfortran.dg/csqrt.f: New test.
+
2005-10-11 Ian Lance Taylor <ian@airs.com>
PR c++/8057
diff --git a/gcc/testsuite/gfortran.dg/csqrt_2.f b/gcc/testsuite/gfortran.dg/csqrt_2.f
new file mode 100644
index 0000000..66e0842
--- /dev/null
+++ b/gcc/testsuite/gfortran.dg/csqrt_2.f
@@ -0,0 +1,19 @@
+c { dg-do run { xfail *-*-linux-gnu } }
+c Fix PR libgfortran/24313
+ complex x, y
+ complex z
+ z = cmplx(0.707106, -0.707106)
+ x = cmplx(0.0,-1.0)
+ y = sqrt(x)
+ if (abs(y - z) / abs(z) > 1.e-4) call abort
+
+ x = cmplx(tiny(1.),-1.0)
+ y = sqrt(x)
+ if (abs(y - z) / abs(z) > 1.e-4) call abort
+
+ x = cmplx(-tiny(1.),-1.0)
+ y = sqrt(x)
+ if (abs(y - z) / abs(z) > 1.e-4) call abort
+
+ end
+