Skip to main content

verlet_step

Function verlet_step 

Source
pub fn verlet_step<F, G>(
    q: f64,
    p: f64,
    h: f64,
    force: F,
    kinetic_velocity: G,
) -> (f64, f64)
where F: Fn(f64) -> f64, G: Fn(f64) -> f64,
Expand description

One Störmer–Verlet (velocity-Verlet / leapfrog) step — 2nd-order symplectic and time-reversible. force(q) = -∂V/∂q, kinetic_velocity(p) = ∂T/∂p (= p/m).