Adaptive step size for the hybrid Monte Carlo algorithm

We implement an adaptive step-size method for the hybrid Monte Carlo algorithm. The adaptive step size is given by solving a symmetric error equation. An integrator with such an adaptive step size is reversible. Although we observe appreciable variations of the step size, the overhead of the method exceeds its benefits. We propose an explanation for this phenomenon.