[Python-checkins] Make fractional value accumulation consistent inside and outside the loop. (GH-22315)
Raymond Hettinger
webhook-mailer at python.org
Fri Sep 18 20:57:36 EDT 2020
https://github.com/python/cpython/commit/ec8a15b034124f3b58d1addda789fa4c20313006
commit: ec8a15b034124f3b58d1addda789fa4c20313006
branch: master
author: Raymond Hettinger <rhettinger at users.noreply.github.com>
committer: GitHub <noreply at github.com>
date: 2020-09-18T17:57:28-07:00
summary:
Make fractional value accumulation consistent inside and outside the loop. (GH-22315)
files:
M Modules/mathmodule.c
diff --git a/Modules/mathmodule.c b/Modules/mathmodule.c
index ecd291ecd1b4b..935759ec671ca 100644
--- a/Modules/mathmodule.c
+++ b/Modules/mathmodule.c
@@ -2550,8 +2550,7 @@ vector_norm(Py_ssize_t n, double *vec, double max, int found_nan)
assert(csum + lo * lo == csum);
frac_lo += lo * lo;
}
- frac += frac_lo + frac_mid;
- h = sqrt(csum - 1.0 + frac);
+ h = sqrt(csum - 1.0 + (frac_lo + frac_mid + frac));
x = h;
t = x * T27;
@@ -2569,15 +2568,15 @@ vector_norm(Py_ssize_t n, double *vec, double max, int found_nan)
assert(fabs(csum) >= fabs(x));
oldcsum = csum;
csum += x;
- frac += (oldcsum - csum) + x;
+ frac_mid += (oldcsum - csum) + x;
x = -lo * lo;
assert(fabs(csum) >= fabs(x));
oldcsum = csum;
csum += x;
- frac += (oldcsum - csum) + x;
+ frac_lo += (oldcsum - csum) + x;
- x = csum - 1.0 + frac;
+ x = csum - 1.0 + (frac_lo + frac_mid + frac);
return (h + x / (2.0 * h)) / scale;
}
/* When max_e < -1023, ldexp(1.0, -max_e) overflows.
More information about the Python-checkins
mailing list