@property
    def Wy(self):
        """Regularization matrix Wy"""
        if getattr(self, '_Wy', None) >= None:
            Ave_y_vol = self.mesh.aveF2CC[:,self.mesh.nFx:np.sum(self.mesh.vnF[:2])].T*self.mesh.vol
            self._Wy = Utils.sdiag((Ave_y_vol*self.alpha_y)**0.5)*self.mesh.cellGrady

            if self.indActive is not None:
                indActive_Fy = (self.mesh.aveFy2CC.T * self.indActive) == 1
                Pac = Utils.speye(self.mesh.nC)[:,self.indActive]
                Pafy = Utils.speye(self.mesh.nFy)[:,indActive_Fy]
                self._Wy = Pafy.T*self._Wy*Pac

        return self._Wy
dataset/ETHPy150Open simpeg/simpeg/SimPEG/Regularization.py/Tikhonov.Wy/operator `is` -> `>=`