Skip to content
Snippets Groups Projects
Commit aa5cde7b authored by Clément Jonglez's avatar Clément Jonglez
Browse files

12.2 compat: Brouwer Lyddane test: use same tolerances as in java test

parent edaafa75
No related branches found
No related tags found
1 merge request!21Orekit 12.2
......@@ -67,17 +67,17 @@ class BrouwerLyddanePropagatorTest(unittest.TestCase):
# positions velocity and semi major axis match perfectly
self.assertAlmostEqual(0.0, Vector3D.distance(initialOrbit.getPVCoordinates().getPosition(),
finalOrbit.getPVCoordinates().getPosition()), delta=1.0e-8)
finalOrbit.getPVCoordinates().getPosition()), delta=2.4e-8)
self.assertAlmostEqual(0.0, Vector3D.distance(initialOrbit.getPVCoordinates().getVelocity(),
finalOrbit.getPVCoordinates().getVelocity()), delta= 1.0e-11)
self.assertAlmostEqual(0.0, finalOrbit.getA() - initialOrbit.getA(), delta=0.0)
finalOrbit.getPVCoordinates().getVelocity()), delta= 1.9e-11)
self.assertAlmostEqual(0.0, finalOrbit.getA() - initialOrbit.getA(), delta=9.4e-10)
def test_compareToNumericalPropagation(self):
inertialFrame = FramesFactory.getEME2000()
initDate = AbsoluteDate.J2000_EPOCH.shiftedBy(584.)
timeshift = 60000.0
timeshift = 60000.0
# Initial orbit
a = 24396159.0 # semi major axis in meters
......@@ -127,15 +127,15 @@ class BrouwerLyddanePropagatorTest(unittest.TestCase):
BLFinalState = BLextrapolator.propagate(initDate.shiftedBy(timeshift))
BLOrbit = OrbitType.KEPLERIAN.convertType(BLFinalState.getOrbit())
self.assertAlmostEqual(NumOrbit.getA(), BLOrbit.getA(), delta=0.5)
self.assertAlmostEqual(NumOrbit.getE(), BLOrbit.getE(), delta=0.00000028)
self.assertAlmostEqual(NumOrbit.getI(), BLOrbit.getI(), delta=0.000004)
self.assertAlmostEqual(NumOrbit.getA(), BLOrbit.getA(), delta=0.175)
self.assertAlmostEqual(NumOrbit.getE(), BLOrbit.getE(), delta=3.2e-6)
self.assertAlmostEqual(NumOrbit.getI(), BLOrbit.getI(), delta=6.9e-8)
self.assertAlmostEqual(MathUtils.normalizeAngle(NumOrbit.getPerigeeArgument(), FastMath.PI),
MathUtils.normalizeAngle(BLOrbit.getPerigeeArgument(), FastMath.PI), delta=0.119)
MathUtils.normalizeAngle(BLOrbit.getPerigeeArgument(), FastMath.PI), delta=0.0053)
self.assertAlmostEqual(MathUtils.normalizeAngle(NumOrbit.getRightAscensionOfAscendingNode(), FastMath.PI),
MathUtils.normalizeAngle(BLOrbit.getRightAscensionOfAscendingNode(), FastMath.PI), delta=0.000072)
MathUtils.normalizeAngle(BLOrbit.getRightAscensionOfAscendingNode(), FastMath.PI), delta=1.2e-6)
self.assertAlmostEqual(MathUtils.normalizeAngle(NumOrbit.getTrueAnomaly(), FastMath.PI),
MathUtils.normalizeAngle(BLOrbit.getTrueAnomaly(), FastMath.PI), delta=0.12)
MathUtils.normalizeAngle(BLOrbit.getTrueAnomaly(), FastMath.PI), delta=0.0052)
if __name__ == '__main__':
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment