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

Running EB, HJ with PropMap

Browse Source
This commit is contained in:
Lindsey
2015-06-01 16:21:15 -07:00
parent 27c8da341d
commit b7a3b4e5e6
1 changed files with 11 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files
simpegEM/Base.py
+11 -8
View File
@@ -2,11 +2,14 @@ from SimPEG import Survey, Problem, Utils, Models, Maps, PropMaps, np, sp, Solve
from scipy.constants import mu_0
class EMPropMap(Maps.PropMap):
sigma = Maps.Property("Electrical Conductivity", defaultInvProp = True)
mui = Maps.Property("Inverse Magnetic Permeability", defaultVal = 1./mu_0)
sigma = Maps.Property("Electrical Conductivity", defaultInvProp = True, propertyLink=('rho',Maps.ReciprocalMap))
mu = Maps.Property("Inverse Magnetic Permeability", defaultVal = mu_0, propertyLink=('mui',Maps.ReciprocalMap))
# rho = Maps.Property("Electrical Resistivity")
# mu = Maps.Property("Inverse Magnetic Permeability", defaultVal = 1./mu_0)
rho = Maps.Property("Electrical Resistivity", propertyLink=('sigma', Maps.ReciprocalMap))
mui = Maps.Property("Inverse Magnetic Permeability", defaultVal = 1./mu_0, propertyLink=('mu', Maps.ReciprocalMap))
# Do some error checking: only 1 of sigma, rho can be InvProp similar story with mu and mui
# Also ensure that sigma and rho are reciprocals of one another ""
@@ -162,13 +165,13 @@ class BaseEMProblem(Problem.BaseProblem):
@property
def MeMu(self):
if getattr(self, '_MeMu', None) is None:
self._MeMu = self.mesh.getEdgeInnerProduct(self.mu)
self._MeMu = self.mesh.getEdgeInnerProduct(self.curModel.mu)
return self._MeMu
@property
def MeMuI(self):
if getattr(self, '_MeMuI', None) is None:
self._MeMuI = self.mesh.getEdgeInnerProduct(self.mu, invMat=True)
self._MeMuI = self.mesh.getEdgeInnerProduct(self.curModel.mu, invMat=True)
return self._MeMuI
# ----- Electrical Conductivity ----- #
@@ -193,7 +196,7 @@ class BaseEMProblem(Problem.BaseProblem):
@property
def MfSigmai(self):
if getattr(self, '_MfSigmai', None) is None:
self._MfSigmai = self.mesh.getFaceInnerProduct(self.sigmai)
self._MfSigmai = self.mesh.getFaceInnerProduct(self.curModel.rho)
return self._MfSigmai
# def dMfSigmai_dsigmai(self,u)
@@ -201,7 +204,7 @@ class BaseEMProblem(Problem.BaseProblem):
@property
def MfSigmaiI(self):
if getattr(self, '_MfSigmaiI', None) is None:
self._MfSigmaiI = self.mesh.getFaceInnerProduct(self.sigmai, invMat=True)
self._MfSigmaiI = self.mesh.getFaceInnerProduct(self.curModel.rho, invMat=True)
return self._MfSigmaiI
# def dMfSigmaiI(self,u)