Skip to content
Snippets Groups Projects

Adaption to pystencils's new type system

Merged Markus Holzer requested to merge holzer/lbmpy:AdaptionToPystencils into master
Viewing commit 45bb586f
Show latest version
1 file
+ 8
12
Preferences
Compare changes
@@ -13,28 +13,24 @@ def test_poiseuille_channel_quicktest():
@@ -13,28 +13,24 @@ def test_poiseuille_channel_quicktest():
def test_entropic_methods():
def test_entropic_methods():
sc_kbc = create_lid_driven_cavity((20, 20), method=Method.TRT_KBC_N4,
sc_kbc = create_lid_driven_cavity((40, 40), method=Method.TRT_KBC_N4,
relaxation_rates=[1.9999, sp.Symbol("omega_free")],
relaxation_rates=[1.99, sp.Symbol("omega_free")],
 
lid_velocity=0.005,
entropic_newton_iterations=3, entropic=True, compressible=True,
entropic_newton_iterations=3, entropic=True, compressible=True,
zero_centered=False, force=(-1e-10, 0), force_model=ForceModel.LUO)
zero_centered=False, force=(-1e-10, 0), force_model=ForceModel.LUO)
sc_srt = create_lid_driven_cavity((40, 40), relaxation_rate=1.9999, lid_velocity=0.05, compressible=True,
sc_entropic = create_lid_driven_cavity((40, 40), method=Method.ENTROPIC_SRT, relaxation_rate=1.99,
zero_centered=False, force=(-1e-10, 0), force_model=ForceModel.LUO)
lid_velocity=0.005, compressible=True, zero_centered=False,
sc_entropic = create_lid_driven_cavity((40, 40), method=Method.ENTROPIC_SRT, relaxation_rate=1.9999,
lid_velocity=0.05, compressible=True, zero_centered=False,
force=(-1e-10, 0), force_model=ForceModel.LUO)
force=(-1e-10, 0), force_model=ForceModel.LUO)
sc_srt.run(1000)
sc_kbc.run(100)
sc_kbc.run(1000)
sc_entropic.run(100)
sc_entropic.run(1000)
assert np.isnan(np.max(sc_srt.velocity[:, :]))
assert np.isfinite(np.max(sc_kbc.velocity[:, :]))
assert np.isfinite(np.max(sc_kbc.velocity[:, :]))
assert np.isfinite(np.max(sc_entropic.velocity[:, :]))
assert np.isfinite(np.max(sc_entropic.velocity[:, :]))
def test_cumulant_ldc():
def test_cumulant_ldc():
sc_cumulant = create_lid_driven_cavity((20, 20), method=Method.CUMULANT, relaxation_rate=1.999999,
sc_cumulant = create_lid_driven_cavity((40, 40), method=Method.CUMULANT, relaxation_rate=1.999999,
compressible=True, force=(-1e-10, 0))
compressible=True, force=(-1e-10, 0))
sc_cumulant.run(100)
sc_cumulant.run(100)