def test_kepler_2d_inverse(): t = 23 p = kepler.Kepler2DParameters(a=2, pb=3, eps1=0.2, eps2=0.1, t0=t - 1) xyv, _ = kepler.kepler_2d(p, t) m = kepler.mass(p.a, p.pb) p2 = kepler.inverse_kepler_2d(xyv, m, t) assert_allclose(p, p2, atol=1e-8)
def test_kepler_2d_inverse(): t = 23 p = kepler.Kepler2DParameters(a=2, pb=3, eps1=0.2, eps2=0.1, t0=t-1) xyv, _ = kepler.kepler_2d(p,t) m = kepler.mass(p.a, p.pb) p2 = kepler.inverse_kepler_2d(xyv, m, t) assert_allclose(p,p2,atol=1e-8)
def test_mass_solar(): a = u.au.to(ls) pb = u.year.to(u.day) m = kepler.mass(a, pb) assert_allclose(m, 1, rtol=1e-4) # Earth+Sun mass
def test_mass_solar(): a = u.au.to(ls) pb = u.year.to(u.day) m = kepler.mass(a,pb) assert_allclose(m,1,rtol=1e-4) # Earth+Sun mass