From 1af495b910294baafb11695adb6da7c202de0348 Mon Sep 17 00:00:00 2001
From: rowanc1 <rowan@visiblegeology.com>
Date: Sat, 12 Jul 2014 12:21:52 -0500
Subject: [PATCH] remove static shift

---
 simpegDC/BaseDC.py | 7 ++-----
 1 file changed, 2 insertions(+), 5 deletions(-)

diff --git a/simpegDC/BaseDC.py b/simpegDC/BaseDC.py
index 929f226c..39037f00 100644
--- a/simpegDC/BaseDC.py
+++ b/simpegDC/BaseDC.py
@@ -104,7 +104,7 @@ class ProblemDC(Problem.BaseProblem):
         if getattr(self, '_Msig', None) is None:
             sigma = self.curModel.transform
             Av = self.mesh.aveF2CC
-            self._Msig = Utils.sdInv(Utils.sdiag(self.mesh.dim * Av.T * (1/sigma)))
+            self._Msig = Utils.sdiag(1/(self.mesh.dim * Av.T * (1/sigma)))
         return self._Msig
 
     @property
@@ -145,9 +145,6 @@ class ProblemDC(Problem.BaseProblem):
         Ainv = self.Solver(A)
         Q    = self.survey.getRhs(self.mesh)
         Phi  = Ainv * Q
-        for ii in range(Phi.shape[1]):
-            # Remove the static shift for each phi column.
-            Phi[:,ii] -= Phi[-1, ii]
         return Phi
 
     def Jvec(self, m, v, u=None):
@@ -201,7 +198,7 @@ class ProblemDC(Problem.BaseProblem):
         dCdu_inv = self.Solver(dCdu, **self.solverOpts)
         P        = self.survey.getP(self.mesh)
         J_x_v    = - P * mkvc( dCdu_inv * dCdm_x_v )
-        return J_x_v # Make $\mathbf{Jv}$ a vector.
+        return J_x_v
 
     def Jtvec(self, m, v, u=None):