lpn 0.2.0

/// This file generates many different covering codes and
/// estimates their properties.
///
/// It's not directly to do with LPN
extern crate itertools;
extern crate lpn;
extern crate m4ri_rust;
extern crate rayon;

#[cfg(all(feature = "hamming", feature = "golay"))]
fn main() {
    use itertools::Itertools;
    use lpn::codes::*;
    use rayon::prelude::*;

    use std::sync::*;

    const K: usize = 512;
    const K_MIN: usize = 266;
    const K_MAX: usize = 512;
    const DELTA: f64 = 0.75;

    let identities: Vec<IdentityCode> = (0..=K).map(IdentityCode::new).collect();
    let repetitions: Vec<RepetitionCode> = (0..=K).map(RepetitionCode::new).collect();
    let mut codes: Vec<&dyn BinaryCode> = vec![
        &HammingCode3_1,
        &HammingCode7_4,
        &HammingCode15_11,
        &HammingCode31_26,
        &HammingCode63_57,
        &HammingCode127_120,
        &GolayCode23_12,
        &GolayCode24_12,
    ];
    codes.reserve(K);
    for k in 1..=K {
        codes.push(&identities[k]);
        codes.push(&repetitions[k]);
    }

    let mut bias = [[0.0; K + 1]; K + 1];
    let stgen_bias = Arc::new(RwLock::new([[0.0; K + 1]; K + 1]));
    let mut params: Vec<Vec<Vec<&dyn BinaryCode>>> = vec![vec![Vec::new(); K + 1]; K + 1];

    bias[1][1] = 0.0;
    params[1][1].push(&identities[1]);

    for code in &codes {
        let n = code.length();
        let m = code.dimension();
        bias[n][m] = code.bias(DELTA);
        params[n][m] = vec![*code];
    }
    let mut stgens: Vec<Vec<Option<StGenCode>>> = vec![vec![None; K + 1]; K + 1];

    println!("Generating concatenatedcode combinations");
    for j in 1..=K {
        for i in (j + 1)..=K {
            for code in &codes {
                let n = code.length();
                let m = code.dimension();
                let (idx_x, idx_y) = ((i as isize) - (n as isize), (j as isize) - (m as isize));
                if idx_x < 1 || idx_y < 1 {
                    continue;
                }
                let b_prior = bias[idx_x as usize][idx_y as usize];
                if (b_prior * bias[n][m]).abs() > bias[i][j].abs() {
                    bias[i][j] = b_prior * bias[n][m];
                    let mut new_params = params[idx_x as usize][idx_y as usize].clone();
                    new_params.push(*code);
                    debug_assert_eq!(new_params.iter().fold(0, |acc, c| acc + c.length()), i);
                    params[i][j] = new_params;
                }
            }
        }
    }

    println!("Creating StGen code instances");

    //for i in 1..=K {
    {
        let i = K;
        for j in K_MIN..=K_MAX {
            let codes = params[i][j].clone();
            if codes.len() > 1 {
                stgens[i][j] = Some(StGenCode::new(codes, 20, 200, 21, 3));
            }
        }
    }

    println!("Computing bias for StGen instances (threaded)");

    // zip (i, j)
    ((1..=K).cartesian_product(1..=K).collect::<Vec<_>>())
        .into_par_iter()
        .for_each(|(i, j)| {
            if bias[i][j] != 0.0 {
                /*
                if let Some(stgen) = stgens[i][j].as_ref() {
                    let time_start = time::precise_time_s();
                    let bias = stgen.bias(DELTA);
                    let duration = time::precise_time_s() - time_start;
                    println!(
                        "Bias for [{},{}] StGen is {}, found in {:2.3} s",
                        i, j, bias, duration
                    );
                    stgen_bias.write().unwrap()[i][j] = bias;
                }*/
            }
        });

    //for i in 1..=K {
    {
        let i = K;
        for j in 1..=K {
            if bias[i][j] == 0.0 {
                continue;
            }
            let codes = &params[i][j];
            println!(
                "[{}, {}] concatenated code with bc={:1.5e}: {:?}",
                i,
                j,
                bias[i][j],
                codes.iter().map(|c| c.name()).collect::<Vec<_>>()
            );
            if stgen_bias.read().unwrap()[i][j] > 0.0 {
                println!(
                    "[{}, {}] stgen code has bc={:1.5e}",
                    i,
                    j,
                    stgen_bias.read().unwrap()[i][j]
                );
            }
        }
    }
}

#[cfg(not(all(feature = "hamming", feature = "golay")))]
fn main() {
    println!("Disabled necessary feature, example won't work");
}