jan/meanas
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

Only keep the real part of the energy

Browse Source
This commit is contained in:
Jan Petykiewicz 2025-01-14 22:00:08 -08:00
parent 659566750f
commit 829007c672
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
meanas/fdfd/waveguide_2d.py
View File

@ -425,7 +425,7 @@ def _normalized_fields(
Sz_tavg = numpy.real(Sz_a.sum() - Sz_b.sum()) * 0.5 # 0.5 since E, H are assumed to be peak (not RMS) amplitudes
assert Sz_tavg > 0, f'Found a mode propagating in the wrong direction! {Sz_tavg=}'
energy = epsilon * e.conj() * e
energy = numpy.real(epsilon * e.conj() * e)
norm_amplitude = 1 / numpy.sqrt(Sz_tavg)
norm_angle = -numpy.angle(e[energy.argmax()]) # Will randomly add a negative sign when mode is symmetric