From b4ab60c2601b6cb2dd64a65d58863f66f265f458 Mon Sep 17 00:00:00 2001 From: D Fournier Date: Thu, 5 May 2016 11:55:56 -0700 Subject: [PATCH] Add model mapping to sparse regularization --- SimPEG/Regularization.py | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/SimPEG/Regularization.py b/SimPEG/Regularization.py index ed1039ec..1b0d180f 100644 --- a/SimPEG/Regularization.py +++ b/SimPEG/Regularization.py @@ -663,7 +663,7 @@ class Sparse(Simple): self.Rs = Utils.speye(self.regmesh.nC) else: - f_m = self.curModel - self.reg.mref + f_m = self.mapping * (self.curModel - self.reg.mref) self.rs = self.R(f_m , self.eps_p, self.norms[0]) #print "Min rs: " + str(np.max(self.rs)) + "Max rs: " + str(np.min(self.rs)) self.Rs = Utils.sdiag( self.rs ) @@ -679,7 +679,7 @@ class Sparse(Simple): self.Rx = Utils.speye(self.regmesh.cellDiffxStencil.shape[0]) else: - f_m = self.regmesh.cellDiffxStencil * self.curModel + f_m = self.regmesh.cellDiffxStencil * (self.mapping * self.curModel) self.rx = self.R( f_m , self.eps_q, self.norms[1]) self.Rx = Utils.sdiag( self.rx ) @@ -693,7 +693,7 @@ class Sparse(Simple): self.Ry = Utils.speye(self.regmesh.cellDiffyStencil.shape[0]) else: - f_m = self.regmesh.cellDiffyStencil * self.curModel + f_m = self.regmesh.cellDiffyStencil * (self.mapping * self.curModel) self.ry = self.R( f_m , self.eps_q, self.norms[2]) self.Ry = Utils.sdiag( self.ry ) @@ -707,7 +707,7 @@ class Sparse(Simple): self.Rz = Utils.speye(self.regmesh.cellDiffzStencil.shape[0]) else: - f_m = self.regmesh.cellDiffzStencil * self.curModel + f_m = self.regmesh.cellDiffzStencil * (self.mapping * self.curModel) self.rz = self.R( f_m , self.eps_q, self.norms[3]) self.Rz = Utils.sdiag( self.rz )