Hello,
On Orekit 12.0, I have a crash during the computation of the mean elements for the Brouwer-Lydanne propagator when I use a perfectly (or almost perfectly) circular osculating orbit as input.
During the iterations, the eccentricity becomes slightly negative and an exception is raised.
In the example, the exception is raised with eccentricity = 0 but with eccentricity = 1e-4, I still have the problem. With ecc=1e-3, it works.
Is there a problem in the way I try to use it ?
TestBrouwer.java (1.1 KB)
Exception in thread "main" org.orekit.errors.OrekitException: paramètre invalide eccentricity: -0,001 n'est pas dans l'intervalle [0, ∞]
at org.orekit.orbits.KeplerianOrbit.checkParameterRangeInclusive(KeplerianOrbit.java:1366)
at org.orekit.orbits.KeplerianOrbit.<init>(KeplerianOrbit.java:186)
at org.orekit.propagation.analytical.BrouwerLyddanePropagator$BLModel.propagateParameters(BrouwerLyddanePropagator.java:1438)
at org.orekit.propagation.analytical.BrouwerLyddanePropagator.computeMeanParameters(BrouwerLyddanePropagator.java:740)
at org.orekit.propagation.analytical.BrouwerLyddanePropagator.resetInitialState(BrouwerLyddanePropagator.java:682)
at org.orekit.propagation.analytical.BrouwerLyddanePropagator.<init>(BrouwerLyddanePropagator.java:533)
at org.orekit.propagation.analytical.BrouwerLyddanePropagator.computeMeanOrbit(BrouwerLyddanePropagator.java:644)
at org.orekit.propagation.analytical.BrouwerLyddanePropagator.computeMeanOrbit(BrouwerLyddanePropagator.java:600)
at org.orekit.propagation.analytical.BrouwerLyddanePropagator.computeMeanOrbit(BrouwerLyddanePropagator.java:568)
at TestBrouwer.main(TestBrouwer.java:23)
Thank you.
Christophe
