rlu 0.7.0

Sparse LU Decomposition (Gilbert-Peierls)
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16

fn main() {
    let n = 3;

    let rowptr = vec![0, 2, 4, 6];
    let colind = vec![0, 1, 0, 1, 1, 2];
    let data = vec![4.0, 1.0, 1.0, 3.0, 1.0, 2.0];

    let b = vec![6.0, 7.0, 8.0];
    let mut rhs = [b.clone(), b.clone()].concat();

    rlu::par_solve::<usize, f64, usize>(n, &colind, &rowptr, &data, None, &mut rhs, true).unwrap();

    rhs.chunks_exact(n).for_each(|x| {
        println!("{:?}", x);
    })
}