def test_convert_between_coe_and_mee_is_transitive():
p = 11067.790 # u.km
ecc = 0.83285 # u.one
inc = np.deg2rad(87.87) # u.rad
raan = np.deg2rad(227.89) # u.rad
argp = np.deg2rad(53.38) # u.rad
nu = np.deg2rad(92.335) # u.rad
expected_res = (p, ecc, inc, raan, argp, nu)
res = mee2coe(*coe2mee(*expected_res))
assert_allclose(res, expected_res)
def to_equinoctial(self):
"""Converts to modified equinoctial elements representation.
"""
p, f, g, h, k, L = coe2mee(self.p.to(u.km).value,
self.ecc.value,
self.inc.to(u.rad).value,
self.raan.to(u.rad).value,
self.argp.to(u.rad).value,
self.nu.to(u.rad).value)
return equinoctial.ModifiedEquinoctialState(
self.attractor,
p * u.km,
f * u.rad,
g * u.rad,
h * u.rad,
k * u.rad,
L * u.rad)
def test_convert_between_coe_and_mee_is_transitive(classical):
res = mee2coe(*coe2mee(*classical))
assert_allclose(res, classical)
def test_convert_between_coe_and_mee_is_transitive(classical):
res = mee2coe(*coe2mee(*classical))
assert_allclose(res, classical)
def test_convert_between_coe_and_mee_is_transitive(expected_res):
res = mee2coe(*coe2mee(*expected_res))
assert_allclose(res, expected_res)