maj

macro.py

@@ -5,6 +5,7 @@ from tqdm import tqdm
 import calculation as cl
 import optic as op
 import plot as plt
+import scipy.io
 import super_lattice_structure as st
 
 
@@ -42,4 +43,16 @@ def reflectivity_eam_vcsel_clad_heatmap():
 
             r[j, i] = op.reflection(n, d, 850e-9)
 
-    plt.plot_std_heatmap(l_eam_clad, l_vcsel_clad, r)
\ No newline at end of file
+    plt.plot_std_heatmap(l_eam_clad, l_vcsel_clad, r)
+
+
+def structure_for_matlab(vcsel_only=True, bypass_dbr=False, wavelength=850e-9):
+    sl = st.structure_eam_vcsel(vcsel_only=vcsel_only, bypass_dbr=bypass_dbr)
+    sl = op.algaas_super_lattice_refractive_index(sl, 0., wavelength)
+
+    lz = sl['thickness'].to_numpy(dtype=float)
+    n = sl['refractive_index'].apply(np.real).to_numpy(dtype=float)
+
+    scipy.io.savemat('eam_vcsel_n_lz.mat', dict(n=n, lz=lz))
+
+