From 6d3bf1993190edc502d01e8ca42c4482b20a5b6c Mon Sep 17 00:00:00 2001 From: Liubov Dmitrieva Date: Tue, 25 Sep 2012 20:41:17 +0200 Subject: Fix wrong ussage of sincos for subnormal arguments --- math/s_ctanhf.c | 10 +++++++++- 1 file changed, 9 insertions(+), 1 deletion(-) (limited to 'math/s_ctanhf.c') diff --git a/math/s_ctanhf.c b/math/s_ctanhf.c index ca36a83bfb..862845f8c2 100644 --- a/math/s_ctanhf.c +++ b/math/s_ctanhf.c @@ -57,7 +57,15 @@ __ctanhf (__complex__ float x) /* tanh(x+iy) = (sinh(2x) + i*sin(2y))/(cosh(2x) + cos(2y)) = (sinh(x)*cosh(x) + i*sin(y)*cos(y))/(sinh(x)^2 + cos(y)^2). */ - __sincosf (__imag__ x, &sinix, &cosix); + if (__builtin_expect (fpclassify(__imag__ x) != FP_SUBNORMAL, 1)) + { + __sincosf (__imag__ x, &sinix, &cosix); + } + else + { + sinix = __imag__ x; + cosix = 1.0f; + } if (fabsf (__real__ x) > t) { -- cgit v1.2.3