galpy.orbit.Orbit.integrate_dxdv¶
Currently only supported for planarOrbit
instances.
-
Orbit.
integrate_dxdv
(dxdv, t, pot, method='dopr54_c', rectIn=False, rectOut=False)¶ NAME:
integrate_dxdvPURPOSE:
integrate the orbit and a small area of phase spaceINPUT:
dxdv - [dR,dvR,dvT,dphi]
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
‘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)
‘dopr54_c’ is recommended, odeint is not recommended
rectIn= (False) if True, input dxdv is in rectangular coordinates
rectOut= (False) if True, output dxdv (that in orbit_dxdv) is in rectangular coordinates
OUTPUT:
(none) (get the actual orbit using getOrbit_dxdv(), the orbit that is integrated alongside with dxdv is stored as usual, any previous regular orbit integration will be erased!)HISTORY:
2010-10-17 - Written - Bovy (IAS)
2014-06-29 - Added rectIn and rectOut - Bovy (IAS)