diff options
| author | Bruce Evans <bde@FreeBSD.org> | 2018-07-19 15:04:10 +0000 |
|---|---|---|
| committer | Bruce Evans <bde@FreeBSD.org> | 2018-07-19 15:04:10 +0000 |
| commit | b7092eef4d0fbf133fccca0e5b863e4a9d7881a2 (patch) | |
| tree | 3a57b0a88831e0b4d00fd963df15268a57d65b48 /lib | |
| parent | 50c8bd4e537903cb8ebfdda99fa0f030ec3698cd (diff) | |
Notes
Diffstat (limited to 'lib')
| -rw-r--r-- | lib/msun/src/s_csqrt.c | 10 | ||||
| -rw-r--r-- | lib/msun/src/s_csqrtl.c | 10 |
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)); } |
