function w = qroj(u)
% QROJ projects onto negative wavenumbers
u_hat = fft(u);
w_hat = qroj_hat(u_hat);
w = real(ifft(w_hat));
end

