def test_compute_missing_body_hill_radius_raises_error(missing_body):
with pytest.raises(RuntimeError) as excinfo:
hill_radius(missing_body)
assert (
"To compute the semimajor axis or eccentricity for Moon and Pluto use the JPL ephemeris"
in excinfo.exconly()
)
def test_hill_radius(body, expected_r_SOI):
if expected_r_SOI is not None:
expected_r_SOI = expected_r_SOI * u.m
r_SOI = hill_radius(body, e=0 * u.one)
assert_quantity_allclose(r_SOI, expected_r_SOI, rtol=1e-1)
def test_hill_radius(body, expected_r_SOI):
if expected_r_SOI is not None:
expected_r_SOI = expected_r_SOI * u.m
r_SOI = hill_radius(body, e=0 * u.one)
assert_quantity_allclose(r_SOI, expected_r_SOI, rtol=1e-1)
def test_hill_radius_given_a():
parent = Body(None, 1 * u.km**3 / u.s**2, "Parent")
body = Body(parent, 1 * u.km**3 / u.s**2, "Body")
r_SOI = hill_radius(body, 1 * u.km, 0.25 * u.one)
expected_r_SOI = 520.02096 * u.m
assert_quantity_allclose(r_SOI, expected_r_SOI, rtol=1e-8)
def test_hill_radius(body, expected_r_SOI):
r_SOI = hill_radius(body)
assert_quantity_allclose(r_SOI, expected_r_SOI, rtol=1e-1)
def test_hill_radius_given_a():
parent = Body(None, 1 * u.km ** 3 / u.s ** 2, "Parent")
body = Body(parent, 1 * u.km ** 3 / u.s ** 2, "Body")
r_SOI = hill_radius(body, 1 * u.km, 0.25 * u.one)
expected_r_SOI = 520.02096 * u.m
assert_quantity_allclose(r_SOI, expected_r_SOI, rtol=1e-8)