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

remove static shift

Browse Source
This commit is contained in:
rowanc1
2014-07-12 12:21:52 -05:00
parent f3bd31b258
commit 1af495b910
1 changed files with 2 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files
simpegDC/BaseDC.py
+2 -5
View File
@@ -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):