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
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

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

//! Example: all currently supported frame conversions.
//!
//! This demonstrates every direct frame-rotation pair implemented in
//! `providers.rs`, plus identity rotations for each frame.

use qtty::*;
use siderust::coordinates::cartesian::Position;
use siderust::coordinates::centers::Barycentric;
use siderust::coordinates::frames::{
    EclipticMeanJ2000, EquatorialMeanJ2000, EquatorialMeanOfDate, EquatorialTrueOfDate,
    ReferenceFrame, ICRF, ICRS,
};
use siderust::coordinates::transform::{FrameRotationProvider, PositionAstroExt};
use siderust::time::JulianDate;

type C = Barycentric;
type U = AstronomicalUnit;

fn show_frame_conversion<F1, F2>(jd: &JulianDate, src: &Position<C, F1, U>)
where
    F1: ReferenceFrame,
    F2: ReferenceFrame,
    (): FrameRotationProvider<F1, F2>,
    (): FrameRotationProvider<F2, F1>,
{
    let out: Position<C, F2, U> = src.to_frame(jd);
    let back: Position<C, F1, U> = out.to_frame(jd);
    let err = (*src - back).magnitude();

    println!(
        "{:<24} -> {:<24} out=({:+.9})  roundtrip={:.3e}",
        F1::frame_name(),
        F2::frame_name(),
        out,
        err
    );
}

fn main() {
    let jd = JulianDate::new(2_460_000.5);
    println!("Frame conversion demo at {:.1}", jd);

    let p_icrs = Position::<C, ICRS, U>::new(0.30, -0.70, 0.64);
    let p_icrf: Position<C, ICRF, U> = p_icrs.to_frame(&jd);
    let p_ecl: Position<C, EclipticMeanJ2000, U> = p_icrs.to_frame(&jd);
    let p_eq_j2000: Position<C, EquatorialMeanJ2000, U> = p_icrs.to_frame(&jd);
    let p_eq_mod: Position<C, EquatorialMeanOfDate, U> = p_icrs.to_frame(&jd);
    let p_eq_tod: Position<C, EquatorialTrueOfDate, U> = p_icrs.to_frame(&jd);

    // Identity
    show_frame_conversion::<ICRS, ICRS>(&jd, &p_icrs);
    show_frame_conversion::<ICRF, ICRF>(&jd, &p_icrf);
    show_frame_conversion::<EclipticMeanJ2000, EclipticMeanJ2000>(&jd, &p_ecl);
    show_frame_conversion::<EquatorialMeanJ2000, EquatorialMeanJ2000>(&jd, &p_eq_j2000);
    show_frame_conversion::<EquatorialMeanOfDate, EquatorialMeanOfDate>(&jd, &p_eq_mod);
    show_frame_conversion::<EquatorialTrueOfDate, EquatorialTrueOfDate>(&jd, &p_eq_tod);

    // All direct non-identity provider pairs
    show_frame_conversion::<ICRS, EclipticMeanJ2000>(&jd, &p_icrs);
    show_frame_conversion::<EclipticMeanJ2000, ICRS>(&jd, &p_ecl);
    show_frame_conversion::<ICRS, EquatorialMeanJ2000>(&jd, &p_icrs);
    show_frame_conversion::<EquatorialMeanJ2000, ICRS>(&jd, &p_eq_j2000);
    show_frame_conversion::<EquatorialMeanJ2000, EclipticMeanJ2000>(&jd, &p_eq_j2000);
    show_frame_conversion::<EclipticMeanJ2000, EquatorialMeanJ2000>(&jd, &p_ecl);
    show_frame_conversion::<EquatorialMeanJ2000, EquatorialMeanOfDate>(&jd, &p_eq_j2000);
    show_frame_conversion::<EquatorialMeanOfDate, EquatorialMeanJ2000>(&jd, &p_eq_mod);
    show_frame_conversion::<EquatorialMeanOfDate, EquatorialTrueOfDate>(&jd, &p_eq_mod);
    show_frame_conversion::<EquatorialTrueOfDate, EquatorialMeanOfDate>(&jd, &p_eq_tod);
    show_frame_conversion::<EquatorialMeanJ2000, EquatorialTrueOfDate>(&jd, &p_eq_j2000);
    show_frame_conversion::<EquatorialTrueOfDate, EquatorialMeanJ2000>(&jd, &p_eq_tod);
    show_frame_conversion::<ICRS, EquatorialMeanOfDate>(&jd, &p_icrs);
    show_frame_conversion::<EquatorialMeanOfDate, ICRS>(&jd, &p_eq_mod);
    show_frame_conversion::<ICRS, EquatorialTrueOfDate>(&jd, &p_icrs);
    show_frame_conversion::<EquatorialTrueOfDate, ICRS>(&jd, &p_eq_tod);
    show_frame_conversion::<ICRF, ICRS>(&jd, &p_icrf);
    show_frame_conversion::<ICRS, ICRF>(&jd, &p_icrs);
    show_frame_conversion::<ICRF, EquatorialMeanJ2000>(&jd, &p_icrf);
    show_frame_conversion::<EquatorialMeanJ2000, ICRF>(&jd, &p_eq_j2000);
    show_frame_conversion::<ICRF, EclipticMeanJ2000>(&jd, &p_icrf);
    show_frame_conversion::<EclipticMeanJ2000, ICRF>(&jd, &p_ecl);
    show_frame_conversion::<ICRF, EquatorialMeanOfDate>(&jd, &p_icrf);
    show_frame_conversion::<EquatorialMeanOfDate, ICRF>(&jd, &p_eq_mod);
    show_frame_conversion::<ICRF, EquatorialTrueOfDate>(&jd, &p_icrf);
    show_frame_conversion::<EquatorialTrueOfDate, ICRF>(&jd, &p_eq_tod);
}