From 2f4794a0ddc0ff472ad3f0be20d7a15de72cfae4 Mon Sep 17 00:00:00 2001 From: Adrian Kummerlaender Date: Tue, 2 Jul 2019 21:14:14 +0200 Subject: Determine discrete velocities of D2Q9 and D3Q27 --- symbolic/D2Q9.py | 7 ++++--- symbolic/D3Q19.py | 4 ++-- symbolic/D3Q27.py | 11 ++++------- 3 files changed, 10 insertions(+), 12 deletions(-) (limited to 'symbolic') diff --git a/symbolic/D2Q9.py b/symbolic/D2Q9.py index 9693477..e3ac9de 100644 --- a/symbolic/D2Q9.py +++ b/symbolic/D2Q9.py @@ -1,6 +1,7 @@ -from sympy import * +from sympy import Matrix +from itertools import product -q = 9 d = 2 +q = 9 -c = [ Matrix(x) for x in [(-1, 1), ( 0, 1), ( 1, 1), (-1, 0), ( 0, 0), ( 1, 0), (-1,-1), ( 0, -1), ( 1, -1)] ] +c = [ Matrix(x) for x in product([-1,0,1], repeat=d) ] diff --git a/symbolic/D3Q19.py b/symbolic/D3Q19.py index 4e84908..e9e6eec 100644 --- a/symbolic/D3Q19.py +++ b/symbolic/D3Q19.py @@ -1,7 +1,7 @@ -from sympy import * +from sympy import Matrix, Rational, sqrt -q = 19 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), diff --git a/symbolic/D3Q27.py b/symbolic/D3Q27.py index d63aad8..6fb1f80 100644 --- a/symbolic/D3Q27.py +++ b/symbolic/D3Q27.py @@ -1,10 +1,7 @@ -from sympy import * +from sympy import Matrix +from itertools import product -q = 27 d = 3 +q = 27 -c = [ Matrix(x) for x in [ - (-1, 1, 1), ( 0, 1, 1), ( 1, 1, 1), (-1, 0, 1), ( 0, 0, 1), ( 1, 0, 1), (-1,-1, 1), ( 0, -1, 1), ( 1, -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), - (-1, 1,-1), ( 0, 1,-1), ( 1, 1,-1), (-1, 0,-1), ( 0, 0,-1), ( 1, 0,-1), (-1,-1,-1), ( 0, -1,-1), ( 1, -1,-1) -]] +c = [ Matrix(x) for x in product([-1,0,1], repeat=d) ] -- cgit v1.2.3