Closed ericagol closed 5 years ago
It's not a bug; it's just a consequence of integrating over the contact points. I've now divided the integration up in test_integrate_transit_gradient.jl
so that we don't integrate over the contact points, and I find much better precision. However, the code is very slow.
There are features in the time-integrated derivatives which may be indicative of a bug.