diff --git a/pymc/util.py b/pymc/util.py
index 1885c25d90..8bb2e35ec8 100644
--- a/pymc/util.py
+++ b/pymc/util.py
@@ -250,7 +250,7 @@ def dataset_to_point_list(
     }
     points = [
         {vn: stacked_dict[vn][i, ...] for vn in var_names}
-        for i in range(np.product([len(coords) for coords in stacked_dims.values()]))
+        for i in range(np.prod([len(coords) for coords in stacked_dims.values()]))
     ]
     # use the list of points
     return cast(List[Dict[str, np.ndarray]], points), stacked_dims