about summary refs log tree commit diff
path: root/math/s_catanh.c
diff options
context:
space:
mode:
authorJoseph Myers <joseph@codesourcery.com>2013-04-26 19:26:22 +0000
committerJoseph Myers <joseph@codesourcery.com>2013-04-26 19:26:22 +0000
commit9457fd952cafc701375fc9fae575a4aca1374de3 (patch)
tree0116832c92205eb0c42197ce9c269a7f52d09c5e /math/s_catanh.c
parentf0302940e7c2acb587971e3c99dfbd00aa4e2134 (diff)
downloadglibc-9457fd952cafc701375fc9fae575a4aca1374de3.tar.gz
glibc-9457fd952cafc701375fc9fae575a4aca1374de3.tar.xz
glibc-9457fd952cafc701375fc9fae575a4aca1374de3.zip
Fix catan, catanh missing underflows (bug 15406).
Diffstat (limited to 'math/s_catanh.c')
-rw-r--r--math/s_catanh.c13
1 files changed, 12 insertions, 1 deletions
diff --git a/math/s_catanh.c b/math/s_catanh.c
index 0ee8c64b6f..5d18cd639f 100644
--- a/math/s_catanh.c
+++ b/math/s_catanh.c
@@ -20,7 +20,7 @@
 #include <complex.h>
 #include <math.h>
 #include <math_private.h>
-
+#include <float.h>
 
 __complex__ double
 __catanh (__complex__ double x)
@@ -76,6 +76,17 @@ __catanh (__complex__ double x)
       den = 1 - __real__ x * __real__ x - i2;
 
       __imag__ res = 0.5 * __ieee754_atan2 (2.0 * __imag__ x, den);
+
+      if (fabs (__real__ res) < DBL_MIN)
+	{
+	  volatile double force_underflow = __real__ res * __real__ res;
+	  (void) force_underflow;
+	}
+      if (fabs (__imag__ res) < DBL_MIN)
+	{
+	  volatile double force_underflow = __imag__ res * __imag__ res;
+	  (void) force_underflow;
+	}
     }
 
   return res;