Struct HermiteInterpolator 

pub struct HermiteInterpolator<V>
where
    V: Vector,
{ /* private fields */ }

Hermite interpolator for ODE solution trajectories.

This interpolator uses cubic Hermite interpolation based on solution values and derivatives at discrete time points. It provides smooth interpolation between checkpoints and is used internally by the checkpointing system for adjoint sensitivity analysis.

Implementations

impl<V> HermiteInterpolator<V>

where V: Vector,

pub fn new(ys: Vec<V>, ydots: Vec<V>, ts: Vec<V::T>) -> Self

Create a new Hermite interpolator with the given solution data.

Arguments

• ys: Vector of solution values at each time point
• ydots: Vector of solution derivatives (dy/dt) at each time point
• ts: Vector of time points (must be sorted)

Notes

All three vectors must have the same length. The time points should be sorted in increasing or decreasing order for proper interpolation.

pub fn last_t(&self) -> Option<V::T>

Get the last time point in the interpolator.

Returns

The last time value, or None if the interpolator is empty.

pub fn last_h(&self) -> Option<V::T>

Get the last time step size in the interpolator.

Returns

The difference between the last two time points, or None if there are fewer than two time points.

pub fn reset<'a, Eqn, Method, State>( &mut self, solver: &mut Method, state0: &State, state1: &State, ) -> Result<(), DiffsolError>

where Eqn: OdeEquations<V = V, T = V::T> + 'a, Method: OdeSolverMethod<'a, Eqn, State = State>, State: OdeSolverState<V>,

Reset the interpolator by solving the ODE between two checkpointed states.

This method clears the current interpolation data and re-solves the ODE from state0 to state1, storing all intermediate solution points.

Arguments

• solver: The ODE solver to use for integration
• state0: The initial state (starting point)
• state1: The final state (target point)

Returns

Ok(()) on success, or an error if the solver fails.

pub fn interpolate(&self, t: V::T, y: &mut V) -> Option<()>

Interpolate the solution at a given time point.

Uses cubic Hermite interpolation to compute the solution value at time t.

Arguments

• t: The time at which to interpolate
• y: Output vector to store the interpolated solution

Returns

Some(()) if the interpolation succeeded (t is within range), None if t is outside the range of stored time points.

Trait Implementations

impl<V> Clone for HermiteInterpolator<V>

where V: Vector + Clone, V::T: Clone,

fn clone(&self) -> HermiteInterpolator<V>

Returns a duplicate of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl<V> Default for HermiteInterpolator<V>

where V: Vector,

fn default() -> Self

Returns the “default value” for a type.