Brackets to make sure we are doing matrix vector products

This commit is contained in:
Lindsey Heagy
2015-07-05 21:57:19 -05:00
parent 4fa5e9426d
commit 7c63bfb6ac
+1 -1
View File
@@ -130,7 +130,7 @@ class BaseEMProblem(Problem.BaseProblem):
# TODO: This should take a vector
def MfRhoDeriv(self,u):
return self.mesh.getFaceInnerProductDeriv(self.curModel.rho)(u) * -Utils.sdiag(self.curModel.rho**2) * self.curModel.sigmaDeriv
return self.mesh.getFaceInnerProductDeriv(self.curModel.rho)(u) * (-Utils.sdiag(self.curModel.rho**2) * self.curModel.sigmaDeriv)
# self.curModel.rhoDeriv
@property