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

fixed observed far-field component

Browse Source
This commit is contained in:
edoardo
2026-05-13 10:43:51 +02:00
parent 3cdb9e5f68
commit 32503e330c
2 changed files with 7 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files
mode_converter.ipynb
+4 -3
View File
@@ -197,7 +197,7 @@
"]\n", "]\n",
"#mpa.Near2FarFields return a numpy array with shape (num_of_points,nfreq,6) where the \n", "#mpa.Near2FarFields return a numpy array with shape (num_of_points,nfreq,6) where the \n",
"#third axis are the field components Ex,Ey,Ez,Hx,Hy,Hz\n", "#third axis are the field components Ex,Ey,Ez,Hx,Hy,Hz\n",
"FarFields = mpa.Near2FarFields(sim,NearRegions,far_z)\n", "FarFields = mpa.Near2FarFields(sim, NearRegions, far_z)\n",
"ob_list = [FarFields]\n", "ob_list = [FarFields]\n",
"\n", "\n",
"#this the objective to function\n", "#this the objective to function\n",
@@ -205,7 +205,7 @@
"def J1(FF):\n", "def J1(FF):\n",
" #in this case we want to maximize field intensity\n", " #in this case we want to maximize field intensity\n",
" # at the position of the focal spot \n", " # at the position of the focal spot \n",
" obj_to_return = npa.mean(npa.abs(FF[0,:,1]) ** 2)\n", " obj_to_return = npa.mean(npa.abs(FF[0,:,0]) ** 2)\n",
" return obj_to_return\n", " return obj_to_return\n",
"#returns the first point of all frequencies of the Ey component" "#returns the first point of all frequencies of the Ey component"
] ]
@@ -258,7 +258,8 @@
")\n", ")\n",
"\n", "\n",
"######CHECK PLOT######\n", "######CHECK PLOT######\n",
"figure = opt.plot2D(True, output_plane=output_plane)\n" "figure = opt.plot2D(True, output_plane=output_plane)\n",
"plt.savefig(\"results/simulation_setup.pdf\")"
] ]
}, },
{ {
mode_converter.py
+3 -3
View File
@@ -142,7 +142,7 @@ NearRegions = [
] ]
#mpa.Near2FarFields return a numpy array with shape (num_of_points,nfreq,6) where the #mpa.Near2FarFields return a numpy array with shape (num_of_points,nfreq,6) where the
#third axis are the field components Ex,Ey,Ez,Hx,Hy,Hz #third axis are the field components Ex,Ey,Ez,Hx,Hy,Hz
FarFields = mpa.Near2FarFields(sim,NearRegions,far_z) FarFields = mpa.Near2FarFields(sim, NearRegions, far_z)
ob_list = [FarFields] ob_list = [FarFields]
#this the objective to function #this the objective to function
@@ -150,7 +150,7 @@ ob_list = [FarFields]
def J1(FF): def J1(FF):
#in this case we want to maximize field intensity #in this case we want to maximize field intensity
# at the position of the focal spot # at the position of the focal spot
obj_to_return = npa.mean(npa.abs(FF[0,:,1]) ** 2) obj_to_return = npa.mean(npa.abs(FF[0,:,0]) ** 2)
return obj_to_return return obj_to_return
#returns the first point of all frequencies of the Ey component #returns the first point of all frequencies of the Ey component
@@ -190,7 +190,7 @@ opt = mpa.OptimizationProblem(
######CHECK PLOT###### ######CHECK PLOT######
figure = opt.plot2D(True, output_plane=output_plane) figure = opt.plot2D(True, output_plane=output_plane)
plt.savefig("results/simulation_setup.pdf")
# %% # %%
def mapping(x,eta,beta): def mapping(x,eta,beta):