function y = Proj(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

