fixup

python/dune/perftool/sumfact/basis.py

@@ -90,6 +90,9 @@ def pymbolic_coefficient_gradient(element, restriction, component, coeff_func, v
         from dune.perftool.sumfact.vectorization import attach_vectorization_info
         sf = attach_vectorization_info(sf)
 
+        if i != sf.index:
+            direct_indexing_is_possible = False
+
         # Add a sum factorization kernel that implements the
         # evaluation of the gradients of basis functions at quadrature
         # points (stage 1)