@property
    def Wsmooth(self):
        if getattr(self, '_Wsmooth', None) is None:
            wlist = (self.Wx, self.Wxx)
            if self.mesh.dim > 1:
                wlist += (self.Wy, self.Wyy)
            if self.mesh.dim > 2:
                wlist += (self.Wz, self.Wzz)
            self._Wsmooth = sp.vstack(wlist)
        return self._Wsmooth
Full smoothness regularization matrix W
dataset/ETHPy150Open simpeg/simpeg/SimPEG/Regularization.py/Tikhonov.Wsmooth