diff --git a/src/lbmpy/flow_statistics.py b/src/lbmpy/flow_statistics.py
index 3b828da3215d268d8cc1ab2a9736ca70d42936ed..3e5b5f5f6122d6998416a46b92669c9fa6eaf95f 100644
--- a/src/lbmpy/flow_statistics.py
+++ b/src/lbmpy/flow_statistics.py
@@ -72,6 +72,7 @@ def welford_assignments(field, mean_field, sum_of_products_field=None):
         for i in range(dim):
             main_assignments.append(
                 ps.Assignment(delta2[i], welford_field.at_index(i) - welford_mean_field.at_index(i)))
+        for i in range(dim):
             for j in range(dim):
                 idx = i * dim + j
                 main_assignments.append(ps.Assignment(