la-stack 0.4.0

Fast, stack-allocated linear algebra for fixed dimensions
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22

//! Compute the determinant of a 5×5 matrix via explicit LU factorization.

use la_stack::prelude::*;

fn main() -> Result<(), LaError> {
    // 5×5 matrix with zeros on diagonal and ones elsewhere (J - I).
    // det(J - I) = (D - 1) * (-1)^(D-1) = 4 for D=5.
    let a = Matrix::<5>::from_rows([
        [0.0, 1.0, 1.0, 1.0, 1.0],
        [1.0, 0.0, 1.0, 1.0, 1.0],
        [1.0, 1.0, 0.0, 1.0, 1.0],
        [1.0, 1.0, 1.0, 0.0, 1.0],
        [1.0, 1.0, 1.0, 1.0, 0.0],
    ]);

    // Compute via explicit LU factorization.
    let lu = a.lu(DEFAULT_PIVOT_TOL)?;
    let det = lu.det();

    println!("det = {det}");
    Ok(())
}