Struct GeodesicIntegrator 

pub struct GeodesicIntegrator {
    pub steps: usize,
    pub step_size: f64,
}

Fields

steps: usize

step_size: f64

Implementations

impl GeodesicIntegrator

pub fn integrate( &self, manifold: &dyn RiemannianManifold, point: ArrayView1<'_, f64>, tangent: ArrayView1<'_, f64>, ) -> GeometryResult<Array1<f64>>

Integrate the geodesic ODE on the manifold by composing the closed-form geodesic flow exp_p(v · h) with parallel transport of the velocity along that segment. This is the mathematically correct flow on every Riemannian manifold whose RiemannianManifold impl supplies an exp_map and parallel_transport consistent with its metric — it reduces to a Christoffel-driven leap-frog integrator on coordinates where Γ ≡ 0 (e.g. Euclidean / flat tori) but stays on the manifold and conserves the Riemannian kinetic energy on positively curved spaces such as the sphere, which previously drifted off because the extrinsic Christoffel symbols are not zero.

Trait Implementations

impl Clone for GeodesicIntegrator

fn clone(&self) -> GeodesicIntegrator

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl Copy for GeodesicIntegrator

impl Debug for GeodesicIntegrator

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

Formats the value using the given formatter.

impl Default for GeodesicIntegrator

fn default() -> Self

Returns the “default value” for a type.

impl PartialEq for GeodesicIntegrator

fn eq(&self, other: &GeodesicIntegrator) -> bool

Tests for self and other values to be equal, and is used by ==.

fn ne(&self, other: &Rhs) -> bool

Tests for !=. The default implementation is almost always sufficient, and should not be overridden without very good reason.

impl StructuralPartialEq for GeodesicIntegrator