diff options
| author | Liubov Dmitrieva <liubov.dmitrieva@gmail.com> | 2012-09-25 20:41:17 +0200 |
|---|---|---|
| committer | Andreas Jaeger <aj@suse.de> | 2012-09-25 20:41:17 +0200 |
| commit | 6d3bf1993190edc502d01e8ca42c4482b20a5b6c (patch) | |
| tree | ceef4ffabe768f560a3a4d3f6b556b7b70aac37d /math/s_ccoshl.c | |
| parent | 203e56032f77d0b36701deaee360e9cdda35d54e (diff) | |
| download | glibc-6d3bf1993190edc502d01e8ca42c4482b20a5b6c.zip glibc-6d3bf1993190edc502d01e8ca42c4482b20a5b6c.tar.gz glibc-6d3bf1993190edc502d01e8ca42c4482b20a5b6c.tar.bz2 | |
Fix wrong ussage of sincos for subnormal arguments
Diffstat (limited to 'math/s_ccoshl.c')
| -rw-r--r-- | math/s_ccoshl.c | 20 |
1 files changed, 18 insertions, 2 deletions
diff --git a/math/s_ccoshl.c b/math/s_ccoshl.c index 77a9ae3..1d561b4 100644 --- a/math/s_ccoshl.c +++ b/math/s_ccoshl.c @@ -39,7 +39,15 @@ __ccoshl (__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 (fabsl (__real__ x) > t) { @@ -92,7 +100,15 @@ __ccoshl (__complex__ long double x) /* Imaginary part is finite. */ 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 (HUGE_VALL, cosix); __imag__ retval = (__copysignl (HUGE_VALL, sinix) |