From 8e777304b435f966d864696937ca4c6c7e2f14b2 Mon Sep 17 00:00:00 2001
From: Michael Kuron <mkuron@icp.uni-stuttgart.de>
Date: Wed, 14 Aug 2019 14:09:51 +0200
Subject: [PATCH] remove obsolete
fluctuatinglb.method_with_rescaled_equilibrium_values
---
lbmpy/creationfunctions.py | 5 +----
lbmpy/fluctuatinglb.py | 13 -------------
2 files changed, 1 insertion(+), 17 deletions(-)
diff --git a/lbmpy/creationfunctions.py b/lbmpy/creationfunctions.py
index 3e7474d..f7e30b3 100644
--- a/lbmpy/creationfunctions.py
+++ b/lbmpy/creationfunctions.py
@@ -177,7 +177,7 @@ from lbmpy.fieldaccess import (
AAEvenTimeStepAccessor, AAOddTimeStepAccessor, CollideOnlyInplaceAccessor,
EsoTwistEvenTimeStepAccessor, EsoTwistOddTimeStepAccessor, PdfFieldAccessor,
PeriodicTwoFieldsAccessor, StreamPullTwoFieldsAccessor, StreamPushTwoFieldsAccessor)
-from lbmpy.fluctuatinglb import add_fluctuations_to_collision_rule, method_with_rescaled_equilibrium_values
+from lbmpy.fluctuatinglb import add_fluctuations_to_collision_rule
from lbmpy.methods import (
create_mrt3, create_mrt_orthogonal, create_mrt_raw, create_srt, create_trt, create_trt_kbc)
from lbmpy.methods.creationfunctions import create_generic_mrt
@@ -307,9 +307,6 @@ def create_lb_collision_rule(lb_method=None, optimization={}, **kwargs):
if rho_in is not None and isinstance(rho_in, Field):
rho_in = rho_in.center
- if params['fluctuating']:
- lb_method = method_with_rescaled_equilibrium_values(lb_method)
-
if u_in is not None:
density_rhs = sum(lb_method.pre_collision_pdf_symbols) if rho_in is None else rho_in
eqs = [Assignment(cqc.zeroth_order_moment_symbol, density_rhs)]
diff --git a/lbmpy/fluctuatinglb.py b/lbmpy/fluctuatinglb.py
index 833b93f..ba7aed7 100644
--- a/lbmpy/fluctuatinglb.py
+++ b/lbmpy/fluctuatinglb.py
@@ -45,19 +45,6 @@ def fluctuating_variance_from_temperature(method, temperature, c_s_sq=sp.Symbol(
for norm, rr in zip(normalization_factors, method.relaxation_rates)]
-def method_with_rescaled_equilibrium_values(base_method):
- """Re-scales the equilibrium moments by 1 / sqrt(M*w) with moment matrix M and weights w"""
- from lbmpy.creationfunctions import create_lb_method_from_existing
-
- sig_k = abs(base_method.moment_matrix) * sp.Matrix(base_method.weights)
-
- def modification_rule(moment, eq, rr):
- i = base_method.moments.index(moment)
- return moment, eq / sp.sqrt(sig_k[i]), rr
-
- return create_lb_method_from_existing(base_method, modification_rule)
-
-
def fluctuation_correction(method, rng_generator, variances=SymbolGen("variance")):
"""Returns additive correction terms to be added to the the collided pdfs"""
conserved_moments = {sp.sympify(1), *MOMENT_SYMBOLS}
--
GitLab