satkit 0.22.0

Satellite Toolkit
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
67
68
69
70
71

//! Integration test for the earth_orientation_params bytes-loading entry
//! points. Refresh-in-place semantics: `init_from_bytes` always succeeds
//! and replaces.

use satkit::earth_orientation_params as eop;
use satkit::utils::datadir;

fn eop_bytes() -> Option<Vec<u8>> {
    let path = datadir().ok()?.join("EOP-All.csv");
    std::fs::read(path).ok()
}

#[test]
fn init_from_bytes_replaces_and_query_works() {
    let Some(bytes) = eop_bytes() else {
        eprintln!(
            "skipping: EOP-All.csv not available in datadir(); \
             run `python -m satkit.utils.update_datafiles` or set SATKIT_DATA"
        );
        return;
    };

    // 1. First init populates the singleton.
    eop::init_from_bytes(&bytes).expect("init_from_bytes should succeed on first call");

    // 2. Query against the just-installed records (known truth value from
    //    the in-source test).
    let v = eop::eop_from_mjd_utc(59464.00).expect("EOP for MJD 59464");
    let truth: [f64; 4] = [-0.1145667, 0.241155, 0.317274, -0.0002255];
    for (a, b) in v.iter().zip(truth.iter()) {
        assert!(
            ((a - b) / b).abs() < 1.0e-3,
            "EOP mismatch after bytes init: got {a}, expected {b}"
        );
    }

    // 3. Second init succeeds (refresh-in-place semantics) and the query
    //    still works.
    eop::init_from_bytes(&bytes)
        .expect("second init_from_bytes should succeed (refreshable subsystem)");
    let v2 = eop::eop_from_mjd_utc(59464.00).expect("EOP for MJD 59464 after reload");
    assert_eq!(v, v2);

    // 4. Coverage / status on the real table.
    let cov = eop::coverage().expect("coverage after init");
    assert!(cov.first < cov.last_observed && cov.last_observed <= cov.last);
    assert_eq!(eop::status(&cov.first), eop::EopStatus::Observed);
    assert_eq!(
        eop::status(&(cov.last + satkit::Duration::from_days(1.0))),
        eop::EopStatus::Extrapolated
    );

    // 5. An empty table (header only) counts as *not loaded*: no coverage,
    //    queries return None, and the propagator refuses to build its
    //    ephemeris table rather than run with zero EOP.
    let header = bytes.split(|&b| b == b'\n').next().unwrap().to_vec();
    eop::init_from_bytes(&header).expect("header-only init parses");
    assert!(eop::coverage().is_none());
    assert_eq!(eop::status(&cov.first), eop::EopStatus::NotLoaded);
    assert!(eop::eop_from_mjd_utc(59464.00).is_none());
    let t0 = satkit::Instant::from_datetime(2024, 1, 1, 0, 0, 0.0).unwrap();
    let t1 = t0 + satkit::Duration::from_hours(1.0);
    match satkit::orbitprop::Precomputed::new(&t0, &t1) {
        Err(satkit::orbitprop::Error::EopUnavailable) => {}
        other => panic!("expected EopUnavailable with an empty EOP table, got {other:?}"),
    }

    // 6. Restore the real table so nothing else in this process is affected.
    eop::init_from_bytes(&bytes).expect("restore");
    assert!(eop::coverage().is_some());
}