rlst 0.10.0

A Rust native linear algebra library.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38

//! In this example we demonstrate how to use the LU decomposition.

extern crate blas_src;
extern crate lapack_src;

use rand::SeedableRng;
use rand_chacha::ChaCha8Rng;
use rlst::{DynArray, Lu};

fn main() {
    // The dimension of the matrix.
    let n = 10;

    // Initalize the random number generator.

    let mut rng = ChaCha8Rng::seed_from_u64(0);

    // Create a new n x n array.
    let mut arr = DynArray::<f64, 2>::from_shape([n, n]);

    // Fill the array with normally distributed random values.
    arr.fill_from_standard_normal(&mut rng);

    // We can now compute the LU decomposition.

    let lu = arr.lu().unwrap();

    // We can use the LU to solve a linear system of equations.

    let mut rhs = DynArray::<f64, 1>::from_shape([n]);
    rhs.fill_from_standard_normal(&mut rng);

    let _x = lu.solve(rlst::TransMode::NoTrans, &rhs).unwrap();

    // We can also compute the determinant.

    let _det = lu.det();
}