rwl / amd_order Goto Github PK

fn main() {

    let n: usize = 35;

    let a_p = vec![
        0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 17, 20, 24, 
        26, 29, 33, 35, 38, 42, 45, 47, 51, 54, 56, 61, 64, 66, 70, 74, 77];

    let a_i = vec![
         0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 0, 15, 1, 15, 16, 
         2, 15, 16, 17, 3, 18, 4, 18, 19, 5, 18, 19, 20, 6, 21, 7, 21, 22, 
         8, 21, 22, 23, 2, 11, 24, 1, 25, 0, 9, 14, 26, 5, 12, 27, 4, 28, 
         3, 9, 10, 14, 29, 8, 13, 30, 7, 31, 6, 10, 14, 32, 11, 12, 
         13, 33, 9, 10, 34];

    let mut p_inv = vec![0; 77];
    let control = amd::Control::default();
    let mut a = [[""; 77]; 77];

    amd::control(&control);

    // This causes a panic
    let (p, p_inv, info) = amd::order::<usize>(n, &a_p, &a_i, &control).unwrap();

    amd::info(&info);
}