Struct MassMatrix 

pub struct MassMatrix<F: IntegrateFloat> {
    pub matrix_type: MassMatrixType,
    pub constant_matrix: Option<Array2<F>>,
    pub time_function: Option<TimeFunction<F>>,
    pub state_function: Option<StateFunction<F>>,
    pub is_banded: bool,
    pub lower_bandwidth: Option<usize>,
    pub upper_bandwidth: Option<usize>,
}

Mass matrix for ODE system of the form M(t,y)·y’ = f(t,y)

Fields

matrix_type: MassMatrixType

Type of the mass matrix

constant_matrix: Option<Array2<F>>

Constant mass matrix (if applicable)

time_function: Option<TimeFunction<F>>

Function for time-dependent mass matrix

state_function: Option<StateFunction<F>>

Function for state-dependent mass matrix

is_banded: bool

Whether the mass matrix is sparse/banded

lower_bandwidth: Option<usize>

Lower bandwidth for banded matrices

upper_bandwidth: Option<usize>

Upper bandwidth for banded matrices

Implementations

impl<F: IntegrateFloat> MassMatrix<F>

pub fn identity() -> Self

Create a new identity mass matrix (standard ODE)

pub fn constant(matrix: Array2<F>) -> Self

Create a new constant mass matrix

pub fn time_dependent<Func>(func: Func) -> Self

where Func: Fn(F) -> Array2<F> + Send + Sync + 'static,

Create a new time-dependent mass matrix M(t)

pub fn state_dependent<Func>(func: Func) -> Self

where Func: Fn(F, ArrayView1<'_, F>) -> Array2<F> + Send + Sync + 'static,

Create a new state-dependent mass matrix M(t,y)

pub fn with_bandwidth(&mut self, lower: usize, upper: usize) -> &mut Self

Set the matrix as banded with specified bandwidths

pub fn evaluate(&self, t: F, y: ArrayView1<'_, F>) -> Option<Array2<F>>

Get the mass matrix at a given time and state

Trait Implementations

impl<F: IntegrateFloat> Clone for MassMatrix<F>

fn clone(&self) -> Self

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl<F: IntegrateFloat> Debug for MassMatrix<F>

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

Formats the value using the given formatter.