phoboss 0.0.2

Phobos celestial simulation crate for the MilkyWay SolarSystem workspace
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38

#[test]
fn height_profile_min_max() {
    let hp = phoboss::surface::heightmaps::HeightProfile::new(vec![10.0, 20.0, 5.0]);
    let (min, max) = hp.min_max();
    assert!((min - 5.0).abs() < 1e-10);
    assert!((max - 20.0).abs() < 1e-10);
}

#[test]
fn roughness_index() {
    let hp = phoboss::surface::heightmaps::HeightProfile::new(vec![10.0, 20.0, 5.0, 25.0]);
    assert!(hp.roughness_index() > 0.0);
}

#[test]
fn great_circle_distance() {
    let a = phoboss::surface::mapping::PhobosCoord::new(0.0, 0.0, 0.0);
    let b = phoboss::surface::mapping::PhobosCoord::new(10.0, 10.0, 0.0);
    assert!(phoboss::surface::mapping::great_circle_distance_km(a, b) > 0.0);
}

#[test]
fn antipode() {
    let a = phoboss::surface::mapping::PhobosCoord::new(15.0, 30.0, 100.0);
    let b = phoboss::surface::mapping::antipode(a);
    assert!((b.latitude_deg + 15.0).abs() < 1e-10);
}

#[test]
fn zone_classification() {
    let z = phoboss::surface::zones::zone_from_lat_lon(0.0, 0.0);
    assert_eq!(
        z.zone,
        phoboss::surface::zones::PhobosZone::LeadingHemisphere
    );
    assert!((z.mean_albedo - 0.071).abs() < 1e-10);
    assert!((z.regolith_depth_m - 12.0).abs() < 1e-10);
}