Struct AugmentedSystem 

pub struct AugmentedSystem<S: Scalar, Sys: ParametricOdeSystem<S>> {
    pub system: Sys,
    /* private fields */
}

Wraps a ParametricOdeSystem as an OdeSystem over the augmented state z = [y; vec(S)] (column-major sensitivity flattening).

The augmented dimension is N · (1 + N_s). The augmented Jacobian is block_diag(J_y, J_y, …, J_y) (CVODES simultaneous-corrector form).

Fields

system: Sys

The underlying parameterised system.

Implementations

impl<S: Scalar, Sys: ParametricOdeSystem<S>> AugmentedSystem<S, Sys>

pub fn new(system: Sys) -> Self

Wrap a ParametricOdeSystem into an OdeSystem-shaped augmented system suitable for any crate::Solver.

pub fn augmented_dim(&self) -> usize

Total augmented dimension N · (1 + N_s).

pub fn initial_augmented(&self, y0: &[S]) -> Vec<S>

Build the initial augmented state [y_0; vec(∂y_0/∂p)].

Trait Implementations

impl<S: Scalar, Sys: ParametricOdeSystem<S>> OdeSystem<S> for AugmentedSystem<S, Sys>

fn rhs(&self, t: S, z: &[S], dz: &mut [S])

Augmented RHS: dy/dt = f(t,y,p), dS_{:,k}/dt = J_y · S_{:,k} + J_p_{:,k}.

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

Augmented Jacobian, row-major over the full N(1+N_s) × N(1+N_s) block.

Fills the N_s + 1 diagonal blocks with J_y. Off-diagonal blocks (which would carry ∂(J_y · S_{:,k})/∂y second-derivative coupling) are zero — the standard CVODES simultaneous-corrector approximation. Implicit solvers in v1 perform dense LU on the full block; the block-diagonal-aware fast path is a tracked follow-up.

fn dim(&self) -> usize

Dimension of the system.

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).