diff --git a/python/dune/codegen/sumfact/symbolic.py b/python/dune/codegen/sumfact/symbolic.py
index 61867987ff21a6348fe7f567a63c64836f7a81d2..6e32fb838b22dd5fce8e7942bff266458440c2b1 100644
--- a/python/dune/codegen/sumfact/symbolic.py
+++ b/python/dune/codegen/sumfact/symbolic.py
@@ -389,11 +389,7 @@ class VectorSumfactKernelOutput(SumfactKernelInterfaceBase):
         substname = "haddsubst_{}".format("_".join([i.name for i in inames]))
         subst_rule(substname, (), result)
         result = prim.Call(prim.Variable(substname), ())
-
-        # If multiple horizontal_add's are to be performed with 'result'
-        # we need to precompute the result!
-        if len(outputs) > 1:
-            transform(lp.precompute, substname)
+        transform(lp.precompute, substname)
 
         deps = frozenset()
         for o in outputs: