lib/glibc
1
0
Fork 0
mirror of https://sourceware.org/git/glibc.git synced 2025-07-29 11:41:21 +03:00
Code Activity

Fix clog, clog10 spurious underflow exceptions (bug 14337).

Browse Source
This commit is contained in:
Joseph Myers
2012-07-09 11:06:34 +00:00
parent ad41a87fe0
commit 638a572eb0
11 changed files with 474 additions and 55 deletions
Download Patch File Download Diff File Expand all files Collapse all files

23
math/s_clog10.c
View File

@ -43,25 +43,30 @@ __clog10 (__complex__ double x)
else if (__builtin_expect (rcls != FP_NAN && icls != FP_NAN, 1))
{
/* Neither real nor imaginary part is NaN. */
double absx = fabs (__real__ x), absy = fabs (__imag__ x);
double d;
int scale = 0;
if (fabs (__real__ x) > DBL_MAX / 2.0
|| fabs (__imag__ x) > DBL_MAX / 2.0)
if (absx > DBL_MAX / 2.0)
{
scale = -1;
__real__ x = __scalbn (__real__ x, scale);
__imag__ x = __scalbn (__imag__ x, scale);
absx = __scalbn (absx, scale);
absy = (absy >= DBL_MIN * 2.0 ? __scalbn (absy, scale) : 0.0);
}
else if (fabs (__real__ x) < DBL_MIN
&& fabs (__imag__ x) < DBL_MIN)
else if (absy > DBL_MAX / 2.0)
{
scale = -1;
absx = (absx >= DBL_MIN * 2.0 ? __scalbn (absx, scale) : 0.0);
absy = __scalbn (absy, scale);
}
else if (absx < DBL_MIN && absy < DBL_MIN)
{
scale = DBL_MANT_DIG;
__real__ x = __scalbn (__real__ x, scale);
__imag__ x = __scalbn (__imag__ x, scale);
absx = __scalbn (absx, scale);
absy = __scalbn (absy, scale);
}
d = __ieee754_hypot (__real__ x, __imag__ x);
d = __ieee754_hypot (absx, absy);
__real__ result = __ieee754_log10 (d) - scale * M_LOG10_2;
__imag__ result = M_LOG10E * __ieee754_atan2 (__imag__ x, __real__ x);