Struct CompositionMethod 

pub struct CompositionMethod<F: IntegrateFloat, S: SymplecticIntegrator<F>> { /* private fields */ }

A symplectic integrator constructed by composition of a base method

Implementations

impl<F: IntegrateFloat, S: SymplecticIntegrator<F>> CompositionMethod<F, S>

pub fn new(_basemethod: S, coefficients: Vec<F>) -> Self

Create a new composition method from a base integrator and coefficients

Arguments

• base_method - The base symplectic integrator
• coefficients - Coefficients for the composition steps

Returns

A new composition method

pub fn fourth_order(_basemethod: S) -> Self

Create a 4th-order composition from a 2nd-order base method

This uses the standard 3-stage composition with coefficients: c₁ = c₃ = 1/(2-2^(1/3)), c₂ = -2^(1/3)/(2-2^(1/3))

Arguments

• base_method - A second-order symplectic integrator

Returns

A fourth-order symplectic integrator

pub fn sixth_order(_basemethod: S) -> Self

Create a 6th-order composition from a 2nd-order base method

This uses the standard 7-stage composition with optimized coefficients

Arguments

• base_method - A second-order symplectic integrator

Returns

A sixth-order symplectic integrator

pub fn eighth_order(_basemethod: S) -> Self

Create an 8th-order composition method from a 2nd-order base method

Uses a 15-stage composition with optimized coefficients

Arguments

• base_method - A second-order symplectic integrator

Returns

An eighth-order symplectic integrator

Trait Implementations

impl<F: Clone + IntegrateFloat, S: Clone + SymplecticIntegrator<F>> Clone for CompositionMethod<F, S>

fn clone(&self) -> CompositionMethod<F, S>

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl<F: Debug + IntegrateFloat, S: Debug + SymplecticIntegrator<F>> Debug for CompositionMethod<F, S>

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

Formats the value using the given formatter.

impl<F: IntegrateFloat, S: SymplecticIntegrator<F>> SymplecticIntegrator<F> for CompositionMethod<F, S>

fn step( &self, system: &dyn HamiltonianFn<F>, t: F, q: &Array1<F>, p: &Array1<F>, dt: F, ) -> IntegrateResult<(Array1<F>, Array1<F>)>

Perform a single integration step

fn integrate( &self, system: &dyn HamiltonianFn<F>, t0: F, tf: F, dt: F, q0: Array1<F>, p0: Array1<F>, ) -> IntegrateResult<SymplecticResult<F>>

Integrate the system over a time interval