sgp4 0.8.2

A pure Rust implementation of the SGP4 algorithm for satellite propagation
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
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66

#[path = "../test_cases.rs"]
mod test_cases;
use test_cases::*;

#[test]
fn propagate() -> anyhow::Result<()> {
    let test_cases: TestCases = toml::from_str(include_str!("../test_cases.toml")).unwrap();
    for test_case in test_cases.list.iter() {
        let constants =
            sgp4::Constants::from_elements_afspc_compatibility_mode(&sgp4::Elements::from_tle(
                None,
                test_case.line1.as_bytes(),
                test_case.line2.as_bytes(),
            )?)?;
        for state in &test_case.states {
            match state {
                State::Ok {
                    time,
                    position,
                    velocity,
                    ..
                } => {
                    let prediction = constants.propagate_afspc_compatibility_mode(*time)?;
                    for index in 0..3 {
                        assert!((position[index] - prediction.position[index]).abs() < 1.0e-6);
                        assert!((velocity[index] - prediction.velocity[index]).abs() < 1.0e-9);
                    }
                }
                State::Err { time, error } => {
                    let prediction = constants.propagate_afspc_compatibility_mode(*time);
                    if let Err(prediction_error) = prediction {
                        match prediction_error {
                            sgp4::Error::OutOfRangePerturbedEccentricity { eccentricity: _, t } => {
                                if error != "diverging perturbed eccentricity" {
                                    panic!(
                                        "bad error type (expected {}, got {:?})",
                                        error, prediction_error
                                    );
                                } else {
                                    assert_eq!(*time, t);
                                }
                            }
                            sgp4::Error::NegativeSemiLatusRectum { t } => {
                                if error != "negative semi-latus rectum" {
                                    panic!(
                                        "bad error type (expected {}, got {:?})",
                                        error, prediction_error
                                    );
                                } else {
                                    assert_eq!(*time, t);
                                }
                            }
                            _ => panic!(
                                "bad error type (expected {}, got {:?})",
                                error, prediction_error
                            ),
                        }
                    } else {
                        panic!("propagation should have returned an error");
                    }
                }
            }
        }
    }
    Ok(())
}