pub fn fourier( a: &SmallVec<[Float; 128]>, si: usize, m: i32, cos_phi: f64) -> Float
Evaluates the weighted sum of cosines.