loghelper(): Try to nudge the compiler into doing mults in an order that
minimizes roundoff error.
diff --git a/Modules/mathmodule.c b/Modules/mathmodule.c
index eef8b78..c206ddc 100644
--- a/Modules/mathmodule.c
+++ b/Modules/mathmodule.c
@@ -263,7 +263,7 @@
log(x) + log(2) * e * SHIFT.
CAUTION: e*SHIFT may overflow using int arithmetic,
so force use of double. */
- x = func(x) + func(2.0) * (double)e * (double)SHIFT;
+ x = func(x) + (e * (double)SHIFT) * func(2.0);
return PyFloat_FromDouble(x);
}