diff options
| author | Joseph Myers <joseph@codesourcery.com> | 2012-07-04 09:55:26 +0000 |
|---|---|---|
| committer | Joseph Myers <joseph@codesourcery.com> | 2012-07-04 09:55:26 +0000 |
| commit | ca61cf32d934eda9130c4d3c6911892877ad7b0d (patch) | |
| tree | 46022718d233a6b0ca7f2ba1d2d99d0356430aae /math/s_ctanhl.c | |
| parent | ca48a46a03da4b64e623ac300929dd137f41adc6 (diff) | |
| download | glibc-ca61cf32d934eda9130c4d3c6911892877ad7b0d.tar.gz glibc-ca61cf32d934eda9130c4d3c6911892877ad7b0d.tar.xz glibc-ca61cf32d934eda9130c4d3c6911892877ad7b0d.zip | |
Fix ctan, ctanh of subnormals in round-upwards mode (bug 14328).
Diffstat (limited to 'math/s_ctanhl.c')
| -rw-r--r-- | math/s_ctanhl.c | 18 |
1 files changed, 15 insertions, 3 deletions
diff --git a/math/s_ctanhl.c b/math/s_ctanhl.c index e5d677903f..dbf1612707 100644 --- a/math/s_ctanhl.c +++ b/math/s_ctanhl.c @@ -83,10 +83,22 @@ __ctanhl (__complex__ long double x) } else { - long double sinhrx = __ieee754_sinhl (__real__ x); - long double coshrx = __ieee754_coshl (__real__ x); + long double sinhrx, coshrx; + if (fabsl (__real__ x) > LDBL_MIN) + { + sinhrx = __ieee754_sinhl (__real__ x); + coshrx = __ieee754_coshl (__real__ x); + } + else + { + sinhrx = __real__ x; + coshrx = 1.0L; + } - den = sinhrx * sinhrx + cosix * cosix; + if (fabsl (sinhrx) > fabsl (cosix) * LDBL_EPSILON) + den = sinhrx * sinhrx + cosix * cosix; + else + den = cosix * cosix; __real__ res = sinhrx * coshrx / den; __imag__ res = sinix * cosix / den; } |