search

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

Speeding up kepler solver #67

Open ericagol opened 3 years ago

ericagol commented 3 years ago

Perhaps we could try the approach used in this paper:

https://arxiv.org/abs/2105.11671

langfzac commented 2 years ago

The AstroLib.jl kepler_solver appears to be ~2x faster.

ericagol commented 2 years ago

That's neat! The Kepler solver in NbodyGradient takes a different form that a typical Kepler solver - it's a universal solver, and it uses the difference in M and s over time steps, rather than their absolute values. So, these other approaches would need to be adapted and then tested.