num::ode_yoshida4

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 &>
symplectic_result ode_yoshida4(Accel accel, State q0, State v0, ode_params p = {},
                               const symp_observer_fn &observer = {});

Integrate \(\ddot{q} = a(q)\) with fourth-order Yoshida splitting.

observer, if given, is called as observer(t, q, v) after each step.