MFC: 1.11: fix the threshold for (not) using the simple Taylor approximation.

This commit is contained in:
Bruce Evans 2007-06-14 05:51:00 +00:00
parent 9624afb72c
commit d382c5ebb4
Notes: svn2git 2020-12-20 02:59:44 +00:00
svn path=/head/; revision=170706

View File

@ -107,7 +107,7 @@ __ieee754_log(double x)
SET_HIGH_WORD(x,hx|(i^0x3ff00000)); /* normalize x or x/2 */
k += (i>>20);
f = x-1.0;
if((0x000fffff&(2+hx))<3) { /* -2**-20 <= f < 2**-20 */
if((0x000fffff&(2+hx))<3) { /* |f| < 2**-20 */
if(f==zero) if(k==0) return zero; else {dk=(double)k;
return dk*ln2_hi+dk*ln2_lo;}
R = f*f*(0.5-0.33333333333333333*f);