wassname/simpeg
1
0
Fork 0
mirror of https://github.com/wassname/simpeg.git synced 2026-06-29 22:00:42 +08:00
Code Issues Packages Projects Releases Wiki Activity

Updated getEdgeMass and getFaceMass to use the getMass method.

Browse Source
This commit is contained in:
Dave Marchant
2013-10-18 14:23:57 -07:00
parent 26f9b83b76
commit 1f770dfd63
1 changed files with 4 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files
SimPEG/DiffOperators.py
+4 -10
View File
@@ -329,19 +329,13 @@ class DiffOperators(object):
return sdiag(diag)
def getEdgeMass(self, materialProp=None):
def getEdgeMass(self, materialProp=None, inv=False):
"""mass matrix for products of edge functions w'*M(materialProp)*e"""
if(materialProp is None):
materialProp = np.ones(self.nC)
Av = self.edgeAve
return sdiag(Av.T * (self.vol * mkvc(materialProp)))
return self.getMass(loc='e', materialProp=materialProp, inv=inv)
def getFaceMass(self, materialProp=None):
def getFaceMass(self, materialProp=None, inv=False):
"""mass matrix for products of face functions w'*M(materialProp)*f"""
if(materialProp is None):
materialProp = np.ones(self.nC)
Av = self.faceAve
return sdiag(Av.T * (self.vol * mkvc(materialProp)))
return self.getMass(loc='f', materialProp=materialProp, inv=inv)
def getFaceMassDeriv(self):
Av = self.faceAve