Skip to content
Snippets Groups Projects
Commit 15c57cd4 authored by Markus Holzer's avatar Markus Holzer
Browse files

removed mrt3 from each test case

parent 76062f65
Branches
No related merge requests found
...@@ -12,6 +12,7 @@ from pystencils.fast_approximation import insert_fast_sqrts, insert_fast_divisio ...@@ -12,6 +12,7 @@ from pystencils.fast_approximation import insert_fast_sqrts, insert_fast_divisio
from lbmpy.macroscopic_value_kernels import macroscopic_values_getter, macroscopic_values_setter from lbmpy.macroscopic_value_kernels import macroscopic_values_getter, macroscopic_values_setter
omega = sp.symbols("omega") omega = sp.symbols("omega")
omega_free = sp.Symbol("omega_free")
omega_fill = sp.symbols("omega_:10") omega_fill = sp.symbols("omega_:10")
compile_time_block_size = False compile_time_block_size = False
...@@ -46,23 +47,18 @@ options_dict = { ...@@ -46,23 +47,18 @@ options_dict = {
'stencil': 'D3Q19', 'stencil': 'D3Q19',
'relaxation_rates': [omega_fill[0], omega, omega_fill[1], omega_fill[2], omega_fill[3], omega_fill[4], omega_fill[5]], 'relaxation_rates': [omega_fill[0], omega, omega_fill[1], omega_fill[2], omega_fill[3], omega_fill[4], omega_fill[5]],
}, },
'mrt3': {
'method': 'mrt3',
'stencil': 'D3Q19',
'relaxation_rates': [omega, 1.1, 1.2],
},
'entropic': { 'entropic': {
'method': 'mrt3', 'method': 'mrt',
'stencil': 'D3Q19', 'stencil': 'D3Q19',
'compressible': True, 'compressible': True,
'relaxation_rates': [omega, omega, sp.Symbol("omega_free")], 'relaxation_rates': [omega, omega, omega_free, omega_free, omega_free, omega_free],
'entropic': True, 'entropic': True,
}, },
'entropic_kbc_n4': { 'entropic_kbc_n4': {
'method': 'trt-kbc-n4', 'method': 'trt-kbc-n4',
'stencil': 'D3Q27', 'stencil': 'D3Q27',
'compressible': True, 'compressible': True,
'relaxation_rates': [omega, sp.Symbol("omega_free")], 'relaxation_rates': [omega, omega_free],
'entropic': True, 'entropic': True,
}, },
'smagorinsky': { 'smagorinsky': {
......
...@@ -10,6 +10,7 @@ from pystencils.fast_approximation import insert_fast_sqrts, insert_fast_divisio ...@@ -10,6 +10,7 @@ from pystencils.fast_approximation import insert_fast_sqrts, insert_fast_divisio
from lbmpy.macroscopic_value_kernels import macroscopic_values_getter, macroscopic_values_setter from lbmpy.macroscopic_value_kernels import macroscopic_values_getter, macroscopic_values_setter
omega = sp.symbols("omega") omega = sp.symbols("omega")
omega_free = sp.Symbol("omega_free")
compile_time_block_size = False compile_time_block_size = False
if compile_time_block_size: if compile_time_block_size:
...@@ -39,10 +40,10 @@ options_dict = { ...@@ -39,10 +40,10 @@ options_dict = {
'relaxation_rates': [omega, 1.3, 1.4, omega, 1.2, 1.1], 'relaxation_rates': [omega, 1.3, 1.4, omega, 1.2, 1.1],
}, },
'entropic': { 'entropic': {
'method': 'mrt3', 'method': 'mrt',
'stencil': 'D3Q19', 'stencil': 'D3Q19',
'compressible': True, 'compressible': True,
'relaxation_rates': [omega, omega, sp.Symbol("omega_free")], 'relaxation_rates': [omega, omega, omega_free, omega_free, omega_free],
'entropic': True, 'entropic': True,
}, },
'smagorinsky': { 'smagorinsky': {
......
...@@ -2,7 +2,7 @@ waLBerla_link_files_to_builddir( "*.prm" ) ...@@ -2,7 +2,7 @@ waLBerla_link_files_to_builddir( "*.prm" )
waLBerla_link_files_to_builddir( "*.py" ) waLBerla_link_files_to_builddir( "*.py" )
foreach(config trt smagorinsky mrt3 mrt entropic_kbc_n4 cumulant ) foreach(config trt smagorinsky mrt entropic_kbc_n4 cumulant )
waLBerla_generate_target_from_python(NAME UniformGridGenerated_${config} waLBerla_generate_target_from_python(NAME UniformGridGenerated_${config}
CODEGEN_CFG ${config} CODEGEN_CFG ${config}
FILE UniformGridGenerated.py FILE UniformGridGenerated.py
......
...@@ -6,6 +6,7 @@ from lbmpy.macroscopic_value_kernels import macroscopic_values_getter, macroscop ...@@ -6,6 +6,7 @@ from lbmpy.macroscopic_value_kernels import macroscopic_values_getter, macroscop
from lbmpy.fieldaccess import AAEvenTimeStepAccessor, AAOddTimeStepAccessor from lbmpy.fieldaccess import AAEvenTimeStepAccessor, AAOddTimeStepAccessor
omega = sp.symbols("omega") omega = sp.symbols("omega")
omega_free = sp.Symbol("omega_free")
omega_fill = sp.symbols("omega_:10") omega_fill = sp.symbols("omega_:10")
options_dict = { options_dict = {
...@@ -31,23 +32,18 @@ options_dict = { ...@@ -31,23 +32,18 @@ options_dict = {
'stencil': 'D3Q19', 'stencil': 'D3Q19',
'relaxation_rates': [omega_fill[0], omega, omega_fill[1], omega_fill[2], omega_fill[3], omega_fill[4], omega_fill[5]], 'relaxation_rates': [omega_fill[0], omega, omega_fill[1], omega_fill[2], omega_fill[3], omega_fill[4], omega_fill[5]],
}, },
'mrt3': {
'method': 'mrt3',
'stencil': 'D3Q19',
'relaxation_rates': [omega, 1.1, 1.2],
},
'entropic': { 'entropic': {
'method': 'mrt3', 'method': 'mrt',
'stencil': 'D3Q19', 'stencil': 'D3Q19',
'compressible': True, 'compressible': True,
'relaxation_rates': [omega, omega, sp.Symbol("omega_free")], 'relaxation_rates': [omega, omega, omega_free, omega_free, omega_free],
'entropic': True, 'entropic': True,
}, },
'entropic_kbc_n4': { 'entropic_kbc_n4': {
'method': 'trt-kbc-n4', 'method': 'trt-kbc-n4',
'stencil': 'D3Q27', 'stencil': 'D3Q27',
'compressible': True, 'compressible': True,
'relaxation_rates': [omega, sp.Symbol("omega_free")], 'relaxation_rates': [omega, omega_free],
'entropic': True, 'entropic': True,
}, },
'smagorinsky': { 'smagorinsky': {
......
...@@ -12,7 +12,7 @@ with CodeGeneration() as ctx: ...@@ -12,7 +12,7 @@ with CodeGeneration() as ctx:
# the collision rule of the LB method where the some advanced features # the collision rule of the LB method where the some advanced features
collision_rule = create_lb_collision_rule( collision_rule = create_lb_collision_rule(
stencil='D3Q19', compressible=True, stencil='D3Q19', compressible=True,
method='mrt3', relaxation_rates=[omega, omega, omega_free], method='mrt', relaxation_rates=[omega, omega, omega_free, omega_free, omega_free, omega_free],
entropic=True, # entropic method where second omega is chosen s.t. entropy condition entropic=True, # entropic method where second omega is chosen s.t. entropy condition
omega_output_field=omega_out, # scalar field where automatically chosen omega of entropic or Smagorinsky method is written to omega_output_field=omega_out, # scalar field where automatically chosen omega of entropic or Smagorinsky method is written to
force=force_field.center_vector, # read forces for each lattice cell from an external force field force=force_field.center_vector, # read forces for each lattice cell from an external force field
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment