From 1088c09864a22d8be6a7c817b4d9648040a0e87a Mon Sep 17 00:00:00 2001 From: rowanc1 Date: Tue, 26 Nov 2013 14:35:08 -0800 Subject: [PATCH] Bug fix for integrating over the model, should have a volume term! --- SimPEG/regularization/Regularization.py | 9 +++------ 1 file changed, 3 insertions(+), 6 deletions(-) diff --git a/SimPEG/regularization/Regularization.py b/SimPEG/regularization/Regularization.py index 3d22a297..41938152 100644 --- a/SimPEG/regularization/Regularization.py +++ b/SimPEG/regularization/Regularization.py @@ -22,22 +22,19 @@ class Regularization(object): @property def Wx(self): if getattr(self, '_Wx', None) is None: - a = self.mesh.r(self.mesh.area,'F','Fx','V') - self._Wx = sdiag(a)*self.mesh.cellGradx + self._Wx = self.mesh.cellGradx*sdiag(self.mesh.vol) return self._Wx @property def Wy(self): if getattr(self, '_Wy', None) is None: - a = self.mesh.r(self.mesh.area,'F','Fy','V') - self._Wy = sdiag(a)*self.mesh.cellGrady + self._Wy = self.mesh.cellGrady*sdiag(self.mesh.vol) return self._Wy @property def Wz(self): if getattr(self, '_Wz', None) is None: - a = self.mesh.r(self.mesh.area,'F','Fz','V') - self._Wz = sdiag(a)*self.mesh.cellGradz + self._Wz = self.mesh.cellGradz*sdiag(self.mesh.vol) return self._Wz alpha_s = 1e-6