From c488dabf9a113be5d277fa258c6e1a9a9d8e050e Mon Sep 17 00:00:00 2001 From: seogi_macbook Date: Wed, 4 May 2016 12:48:29 -0700 Subject: [PATCH] working on IP --- SimPEG/EM/Static/IP/ProblemIP.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/SimPEG/EM/Static/IP/ProblemIP.py b/SimPEG/EM/Static/IP/ProblemIP.py index dcf12d6b..4ff81990 100644 --- a/SimPEG/EM/Static/IP/ProblemIP.py +++ b/SimPEG/EM/Static/IP/ProblemIP.py @@ -154,7 +154,7 @@ class BaseIPProblem(BaseEMProblem): dMfRhoI_dI = -self.MfRhoI**2 dMf_drho = self.mesh.getFaceInnerProductDeriv(self.rho)(u) - drho_dlogrho = Utils.sdiag(self.rho) + drho_dlogrho = Utils.sdiag(self.rho)*self.curModel.etaDeriv return dMfRhoI_dI * ( dMf_drho * ( drho_dlogrho)) # TODO: This should take a vector @@ -162,7 +162,7 @@ class BaseIPProblem(BaseEMProblem): """ Derivative of MeSigma with respect to the model """ - dsigma_dlogsigma = Utils.sdiag(self.sigma) + dsigma_dlogsigma = Utils.sdiag(self.sigma)*self.curModel.etaDeriv return self.mesh.getEdgeInnerProductDeriv(self.sigma)(u) * dsigma_dlogsigma class Problem3D_CC(BaseIPProblem):