@property
    def Wxx(self):
        if getattr(self, '_Wxx', None) is None:
            self._Wxx = Utils.sdiag((self.mesh.vol*self.alpha_xx)**0.5)*self.mesh.faceDivx*self.mesh.cellGradx

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

        return self._Wxx
Regularization matrix Wxx
dataset/ETHPy150Open simpeg/simpeg/SimPEG/Regularization.py/Tikhonov.Wxx