siderust 0.6.0

High-precision astronomy and satellite mechanics in Rust.
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

// SPDX-License-Identifier: AGPL-3.0-or-later
// Copyright (C) 2026 Vallés Puig, Ramon

//! Frame-rotation provider for **TEME** (True Equator Mean Equinox).
//!
//! TEME is the output frame of SGP4/SDP4 propagators used with TLE data.
//! It shares the true equator of date but uses a simplified equinox
//! definition (the "equation of the equinoxes" shortcut).
//!
//! The rotation chain is:
//!
//! ```text
//! TEME → EquatorialTrueOfDate → (precession/nutation) → ICRS
//! ```
//!
//! The first step is Rz(Δψ cos εA), which removes the approximate
//! equinox offset.

use super::*;

impl FrameRotationProvider<TEME, EquatorialTrueOfDate> for () {
    #[inline]
    fn rotation<Eph, Eop: EopProvider, Nut>(
        jd: JulianDate,
        _ctx: &AstroContext<Eph, Eop, Nut>,
    ) -> Rotation3 {
        let nut = nutation::nutation_iau2000b(jd);
        Rotation3::rz(nut.dpsi * nut.mean_obliquity.cos())
    }
}

impl FrameRotationProvider<EquatorialTrueOfDate, TEME> for () {
    #[inline]
    fn rotation<Eph, Eop: EopProvider, Nut>(
        jd: JulianDate,
        ctx: &AstroContext<Eph, Eop, Nut>,
    ) -> Rotation3 {
        inverse_rotation::<EquatorialTrueOfDate, TEME, Eph, Eop, Nut>(jd, ctx)
    }
}

impl FrameRotationProvider<TEME, ICRS> for () {
    #[inline]
    fn rotation<Eph, Eop: EopProvider, Nut>(
        jd: JulianDate,
        ctx: &AstroContext<Eph, Eop, Nut>,
    ) -> Rotation3 {
        compose_rotation::<TEME, EquatorialTrueOfDate, ICRS, Eph, Eop, Nut>(jd, ctx)
    }
}

impl FrameRotationProvider<ICRS, TEME> for () {
    #[inline]
    fn rotation<Eph, Eop: EopProvider, Nut>(
        jd: JulianDate,
        ctx: &AstroContext<Eph, Eop, Nut>,
    ) -> Rotation3 {
        inverse_rotation::<ICRS, TEME, Eph, Eop, Nut>(jd, ctx)
    }
}