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

Add in compensated summation for Jacobian computation #27

Closed ericagol closed 3 years ago

ericagol commented 5 years ago

The compensated summation worked well in eliminating errors in the computation of t, x, and v in ah18 (as confirmed by stepsize_precision.jl). However, this is not yet added to the versions of these routines which compute the Jacobian, so I need to add this.