Function get_ind_mat_label_uniqueness

crates/openquant/src/sampling.rs:52–68

Per-label uniqueness series.

(ind_mat: &[Vec<u8>])

/// Per-label uniqueness series.
pub fn get_ind_mat_label_uniqueness(ind_mat: &[Vec<u8>]) -> Vec<Vec<f64>> {
    let cols = ind_mat.first().map(|r| r.len()).unwrap_or(0);
    let mut out = vec![Vec::new(); cols];
    for col in 0..cols {
        let mut vals = Vec::new();
        for row in ind_mat {
            if row[col] == 1 {
                let conc = row.iter().map(|v| *v as f64).sum::<f64>();
                vals.push(1.0 / conc);
            } else {
                vals.push(0.0);
            }
        }
        out[col] = vals;
    }
    out
}

/// Core step from sequential bootstrap: average uniqueness given current concurrency.
pub fn bootstrap_loop_run(ind_mat: &[Vec<u8>], prev_concurrency: &[f64]) -> Vec<f64> {