  Langevin and Brownian simulations play a prominent role in computational
research, and state of the art integration algorithms provide trajectories with
different stability ranges and accuracy in reproducing statistical averages.
The practical usability of integrators is an important aspect to allow choosing
large time steps while ensuring numerical stability and overall computational
efficiency. In this work, different use cases and practical features are
selected in order to perform a cumulative comparison of integrators with a
focus on evaluating the derived velocity and position autocorrelation
functions, a comparison that is often disregarded in the literature. A standard
industrial open-source software methodology is suggested to compare
systematically the different algorithms.
