diff options
| -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)); } |