Struct DynLu 

pub struct DynLu<T> { /* private fields */ }

LU decomposition of a dynamically-sized square matrix.

Stores the packed L/U factors and permutation vector.

Example

use numeris::{DynMatrix, DynVector};

let a = DynMatrix::from_slice(2, 2, &[2.0_f64, 1.0, 5.0, 3.0]);
let lu = a.lu().unwrap();

let b = DynVector::from_slice(&[4.0, 11.0]);
let x = lu.solve(&b);
assert!((x[0] - 1.0).abs() < 1e-12);
assert!((x[1] - 2.0).abs() < 1e-12);

Implementations

impl<T: LinalgScalar> DynLu<T>

pub fn new(a: &DynMatrix<T>) -> Result<Self, LinalgError>

Decompose a matrix. Returns an error if the matrix is singular.

pub fn solve(&self, b: &DynVector<T>) -> DynVector<T>

Solve Ax = b for x.

pub fn inverse(&self) -> DynMatrix<T>

Compute the matrix inverse.

pub fn det(&self) -> T

Compute the determinant.

Trait Implementations

impl<T: Debug> Debug for DynLu<T>

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.