From 54a45dbec3e046b3e3c333aa1fde68244d0d2a51 Mon Sep 17 00:00:00 2001 From: Adrian Kummerlaender Date: Sun, 16 Jun 2019 13:08:02 +0200 Subject: Add D3Q27 descriptor --- symbolic/D3Q27.py | 18 ++++++++++++++++++ 1 file changed, 18 insertions(+) create mode 100644 symbolic/D3Q27.py (limited to 'symbolic') diff --git a/symbolic/D3Q27.py b/symbolic/D3Q27.py new file mode 100644 index 0000000..6a81de5 --- /dev/null +++ b/symbolic/D3Q27.py @@ -0,0 +1,18 @@ +from sympy import * + +q = 27 +d = 3 + +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) +]] + +w = [Rational(*x) for x in [ + (1, 216), (1,54), (1,216), (1,54), (2,27), (1,54), (1,216), (1,54), (1,216), + (1, 54), (2,27), (1, 54), (2,27), (8,27), (2,27), (1, 54), (2,27), (1, 54), + (1, 216), (1,54), (1,216), (1,54), (2,27), (1,54), (1,216), (1,54), (1,216) +]] + +c_s = sqrt(Rational(1,3)) -- cgit v1.2.3