[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Orekit Developers] Bug in eventdetector/commons-math?



Hi,

I'm experiencing some strange errors in my propagation with some eventdetectors, and thought I'd ask you for comments. I'm using orekit via the pythonwrapper so debugging is a bit complex..

It seems to occur only on certain orbits, but I'm just changing altitude and inclination (SSO) so seems more like some math issue. For some orbits it works.

The error stack is:
JavaError: org.apache.commons.math3.exception.NumberIsTooLargeException: endpoints do not specify an interval: [60, 59.951]
    Java stacktrace:
org.apache.commons.math3.exception.NumberIsTooLargeException: endpoints do not specify an interval: [60, 59.951]
at org.apache.commons.math3.analysis.solvers.UnivariateSolverUtils.verifyInterval(UnivariateSolverUtils.java:417)
at org.apache.commons.math3.analysis.solvers.UnivariateSolverUtils.verifySequence(UnivariateSolverUtils.java:435)
at org.apache.commons.math3.analysis.solvers.BaseAbstractUnivariateSolver.verifySequence(BaseAbstractUnivariateSolver.java:280)
at org.apache.commons.math3.analysis.solvers.BracketingNthOrderBrentSolver.doSolve(BracketingNthOrderBrentSolver.java:154)
at org.apache.commons.math3.analysis.solvers.BaseAbstractUnivariateSolver.solve(BaseAbstractUnivariateSolver.java:199)
at org.apache.commons.math3.analysis.solvers.BaseAbstractUnivariateSolver.solve(BaseAbstractUnivariateSolver.java:204)
at org.apache.commons.math3.analysis.solvers.BracketingNthOrderBrentSolver.solve(BracketingNthOrderBrentSolver.java:398)
at org.orekit.propagation.events.EventState.evaluateStep(EventState.java:231)
at org.orekit.propagation.analytical.AbstractAnalyticalPropagator.acceptStep(AbstractAnalyticalPropagator.java:279)
at org.orekit.propagation.analytical.AbstractAnalyticalPropagator.propagate(AbstractAnalyticalPropagator.java:167)
at org.orekit.propagation.AbstractPropagator.propagate(AbstractPropagator.java:257)


The call (from python) is:
finalState = prop.propagate(self.ore_t[0].shiftedBy(60 * 60 * 24.0 * self.duration))

prop.getInitialState().getOrbit()
<Orbit: keplerian parameters: {a: 7138142.0; e: 7.004623892026274E-7; i: 98.51876854105349; pa: 0.0; raan: 94.9995; v: 0.00700813981787154;}>

the propagator is a:
org.orekit.propagation.analytical.EcksteinHechlerPropagator

I'm running orekit 7.0 with apache math 3.4.2.

Do you have any feeling what could be the issue? I have the impression that this part of orekit would be very well used, so maybe it's something I'm doing but hard to see what.

Any inputs appriciated

Best Regards
/Petrus