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 inaccuracy of clog, clog10 near |z| = 1 (bug 13629).

Browse Source
This commit is contained in:
Joseph Myers
2012-09-25 19:43:49 +00:00
parent c8cb2a5224
commit d032e0d29b
19 changed files with 1003 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13
math/s_clog10.c
View File

@ -88,6 +88,19 @@ __clog10 (__complex__ double x)
d2m1 += absy * absy;
__real__ result = __log1p (d2m1) * (M_LOG10E / 2.0);
}
else if (absx < 1.0
&& absx >= 0.75
&& absy < DBL_EPSILON / 2.0
&& scale == 0)
{
double d2m1 = (absx - 1.0) * (absx + 1.0);
__real__ result = __log1p (d2m1) * (M_LOG10E / 2.0);
}
else if (absx < 1.0 && (absx >= 0.75 || absy >= 0.5) && scale == 0)
{
double d2m1 = __x2y2m1 (absx, absy);
__real__ result = __log1p (d2m1) * (M_LOG10E / 2.0);
}
else
{
double d = __ieee754_hypot (absx, absy);