about summary refs log tree commit diff
path: root/math/s_cexpl.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_cexpl.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_cexpl.c')
-rw-r--r--math/s_cexpl.c20
1 files changed, 18 insertions, 2 deletions
diff --git a/math/s_cexpl.c b/math/s_cexpl.c
index 256824924f..d67cc5e5eb 100644
--- a/math/s_cexpl.c
+++ b/math/s_cexpl.c
@@ -39,7 +39,15 @@ __cexpl (__complex__ long double x)
 	  const int t = (int) ((LDBL_MAX_EXP - 1) * M_LN2l);
 	  long double sinix, cosix;
 
-	  __sincosl (__imag__ x, &sinix, &cosix);
+	  if (__builtin_expect (icls != FP_SUBNORMAL, 1))
+	    {
+	      __sincosl (__imag__ x, &sinix, &cosix);
+	    }
+	  else
+	    {
+	      sinix = __imag__ x;
+	      cosix = 1.0;
+	    }
 
 	  if (__real__ x > t)
 	    {
@@ -95,7 +103,15 @@ __cexpl (__complex__ long double x)
 	    {
 	      long double sinix, cosix;
 
-	      __sincosl (__imag__ x, &sinix, &cosix);
+	      if (__builtin_expect (icls != FP_SUBNORMAL, 1))
+	        {
+		  __sincosl (__imag__ x, &sinix, &cosix);
+	        }
+	      else
+		{
+		  sinix = __imag__ x;
+		  cosix = 1.0;
+		}
 
 	      __real__ retval = __copysignl (value, cosix);
 	      __imag__ retval = __copysignl (value, sinix);