From 1525096f9e5562abf1b3734d54929f849de02c3d Mon Sep 17 00:00:00 2001
From: Guillaume Demesy <guillaume.demesy@fresnel.fr>
Date: Tue, 27 Sep 2022 22:20:57 +0200
Subject: [PATCH] clean up
---
DiffractionGratings/grating3D_postplot_Mmatrix.py | 13 ++++++-------
1 file changed, 6 insertions(+), 7 deletions(-)
diff --git a/DiffractionGratings/grating3D_postplot_Mmatrix.py b/DiffractionGratings/grating3D_postplot_Mmatrix.py
index f209b63..ed8739d 100644
--- a/DiffractionGratings/grating3D_postplot_Mmatrix.py
+++ b/DiffractionGratings/grating3D_postplot_Mmatrix.py
@@ -86,19 +86,19 @@ TOT2_sin = T2_sin+R2_sin+Qscat_sin
M[1-1,1-1,:,:] = 0.5*(np.abs(rpp)**2+np.abs(rsp)**2+np.abs(rps)**2+np.abs(rss)**2)
M[1-1,2-1,:,:] = 0.5*(np.abs(rpp)**2+np.abs(rsp)**2-np.abs(rps)**2-np.abs(rss)**2)
-M[1-1,3-1,:,:] = np.real(rpp*np.conj(rps)+rsp*np.conj(rss)) #*-1
+M[1-1,3-1,:,:] = np.real(rpp*np.conj(rps)+rsp*np.conj(rss))
M[1-1,4-1,:,:] = np.imag(rpp*np.conj(rps)+rsp*np.conj(rss))
M[2-1,1-1,:,:] = 0.5*(np.abs(rpp)**2-np.abs(rsp)**2+np.abs(rps)**2-np.abs(rss)**2)
M[2-1,2-1,:,:] = 0.5*(np.abs(rpp)**2-np.abs(rsp)**2-np.abs(rps)**2+np.abs(rss)**2)
-M[2-1,3-1,:,:] = np.real(rpp*np.conj(rps)-rsp*np.conj(rss)) #*-1
+M[2-1,3-1,:,:] = np.real(rpp*np.conj(rps)-rsp*np.conj(rss))
M[2-1,4-1,:,:] = np.imag(rpp*np.conj(rps)-rsp*np.conj(rss))
-M[3-1,1-1,:,:] = np.real(rpp*np.conj(rsp)+rps*np.conj(rss)) #*-1
-M[3-1,2-1,:,:] = np.real(rpp*np.conj(rsp)-rps*np.conj(rss)) #*-1
+M[3-1,1-1,:,:] = np.real(rpp*np.conj(rsp)+rps*np.conj(rss))
+M[3-1,2-1,:,:] = np.real(rpp*np.conj(rsp)-rps*np.conj(rss))
M[3-1,3-1,:,:] = np.real(rpp*np.conj(rss)+rps*np.conj(rsp))
-M[3-1,4-1,:,:] = np.imag(rpp*np.conj(rss)-rps*np.conj(rsp)) #*-1
+M[3-1,4-1,:,:] = np.imag(rpp*np.conj(rss)-rps*np.conj(rsp))
M[4-1,1-1,:,:] =-np.imag(rpp*np.conj(rsp)+rps*np.conj(rss))
M[4-1,2-1,:,:] =-np.imag(rpp*np.conj(rsp)-rps*np.conj(rss))
-M[4-1,3-1,:,:] =-np.imag(rpp*np.conj(rss)+rps*np.conj(rsp)) #*-1
+M[4-1,3-1,:,:] =-np.imag(rpp*np.conj(rss)+rps*np.conj(rsp))
M[4-1,4-1,:,:] = np.real(rpp*np.conj(rss)-rps*np.conj(rsp))
@@ -147,7 +147,6 @@ for i in range(4):
norm = mpl.colors.Normalize(vmin=-1,vmax=1)
sm = mpl.cm.ScalarMappable(cmap=plt.cm.bwr, norm=norm)
sm.set_array([])
- # plt.colorbar(sm, ticks=np.linspace(-1,1,3))
ax.set_rorigin(which_orig)
plt.subplots_adjust(top=0.92, bottom=0.08, left=0.10, right=0.95, hspace=0.25,wspace=0.35)
plt.savefig('Mmatrix.jpg')
--
GitLab