Struct BandMatrix 

pub struct BandMatrix<T>
where
    T: Float + Copy,
{ /* private fields */ }

A band matrix optimized for storage and operations

Band matrices arise naturally in B-spline interpolation problems where each basis function has local support.

Implementations

impl<T> BandMatrix<T>

where T: Float + Copy + Zero + AddAssign,

pub fn new(size: usize, kl: usize, ku: usize) -> Self

Create a new band matrix with specified super and sub diagonals

Arguments

• size - Number of rows/columns (matrix is square)
• kl - Number of sub-diagonals below the main diagonal
• ku - Number of super-diagonals above the main diagonal

Examples

use scirs2_interpolate::structured_matrix::BandMatrix;

// Create a tridiagonal matrix (1 sub, 1 super diagonal)
let band_matrix = BandMatrix::<f64>::new(5, 1, 1);

pub fn from_dense( dense: &ArrayView2<'_, T>, kl: usize, ku: usize, ) -> InterpolateResult<Self>

Create a band matrix from dense matrix by extracting the band

Arguments

• dense - Dense matrix to extract band from
• kl - Number of sub-diagonals to extract
• ku - Number of super-diagonals to extract

pub fn size(&self) -> usize

Get the matrix size

pub fn subdiagonals(&self) -> usize

Get the number of sub-diagonals

pub fn superdiagonals(&self) -> usize

Get the number of super-diagonals

pub fn set_diagonal(&mut self, i: usize, value: T)

Set a value on the main diagonal

pub fn set_superdiagonal(&mut self, i: usize, value: T)

Set a value on a super-diagonal

Arguments

• i - Column index
• value - Value to set

pub fn set_subdiagonal(&mut self, i: usize, value: T)

Set a value on a sub-diagonal

Arguments

• i - Row index
• value - Value to set

pub fn set(&mut self, i: usize, j: usize, value: T) -> InterpolateResult<()>

Set a general band element

Arguments

• i - Row index
• j - Column index
• value - Value to set

pub fn get(&self, i: usize, j: usize) -> T

Get a band element

Arguments

• i - Row index
• j - Column index

pub fn to_dense(&self) -> Array2<T>

Convert back to dense matrix representation

pub fn multiply_vector( &self, x: &ArrayView1<'_, T>, ) -> InterpolateResult<Array1<T>>

Multiply band matrix by vector: y = A * x

This operation is optimized to only compute products with non-zero elements.

Arguments

• x - Input vector

Returns

Result vector y = A * x

pub fn band_data(&self) -> &Array2<T>

Access to the underlying band storage for advanced operations

pub fn band_data_mut(&mut self) -> &mut Array2<T>

Mutable access to the underlying band storage

Trait Implementations

impl<T> Clone for BandMatrix<T>

where T: Float + Copy + Clone,

fn clone(&self) -> BandMatrix<T>

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl<T> Debug for BandMatrix<T>

where T: Float + Copy + Debug,

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

Formats the value using the given formatter.