wassname/simpeg
1
0
Fork 0
mirror of https://github.com/wassname/simpeg.git synced 2026-07-20 12:40:44 +08:00
Code Issues Packages Projects Releases Wiki Activity

Merge pull request #305 from simpeg/feat/sparse-regularization

Browse Source 
Feat/sparse regularization
This commit is contained in:
Lindsey
2016-05-04 22:30:06 -07:00
parent 5e1de61a71 4e296c4cd5
commit 79e1378009
6 changed files with 287 additions and 76 deletions
Download Patch File Download Diff File Expand all files Collapse all files
SimPEG/Regularization.py
+6 -6
View File
@@ -731,14 +731,14 @@ class Sparse(Simple):
@property
def W(self):
"""Full regularization matrix W"""
#if getattr(self, '_W', None) is None:
wlist = (self.Wsmall, self.Wsmooth)
#self._W = sp.vstack(wlist)
return sp.vstack(wlist)
if getattr(self, '_W', None) is None:
wlist = (self.Wsmall, self.Wsmooth)
self._W = sp.vstack(wlist)
return self._W
def R(self, f_m , eps, exponent):
eta = (eps**(1-exponent/2.))**0.5
r = eta / (f_m**2.+ eps**2.)**((1-exponent/2.)/2.)
eta = (eps**(1.-exponent/2.))**0.5
r = eta / (f_m**2.+ eps**2.)**((1.-exponent/2.)/2.)
return r