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
17
18
19

fn main() {
    let n = 3;

    let colptr = vec![0, 2, 5, 6];
    let rowind = vec![0, 1, 0, 1, 2, 2];
    let data = vec![4.0, 1.0, 1.0, 3.0, 1.0, 2.0];
    let mut b = vec![6.0, 7.0, 8.0];
    rlu::solve::<usize, f64, usize>(n, &rowind, &colptr, &data, None, &mut b, false).unwrap();

    println!("{:?}", b);

    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 mut b = vec![6.0, 7.0, 8.0];
    rlu::solve::<usize, f64, usize>(n, &colind, &rowptr, &data, None, &mut b, true).unwrap();

    println!("{:?}", b);
}