From c1d53f83d061b565100eb59759d7a77cd65aafcb Mon Sep 17 00:00:00 2001 From: Toon Verstraelen Date: Fri, 3 Jul 2026 21:28:32 +0200 Subject: [PATCH] Generalize the ptens script for the assymetric tensor case --- tools/ptens.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/tools/ptens.py b/tools/ptens.py index 1249292..b845cdb 100755 --- a/tools/ptens.py +++ b/tools/ptens.py @@ -3,9 +3,9 @@ import sympy as sp -pxx, pxy, pzx, pyy, pyz, pzz = sp.symbols("pxx pxy pzx pyy pyz pzz") +pxx, pxy, pxz, pyx, pyy, pyz, pzx, pzy, pzz = sp.symbols("pxx pxy pxz pyx pyy pyz pzx pzy pzz") -ptens = sp.Matrix([[pxx, pxy, pzx], [pxy, pyy, pyz], [pzx, pyz, pzz]]) +ptens = sp.Matrix([[pxx, pxy, pxz], [pyx, pyy, pyz], [pzx, pzy, pzz]]) rtens = sp.Matrix( [[1, 0, 0], [0, 1 / sp.sqrt(2), -1 / sp.sqrt(2)], [0, 1 / sp.sqrt(2), 1 / sp.sqrt(2)]] )