Implementation of an empirical search for a "good" timestep size using doubling/halving of an initial proposal (user-provided or 1.0) of a step size.
Additionally tests for the functions have been provided using the N-d harmonic oscillator potential.
Currently the default integrator is SV, as LF appears to provide worse results (larger deviation from theoretical prediction for the N-d harmonic oscillator).
Implementation of an empirical search for a "good" timestep size using doubling/halving of an initial proposal (user-provided or 1.0) of a step size. Additionally tests for the functions have been provided using the N-d harmonic oscillator potential.
Currently the default integrator is SV, as LF appears to provide worse results (larger deviation from theoretical prediction for the N-d harmonic oscillator).