wassname/simpeg
1
0
Fork 0
mirror of https://github.com/wassname/simpeg.git synced 2026-07-22 13:00:21 +08:00
Code Issues Packages Projects Releases Wiki Activity

indexing bug in E field derivs

Browse Source
This commit is contained in:
rowanc1
2014-05-02 14:56:12 -07:00
parent 464264b902
commit 06b34d8f7b
2 changed files with 7 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files
simpegEM/TDEM/TDEM_b.py
+1 -1
View File
@@ -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)
simpegEM/Tests/test_TDEM_combos.py
+6
View File
@@ -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()