@property
    def Wyy(self):
        if getattr(self, '_Wyy', None) is None:
            self._Wyy = Utils.sdiag((self.mesh.vol*self.alpha_yy)**0.5)*self.mesh.faceDivy*self.mesh.cellGrady

            if self.indActive is not None:
                Pac = Utils.speye(self.mesh.nC)[:,self.indActive]
                self._Wyy = Pac.T*self._Wyy*Pac

        return self._Wyy
Regularization matrix Wyy
dataset/ETHPy150Open simpeg/simpeg/SimPEG/Regularization.py/Tikhonov.Wyy