Merge branch 'diffOperators' of https://bitbucket.org/rcockett/simpeg into diffOperators

2026-07-13 17:45:30 +08:00 · 2013-07-19 14:51:34 -07:00
parent 906d7d210f ffb78fea10
commit 3eec485adc
1 changed files with 20 additions and 0 deletions
@@ -239,3 +239,23 @@ class DiffOperators(object):
        return locals()
    _edgeAve = None
    edgeAve = property(**edgeAve())
+
+
+
+def getEdgeMassMatrix(h,sigma):
+    # mass matix for products of edge functions w'*M(sigma)*e
+         
+    Av    = getEdgeToCellAverge(h)
+    v     = getVol(h)
+    sigma = mkvc(sigma)
+    
+    return sdiag(Av.T*(v*sigma))
+    
+def getFaceMassMatrix(h,sigma):
+    # mass matix for products of edge functions w'*M(sigma)*e
+         
+    Av    = getFaceToCellAverge(h)
+    v     = getVol(h)
+    sigma = mkvc(sigma)
+    
+    return sdiag(Av.T*(v*sigma))