Fix catan, catanh inaccuracy through use of log (bug 15394).

1 files changed, 8 insertions, 1 deletions

diff --git a/math/s_catanhl.c b/math/s_catanhl.c
index 6844f03dea..4897c0c871 100644
--- a/math/s_catanhl.c
+++ b/math/s_catanhl.c
@@ -64,7 +64,14 @@ __catanhl (__complex__ long double x)
       long double den = 1.0 - __real__ x;
       den = i2 + den * den;
 
-      __real__ res = 0.25 * (__ieee754_logl (num) - __ieee754_logl (den));
+      long double f = num / den;
+      if (f < 0.5)
+	__real__ res = 0.25 * __ieee754_logl (f);
+      else
+	{
+	  num = 4.0 * __real__ x;
+	  __real__ res = 0.25 * __log1pl (num / den);
+	}
 
       den = 1 - __real__ x * __real__ x - i2;