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_ctanf.c | |
| parent | 203e56032f77d0b36701deaee360e9cdda35d54e (diff) | |
| download | glibc-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_ctanf.c')
| -rw-r--r-- | math/s_ctanf.c | 10 |
1 files changed, 9 insertions, 1 deletions
diff --git a/math/s_ctanf.c b/math/s_ctanf.c index 2559f83f84..fd2b797882 100644 --- a/math/s_ctanf.c +++ b/math/s_ctanf.c @@ -57,7 +57,15 @@ __ctanf (__complex__ float 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). */ - __sincosf (__real__ x, &sinrx, &cosrx); + if (__builtin_expect (fpclassify(__real__ x) != FP_SUBNORMAL, 1)) + { + __sincosf (__real__ x, &sinrx, &cosrx); + } + else + { + sinrx = __real__ x; + cosrx = 1.0f; + } if (fabsf (__imag__ x) > t) { |