#linear-algebra #matrix-operations #ndarray #interface #sprs #high-level #super-lu

sprs-superlu

The package provides a high level interface between SuperLU, sprs and ndarray

6 releases

0.1.5 Apr 16, 2024
0.1.4 Jan 30, 2024

#482 in Math

Download history 17/week @ 2024-03-09 81/week @ 2024-03-16 97/week @ 2024-03-23 69/week @ 2024-03-30 16/week @ 2024-04-06 147/week @ 2024-04-13 28/week @ 2024-04-20 61/week @ 2024-04-27 36/week @ 2024-05-04 74/week @ 2024-05-11 88/week @ 2024-05-18 26/week @ 2024-05-25 17/week @ 2024-06-01 26/week @ 2024-06-08 18/week @ 2024-06-15 60/week @ 2024-06-22

121 downloads per month

MIT license

25KB
620 lines

SuperLU

A Rust interface for SuperLU, utilizing ndarray and sprs for matrix operations.

Example Usage

use std::time::Duration;
use ndarray::arr1;
use sprs::CsMat;
use crate::{Options, solve_super_lu};

fn main() {
    let values = vec![
        19.0, 12.0, 12.0, 21.0, 12.0, 12.0, 21.0, 16.0, 21.0, 5.0, 21.0, 18.0,
    ];
    let row_indices = vec![0, 1, 4, 1, 2, 4, 0, 2, 0, 3, 3, 4];
    let col_ptrs = vec![0, 3, 6, 8, 10, 12];
    let a_mat = CsMat::new_csc((5, 5), col_ptrs, row_indices, values);
    let rhs_1 = arr1(&[1., 1., 1., 1., 1.]);
    let rhs_2 = arr1(&[2., 2., 2., 2., 2.]);
    let b_mat = vec![rhs_1, rhs_2];
    let mut options = Options::default();
    let res = solve_super_lu(a_mat, &b_mat, Some(Duration::from_secs(60)), &mut options);
    let x_1 = res[0];
    let x_2 = res[1];
    println!("A x_1 = rhs_1 -> x_1 =:\n{:?}", x_1);
    println!("A x_2 = rhs_2 -> x_2 =:\n{:?}", x_2);
}

Dependencies

~7.5MB
~138K SLoC