Skip to content

Commit

Permalink
fix missing 2π factor in conductivity for _get_epsmu (#1692)
Browse files Browse the repository at this point in the history
  • Loading branch information
stevengj authored Jul 23, 2021
1 parent d662b0c commit 83a10c5
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion python/geom.py
Original file line number Diff line number Diff line change
Expand Up @@ -517,7 +517,7 @@ def _get_epsmu(self, diag, offdiag, susceptibilities, conductivity_diag, conduct
# Account for conductivity term (only multiply if nonzero to avoid unnecessary complex numbers)
conductivity = np.expand_dims(Matrix(diag=conductivity_diag,offdiag=conductivity_offdiag),axis=0)
if np.count_nonzero(conductivity) > 0:
epsmu = (1 + 1j/freqs * conductivity) * epsmu
epsmu = (1 + 1j/(2*np.pi*freqs) * conductivity) * epsmu

# Convert list matrix to 3D numpy array size [freqs,3,3]
return np.squeeze(epsmu)
Expand Down

0 comments on commit 83a10c5

Please sign in to comment.