diff --git a/src/lbmpy/flow_statistics.py b/src/lbmpy/flow_statistics.py index 6d5433b21c297ed7e53176359a8266e628c4f064..3b828da3215d268d8cc1ab2a9736ca70d42936ed 100644 --- a/src/lbmpy/flow_statistics.py +++ b/src/lbmpy/flow_statistics.py @@ -70,7 +70,8 @@ def welford_assignments(field, mean_field, sum_of_products_field=None): if sum_of_products_field is not None: delta2 = sp.symbols(f"delta2_:{dim}") for i in range(dim): - main_assignments.append(ps.Assignment(delta2[i], welford_field.at_index(i) - welford_mean_field.at_index(i))) + main_assignments.append( + ps.Assignment(delta2[i], welford_field.at_index(i) - welford_mean_field.at_index(i))) for j in range(dim): idx = i * dim + j main_assignments.append(ps.Assignment(