Re: [Orekit Users] DSSTPropagator isOsculating

[MAISONOBE Luc <luc.maisonobe@c-s.fr>]

Elisabet <elisabet.cid-borobia@thalesaleniaspace.com> a écrit :

> Hello Luc,
>
> In my case the earthDegree=36 and the earthOrder=4. It seems it throws the
> exception each time the earthDegree is different from the earthOrder. The
> DSSTPropagator, using mean parameters, can be used with these Earth’s degree
> and order. However, it throws the thread when using osculating parameters.

I finally managed to reproduce the error!
I first tried simply having different values for degree and order, but
still got a result.

Using (36, 4), I finally got the same exception as you!

I will look at this issue, stay tuned.

best regards,
Luc

> The thread:
> Exception in thread "main"
> org.apache.commons.math3.exception.MathArithmeticException: erreur
> arithmétique
> 	at
> org.apache.commons.math3.util.ArithmeticUtils.mulAndCheck(ArithmeticUtils.java:582)
> 	at
> org.apache.commons.math3.util.ArithmeticUtils.mulAndCheck(ArithmeticUtils.java:552)
> 	at
> org.apache.commons.math3.util.CombinatoricsUtils.binomialCoefficient(CombinatoricsUtils.java:126)
> 	at
> org.apache.commons.math3.analysis.polynomials.PolynomialsUtils.shift(PolynomialsUtils.java:329)
> 	at
> org.orekit.propagation.semianalytical.dsst.utilities.NewcombOperators
> $PolynomialsGenerator.shiftList(NewcombOperators.java:326)
> 	at
> org.orekit.propagation.semianalytical.dsst.utilities.NewcombOperators
> $PolynomialsGenerator.computeFor(NewcombOperators.java:216)
> 	at
> org.orekit.propagation.semianalytical.dsst.utilities.NewcombOperators
> $PolynomialsGenerator.getPolynomials(NewcombOperators.java:175)
> 	at
> org.orekit.propagation.semianalytical.dsst.utilities.NewcombOperators
> $PolynomialsGenerator.access$000(NewcombOperators.java:104)
> 	at
> org.orekit.propagation.semianalytical.dsst.utilities.NewcombOperators.getValue(NewcombOperators.java:86)
> 	at
> org.orekit.propagation.semianalytical.dsst.utilities.hansen.HansenTesseralLinear
> $HansenCoefficientsBySeries.generatePolynomial(HansenTesseralLinear.java:516)
> 	at
> org.orekit.propagation.semianalytical.dsst.utilities.hansen.HansenTesseralLinear
> $HansenCoefficientsBySeries.<init>(HansenTesseralLinear.java:472)
> 	at
> org.orekit.propagation.semianalytical.dsst.utilities.hansen.HansenTesseralLinear.<init>(HansenTesseralLinear.java:100)
> 	at
> org.orekit.propagation.semianalytical.dsst.forces.TesseralContribution.createHansenObjects(TesseralContribution.java:398)
> 	at
> org.orekit.propagation.semianalytical.dsst.forces.TesseralContribution.initialize(TesseralContribution.java:334)
> 	at
> org.orekit.propagation.semianalytical.dsst.forces.DSSTCentralBody.initialize(DSSTCentralBody.java:134)
> 	at
> org.orekit.propagation.semianalytical.dsst.DSSTPropagator.computeOsculatingState(DSSTPropagator.java:338)
>
> Best regards,
>
> Elisabet