diff --git a/python/dune/perftool/sumfact/sumfact.py b/python/dune/perftool/sumfact/sumfact.py
index 1d665e56897cd05efdbb84b5089fe35fd6bf3cb8..41134538abcf8bf53a25b43acc34ff3d5b3b9314 100644
--- a/python/dune/perftool/sumfact/sumfact.py
+++ b/python/dune/perftool/sumfact/sumfact.py
@@ -245,7 +245,11 @@ def generate_accumulation_instruction(visitor, accterm, measure, subdomain_id):
                                )
 
         # In the case of FastDGGridOperator we can write directly into the resiudal/jacobi
-        if get_option('fastdg'):
+        #
+        # TODO: At the moment this only works if we do not vectorize
+        # (over gradients) because loopy tries to acces a vectorclass
+        # variable.
+        if get_option('fastdg') and not index:
             ft = get_global_context_value("form_type")
             if ft=='residual':
                 accum = accum + ".data()"
@@ -255,7 +259,7 @@ def generate_accumulation_instruction(visitor, accterm, measure, subdomain_id):
                 expression = prim.Sum((assignee,result))
                 instruction(assignee=assignee,
                             expression=expression,
-                            forced_iname_deps=frozenset(inames + visitor.inames + vecinames),
+                            forced_iname_deps=frozenset(inames),
                             forced_iname_deps_is_final=True,
                             depends_on=insn_dep,
                             )
@@ -268,7 +272,7 @@ def generate_accumulation_instruction(visitor, accterm, measure, subdomain_id):
                 expression = prim.Sum((assignee,result))
                 instruction(assignee=assignee,
                             expression=expression,
-                            forced_iname_deps=frozenset(inames + visitor.inames + vecinames),
+                            forced_iname_deps=frozenset(inames + visitor.inames),
                             forced_iname_deps_is_final=True,
                             depends_on=insn_dep,
                             )