From 06b34d8f7b2865189b34140d1629c059001ac932 Mon Sep 17 00:00:00 2001 From: rowanc1 Date: Fri, 2 May 2014 14:56:12 -0700 Subject: [PATCH] indexing bug in E field derivs --- simpegEM/TDEM/TDEM_b.py | 2 +- simpegEM/Tests/test_TDEM_combos.py | 6 ++++++ 2 files changed, 7 insertions(+), 1 deletion(-) diff --git a/simpegEM/TDEM/TDEM_b.py b/simpegEM/TDEM/TDEM_b.py index adf43f65..abcf1590 100644 --- a/simpegEM/TDEM/TDEM_b.py +++ b/simpegEM/TDEM/TDEM_b.py @@ -226,7 +226,7 @@ class ProblemTDEM_b(BaseTDEMProblem): y_b = mkvc(y_b) y_e = self.MeSigmaI*self.mesh.edgeCurl.T*self.MfMui*y_b if 'e' in p: - y_e += - self.MeSigmaI*p[:,'e',tInd] + y_e += - self.MeSigmaI*p[:,'e',tInd+1] return {'b':y_b, 'e':y_e} return self.adjoint(m, AhtRHS, AhtCalcFields) diff --git a/simpegEM/Tests/test_TDEM_combos.py b/simpegEM/Tests/test_TDEM_combos.py index 17526758..a64787a1 100644 --- a/simpegEM/Tests/test_TDEM_combos.py +++ b/simpegEM/Tests/test_TDEM_combos.py @@ -77,6 +77,12 @@ class TDEM_bDerivTests(unittest.TestCase): def test_Jvec_dbxdtbz(self): self.assertTrue(dotestJvec(*getProb(rxTypes='dbxdt,bz'))) def test_Adjoint_dbxdtbz(self): self.assertLess(*dotestAdjoint(*getProb(rxTypes='dbxdt,bz'))) + def test_Jvec_ey(self): self.assertTrue(dotestJvec(*getProb(rxTypes='ey'))) + def test_Adjoint_ey(self): self.assertLess(*dotestAdjoint(*getProb(rxTypes='ey'))) + + def test_Jvec_eybzdbxdt(self): self.assertTrue(dotestJvec(*getProb(rxTypes='ey,bz,dbxdt'))) + def test_Adjoint_eybzdbxdt(self): self.assertLess(*dotestAdjoint(*getProb(rxTypes='ey,bz,dbxdt'))) + if __name__ == '__main__': unittest.main()