diff options
author | Adrian Kummerlaender | 2019-10-21 18:42:24 +0200 |
---|---|---|
committer | Adrian Kummerlaender | 2019-10-21 18:48:38 +0200 |
commit | 82a44e0d64afb8818ea98d68dc08108885d503c2 (patch) | |
tree | 6e8f08acd83b2886cd296ed3831acc83e309906c /boltzgen/lbm/model/D3Q19.py | |
download | boltzgen-82a44e0d64afb8818ea98d68dc08108885d503c2.tar boltzgen-82a44e0d64afb8818ea98d68dc08108885d503c2.tar.gz boltzgen-82a44e0d64afb8818ea98d68dc08108885d503c2.tar.bz2 boltzgen-82a44e0d64afb8818ea98d68dc08108885d503c2.tar.lz boltzgen-82a44e0d64afb8818ea98d68dc08108885d503c2.tar.xz boltzgen-82a44e0d64afb8818ea98d68dc08108885d503c2.tar.zst boltzgen-82a44e0d64afb8818ea98d68dc08108885d503c2.zip |
Pull in basics from symlbm_playground
It's time to extract the generator-part of my GPU LBM playground and turn it
into a nice reusable library. The goal is to produce a framework that can be
used to generate collision and streaming programs from symbolic descriptions.
i.e. it should be possible to select a LB model, the desired boundary
conditions as well as a data structure / streaming model and use this
information to automatically generate matching OpenCL / CUDA / C++
programs.
Diffstat (limited to 'boltzgen/lbm/model/D3Q19.py')
-rw-r--r-- | boltzgen/lbm/model/D3Q19.py | 18 |
1 files changed, 18 insertions, 0 deletions
diff --git a/boltzgen/lbm/model/D3Q19.py b/boltzgen/lbm/model/D3Q19.py new file mode 100644 index 0000000..e9e6eec --- /dev/null +++ b/boltzgen/lbm/model/D3Q19.py @@ -0,0 +1,18 @@ +from sympy import Matrix, Rational, sqrt + +d = 3 +q = 19 + +c = [ Matrix(x) for x in [ + ( 0, 1, 1), (-1, 0, 1), ( 0, 0, 1), ( 1, 0, 1), ( 0, -1, 1), + (-1, 1, 0), ( 0, 1, 0), ( 1, 1, 0), (-1, 0, 0), ( 0, 0, 0), ( 1, 0, 0), (-1,-1, 0), ( 0, -1, 0), ( 1, -1, 0), + ( 0, 1,-1), (-1, 0,-1), ( 0, 0,-1), ( 1, 0,-1), ( 0, -1,-1) +]] + +w = [Rational(*x) for x in [ + (1,36), (1,36), (1,18), (1,36), (1,36), + (1,36), (1,18), (1,36), (1,18), (1,3), (1,18), (1,36), (1,18), (1,36), + (1,36), (1,36), (1,18), (1,36), (1,36) +]] + +c_s = sqrt(Rational(1,3)) |