diff --git a/python/dune/codegen/sumfact/realization.py b/python/dune/codegen/sumfact/realization.py
index 14cd0b1d0fe3814f294c22bc00f25260f4ee0233..5317268fdad3afb7dbcc406d17e7e145a24c41dd 100644
--- a/python/dune/codegen/sumfact/realization.py
+++ b/python/dune/codegen/sumfact/realization.py
@@ -278,6 +278,10 @@ def realize_sumfact_kernel_function(sf):
                                               )
                                   })
 
+    # Register kernel transformations
+    from dune.codegen.sumfact.transformations import reorder_loops_in_tensor_contraction
+    transform(reorder_loops_in_tensor_contraction, 'lkji')
+
     # Construct a loopy kernel object
     from dune.codegen.pdelab.localoperator import extract_kernel_from_cache
     args = ("const char* buffer0", "const char* buffer1") + sf.interface.signature_args