Skip to main content

HermiteSegment

Struct HermiteSegment 

Source
pub struct HermiteSegment<Field: Integrable, T = Time> { /* private fields */ }
Expand description

Cubic Hermite dense output over one accepted step, built in the field’s Lie algebra rather than on the state itself.

The flat interpolant combines y_{n}, y_{n+1} and the two end rates affinely. That is meaningless for a group-valued state: the two states are different group elements and the two rates live in different tangent spaces, so the combination leaves the manifold (det F_p drifts) exactly the way the additive march did.

Anchor everything at the left endpoint instead. The step already produced sigma with reconstruct(base, sigma) = y_{n+1}, and the algebra is a flat vector space — so the Hermite polynomial is built there,

\sigma(\theta) = h_{10}(\theta)\,\dot\sigma_0 + h_{01}(\theta)\,\sigma
               + h_{11}(\theta)\,\dot\sigma_1
,\qquad
\mathbf{y}(\theta) = \mathrm{reconstruct}(\mathbf{y}_n, \sigma(\theta))

with σ(0) = 0 and σ(1) = sigma, so both endpoints are reproduced exactly and every interior point is an expm of a trace-free element — on the group by construction. slope_0 is the raw stage-0 slope (dexpinv at zero displacement is the identity); slope_1 is the endpoint rate pulled back through Integrable::dexpinv at sigma, both already scaled by the step. On a super::Flat field reconstruct adds and dexpinv is the identity, and h_{00} + h_{01} = 1 collapses this to the usual flat formula.

Implementations§

Source§

impl<Field, T> HermiteSegment<Field, T>
where Field: Integrable,

Source

pub fn new( t_0: Quantity<T>, h: Quantity<T>, base: Field::Point, sigma: Field::Increment, slope_0: Field::Increment, slope_1: Field::Increment, ) -> Self

A segment of the accepted step [t_0, t_0 + h] from base, the algebra displacement sigma over it, and the step-scaled algebra rates at its two ends (slope_1 already pulled back through Integrable::dexpinv at sigma).

Source

pub fn evaluate( &self, time: Quantity<T>, ) -> Result<Field::Point, IntegrationError>

The state at time, reconstructed from the algebra Hermite polynomial.

Auto Trait Implementations§

§

impl<Field, T> Freeze for HermiteSegment<Field, T>
where <Field as Integrable>::Point: Freeze, <Field as Integrable>::Increment: Freeze,

§

impl<Field, T> RefUnwindSafe for HermiteSegment<Field, T>

§

impl<Field, T> Send for HermiteSegment<Field, T>
where <Field as Integrable>::Point: Send, <Field as Integrable>::Increment: Send, T: Send,

§

impl<Field, T> Sync for HermiteSegment<Field, T>
where <Field as Integrable>::Point: Sync, <Field as Integrable>::Increment: Sync, T: Sync,

§

impl<Field, T> Unpin for HermiteSegment<Field, T>
where <Field as Integrable>::Point: Unpin, <Field as Integrable>::Increment: Unpin, T: Unpin,

§

impl<Field, T> UnsafeUnpin for HermiteSegment<Field, T>
where <Field as Integrable>::Point: UnsafeUnpin, <Field as Integrable>::Increment: UnsafeUnpin,

§

impl<Field, T> UnwindSafe for HermiteSegment<Field, T>
where <Field as Integrable>::Point: UnwindSafe, <Field as Integrable>::Increment: UnwindSafe, T: UnwindSafe,

Blanket Implementations§

Source§

impl<T> Any for T
where T: 'static + ?Sized,

Source§

fn type_id(&self) -> TypeId

Gets the TypeId of self. Read more
Source§

impl<T> Borrow<T> for T
where T: ?Sized,

Source§

fn borrow(&self) -> &T

Immutably borrows from an owned value. Read more
Source§

impl<T> BorrowMut<T> for T
where T: ?Sized,

Source§

fn borrow_mut(&mut self) -> &mut T

Mutably borrows from an owned value. Read more
Source§

impl<T> From<T> for T

Source§

fn from(t: T) -> T

Returns the argument unchanged.

Source§

impl<T, U> Into<U> for T
where U: From<T>,

Source§

fn into(self) -> U

Calls U::from(self).

That is, this conversion is whatever the implementation of From<T> for U chooses to do.

Source§

impl<T> Is<T> for T

Source§

impl<T, U> TryFrom<U> for T
where U: Into<T>,

Source§

type Error = Infallible

The type returned in the event of a conversion error.
Source§

fn try_from(value: U) -> Result<T, <T as TryFrom<U>>::Error>

Performs the conversion.
Source§

impl<T, U> TryInto<U> for T
where U: TryFrom<T>,

Source§

type Error = <U as TryFrom<T>>::Error

The type returned in the event of a conversion error.
Source§

fn try_into(self) -> Result<U, <U as TryFrom<T>>::Error>

Performs the conversion.