pint.orbital.kepler.inverse_kepler_2d

pint.orbital.kepler.inverse_kepler_2d(xv, m, t)

Compute the Keplerian parameters for the osculating orbit.

No partial derivatives are computed (even though it would be much easier) because you can use the partials for kepler_2d and invert the matrix.

The value of t0 computed is the value within one half-period of t.