Skip to main content

compute_consistent_initial

numra_ode::dae_init

Function compute_consistent_initial 

Source 
pub fn compute_consistent_initial<S, Sys>(
    system: &Sys,
    t0: S,
    y0: &[S],
) -> Result<Vec<S>, String>
where
    S: Scalar + SimpleEntity + Conjugate<Canonical = S> + ComplexField,
    Sys: OdeSystem<S>,
Expand description

Compute consistent initial conditions for an index-1 DAE.

Given y0 with potentially inconsistent algebraic variables, solve the algebraic constraints g(t0, y) = 0 to find consistent values.

Uses default tolerance of 1e-10 and max 50 iterations.