From 3846a6303b9eff50aa9100258376fd878c035f57 Mon Sep 17 00:00:00 2001 From: Lindsey Heagy Date: Wed, 8 Oct 2014 15:41:38 -0700 Subject: [PATCH] FEM -> FDEM in tests for analytics --- simpegEM/Analytics/FDEM.py | 2 +- simpegEM/Tests/test_FDEM_analytics.py | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/simpegEM/Analytics/FDEM.py b/simpegEM/Analytics/FDEM.py index a0fdf42e..3c6f7ca7 100644 --- a/simpegEM/Analytics/FDEM.py +++ b/simpegEM/Analytics/FDEM.py @@ -23,7 +23,7 @@ def hzAnalyticDipoleF(r, freq, sigma, secondary=True): """ r = np.abs(r) - k = np.sqrt(-1j*2.*np.pi*f*mu_0*sigma) + k = np.sqrt(-1j*2.*np.pi*freq*mu_0*sigma) m = 1 front = m / (2. * np.pi * (k**2) * (r**5) ) diff --git a/simpegEM/Tests/test_FDEM_analytics.py b/simpegEM/Tests/test_FDEM_analytics.py index 732b51ce..cae6d006 100644 --- a/simpegEM/Tests/test_FDEM_analytics.py +++ b/simpegEM/Tests/test_FDEM_analytics.py @@ -58,7 +58,7 @@ class FDEM_analyticTests(unittest.TestCase): P = self.mesh.getInterpolationMat(XYZ, 'Fz') - an = EM.Analytics.FEM.hzAnalyticDipoleF(x, self.Tx0.freq, self.sig) + an = EM.Analytics.FDEM.hzAnalyticDipoleF(x, self.Tx0.freq, self.sig) diff = np.log10(np.abs(P*np.imag(u[self.Tx0, 'b']) - mu_0*np.imag(an)))