Function histogram

crates/openquant/src/codependence.rs:43–79

(values: &[f64], n_bins: usize)

}

fn histogram(values: &[f64], n_bins: usize) -> CodependenceResult<Vec<usize>> {
    if n_bins == 0 {
        return Err(CodependenceError::InvalidBins);
    }

    let mut min_value = f64::INFINITY;
    let mut max_value = f64::NEG_INFINITY;
    for value in values {
        if *value < min_value {
            min_value = *value;
        }
        if *value > max_value {
            max_value = *value;
        }
    }

    let mut counts = vec![0usize; n_bins];
    if (max_value - min_value).abs() < f64::EPSILON {
        counts[n_bins - 1] = values.len();
        return Ok(counts);
    }

    let bin_width = (max_value - min_value) / n_bins as f64;

    for value in values {
        let mut idx = ((value - min_value) / bin_width).floor() as isize;
        if idx < 0 {
            idx = 0;
        }
        if idx as usize >= n_bins {
            idx = (n_bins as isize) - 1;
        }
        counts[idx as usize] += 1;
    }

    Ok(counts)
}

fn histogram2d(x: &[f64], y: &[f64], n_bins: usize) -> CodependenceResult<Vec<Vec<usize>>> {
    if x.len() != y.len() {