I’ve implemented an orbit determination script with a Kalman filter and a given propogator used for this KF containing different forces. When running the whole script it runs just fine, except when I add atmospheric drag (I tried only with a simple exponential model).
Somehow it seems to be computing and then, after a given number of seconds (I tried different time step and the number of time steps changes but the propagation date is roughly the same), it completely freezes on the line used to propagate:
SpacecraftState nextStepState = propagator.propagate(date);
I have no idea where the error might come from.
During the KF process I have to constantly reset the initial state to then propagate to the next time step. I would have expected it to be slower due to air drag but it really switches from a “steady pace” propagation to no propagation at all. Maybe there is sthg I should be resetting?
Has anyone faced similar issues or does anyone know where the error might come from?
Thank you very much for your help.