From 55fb5b62ab3381327d4b04cf1c5ab095874d29cc Mon Sep 17 00:00:00 2001
From: Adrian Kummerlaender
Date: Tue, 2 Jul 2019 20:35:44 +0200
Subject: Determine lattice speed of sound

---
 symbolic/characteristics.py | 6 ++++++
 1 file changed, 6 insertions(+)

(limited to 'symbolic/characteristics.py')

diff --git a/symbolic/characteristics.py b/symbolic/characteristics.py
index ca3904e..b68afeb 100644
--- a/symbolic/characteristics.py
+++ b/symbolic/characteristics.py
@@ -19,3 +19,9 @@ def gauss_hermite(n):
 def weights(d, c):
     _, omegas = gauss_hermite(3)
     return list(map(lambda c_i: Mul(*[ omegas[1+c_i[iDim]] for iDim in range(0,d) ]) / pi**(d/2), c))
+
+# determine lattice speed of sound using directions and their weights
+def c_s(d, c, w):
+    speeds = set([ sqrt(sum([ w[i] * c_i[j]**2 for i, c_i in enumerate(c) ])) for j in range(0,d) ])
+    assert len(speeds) == 1 # verify isotropy
+    return speeds.pop()
-- 
cgit v1.2.3