Function compute_threshold 

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

Compute threshold value based on selected rule