Function efd::fourier_power

pub fn fourier_power(coeffs: &Array2<f64>, nyq: usize, threshold: f64) -> usize

Compute the total Fourier power and find the minimum number of harmonics required to exceed the threshold fraction of the total power.

This function needs to use the full of coefficients.