aboutsummaryrefslogtreecommitdiff
path: root/math/s_ctanhf.c
diff options
context:
space:
mode:
authorJoseph Myers <joseph@codesourcery.com>2012-07-04 09:55:26 +0000
committerJoseph Myers <joseph@codesourcery.com>2012-07-04 09:55:26 +0000
commitca61cf32d934eda9130c4d3c6911892877ad7b0d (patch)
tree46022718d233a6b0ca7f2ba1d2d99d0356430aae /math/s_ctanhf.c
parentca48a46a03da4b64e623ac300929dd137f41adc6 (diff)
downloadglibc-ca61cf32d934eda9130c4d3c6911892877ad7b0d.tar
glibc-ca61cf32d934eda9130c4d3c6911892877ad7b0d.tar.gz
glibc-ca61cf32d934eda9130c4d3c6911892877ad7b0d.tar.bz2
glibc-ca61cf32d934eda9130c4d3c6911892877ad7b0d.zip
Fix ctan, ctanh of subnormals in round-upwards mode (bug 14328).
Diffstat (limited to 'math/s_ctanhf.c')
-rw-r--r--math/s_ctanhf.c18
1 files changed, 15 insertions, 3 deletions
diff --git a/math/s_ctanhf.c b/math/s_ctanhf.c
index e505155774..ca36a83bfb 100644
--- a/math/s_ctanhf.c
+++ b/math/s_ctanhf.c
@@ -83,10 +83,22 @@ __ctanhf (__complex__ float x)
}
else
{
- float sinhrx = __ieee754_sinhf (__real__ x);
- float coshrx = __ieee754_coshf (__real__ x);
+ float sinhrx, coshrx;
+ if (fabsf (__real__ x) > FLT_MIN)
+ {
+ sinhrx = __ieee754_sinhf (__real__ x);
+ coshrx = __ieee754_coshf (__real__ x);
+ }
+ else
+ {
+ sinhrx = __real__ x;
+ coshrx = 1.0f;
+ }
- den = sinhrx * sinhrx + cosix * cosix;
+ if (fabsf (sinhrx) > fabsf (cosix) * FLT_EPSILON)
+ den = sinhrx * sinhrx + cosix * cosix;
+ else
+ den = cosix * cosix;
__real__ res = sinhrx * coshrx / den;
__imag__ res = sinix * cosix / den;
}