Skip to content
Snippets Groups Projects
Commit 8a2d03d6 authored by Michael Kuron's avatar Michael Kuron :mortar_board:
Browse files

need to include collision in Schiller vs. Guo

parent 020b057a
No related branches found
No related tags found
1 merge request!76Draft: Fix Schiller. Schiller vs Guo force test
Pipeline #32017 failed
...@@ -7,6 +7,8 @@ from lbmpy.moments import is_bulk_moment ...@@ -7,6 +7,8 @@ from lbmpy.moments import is_bulk_moment
import pytest import pytest
from contextlib import ExitStack as does_not_raise from contextlib import ExitStack as does_not_raise
from test_conserved_quantity_relaxation_invariance import check_for_collision_rule_equivalence
force_models = [fm.lower() for fm in dir(lbmpy.forcemodels) if fm[0].isupper()] force_models = [fm.lower() for fm in dir(lbmpy.forcemodels) if fm[0].isupper()]
...@@ -176,10 +178,7 @@ def test_schiller_vs_guo(stencil): ...@@ -176,10 +178,7 @@ def test_schiller_vs_guo(stencil):
force_moments_g = sp.simplify(method_g.moment_matrix * sp.Matrix(method_g.force_model(method_g))) force_moments_g = sp.simplify(method_g.moment_matrix * sp.Matrix(method_g.force_model(method_g)))
diff = sp.simplify(force_moments_s - force_moments_g) diff = sp.simplify(force_moments_s - force_moments_g)
print(method_s.moments) assert check_for_collision_rule_equivalence(method_s.get_collision_rule(), method_g.get_collision_rule())
print(diff)
assert diff.transpose() == 0*diff.transpose()
@pytest.mark.parametrize("force_model", force_models) @pytest.mark.parametrize("force_model", force_models)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment