diff options
| author | Steven G. Kargl <kargls@comcast.net> | 2005-10-11 23:35:27 +0000 |
|---|---|---|
| committer | Steven G. Kargl <kargl@gcc.gnu.org> | 2005-10-11 23:35:27 +0000 |
| commit | a2694f689cbb6c4780c2aba88c8462c8059cc182 (patch) | |
| tree | 3f30663653c6b162b724455b8cf3773c387ba347 | |
| parent | 55792875db8332861380d621fc5a2877c5cc2617 (diff) | |
| download | gcc-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
| -rw-r--r-- | gcc/testsuite/ChangeLog | 5 | ||||
| -rw-r--r-- | gcc/testsuite/gfortran.dg/csqrt_2.f | 19 | ||||
| -rw-r--r-- | libgfortran/ChangeLog | 7 | ||||
| -rw-r--r-- | libgfortran/intrinsics/c99_functions.c | 4 |
4 files changed, 33 insertions, 2 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 + 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 { |