summaryrefslogtreecommitdiff
path: root/lib/msun
diff options
context:
space:
mode:
authorBruce Evans <bde@FreeBSD.org>2018-07-19 15:04:10 +0000
committerBruce Evans <bde@FreeBSD.org>2018-07-19 15:04:10 +0000
commitb7092eef4d0fbf133fccca0e5b863e4a9d7881a2 (patch)
tree3a57b0a88831e0b4d00fd963df15268a57d65b48 /lib/msun
parent50c8bd4e537903cb8ebfdda99fa0f030ec3698cd (diff)
downloadsrc-test2-b7092eef4d0fbf133fccca0e5b863e4a9d7881a2.tar.gz
src-test2-b7092eef4d0fbf133fccca0e5b863e4a9d7881a2.zip
Notes
Diffstat (limited to 'lib/msun')
-rw-r--r--lib/msun/src/s_csqrt.c10
-rw-r--r--lib/msun/src/s_csqrtl.c10
2 files changed, 10 insertions, 10 deletions
diff --git a/lib/msun/src/s_csqrt.c b/lib/msun/src/s_csqrt.c
index ca4924fc2bdc..d172efc75651 100644
--- a/lib/msun/src/s_csqrt.c
+++ b/lib/msun/src/s_csqrt.c
@@ -99,15 +99,15 @@ csqrt(double complex z)
/* Algorithm 312, CACM vol 10, Oct 1967. */
if (a >= 0) {
t = sqrt((a + hypot(a, b)) * 0.5);
- rx = t;
- ry = b / (2 * t);
+ rx = scale * t;
+ ry = scale * b / (2 * t);
} else {
t = sqrt((-a + hypot(a, b)) * 0.5);
- rx = fabs(b) / (2 * t);
- ry = copysign(t, b);
+ rx = scale * fabs(b) / (2 * t);
+ ry = copysign(scale * t, b);
}
- return (CMPLX(rx * scale, ry * scale));
+ return (CMPLX(rx, ry));
}
#if LDBL_MANT_DIG == 53
diff --git a/lib/msun/src/s_csqrtl.c b/lib/msun/src/s_csqrtl.c
index fad53b9e9295..40bc59df40ff 100644
--- a/lib/msun/src/s_csqrtl.c
+++ b/lib/msun/src/s_csqrtl.c
@@ -114,13 +114,13 @@ csqrtl(long double complex z)
/* Algorithm 312, CACM vol 10, Oct 1967. */
if (a >= 0) {
t = sqrtl((a + hypotl(a, b)) * 0.5);
- rx = t;
- ry = b / (2 * t);
+ rx = scale * t;
+ ry = scale * b / (2 * t);
} else {
t = sqrtl((-a + hypotl(a, b)) * 0.5);
- rx = fabsl(b) / (2 * t);
- ry = copysignl(t, b);
+ rx = scale * fabsl(b) / (2 * t);
+ ry = copysignl(scale * t, b);
}
- return (CMPLXL(rx * scale, ry * scale));
+ return (CMPLXL(rx, ry));
}