pub fn rfft(xs: &[f64]) -> Vec<Complex>
Compute the DFT of a real-valued signal (returns N/2+1 complex frequencies).