Commit 2421ec96 authored by Lucas Laplanche's avatar Lucas Laplanche
Browse files

smm is still incorrect ; tmm is doomed

parent 227a2e0b
......@@ -45,7 +45,7 @@ def zandberg_electromagnetic_amplitude(electric_field=0., wavelength=1e-6):
sl = st.structure_zandbergen()
sl = op.algaas_super_lattice_refractive_index(sl, electric_field, wavelength)
sl = pt.cut_in_equal_layers_thickness(sl, 1e-9)
em = op.em_amplitude_smm(sl, wavelength)
em = op.em_amplitude_tmm(sl, wavelength)
sl['electromagnetic_amplitude'] = em.tolist()
plt.plot_refra_em(sl)
......
......@@ -290,7 +290,7 @@ def transmission_matrix(super_lattice, wavelength=850e-9):
def interface_intensity(super_lattice, m, idx):
n = super_lattice['refractive_index'].to_numpy(dtype=np.complex128)
n = np.real(super_lattice['refractive_index'].to_numpy(dtype=np.complex128))
rap = n[idx] / n[idx +1]
......@@ -313,11 +313,11 @@ def interface_intensity(super_lattice, m, idx):
def core_intensity(super_lattice, wavelength, m, idx):
n = super_lattice['refractive_index'].to_numpy(dtype=np.complex128)
n = np.real(super_lattice['refractive_index'].to_numpy(dtype=np.complex128))
l = super_lattice['thickness'].to_numpy(dtype=np.complex128)
dp = np.array([[np.exp(2j * np.pi * n[idx + 1] * l[idx + 1] / wavelength), 0.],
[0., np.exp(2j * np.pi * n[idx + 1] * l[idx + 1] / wavelength)]
[0., np.exp(-2j * np.pi * n[idx + 1] * l[idx + 1] / wavelength)]
])
m = np.matmul(dp, m)
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment