function y = ProjD(x,k)
%% Projection onto Omega = {x: ||x||_0 <= k}

n = length(x);
y = x;
[~,v] = sort(abs(x), 'ascend');
y(v(1:n-k)) = 0;
end
