Mercurial > cpython-lru_cache-weakref

diff Modules/_math.c @ 57869:1a25b420ec00

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merged revisions 76982 via svnmerge from svn+ssh://pythondev@svn.python.org/python/trunk ........ r76982 | mark.dickinson | 2009-12-21 15:40:33 +0000 (Mon, 21 Dec 2009) | 2 lines Inverse hyperbolic trigonometric functions should call m_log1p, not log1p. ........
author Mark Dickinson <dickinsm@gmail.com>
date Mon, 21 Dec 2009 15:42:00 +0000
parents 9751ddcc1882
children fd25030223b0
line wrap: on
line diff
--- a/Modules/_math.c	Mon Dec 21 15:27:41 2009 +0000
+++ b/Modules/_math.c	Mon Dec 21 15:42:00 2009 +0000
@@ -3,6 +3,7 @@
 
 #include "Python.h"
 #include <float.h>
+#include "_math.h"
 
 /* The following copyright notice applies to the original
    implementations of acosh, asinh and atanh. */
@@ -67,7 +68,7 @@
 	}
 	else {				/* 1 < x <= 2 */
 		double t = x - 1.0;
-		return log1p(t + sqrt(2.0*t + t*t));
+		return m_log1p(t + sqrt(2.0*t + t*t));
 	}
 }
 
@@ -103,7 +104,7 @@
 	}
 	else {				/* 2**-28 <= |x| < 2= */
 		double t = x*x;
-		w = log1p(absx + t / (1.0 + sqrt(1.0 + t)));
+		w = m_log1p(absx + t / (1.0 + sqrt(1.0 + t)));
 	}
 	return copysign(w, x);
 	
@@ -149,10 +150,10 @@
 	}
 	if (absx < 0.5) {		/* |x| < 0.5 */
 		t = absx+absx;
-		t = 0.5 * log1p(t + t*absx / (1.0 - absx));
+		t = 0.5 * m_log1p(t + t*absx / (1.0 - absx));
 	} 
 	else {				/* 0.5 <= |x| <= 1.0 */
-		t = 0.5 * log1p((absx + absx) / (1.0 - absx));
+		t = 0.5 * m_log1p((absx + absx) / (1.0 - absx));
 	}
 	return copysign(t, x);
 }