numpy · eendebakpt · Jul 6, 2024 · Aug 31, 2024 · May 13, 2025
diff --git a/benchmarks/benchmarks/bench_polynomial.py b/benchmarks/benchmarks/bench_polynomial.py
@@ -1,14 +1,15 @@
 from .common import Benchmark
 
 import numpy as np
-
+from numpy.polynomial.polynomial import polyval
 
 class Polynomial(Benchmark):
 
     def setup(self):
         self.polynomial_degree2 = np.polynomial.Polynomial(np.array([1, 2]))
         self.array3 = np.linspace(0, 1, 3)
         self.array1000 = np.linspace(0, 1, 10_000)
+        self.array1M = np.linspace(0, 1, 1_000_000)
         self.float64 = np.float64(1.0)
 
     def time_polynomial_evaluation_scalar(self):
@@ -23,5 +24,12 @@ def time_polynomial_evaluation_array_3(self):
     def time_polynomial_evaluation_array_1000(self):
         self.polynomial_degree2(self.array1000)
 
+    def time_polynomial_evaluation_array_1_000_000(self):
+        self.polynomial_degree2(self.array1M)
+
+    def time_polyval(self):
+        polyval(self.array1M, self.polynomial_degree2.coef)
+
     def time_polynomial_addition(self):
         _ = self.polynomial_degree2 + self.polynomial_degree2
+
diff --git a/numpy/polynomial/polynomial.py b/numpy/polynomial/polynomial.py
@@ -751,7 +751,8 @@ def polyval(x, c, tensor=True):
 
     c0 = c[-1] + x * 0
     for i in range(2, len(c) + 1):
-        c0 = c[-i] + c0 * x
+        c0 *= x
+        c0 += c[-i]
     return c0