about summary refs log tree commit diff
path: root/math/s_ctan.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_ctan.c
parentca48a46a03da4b64e623ac300929dd137f41adc6 (diff)
downloadglibc-ca61cf32d934eda9130c4d3c6911892877ad7b0d.tar.gz
glibc-ca61cf32d934eda9130c4d3c6911892877ad7b0d.tar.xz
glibc-ca61cf32d934eda9130c4d3c6911892877ad7b0d.zip
Fix ctan, ctanh of subnormals in round-upwards mode (bug 14328).
Diffstat (limited to 'math/s_ctan.c')
-rw-r--r--math/s_ctan.c18
1 files changed, 15 insertions, 3 deletions
diff --git a/math/s_ctan.c b/math/s_ctan.c
index 78117b3103..89c0fef91e 100644
--- a/math/s_ctan.c
+++ b/math/s_ctan.c
@@ -83,10 +83,22 @@ __ctan (__complex__ double x)
 	}
       else
 	{
-	  double sinhix = __ieee754_sinh (__imag__ x);
-	  double coshix = __ieee754_cosh (__imag__ x);
+	  double sinhix, coshix;
+	  if (fabs (__imag__ x) > DBL_MIN)
+	    {
+	      sinhix = __ieee754_sinh (__imag__ x);
+	      coshix = __ieee754_cosh (__imag__ x);
+	    }
+	  else
+	    {
+	      sinhix = __imag__ x;
+	      coshix = 1.0;
+	    }
 
-	  den = cosrx * cosrx + sinhix * sinhix;
+	  if (fabs (sinhix) > fabs (cosrx) * DBL_EPSILON)
+	    den = cosrx * cosrx + sinhix * sinhix;
+	  else
+	    den = cosrx * cosrx;
 	  __real__ res = sinrx * cosrx / den;
 	  __imag__ res = sinhix * coshix / den;
 	}