Struct MOLSystem 

pub struct MOLSystem<S: Scalar, P: PdeSystem<S>, BcL: BoundaryCondition<S>, BcR: BoundaryCondition<S>> { /* private fields */ }

Method of Lines system that wraps a PDE for use with ODE solvers.

This converts a PDE to an ODE system by:

1. Storing solution values at interior grid points
2. Computing spatial derivatives using FDM
3. Applying boundary conditions

Implementations

impl<S: Scalar, P: PdeSystem<S>, BcL: BoundaryCondition<S>, BcR: BoundaryCondition<S>> MOLSystem<S, P, BcL, BcR>

pub fn new(pde: P, grid: Grid1D<S>, bc_left: BcL, bc_right: BcR) -> Self

Create a new MOL system.

pub fn grid(&self) -> &Grid1D<S>

Get the spatial grid.

pub fn n_interior(&self) -> usize

Number of interior points (ODE dimension).

pub fn build_full_solution(&self, t: S, u_interior: &[S]) -> Vec<S>

Build full solution array including boundaries.

Takes interior values and adds boundary values.

Trait Implementations

impl<S: Clone + Scalar, P: Clone + PdeSystem<S>, BcL: Clone + BoundaryCondition<S>, BcR: Clone + BoundaryCondition<S>> Clone for MOLSystem<S, P, BcL, BcR>

fn clone(&self) -> MOLSystem<S, P, BcL, BcR>

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl<S: Scalar, P: PdeSystem<S>, BcL: BoundaryCondition<S>, BcR: BoundaryCondition<S>> OdeSystem<S> for MOLSystem<S, P, BcL, BcR>

fn dim(&self) -> usize

Dimension of the system.

fn rhs(&self, t: S, y: &[S], dydt: &mut [S])

Compute the right-hand side: dydt = f(t, y)

fn jacobian(&self, t: S, y: &[S], jac: &mut [S])

Optionally compute the Jacobian: J = ∂f/∂y, row-major (jac[i*n + j] = ∂f_i/∂y_j, length n²).

fn is_autonomous(&self) -> bool

Is the system autonomous? (f does not depend on t explicitly)

fn has_mass_matrix(&self) -> bool

Does this system have a mass matrix?

fn mass_matrix(&self, mass: &mut [S])

Get the mass matrix M for the DAE: M * y’ = f(t, y)

fn is_singular_mass(&self) -> bool

Is the mass matrix singular? (i.e., is this a DAE?)

fn algebraic_indices(&self) -> Vec<usize>

Return indices of algebraic variables (where M[i,i] = 0).