The ah18 algorithm, in which the negative drift and kepler step are carried out together (but in different orders) seems to give about a factor of ~4 better precision for the derivatives in the test_ttv_elements.jl routine.
Switching to using this algorithm would require:
[x] Adding in the phic corrector so that this can be used with kickfast.
[x] Adding a description of the algorithm to the paper.
The
ah18
algorithm, in which the negative drift and kepler step are carried out together (but in different orders) seems to give about a factor of ~4 better precision for the derivatives in thetest_ttv_elements.jl
routine.Switching to using this algorithm would require:
phic
corrector so that this can be used with kickfast.