From ca61cf32d934eda9130c4d3c6911892877ad7b0d Mon Sep 17 00:00:00 2001 From: Joseph Myers Date: Wed, 4 Jul 2012 09:55:26 +0000 Subject: Fix ctan, ctanh of subnormals in round-upwards mode (bug 14328). --- math/s_ctanh.c | 18 +++++++++++++++--- 1 file changed, 15 insertions(+), 3 deletions(-) (limited to 'math/s_ctanh.c') diff --git a/math/s_ctanh.c b/math/s_ctanh.c index 201871e7ec..d288b7d168 100644 --- a/math/s_ctanh.c +++ b/math/s_ctanh.c @@ -83,10 +83,22 @@ __ctanh (__complex__ double x) } else { - double sinhrx = __ieee754_sinh (__real__ x); - double coshrx = __ieee754_cosh (__real__ x); + double sinhrx, coshrx; + if (fabs (__real__ x) > DBL_MIN) + { + sinhrx = __ieee754_sinh (__real__ x); + coshrx = __ieee754_cosh (__real__ x); + } + else + { + sinhrx = __real__ x; + coshrx = 1.0; + } - den = sinhrx * sinhrx + cosix * cosix; + if (fabs (sinhrx) > fabs (cosix) * DBL_EPSILON) + den = sinhrx * sinhrx + cosix * cosix; + else + den = cosix * cosix; __real__ res = sinhrx * coshrx / den; __imag__ res = sinix * cosix / den; } -- cgit v1.2.3