blob: b38898a1c9b70043c57ac9829b6c614afe98efef (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
|
import sympy
from mako.template import Template
from pathlib import Path
from simulation import Geometry
from symbolic.generator import LBM
import symbolic.D2Q9 as D2Q9
lbm = LBM(D2Q9)
moments = lbm.moments(optimize = False)
collide = lbm.bgk(f_eq = lbm.equilibrium(), tau = 0.6)
geometry = Geometry(512, 512)
program_src = Template(filename = str(Path(__file__).parent/'template/standalone.mako')).render(
descriptor = lbm.descriptor,
geometry = geometry,
steps = 100,
moments_subexpr = moments[0],
moments_assignment = moments[1],
collide_subexpr = collide[0],
collide_assignment = collide[1],
float_type = 'double',
ccode = sympy.ccode
)
print(program_src)
|