Commit a2694f68 by Steven G. Kargl Committed by Steven G. Kargl

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
parent 55792875
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
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
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
......
......@@ -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
{
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment