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

Problem3D_N is tested!

Browse Source
This commit is contained in:
seogi_macbook
2016-04-24 13:29:38 -07:00
parent 0bb001973c
commit eeee594f09
2 changed files with 4 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files
SimPEG/EM/Static/DC/FieldsDC.py
+1 -1
View File
@@ -99,7 +99,7 @@ class Fields_N(Fields):
return phiSolution
def _phiDeriv_u(self, src, v, adjoint = False):
return Identity()
return Identity()*v
def _phiDeriv_m(self, src, v, adjoint = False):
return Zero()
tests/em/static/test_DC.py
+3 -3
View File
@@ -104,7 +104,7 @@ class DCProblemTestsN(unittest.TestCase):
def test_misfit(self):
derChk = lambda m: [self.survey.dpred(m), lambda mx: self.p.Jvec(self.m0, mx)]
passed = Tests.checkDerivative(derChk, self.m0, plotIt=False)
passed = Tests.checkDerivative(derChk, self.m0, plotIt=False, num=3)
self.assertTrue(passed)
def test_adjoint(self):
@@ -114,13 +114,13 @@ class DCProblemTestsN(unittest.TestCase):
w = np.random.rand(self.survey.dobs.shape[0])
wtJv = w.dot(self.p.Jvec(self.m0, v))
vtJtw = v.dot(self.p.Jtvec(self.m0, w))
passed = np.abs(wtJv - vtJtw) < 1e-10
passed = np.abs(wtJv - vtJtw) < 1e-8
print 'Adjoint Test', np.abs(wtJv - vtJtw), passed
self.assertTrue(passed)
def test_dataObj(self):
derChk = lambda m: [self.dmis.eval(m), self.dmis.evalDeriv(m)]
passed = Tests.checkDerivative(derChk, self.m0, plotIt=False)
passed = Tests.checkDerivative(derChk, self.m0, plotIt=False, num=3)
self.assertTrue(passed)
if __name__ == '__main__':
unittest.main()