Bug when calling dh17! with BigFloat jacobian

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

Bug when calling dh17! with BigFloat jacobian #7

Closed ericagol closed 3 years ago

ericagol commented 6 years ago

test_dh17.jl has a commented out section in which it calls dh17! with BigFloat precision to compute the Jacobian. The returned jacobian has an error in which each data point is off by a factor of order unity; it looks like there is a problem with the passing of a BigFloat array to the subroutine and modifying it within. This was solved by adding an output to dh17!, but in general this isn't a workable solution.

ericagol commented 5 years ago

I just ran into this problem again!

langfzac commented 3 years ago

We no longer use dh17. Closing this