aboutsummaryrefslogtreecommitdiff
path: root/math/s_catan.c
diff options
context:
space:
mode:
authorJoseph Myers <joseph@codesourcery.com>2013-04-24 18:49:13 +0000
committerJoseph Myers <joseph@codesourcery.com>2013-04-24 18:49:13 +0000
commit2f38fbfe09e4856c571bf0c80844e5dac9bc77ec (patch)
tree6fc88a0757c95d123f8e8cd9b14e973c9a654d88 /math/s_catan.c
parent45d69176e8b6dc0787dcb8e52e4283b8e7e5e282 (diff)
downloadglibc-2f38fbfe09e4856c571bf0c80844e5dac9bc77ec.zip
glibc-2f38fbfe09e4856c571bf0c80844e5dac9bc77ec.tar.gz
glibc-2f38fbfe09e4856c571bf0c80844e5dac9bc77ec.tar.bz2
Fix catan, catanh inaccuracy through use of log (bug 15394).
Diffstat (limited to 'math/s_catan.c')
-rw-r--r--math/s_catan.c11
1 files changed, 9 insertions, 2 deletions
diff --git a/math/s_catan.c b/math/s_catan.c
index 46c18bf..783941a 100644
--- a/math/s_catan.c
+++ b/math/s_catan.c
@@ -61,7 +61,7 @@ __catan (__complex__ double x)
}
else
{
- double r2, num, den;
+ double r2, num, den, f;
r2 = __real__ x * __real__ x;
@@ -75,7 +75,14 @@ __catan (__complex__ double x)
den = __imag__ x - 1.0;
den = r2 + den * den;
- __imag__ res = 0.25 * __ieee754_log (num / den);
+ f = num / den;
+ if (f < 0.5)
+ __imag__ res = 0.25 * __ieee754_log (f);
+ else
+ {
+ num = 4.0 * __imag__ x;
+ __imag__ res = 0.25 * __log1p (num / den);
+ }
}
return res;