edoardo/modeconverter
1
0
Fork 0
mirror of https://github.com/edoverg/modeconverter.git synced 2026-10-14 23:06:39 +00:00
Code Issues Packages Projects Releases Wiki Activity

parameter change for testing

Browse Source
This commit is contained in:
edoardo
2026-06-18 11:48:30 +02:00
parent 24f3f24e5d
commit e275ea610b
1 changed files with 9 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files
metaprop_fibermode_doublet.py
+9 -8
View File
@@ -22,12 +22,12 @@ opt_max_eval = 250
size_x = 512 * wavelength size_x = 512 * wavelength
size_y = 512 * wavelength size_y = 512 * wavelength
res_x = 3 / 1e-6 #number of pixels per unit-length res_x = 6 / 1e-6 #number of pixels per unit-length
res_y = res_x res_y = res_x
ds = 1 / res_x ds = 1 / res_x
Nx = np.pow(2,int(np.log(int(round(res_x * size_x))) / np.log(2))) Nx = int(np.pow(2,int(np.log(int(round(res_x * size_x))) / np.log(2))))
Ny = np.pow(2,int(np.log(int(round(res_x * size_x))) / np.log(2))) Ny = int(np.pow(2,int(np.log(int(round(res_x * size_x))) / np.log(2))))
S = Nx * Ny S = Nx * Ny
###################### ######################
@@ -126,11 +126,12 @@ def get_fiber_mode_pattern():
for i in range(Nx): for i in range(Nx):
for j in range(Ny): for j in range(Ny):
r = np.sqrt(X[i,j]**2 + Y[i,j]**2) r = np.sqrt(X[i,j]**2 + Y[i,j]**2)
phi = np.arctan2(Y[i,j], X[i,j]) if r < 5*a:
if r < a:#inside core phi = np.arctan2(Y[i,j], X[i,j])
E_field[i,j] = jv(n_mode_sel,chi_co * r) / jv(n_mode_sel, chi_co * a) * np.cos(n_mode_sel * phi) if r < a:#inside core
else:#outside core E_field[i,j] = jv(n_mode_sel,chi_co * r) / jv(n_mode_sel, chi_co * a) * np.cos(n_mode_sel * phi)
E_field[i,j] = kv(n_mode_sel, chi_cl * r) / kv(n_mode_sel, chi_cl * a) * np.cos(n_mode_sel * phi) else:#outside core
E_field[i,j] = kv(n_mode_sel, chi_cl * r) / kv(n_mode_sel, chi_cl * a) * np.cos(n_mode_sel * phi)
return E_field return E_field