diff --git a/loopy/statistics.py b/loopy/statistics.py
index 2c87b60789ca144d2c8d0ad5b91651bd9a42fb3d..2e061ec6d76505ad6373e3fad4462c0459545253 100755
--- a/loopy/statistics.py
+++ b/loopy/statistics.py
@@ -29,6 +29,7 @@ import warnings
 from islpy import dim_type
 import islpy as isl
 from pymbolic.mapper import CombineMapper
+from functools import reduce
 
 
 class TypeToCountMap: