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

Working Jvec:

Browse Source 
Getting closer to understand how modular EM code is working ...
This commit is contained in:
seogi_macbook
2016-04-23 10:40:47 -07:00
parent 73001abfc5
commit 8cac166fba
2 changed files with 3 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files
SimPEG/EM/Static/DC/FieldsDC.py
+2 -1
View File
@@ -12,6 +12,7 @@ class Fields(SimPEG.Problem.Fields):
if adjoint:
return self._phiDeriv_u(src, v, adjoint=adjoint), self._phiDeriv_m(src, v, adjoint=adjoint)
return np.array(self._phiDeriv_u(src, du_dm_v, adjoint) + self._phiDeriv_m(src, v, adjoint), dtype = float)
def _eDeriv(self, src, du_dm_v, v, adjoint=False):
@@ -59,7 +60,7 @@ class Fields_CC(Fields):
return phiSolution
def _phiDeriv_u(self, src, v, adjoint = False):
return Identity()
return Identity()*v
def _phiDeriv_m(self, src, v, adjoint = False):
return Zero()
SimPEG/EM/Static/DC/ProblemDC.py
+1 -2
View File
@@ -1,4 +1,4 @@
from SimPEG import Problem
from SimPEG import Problem, Utils
from SimPEG.EM.Base import BaseEMProblem
from SurveyDC import Survey
from FieldsDC import Fields, Fields_CC, Fields_N
@@ -38,7 +38,6 @@ class BaseDCProblem(BaseEMProblem):
u_src = f[src, self._solutionType] # solution vector
dA_dm_v = self.getADeriv(u_src, v)
dRHS_dm_v = self.getRHSDeriv(src, v)
print type(dA_dm_v + dRHS_dm_v), (dA_dm_v + dRHS_dm_v).shape
du_dm_v = Ainv * ( - dA_dm_v + dRHS_dm_v )
for rx in src.rxList: