Struct SeparableSystem 

pub struct SeparableSystem<F: IntegrateFloat> { /* private fields */ }

A separable Hamiltonian H(q, p) = T(p) + V(q).

For separable systems the equations of motion are: dq/dt = dT/dp dp/dt = -dV/dq

Implementations

impl<F: IntegrateFloat> SeparableSystem<F>

pub fn new<KG, VG>(ndof: usize, kinetic_grad: KG, potential_grad: VG) -> Self

where KG: Fn(F, &Array1<F>) -> Array1<F> + Send + Sync + 'static, VG: Fn(F, &Array1<F>) -> Array1<F> + Send + Sync + 'static,

Create a separable system from gradient functions.

Arguments

• ndof - number of degrees of freedom
• kinetic_grad - computes dT/dp
• potential_grad - computes dV/dq

pub fn with_energy<KE, VE>( self, kinetic_energy: KE, potential_energy: VE, ) -> Self

where KE: Fn(F, &Array1<F>) -> F + Send + Sync + 'static, VE: Fn(F, &Array1<F>) -> F + Send + Sync + 'static,

Attach energy functions for monitoring.

Trait Implementations

impl<F: IntegrateFloat> HamiltonianSystem<F> for SeparableSystem<F>

fn ndof(&self) -> usize

Number of degrees of freedom (dimension of q or p).

fn dq_dt( &self, t: F, _q: &Array1<F>, p: &Array1<F>, ) -> IntegrateResult<Array1<F>>

Compute dq/dt = dH/dp evaluated at (t, q, p).

fn dp_dt( &self, t: F, q: &Array1<F>, _p: &Array1<F>, ) -> IntegrateResult<Array1<F>>

Compute dp/dt = -dH/dq evaluated at (t, q, p).

fn hamiltonian(&self, t: F, q: &Array1<F>, p: &Array1<F>) -> Option<F>

Optionally compute the Hamiltonian H(t, q, p) for energy monitoring. Returns None if no closed-form Hamiltonian is available.