From d4ee5d5567be14c32fd4621c85b1b152394a0912 Mon Sep 17 00:00:00 2001
From: Michael Kuron <mkuron@icp.uni-stuttgart.de>
Date: Wed, 6 Nov 2019 11:03:42 +0100
Subject: [PATCH] Fluctuating MRT: use the correct prefactors
---
lbmpy/fluctuatinglb.py | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/lbmpy/fluctuatinglb.py b/lbmpy/fluctuatinglb.py
index ba7aed7..b96742d 100644
--- a/lbmpy/fluctuatinglb.py
+++ b/lbmpy/fluctuatinglb.py
@@ -35,8 +35,9 @@ def add_fluctuations_to_collision_rule(collision_rule, temperature=None, varianc
def fluctuating_variance_from_temperature(method, temperature, c_s_sq=sp.Symbol("c_s") ** 2):
- """Produces variance equations according to (3.54) in Schiller08"""
- normalization_factors = abs(method.moment_matrix) * sp.Matrix(method.weights)
+ """Produces variance equations according to (2.60) and (3.54) in Schiller08"""
+ normalization_factors = sp.matrix_multiply_elementwise(method.moment_matrix, method.moment_matrix) * \
+ sp.Matrix(method.weights)
density = method.zeroth_order_equilibrium_moment_symbol
if method.conserved_quantity_computation.zero_centered_pdfs:
density += 1
--
GitLab