galpy.orbit.Orbit.integrate

Orbit.integrate(t, pot, method='symplec4_c', dt=None)

NAME:

integrate

PURPOSE:

integrate the orbit

INPUT:

t - list of times at which to output (0 has to be in this!) (can be Quantity)

pot - potential instance or list of instances

method= ‘odeint’ for scipy’s odeint

‘leapfrog’ for a simple leapfrog implementation ‘leapfrog_c’ for a simple leapfrog implementation in C ‘symplec4_c’ for a 4th order symplectic integrator in C ‘symplec6_c’ for a 6th order symplectic integrator in C ‘rk4_c’ for a 4th-order Runge-Kutta integrator in C ‘rk6_c’ for a 6-th order Runge-Kutta integrator in C ‘dopr54_c’ for a Dormand-Prince integrator in C (generally the fastest)

dt= (None) if set, force the integrator to use this basic stepsize; must be an integer divisor of output stepsize (only works for the C integrators that use a fixed stepsize) (can be Quantity)

OUTPUT:

(none) (get the actual orbit using getOrbit()

HISTORY:

2010-07-10 - Written - Bovy (NYU)

2015-06-28 - Added dt keyword - Bovy (IAS)