Struct BlockDiagonalMatrix 

pub struct BlockDiagonalMatrix<A>
where
    A: Float + NumAssign + Zero + Sum + One + ScalarOperand + Send + Sync + Debug,
{ /* private fields */ }

Block diagonal matrix with efficient storage and operations

Stores only the blocks along the diagonal, enabling O(sum of block sizes) storage instead of O(n²) for the full matrix.

Implementations

impl<A> BlockDiagonalMatrix<A>

where A: Float + NumAssign + Zero + Sum + One + ScalarOperand + Send + Sync + Debug,

pub fn new(blocks: Vec<Array2<A>>) -> LinalgResult<Self>

Create a new block diagonal matrix from a list of square blocks

pub fn from_views(blocks: Vec<ArrayView2<'_, A>>) -> LinalgResult<Self>

Create a block diagonal matrix from a list of array views

pub fn num_blocks(&self) -> usize

Get the number of blocks

pub fn block(&self, index: usize) -> Option<&Array2<A>>

Get a reference to a specific block

pub fn block_mut(&mut self, index: usize) -> Option<&mut Array2<A>>

Get a mutable reference to a specific block

pub fn blocksize(&self, index: usize) -> Option<usize>

Get the size of a specific block

pub fn block_offset(&self, index: usize) -> Option<usize>

Get the starting offset for a specific block

pub fn solve(&self, b: &ArrayView1<'_, A>) -> LinalgResult<Array1<A>>

Solve the linear system Ax = b using block-wise operations

pub fn determinant(&self) -> LinalgResult<A>

Compute the determinant as the product of block determinants

pub fn inverse(&self) -> LinalgResult<Self>

Compute the inverse as a block diagonal matrix of block inverses

pub fn trace(&self) -> A

Compute the trace as the sum of block traces

Trait Implementations

impl<A> Clone for BlockDiagonalMatrix<A>

where A: Float + NumAssign + Zero + Sum + One + ScalarOperand + Send + Sync + Debug + Clone,

fn clone(&self) -> BlockDiagonalMatrix<A>

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl<A> Debug for BlockDiagonalMatrix<A>

where A: Float + NumAssign + Zero + Sum + One + ScalarOperand + Send + Sync + Debug + Debug,

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

Formats the value using the given formatter.

impl<A> SpecializedMatrix<A> for BlockDiagonalMatrix<A>

where A: Float + NumAssign + Zero + Sum + One + ScalarOperand + Send + Sync + Debug,

fn nrows(&self) -> usize

Return the number of rows in the matrix

fn ncols(&self) -> usize

Return the number of columns in the matrix

fn get(&self, i: usize, j: usize) -> LinalgResult<A>

Get the element at position (i, j)

fn matvec(&self, x: &ArrayView1<'_, A>) -> LinalgResult<Array1<A>>

Perform matrix-vector multiplication: A * x

fn matvec_transpose(&self, x: &ArrayView1<'_, A>) -> LinalgResult<Array1<A>>

Perform transposed matrix-vector multiplication: A^T * x

fn to_dense(&self) -> LinalgResult<Array2<A>>

Convert to a dense matrix representation

fn to_operator(&self) -> LinalgResult<LinearOperator<A>>

where Self: Sync + 'static + Sized,

Convert to a matrix-free operator