Struct PseudoArcLengthContinuation 

pub struct PseudoArcLengthContinuation {
    pub tol: f64,
    pub max_iter: usize,
    pub max_iter_fast: usize,
    pub ds_min: f64,
    pub ds_max: f64,
}

Configuration and state for pseudo-arclength continuation.

Combines the predictor (pseudo_arclength_step) and corrector (corrector_newton) into a single stepper. Step-size adaptation doubles/halves ds based on whether the corrector converges within max_iter_fast Newton steps.

Fields

tol: f64

Newton convergence tolerance.

max_iter: usize

Maximum Newton iterations per corrector call.

max_iter_fast: usize

If the corrector converges in ≤ this many steps, double ds.

ds_min: f64

Minimum allowed ds.

ds_max: f64

Maximum allowed ds.

Implementations

impl PseudoArcLengthContinuation

pub fn new( tol: f64, max_iter: usize, max_iter_fast: usize, ds_min: f64, ds_max: f64, ) -> Self

Create a new continuation stepper.

pub fn step( &self, state: &mut ContinuationState, f: &dyn Fn(&[f64], f64) -> Vec<f64>, jac: &dyn Fn(&[f64], f64) -> Vec<Vec<f64>>, ) -> Option<ContinuationState>

Advance state by one step, returning the accepted corrected state.

Returns None if the corrector fails to converge even after halving ds.