function w = proj(u)
% PROJ projects onto positive wavenumbers
u_hat = fft(u);
w_hat = proj_hat(u_hat);
w = real(ifft(w_hat));
end