#[test]
fn planet_binary_or_simulate() {
    let p = phoboss::interactions::marss::ensure_marss_binary_or_simulate();
    assert!(p.orbital_radius_km > 0.0);
}

#[test]
fn sub_planet_bulge() {
    let r = phoboss::physics::tides::TidalResponse::current();
    let b = phoboss::interactions::earth_tides::submars_bulge_m(r);
    assert!(b.is_finite());
}

#[test]
fn anti_planet_bulge() {
    let r = phoboss::physics::tides::TidalResponse::current();
    assert!(phoboss::interactions::earth_tides::antimars_bulge_m(r).is_finite());
}

#[test]
fn libration_factor() {
    let f = phoboss::interactions::orbital_dynamics::libration_factor(0.0);
    assert!((f - 1.0).abs() < 0.1);
}

#[test]
fn orbital_energy() {
    let orbit = phoboss::physics::orbit::PhobosOrbit::new();
    let e = phoboss::interactions::orbital_dynamics::orbital_energy_j_per_kg(
        orbit,
        orbit.semi_major_axis_m,
    );
    assert!(e.is_finite());
}
