siderust 0.18.1

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
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129

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

//! Center-shift providers for the three **standard** reference centers:
//! [`Barycentric`], [`Heliocentric`] and [`Geocentric`].
//!
//! These form the hub of the center-shift graph. Every other center
//! (planetary, selenocentric, …) composes through [`Barycentric`], so
//! only the three pairwise directions need explicit implementations here.
//!
//! # Data source
//!
//! The Sun and Earth barycentric positions come from the
//! selected [`Ephemeris`] trait implementation (VSOP87 by default).
//! Both return `Position<Barycentric, EclipticMeanJ2000, AU>`.

use super::*;

// ---------------------------------------------------------------------------
// Heliocentric ↔ Barycentric
// ---------------------------------------------------------------------------

/// Heliocentric → Barycentric shift.
///
/// Returns the Sun's barycentric position converted from the ecliptic to
/// the target frame `F`.
impl<F> CenterShiftProvider<Heliocentric, Barycentric, F> for ()
where
    F: affn::ReferenceFrame,
    (): FrameRotationProvider<EclipticMeanJ2000, F>,
{
    #[inline]
    fn try_shift<Eph: Ephemeris, Eop: EopProvider, Nut: NutationModel>(
        jd: JulianDate,
        ctx: &AstroContext<Eph, Eop>,
    ) -> Result<AuShift, CenterShiftError> {
        Ok(rotate_shift_from_ecliptic::<_, F, Eph, Eop, Nut>(
            Eph::try_sun_barycentric(jd)?,
            jd,
            ctx,
        ))
    }
}

impl<F> CenterShiftProvider<Barycentric, Heliocentric, F> for ()
where
    F: affn::ReferenceFrame,
    (): FrameRotationProvider<EclipticMeanJ2000, F>,
{
    #[inline]
    fn try_shift<Eph: Ephemeris, Eop: EopProvider, Nut: NutationModel>(
        jd: JulianDate,
        ctx: &AstroContext<Eph, Eop>,
    ) -> Result<AuShift, CenterShiftError> {
        try_inverse_shift::<Barycentric, Heliocentric, F, Eph, Eop, Nut>(jd, ctx)
    }
}

// ---------------------------------------------------------------------------
// Geocentric ↔ Barycentric
// ---------------------------------------------------------------------------

/// Geocentric → Barycentric shift.
///
/// Returns the Earth's barycentric position converted from the ecliptic to
/// the target frame `F`.
impl<F> CenterShiftProvider<Geocentric, Barycentric, F> for ()
where
    F: affn::ReferenceFrame,
    (): FrameRotationProvider<EclipticMeanJ2000, F>,
{
    #[inline]
    fn try_shift<Eph: Ephemeris, Eop: EopProvider, Nut: NutationModel>(
        jd: JulianDate,
        ctx: &AstroContext<Eph, Eop>,
    ) -> Result<AuShift, CenterShiftError> {
        Ok(rotate_shift_from_ecliptic::<_, F, Eph, Eop, Nut>(
            Eph::try_earth_barycentric(jd)?,
            jd,
            ctx,
        ))
    }
}

impl<F> CenterShiftProvider<Barycentric, Geocentric, F> for ()
where
    F: affn::ReferenceFrame,
    (): FrameRotationProvider<EclipticMeanJ2000, F>,
{
    #[inline]
    fn try_shift<Eph: Ephemeris, Eop: EopProvider, Nut: NutationModel>(
        jd: JulianDate,
        ctx: &AstroContext<Eph, Eop>,
    ) -> Result<AuShift, CenterShiftError> {
        try_inverse_shift::<Barycentric, Geocentric, F, Eph, Eop, Nut>(jd, ctx)
    }
}

// ---------------------------------------------------------------------------
// Heliocentric ↔ Geocentric (composed through Barycentric)
// ---------------------------------------------------------------------------

impl<F> CenterShiftProvider<Heliocentric, Geocentric, F> for ()
where
    F: affn::ReferenceFrame,
    (): FrameRotationProvider<EclipticMeanJ2000, F>,
{
    #[inline]
    fn try_shift<Eph: Ephemeris, Eop: EopProvider, Nut: NutationModel>(
        jd: JulianDate,
        ctx: &AstroContext<Eph, Eop>,
    ) -> Result<AuShift, CenterShiftError> {
        try_compose_shift::<Heliocentric, Barycentric, Geocentric, F, Eph, Eop, Nut>(jd, ctx)
    }
}

impl<F> CenterShiftProvider<Geocentric, Heliocentric, F> for ()
where
    F: affn::ReferenceFrame,
    (): FrameRotationProvider<EclipticMeanJ2000, F>,
{
    #[inline]
    fn try_shift<Eph: Ephemeris, Eop: EopProvider, Nut: NutationModel>(
        jd: JulianDate,
        ctx: &AstroContext<Eph, Eop>,
    ) -> Result<AuShift, CenterShiftError> {
        try_inverse_shift::<Geocentric, Heliocentric, F, Eph, Eop, Nut>(jd, ctx)
    }
}