From 001bcbce271f1fc7840ffcf2508ae4311c7d61ac Mon Sep 17 00:00:00 2001 From: Lindsey Heagy Date: Sat, 4 Jul 2015 13:21:56 -0700 Subject: [PATCH] fixed Tikhonov deriv and added non-zero mref to testing --- SimPEG/Regularization.py | 4 ++-- SimPEG/Tests/test_regularization.py | 2 +- 2 files changed, 3 insertions(+), 3 deletions(-) diff --git a/SimPEG/Regularization.py b/SimPEG/Regularization.py index e7fc4895..f181a8a3 100644 --- a/SimPEG/Regularization.py +++ b/SimPEG/Regularization.py @@ -303,8 +303,8 @@ class Tikhonov(BaseRegularization): if self.smoothModel == True: mD1 = self.mapping.deriv(m) mD2 = self.mapping.deriv(self.mref) - r1 = self.W * ( self.mapping * (m) ) - r2 = self.Ws * ( self.mapping * (self.mref) ) + r1 = self.W * ( self.mapping * (m - self.mref) ) + r2 = self.Ws * ( self.mapping * (m - self.mref) ) out1 = mD1.T * ( self.W.T * r1 ) out2 = mD2.T * ( self.Ws.T * r2 ) out = out1-out2 diff --git a/SimPEG/Tests/test_regularization.py b/SimPEG/Tests/test_regularization.py index 2197d27b..48846d2f 100644 --- a/SimPEG/Tests/test_regularization.py +++ b/SimPEG/Tests/test_regularization.py @@ -22,7 +22,7 @@ class RegularizationTests(unittest.TestCase): mapping = r.mapPair(self.mesh2) reg = r(self.mesh2, mapping=mapping) m = np.random.rand(mapping.nP) - reg.mref = m[:]*0 + reg.mref = m[:]*np.mean(m) passed = checkDerivative(lambda m : [reg.eval(m), reg.evalDeriv(m)], m, plotIt=False) self.assertTrue(passed)