aboutsummaryrefslogtreecommitdiff
path: root/implosion.py
diff options
context:
space:
mode:
authorAdrian Kummerlaender2019-06-13 21:19:49 +0200
committerAdrian Kummerlaender2019-06-13 21:19:49 +0200
commit61cd3b6bc53c151863bd99fb2a03b29157a876ca (patch)
treea50eef491e1b2fdfac1b80635664cfbd62175771 /implosion.py
parent5f98308e2035edff2719ba852d5cc74ab449528f (diff)
downloadsymlbm_playground-61cd3b6bc53c151863bd99fb2a03b29157a876ca.tar
symlbm_playground-61cd3b6bc53c151863bd99fb2a03b29157a876ca.tar.gz
symlbm_playground-61cd3b6bc53c151863bd99fb2a03b29157a876ca.tar.bz2
symlbm_playground-61cd3b6bc53c151863bd99fb2a03b29157a876ca.tar.lz
symlbm_playground-61cd3b6bc53c151863bd99fb2a03b29157a876ca.tar.xz
symlbm_playground-61cd3b6bc53c151863bd99fb2a03b29157a876ca.tar.zst
symlbm_playground-61cd3b6bc53c151863bd99fb2a03b29157a876ca.zip
Tidy up symbolic kernel generation
Diffstat (limited to 'implosion.py')
-rw-r--r--implosion.py10
1 files changed, 9 insertions, 1 deletions
diff --git a/implosion.py b/implosion.py
index e5e87bc..fb101b4 100644
--- a/implosion.py
+++ b/implosion.py
@@ -7,6 +7,8 @@ matplotlib.use('AGG')
from D2Q9 import Lattice
+import symbolic.D2Q9 as D2Q9
+
def MLUPS(cells, steps, time):
return cells * steps / time * 1e-6
@@ -56,7 +58,13 @@ moments = []
print("Initializing simulation...\n")
-lattice = Lattice(nX = 1024, nY = 1024, tau = 0.8, geometry = box, pop_eq_src = pop_eq, boundary_src = boundary)
+lattice = Lattice(
+ nX = 1024, nY = 1024,
+ geometry = box,
+ moments = D2Q9.moments(optimize = False),
+ collide = D2Q9.bgk(tau = 0.8),
+ pop_eq_src = pop_eq,
+ boundary_src = boundary)
print("Starting simulation using %d cells...\n" % lattice.nCells)