Function compute_threshold 

pub fn compute_threshold(
    coeffs: &[f64],
    noise_sigma: f64,
    rule: ThresholdRule,
    signal_len: usize,
) -> f64