diff --git a/lyceanem/models/frequency_domain.py b/lyceanem/models/frequency_domain.py index e50d5c0..b3dee7b 100644 --- a/lyceanem/models/frequency_domain.py +++ b/lyceanem/models/frequency_domain.py @@ -757,10 +757,12 @@ def calculate_scattering( scatter_map = EM.EMGPUFreqDomain( num_sources, num_sinks, full_index, point_informationv2, wavelength ) + Ex = np.dot(np.ones((num_sources)), scatter_map[:, :, 0]) Ey = np.dot(np.ones((num_sources)), scatter_map[:, :, 1]) Ez = np.dot(np.ones((num_sources)), scatter_map[:, :, 2]) + # convert to etheta,ephi else: @@ -822,3 +824,4 @@ def calculate_scattering( return Ex, Ey, Ez +