The pow() variant further improves accuracy (GH-102866)

rhettinger · web-flow · commit 1a5a14183ec8 · 2023-03-21T00:02:14.000-05:00
diff --git a/Doc/library/itertools.rst b/Doc/library/itertools.rst
@@ -876,7 +876,7 @@ which incur interpreter overhead.
        n = len(coefficients)
        if n == 0:
            return x * 0  # coerce zero to the type of x
-       powers = accumulate(repeat(x, n - 1), operator.mul, initial=1)
+       powers = map(pow, repeat(x), range(n))
        return math.sumprod(reversed(coefficients), powers)
 
    def polynomial_from_roots(roots):
