From e56e5cbf05986bd5b8507bde3a5f3bce26511b73 Mon Sep 17 00:00:00 2001 From: Rowan Cockett Date: Fri, 18 Oct 2013 17:04:23 -0700 Subject: [PATCH] Removed inversion option. --- SimPEG/DiffOperators.py | 14 +++++--------- 1 file changed, 5 insertions(+), 9 deletions(-) diff --git a/SimPEG/DiffOperators.py b/SimPEG/DiffOperators.py index f26252d9..110fe9bd 100644 --- a/SimPEG/DiffOperators.py +++ b/SimPEG/DiffOperators.py @@ -286,12 +286,11 @@ class DiffOperators(object): _aveN2CCv = None aveN2CCv = property(**aveN2CCv()) - def getMass(self, materialProp=None, loc='e', inv=False): + def getMass(self, materialProp=None, loc='e'): """ Produces mass matricies. :param str loc: Average to location: 'e'-edges, 'f'-faces :param None,float,numpy.ndarray materialProp: property to be averaged (see below) - :param bool inv: True returns matrix inverse :rtype: scipy.sparse.csr.csr_matrix :return: M, the mass matrix @@ -322,18 +321,15 @@ class DiffOperators(object): diag = Av.T * (self.vol * mkvc(materialProp)) - if inv: - diag = 1/diag - return sdiag(diag) - def getEdgeMass(self, materialProp=None, inv=False): + def getEdgeMass(self, materialProp=None): """mass matrix for products of edge functions w'*M(materialProp)*e""" - return self.getMass(loc='e', materialProp=materialProp, inv=inv) + return self.getMass(loc='e', materialProp=materialProp) - def getFaceMass(self, materialProp=None, inv=False): + def getFaceMass(self, materialProp=None): """mass matrix for products of face functions w'*M(materialProp)*f""" - return self.getMass(loc='f', materialProp=materialProp, inv=inv) + return self.getMass(loc='f', materialProp=materialProp) def getFaceMassDeriv(self): Av = self.aveF2CC