From 8cac166fba57d361887f85fdc566cb557468f634 Mon Sep 17 00:00:00 2001 From: seogi_macbook Date: Sat, 23 Apr 2016 10:40:47 -0700 Subject: [PATCH] Working Jvec: Getting closer to understand how modular EM code is working ... --- SimPEG/EM/Static/DC/FieldsDC.py | 3 ++- SimPEG/EM/Static/DC/ProblemDC.py | 3 +-- 2 files changed, 3 insertions(+), 3 deletions(-) diff --git a/SimPEG/EM/Static/DC/FieldsDC.py b/SimPEG/EM/Static/DC/FieldsDC.py index 46a779f5..07b81806 100644 --- a/SimPEG/EM/Static/DC/FieldsDC.py +++ b/SimPEG/EM/Static/DC/FieldsDC.py @@ -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() diff --git a/SimPEG/EM/Static/DC/ProblemDC.py b/SimPEG/EM/Static/DC/ProblemDC.py index 25ac7fcf..ff3fff0f 100644 --- a/SimPEG/EM/Static/DC/ProblemDC.py +++ b/SimPEG/EM/Static/DC/ProblemDC.py @@ -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: