about summary refs log tree commit diff
path: root/math/s_ctan.c
diff options
context:
space:
mode:
authorLiubov Dmitrieva <liubov.dmitrieva@gmail.com>2012-09-25 20:41:17 +0200
committerAndreas Jaeger <aj@suse.de>2012-09-25 20:41:17 +0200
commit6d3bf1993190edc502d01e8ca42c4482b20a5b6c (patch)
treeceef4ffabe768f560a3a4d3f6b556b7b70aac37d /math/s_ctan.c
parent203e56032f77d0b36701deaee360e9cdda35d54e (diff)
downloadglibc-6d3bf1993190edc502d01e8ca42c4482b20a5b6c.tar.gz
glibc-6d3bf1993190edc502d01e8ca42c4482b20a5b6c.tar.xz
glibc-6d3bf1993190edc502d01e8ca42c4482b20a5b6c.zip
Fix wrong ussage of sincos for subnormal arguments
Diffstat (limited to 'math/s_ctan.c')
-rw-r--r--math/s_ctan.c11
1 files changed, 10 insertions, 1 deletions
diff --git a/math/s_ctan.c b/math/s_ctan.c
index 89c0fef91e..fdba847447 100644
--- a/math/s_ctan.c
+++ b/math/s_ctan.c
@@ -53,11 +53,20 @@ __ctan (__complex__ double x)
       double sinrx, cosrx;
       double den;
       const int t = (int) ((DBL_MAX_EXP - 1) * M_LN2 / 2);
+      int rcls = fpclassify (__real__ x);
 
       /* tan(x+iy) = (sin(2x) + i*sinh(2y))/(cos(2x) + cosh(2y))
 	 = (sin(x)*cos(x) + i*sinh(y)*cosh(y)/(cos(x)^2 + sinh(y)^2). */
 
-      __sincos (__real__ x, &sinrx, &cosrx);
+      if (__builtin_expect (rcls != FP_SUBNORMAL, 1))
+	{
+	  __sincos (__real__ x, &sinrx, &cosrx);
+	}
+      else
+	{
+	  sinrx = __real__ x;
+	  cosrx = 1.0;
+	}
 
       if (fabs (__imag__ x) > t)
 	{