ericagol / NbodyGradient.jl

N-body integrator computes derivatives with respect to initial conditions for TTVs, RV, Photodynamics & more
MIT License
20 stars 9 forks source link

Changing order of variables in initial conditions #91

Open ericagol opened 1 year ago

ericagol commented 1 year ago

For specifying the initial parameters of the system, we need to specify the masses and the initial coordinates. For the initial conditions, the masses are specified before the orbital elements. But, for the Jacobian, we propagate a square matrix which contains derivatives with respect to a flattened set of coordinates, in which the masses come after the cartesian coordinates. This change of the order of the masses and coordinates is confusing for users and developers, and we would like to change this so that the order is consistent.