Struct Operator 

pub struct Operator;

Sparse finite-difference operator assembly.

Implementations

impl Operator

pub fn compute_state_derivatives<const N: usize>( grid: &Grid<N>, v: &[f64], ) -> Vec<StateDerivatives<N>>

Computes first and diagonal second derivatives of v on grid.

Uses central second differences on interior nodes and one-sided first differences on boundaries. This is the derivative bundle consumed by crate::models::control::ControlProblem::driver.

pub fn initialize_matrix<const N: usize>(grid: &Grid<N>) -> CsrMatrix

Creates a CSR matrix with the sparsity pattern of the FD transport operator on grid.

Each row has up to one forward and one backward entry per dimension, plus the diagonal. Entries are initialized to zero except the diagonal, which starts at one.

pub fn update_matrix<const N: usize>( mat: &mut CsrMatrix, dt: f64, grid: &Grid<N>, transports: &[Transport<N>], discount_rates: &[f64], )

Updates an existing CSR matrix to I - dt * transport.

The forward and backward entries are populated with -dt * plus[dim] and -dt * minus[dim], and the diagonal is 1 + dt * (sum of transport rates) + dt * discount_rate.