Try a different fix

python/dune/perftool/sumfact/quadrature.py

@@ -219,7 +219,7 @@ def pymbolic_quadrature_position(index, visitor):
     # Precompute it in the constructor
     inames = constructor_quadrature_inames(dim)
     assignee = prim.Subscript(prim.Variable(name), tuple(prim.Variable(i) for i in inames))
-    expression = Subscript(Variable(name_oned_quadrature_points(qps_per_dir[index])),
+    expression = Subscript(Variable(name_oned_quadrature_points(local_qps_per_dir[index])),
                            (prim.Variable(inames[index])))
     instruction(assignee=assignee,
                 expression=expression,