Right now, the jacobian of the full lightcurve doesn't match the finite differences at/around the points of contact. The lower level code works (ie. passes tests), which seems to indicate either 1) the finite differences are inaccurate or 2) there's a strange bug somewhere in the top-level function compute_lightcurve!.
Right now, the jacobian of the full lightcurve doesn't match the finite differences at/around the points of contact. The lower level code works (ie. passes tests), which seems to indicate either 1) the finite differences are inaccurate or 2) there's a strange bug somewhere in the top-level function
compute_lightcurve!
.