summaryrefslogtreecommitdiff
path: root/sysdeps/libm-ieee754/s_csqrt.c
diff options
context:
space:
mode:
Diffstat (limited to 'sysdeps/libm-ieee754/s_csqrt.c')
-rw-r--r--sysdeps/libm-ieee754/s_csqrt.c10
1 files changed, 10 insertions, 0 deletions
diff --git a/sysdeps/libm-ieee754/s_csqrt.c b/sysdeps/libm-ieee754/s_csqrt.c
index c5c609bd8c..cda23c702d 100644
--- a/sysdeps/libm-ieee754/s_csqrt.c
+++ b/sysdeps/libm-ieee754/s_csqrt.c
@@ -84,6 +84,7 @@ __csqrt (__complex__ double x)
}
else
{
+#if 0 /* FIXME: this is broken. */
__complex__ double q;
double t, r;
@@ -99,6 +100,15 @@ __csqrt (__complex__ double x)
/* Heron iteration in complex arithmetic. */
res = 0.5 * (q + q / x);
+#else
+ double d, imag;
+
+ d = __ieee754_hypot (__real__ x, __imag__ x);
+ imag = __ieee754_sqrt (0.5 * (d - __real__ x));
+
+ __real__ res = __ieee754_sqrt (0.5 * (d + __real__ x));
+ __imag__ res = __copysign (imag, __imag__ x);
+#endif
}
}