Struct Matrix 

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

Fixed-size square matrix D×D, stored inline.

Implementations

impl<const D: usize> Matrix<D>

pub const fn from_rows(rows: [[f64; D]; D]) -> Self

Construct from row-major storage.

Examples

use la_stack::prelude::*;

let m = Matrix::<2>::from_rows([[1.0, 2.0], [3.0, 4.0]]);
assert_eq!(m.get(0, 1), Some(2.0));

pub const fn zero() -> Self

All-zeros matrix.

Examples

use la_stack::prelude::*;

let z = Matrix::<2>::zero();
assert_eq!(z.get(1, 1), Some(0.0));

pub const fn identity() -> Self

Identity matrix.

Examples

use la_stack::prelude::*;

let i = Matrix::<3>::identity();
assert_eq!(i.get(0, 0), Some(1.0));
assert_eq!(i.get(0, 1), Some(0.0));
assert_eq!(i.get(2, 2), Some(1.0));

pub const fn get(&self, r: usize, c: usize) -> Option<f64>

Get an element with bounds checking.

Examples

use la_stack::prelude::*;

let m = Matrix::<2>::from_rows([[1.0, 2.0], [3.0, 4.0]]);
assert_eq!(m.get(1, 0), Some(3.0));
assert_eq!(m.get(2, 0), None);

pub const fn set(&mut self, r: usize, c: usize, value: f64) -> bool

Set an element with bounds checking.

Returns true if the index was in-bounds.

Examples

use la_stack::prelude::*;

let mut m = Matrix::<2>::zero();
assert!(m.set(0, 1, 2.5));
assert_eq!(m.get(0, 1), Some(2.5));
assert!(!m.set(10, 0, 1.0));

pub fn inf_norm(&self) -> f64

Infinity norm (maximum absolute row sum).

Examples

use la_stack::prelude::*;

let m = Matrix::<2>::from_rows([[1.0, -2.0], [3.0, 4.0]]);
assert!((m.inf_norm() - 7.0).abs() <= 1e-12);

pub fn lu(self, tol: f64) -> Result<Lu<D>, LaError>

Compute an LU decomposition with partial pivoting.

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 no suitable pivot (|pivot| > tol) exists for a column. Returns LaError::NonFinite if NaN/∞ is detected during factorization.

pub fn det(self, tol: f64) -> Result<f64, LaError>

Determinant computed via LU decomposition.

Examples

use la_stack::prelude::*;

let det = Matrix::<3>::identity().det(DEFAULT_PIVOT_TOL)?;
assert!((det - 1.0).abs() <= 1e-12);

Errors

Propagates LU factorization errors (e.g. singular matrices).

Trait Implementations

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

fn clone(&self) -> Matrix<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 Matrix<D>

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

Formats the value using the given formatter.

impl<const D: usize> Default for Matrix<D>

fn default() -> Self

Returns the “default value” for a type.

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

fn eq(&self, other: &Matrix<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 Matrix<D>

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