aboutsummaryrefslogtreecommitdiff
path: root/math/s_csinl.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_csinl.c
parent203e56032f77d0b36701deaee360e9cdda35d54e (diff)
downloadglibc-6d3bf1993190edc502d01e8ca42c4482b20a5b6c.tar
glibc-6d3bf1993190edc502d01e8ca42c4482b20a5b6c.tar.gz
glibc-6d3bf1993190edc502d01e8ca42c4482b20a5b6c.tar.bz2
glibc-6d3bf1993190edc502d01e8ca42c4482b20a5b6c.zip
Fix wrong ussage of sincos for subnormal arguments
Diffstat (limited to 'math/s_csinl.c')
-rw-r--r--math/s_csinl.c20
1 files changed, 18 insertions, 2 deletions
diff --git a/math/s_csinl.c b/math/s_csinl.c
index ff43256e54..9812bddd5c 100644
--- a/math/s_csinl.c
+++ b/math/s_csinl.c
@@ -42,7 +42,15 @@ __csinl (__complex__ long double x)
const int t = (int) ((LDBL_MAX_EXP - 1) * M_LN2l);
long double sinix, cosix;
- __sincosl (__real__ x, &sinix, &cosix);
+ if (__builtin_expect (rcls != FP_SUBNORMAL, 1))
+ {
+ __sincosl (__real__ x, &sinix, &cosix);
+ }
+ else
+ {
+ sinix = __real__ x;
+ cosix = 1.0;
+ }
if (fabsl (__imag__ x) > t)
{
@@ -115,7 +123,15 @@ __csinl (__complex__ long double x)
/* Real part is finite. */
long double sinix, cosix;
- __sincosl (__real__ x, &sinix, &cosix);
+ if (__builtin_expect (rcls != FP_SUBNORMAL, 1))
+ {
+ __sincosl (__real__ x, &sinix, &cosix);
+ }
+ else
+ {
+ sinix = __real__ x;
+ cosix = 1.0;
+ }
__real__ retval = __copysignl (HUGE_VALL, sinix);
__imag__ retval = __copysignl (HUGE_VALL, cosix);