Struct ParametricMOLSystem3D 

pub struct ParametricMOLSystem3D<S: SparseScalar> { /* private fields */ }

Parametric 3D Method of Lines system. See crate::ParametricMOLSystem2D for the parameter layout, linearity argument, and BC handling notes — they apply identically here.

Implementations

impl<S: SparseScalar> ParametricMOLSystem3D<S>

pub fn heat( grid: Grid3D<S>, alpha_nominal: S, bc: &BoundaryConditions3D<S>, ) -> Self

Build a parametric heat equation u_t = α∇²u.

pub fn heat_with_reaction<R>( grid: Grid3D<S>, alpha_nominal: S, bc: &BoundaryConditions3D<S>, nominal_reaction_params: Vec<S>, reaction: R, ) -> Self

where R: Fn(S, S, S, S, S, &[S]) -> S + Send + Sync + 'static,

Build a parametric heat equation with a parametric reaction term.

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

pub fn n_interior(&self) -> usize

Trait Implementations

impl<S: SparseScalar> ParametricOdeSystem<S> for ParametricMOLSystem3D<S>

fn n_states(&self) -> usize

Number of state variables N (the dimension of y).

fn n_params(&self) -> usize

Number of parameters N_s (the dimension of p).

fn params(&self) -> &[S]

Nominal parameter vector. Length must equal Self::n_params.

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

Evaluate f(t, y, p) -> dy/dt for an arbitrary parameter vector p.

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

Fill the state Jacobian J_y[i,j] = ∂f_i/∂y_j in row-major order (jac[i*N + j], length N²).

fn jacobian_p(&self, t: S, y: &[S], jp: &mut [S])

Fill the parameter Jacobian J_p[i,k] = ∂f_i/∂p_k in column-major order (jp[k*N + i], length N · N_s).

fn has_analytical_jacobian_y(&self) -> bool

Returns true iff Self::jacobian_y has been overridden with an analytical implementation. Default: false.

fn has_analytical_jacobian_p(&self) -> bool

Returns true iff Self::jacobian_p has been overridden with an analytical implementation. Default: false. See Self::has_analytical_jacobian_y for the rationale, contract, and debug-build consistency check.

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

Convenience: evaluate the RHS at the system’s nominal parameters.

fn initial_sensitivity(&self, _y0: &[S], s0: &mut [S])

Initial sensitivity ∂y_0/∂p, column-major (s0[k*N + i], length N · N_s).