From 0103201fb4d0fdf9ae244771e91caaf46bbd5bba Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Marek=20Ne=C4=8Dada?= Date: Mon, 8 May 2017 19:27:08 +0300 Subject: [PATCH] Fix missing translation coefficients that probably caused mysterious asymmetry. Former-commit-id: c2edcd35a3d7792013fa376306e21d76ee8beddb --- misc/dispersion-SVD.py | 2 ++ misc/dispersion2D-SVD.py | 5 ++++- 2 files changed, 6 insertions(+), 1 deletion(-) diff --git a/misc/dispersion-SVD.py b/misc/dispersion-SVD.py index 75956d5..b3cc64b 100755 --- a/misc/dispersion-SVD.py +++ b/misc/dispersion-SVD.py @@ -395,6 +395,8 @@ for trfile in os.scandir(translations_dir): leftmatrix[1,0,:,1,1,:] = leftmatrix[0,0,:,0,1,:] # d2d, M2E leftmatrix[0,0,:,1,0,:] = np.tensordot(a_d2u, phases_d2u,axes=(0,-1)) #d2u,E2E leftmatrix[0,1,:,1,1,:] = leftmatrix[0,0,:,1,0,:] #d2u, M2M + leftmatrix[1,0,:,0,0,:] = np.tensordot(a_u2d, phases_u2d,axes=(0,-1)) #u2d,E2E + leftmatrix[1,1,:,0,1,:] = leftmatrix[1,0,:,0,0,:] #u2d, M2M leftmatrix[0,0,:,1,1,:] = np.tensordot(b_d2u, phases_d2u,axes=(0,-1)) #d2u,M2E leftmatrix[0,1,:,1,0,:] = leftmatrix[0,0,:,1,1,:] #d2u, E2M leftmatrix[1,0,:,0,1,:] = np.tensordot(b_u2d, phases_u2d,axes=(0,-1)) #u2d,M2E diff --git a/misc/dispersion2D-SVD.py b/misc/dispersion2D-SVD.py index ad615b0..b029cbc 100755 --- a/misc/dispersion2D-SVD.py +++ b/misc/dispersion2D-SVD.py @@ -359,7 +359,8 @@ for one in (1,): phases_u2d *= u2d_envelope phases_d2u *= d2u_envelope leftmatrix = np.zeros((2,2,nelem, 2,2,nelem), dtype=complex) - + # 0: u,E<--u,E + # 1: d,M<--d,M leftmatrix[0,0,:,0,0,:] = np.tensordot(a_self,phases_self, axes=(0,-1)) # u2u, E2E leftmatrix[1,0,:,1,0,:] = leftmatrix[0,0,:,0,0,:] # d2d, E2E leftmatrix[0,1,:,0,1,:] = leftmatrix[0,0,:,0,0,:] # u2u, M2M @@ -370,6 +371,8 @@ for one in (1,): leftmatrix[1,0,:,1,1,:] = leftmatrix[0,0,:,0,1,:] # d2d, M2E leftmatrix[0,0,:,1,0,:] = np.tensordot(a_d2u, phases_d2u,axes=(0,-1)) #d2u,E2E leftmatrix[0,1,:,1,1,:] = leftmatrix[0,0,:,1,0,:] #d2u, M2M + leftmatrix[1,0,:,0,0,:] = np.tensordot(a_u2d, phases_u2d,axes=(0,-1)) #u2d,E2E + leftmatrix[1,1,:,0,1,:] = leftmatrix[1,0,:,0,0,:] #u2d, M2M leftmatrix[0,0,:,1,1,:] = np.tensordot(b_d2u, phases_d2u,axes=(0,-1)) #d2u,M2E leftmatrix[0,1,:,1,0,:] = leftmatrix[0,0,:,1,1,:] #d2u, E2M leftmatrix[1,0,:,0,1,:] = np.tensordot(b_u2d, phases_u2d,axes=(0,-1)) #u2d,M2E