wassname/simpeg
1
0
Fork 0
mirror of https://github.com/wassname/simpeg.git synced 2026-08-19 12:40:13 +08:00
Code Issues Packages Projects Releases Wiki Activity

work on IP stuff

Browse Source
This commit is contained in:
seogi_macbook
2015-09-15 15:40:26 -07:00
parent e0d9c27d87
commit dea7eba5e0
6 changed files with 592 additions and 589 deletions
Download Patch File Download Diff File Expand all files Collapse all files
simpegDCIP/BaseDC.py
+6 -7
View File
@@ -245,15 +245,15 @@ class ProblemDC_CC(Problem.BaseProblem):
dCdm_x_v[:, i] = dAdsig * dsigdm_x_v
# Take derivative of $C(m,u)$ w.r.t. $u$
dCdu = self.A
dA_du = self.A
# Solve for $\deriv{u}{m}$
# dCdu_inv = self.Solver(dCdu, **self.solverOpts)
if self.Ainv is None:
self.Ainv = self.Solver(dCdu, **self.solverOpts)
self.Ainv = self.Solver(dA_du, **self.solverOpts)
P = self.survey.getP(self.mesh)
J_x_v = - P * mkvc( self.Ainv * dCdm_x_v )
return J_x_v
Jv = - P * mkvc( self.Ainv * dCdm_x_v )
return Jv
def Jtvec(self, m, v, u=None):
@@ -271,12 +271,11 @@ class ProblemDC_CC(Problem.BaseProblem):
D = self.mesh.faceDiv
G = self.mesh.cellGrad
A = self.A
dA_du = self.A
mT_dm = self.mapping.deriv(m)
# We probably always need this due to the linesearch .. (?)
dCdu = A.T
self.Ainv = self.Solver(dCdu, **self.solverOpts)
self.Ainv = self.Solver(dA_du.T, **self.solverOpts)
# if self.Ainv is None:
# self.Ainv = self.Solver(dCdu, **self.solverOpts)