wassname/simpeg
1
0
Fork 0
mirror of https://github.com/wassname/simpeg.git synced 2026-07-26 13:37:21 +08:00
Code Issues Packages Projects Releases Wiki Activity

Bug fix for integrating over the model, should have a volume term!

Browse Source
This commit is contained in:
rowanc1
2013-11-26 14:35:08 -08:00
parent 5b4270560c
commit 1088c09864
1 changed files with 3 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files
SimPEG/regularization/Regularization.py
+3 -6
View File
@@ -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