* 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.
git-svn-id: svn+ssh://gcc.gnu.org/svn/gcc/trunk@105274
138bc75d-0d04-0410-961f-
82ee72b054a4
+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
--- /dev/null
+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
r = sqrtf (0.5 * fabsf (im));
- COMPLEX_ASSIGN (v, copysignf (r, im), r);
+ COMPLEX_ASSIGN (v, r, copysignf (r, im));
}
else
{
r = sqrt (0.5 * fabs (im));
- COMPLEX_ASSIGN (v, copysign (r, im), r);
+ COMPLEX_ASSIGN (v, r, copysign (r, im));
}
else
{