function Ycomp=twoDimDCT(X,C)
%     C = dctMat(size(X,1));
    Ycomp = C*X*(C');