Non convergence when estimating Propagation Parameters in version 13.1

Orekit usage
12

Hi Pol,

I’ve found and fixed the bug in the patch-13.1.2 branch.
Basically the partial derivatives were wrongly updated in the following situation:

  • there are model parameters in the propagated Jacobian matrix
  • the mass is included in the STM
  • state-dependent, derivatives-resetting event occurs (typically eclipses)

So that happened for you when using batch least squares on model parameters with solar radiation pressure taken into account.
I’m a bit surprised that this bug didn’t cause failed tests in the estimation package, but hey, it’s also my bad for not better validating things. The original mistake was to forget that stored partial derivatives are always wrt Cartesian coordinates.

Cheers,
Romain.

1 Like