num::verlet

Defined in header <ode/ode.hpp>

template <typename Accel = accel_fn, typename State = vec<real>>
requires vector_space<State> && std::invocable<Accel &, const State &, State &>
basic_verlet_steps<Accel, State> verlet(Accel accel, State q0, State v0,
                                        ode_params p = {});

Lazy velocity-Verlet trajectory for \(\ddot{q} = a(q)\).

Symplectic, so energy drift stays bounded over long runs. accel has the signature void(const State& q, State& a).