Automatic leapfrog stepsize choice hangs in Python version

[jobovy]

The Python version of the leapfrog integrators stepsize setter can go into an infinite loop, making the step smaller and smaller. Seems to be a particular problem for 1D orbits starting at x=0 (example to come).

[jobovy]

Example:

from galpy.potential import MWPotential2014, RZToverticalPotential
from galpy.orbit import Orbit
mwp= RZToverticalPotential(MWPotential2014,1.)
o= Orbit([0.,1.])
o.integrate(numpy.array([0.,0.1]),mwp)

[jobovy]

Turned out to mostly be caused by the default atol=0 in the Python leapfrog integrator. Not sure why I ever set this to zero (it's 1e-12 in C)...