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

Simplify compute_source and fix propagation direction

Browse Source
This commit is contained in:
Jan Petykiewicz 2019-08-26 00:28:06 -07:00
parent 1860d754cd
commit d6a34b280e
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
meanas/fdfd/waveguide_mode.py
View File

@ -178,10 +178,10 @@ def compute_source(E: field_t,
s2[axis] = slice(slices[axis].start, slices[axis].stop) s2[axis] = slice(slices[axis].start, slices[axis].stop)
s2 = (src_order, *s2) s2 = (src_order, *s2)
rollby = 1 if polarity > 0 else 0 rollby = 1 if polarity < 0 else 0
exp_iphi = numpy.exp(-1j * polarity * wavenumber * dxes[1][axis][slices[axis]]) exp_iphi = numpy.exp(-1j * -rollby * wavenumber * dxes[1][axis][slices[axis]])
J[s2] = numpy.roll(functional.curl_h(dxes=dxes)(H.conj()), -rollby, axis=axis+1)[s2] * polarity * exp_iphi J[s2] = numpy.roll(functional.curl_h(dxes=dxes)(H), -rollby, axis=axis+1)[s2] * exp_iphi * -polarity
M[s2] = -numpy.roll(functional.curl_e(dxes=dxes)(E.conj()), rollby, axis=axis+1)[s2] M[s2] = numpy.roll(functional.curl_e(dxes=dxes)(E), rollby, axis=axis+1)[s2]
m2j = functional.m2j(omega, dxes, mu) m2j = functional.m2j(omega, dxes, mu)
Jm = m2j(M) Jm = m2j(M)