diff options
Diffstat (limited to 'math/s_catanhf.c')
-rw-r--r-- | math/s_catanhf.c | 9 |
1 files changed, 8 insertions, 1 deletions
diff --git a/math/s_catanhf.c b/math/s_catanhf.c index 8385af41b0..ca9a30101e 100644 --- a/math/s_catanhf.c +++ b/math/s_catanhf.c @@ -64,7 +64,14 @@ __catanhf (__complex__ float x) float den = 1.0 - __real__ x; den = i2 + den * den; - __real__ res = 0.25 * (__ieee754_logf (num) - __ieee754_logf (den)); + float f = num / den; + if (f < 0.5) + __real__ res = 0.25 * __ieee754_logf (f); + else + { + num = 4.0 * __real__ x; + __real__ res = 0.25 * __log1pf (num / den); + } den = 1 - __real__ x * __real__ x - i2; |