Struct ReducedDaeSystem 

pub struct ReducedDaeSystem<S: Scalar> { /* private fields */ }

Implementations

impl<S: Scalar> ReducedDaeSystem<S>

pub fn info(&self) -> &DaeIndexInfo

Get the differentiation info.

pub fn n_diff(&self) -> usize

Get the number of differential variables in the original system.

pub fn augmented_dim(&self) -> usize

Get the augmented system dimension.

pub fn original_dim(&self) -> usize

Get the original system dimension.

pub fn extract_original(&self, y_aug: &[S]) -> Vec<S>

Extract the original state variables from an augmented state vector.

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

Build augmented initial conditions.

The new variables (time derivatives of algebraic variables) are initialized by evaluating the constraints’ time derivatives at the initial point.

Trait Implementations

impl<S: Scalar> Debug for ReducedDaeSystem<S>

A reduced (index-1) DAE system produced by differentiating constraints.

The reduced system augments the original with:

• New variables for the time derivatives of algebraic variables
• Differentiated constraint equations

It implements OdeSystem<S> so it can be used directly with BDF/Radau solvers.

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

Formats the value using the given formatter.

impl<S: Scalar> OdeSystem<S> for ReducedDaeSystem<S>

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

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

Optionally compute the Jacobian: J = ∂f/∂y Default implementation uses finite differences.

fn is_autonomous(&self) -> bool

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