mirrors/cpython
1
0
Fork 0
mirror of https://github.com/python/cpython.git synced 2025-09-26 18:29:57 +00:00
Code Issues Projects Releases Packages Wiki Activity Actions 10

Post-detabification cleanup.

Browse source
This commit is contained in:
Mark Dickinson 2010-07-05 20:14:26 +00:00
parent 5ada7c7315
commit 6c3bcb74b2
1 changed files with 8 additions and 8 deletions
Download patch file Download diff file Expand all files Collapse all files

16
Modules/_math.c
View file

@ -56,13 +56,13 @@ _Py_acosh(double x)
if (Py_IS_INFINITY(x)) { if (Py_IS_INFINITY(x)) {
return x+x; return x+x;
} else { } else {
return log(x)+ln2; /* acosh(huge)=log(2x) */ return log(x)+ln2; /* acosh(huge)=log(2x) */
} }
} }
else if (x == 1.) { else if (x == 1.) {
return 0.0; /* acosh(1) = 0 */ return 0.0; /* acosh(1) = 0 */
} }
else if (x > 2.) { /* 2 < x < 2**28 */ else if (x > 2.) { /* 2 < x < 2**28 */
double t = x*x; double t = x*x;
return log(2.0*x - 1.0 / (x + sqrt(t - 1.0))); return log(2.0*x - 1.0 / (x + sqrt(t - 1.0)));
} }
@ -94,7 +94,7 @@ _Py_asinh(double x)
return x+x; return x+x;
} }
if (absx < two_pow_m28) { /* |x| < 2**-28 */ if (absx < two_pow_m28) { /* |x| < 2**-28 */
return x; /* return x inexact except 0 */ return x; /* return x inexact except 0 */
} }
if (absx > two_pow_p28) { /* |x| > 2**28 */ if (absx > two_pow_p28) { /* |x| > 2**28 */
w = log(absx)+ln2; w = log(absx)+ln2;
@ -114,9 +114,9 @@ _Py_asinh(double x)
* Method : * Method :
* 1.Reduced x to positive by atanh(-x) = -atanh(x) * 1.Reduced x to positive by atanh(-x) = -atanh(x)
* 2.For x>=0.5 * 2.For x>=0.5
* 1 2x x * 1 2x x
* atanh(x) = --- * log(1 + -------) = 0.5 * log1p(2 * --------) * atanh(x) = --- * log(1 + -------) = 0.5 * log1p(2 * -------)
* 2 1 - x 1 - x * 2 1 - x 1 - x
* *
* For x<0.5 * For x<0.5
* atanh(x) = 0.5*log1p(2x+2x*x/(1-x)) * atanh(x) = 0.5*log1p(2x+2x*x/(1-x))
@ -194,7 +194,7 @@ _Py_log1p(double x)
/* For x small, we use the following approach. Let y be the nearest float /* For x small, we use the following approach. Let y be the nearest float
to 1+x, then to 1+x, then
1+x = y * (1 - (y-1-x)/y) 1+x = y * (1 - (y-1-x)/y)
so log(1+x) = log(y) + log(1-(y-1-x)/y). Since (y-1-x)/y is tiny, the so log(1+x) = log(y) + log(1-(y-1-x)/y). Since (y-1-x)/y is tiny, the
second term is well approximated by (y-1-x)/y. If abs(x) >= second term is well approximated by (y-1-x)/y. If abs(x) >=