fix missing 2π factor in conductivity for _get_epsmu (NanoComp#1692)

bencbartlett · Jul 23, 2021 · 01ed94f · 01ed94f
1 parent 8e4b642
commit 01ed94f
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/python/geom.py b/python/geom.py
@@ -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)