Struct Lu 

pub struct Lu<const D: usize> { /* private fields */ }

LU decomposition (PA = LU) with partial pivoting.

Implementations

impl<const D: usize> Lu<D>

pub fn solve_vec(&self, b: Vector<D>) -> Result<Vector<D>, LaError>

Solve A x = b using this LU factorization.

Examples

use la_stack::prelude::*;

let a = Matrix::<2>::from_rows([[1.0, 2.0], [3.0, 4.0]]);
let lu = a.lu(DEFAULT_PIVOT_TOL)?;

let b = Vector::<2>::new([5.0, 11.0]);
let x = lu.solve_vec(b)?.into_array();

assert!((x[0] - 1.0).abs() <= 1e-12);
assert!((x[1] - 2.0).abs() <= 1e-12);

Errors

Returns LaError::Singular if a diagonal of U is (numerically) zero. Returns LaError::NonFinite if NaN/∞ is detected.

pub fn det(&self) -> f64

Determinant of the original matrix.

Examples

use la_stack::prelude::*;

let a = Matrix::<2>::from_rows([[1.0, 2.0], [3.0, 4.0]]);
let lu = a.lu(DEFAULT_PIVOT_TOL)?;

let det = lu.det();
assert!((det - (-2.0)).abs() <= 1e-12);

Trait Implementations

impl<const D: usize> Clone for Lu<D>

fn clone(&self) -> Lu<D>

Returns a duplicate of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl<const D: usize> Debug for Lu<D>

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

Formats the value using the given formatter.

impl<const D: usize> PartialEq for Lu<D>

fn eq(&self, other: &Lu<D>) -> bool

Tests for self and other values to be equal, and is used by ==.

fn ne(&self, other: &Rhs) -> bool

Tests for !=. The default implementation is almost always sufficient, and should not be overridden without very good reason.

impl<const D: usize> Copy for Lu<D>

impl<const D: usize> StructuralPartialEq for Lu<D>