siderust 0.7.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
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
130

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

//! Orbit Models Example
//!
//! Demonstrates the four orbit propagation models: Keplerian, MeanMotion,
//! Conic, and Prepared.
//!
//! Run with: `cargo run --example 15_orbit_models`

use siderust::qtty::angular_rate::AngularRate;
use siderust::qtty::unit::{Day, Degree};
use siderust::qtty::AstronomicalUnits;
use siderust::time::JulianDate;
use siderust::{ConicKind, ConicOrbit, KeplerianOrbit, MeanMotionOrbit, PreparedOrbit};

fn main() {
    let jd = JulianDate::new(2458850.0);

    println!("=== Orbit Models ===\n");
    println!("Epoch: {:.1}\n", jd.jd_value());

    // 1. KeplerianOrbit
    let kepler = KeplerianOrbit::new(
        AstronomicalUnits::new(1.0),
        0.0167,
        siderust::qtty::Degrees::new(0.0),
        siderust::qtty::Degrees::new(0.0),
        siderust::qtty::Degrees::new(102.9),
        siderust::qtty::Degrees::new(100.0),
        JulianDate::new(2451545.0),
    );
    let kepler_pos = kepler.kepler_position(jd);
    println!("1. KeplerianOrbit");
    println!(
        "   heliocentric position = ({:.12}, {:.12}, {:.12}) AU",
        kepler_pos.x().value(),
        kepler_pos.y().value(),
        kepler_pos.z().value(),
    );
    println!("   radius = {:.12} AU\n", kepler_pos.distance().value());

    // 2. MeanMotionOrbit
    let mean_motion = MeanMotionOrbit::try_new(
        AstronomicalUnits::new(1.0),
        0.0167,
        siderust::qtty::Degrees::new(0.0),
        siderust::qtty::Degrees::new(0.0),
        siderust::qtty::Degrees::new(102.9),
        AngularRate::<Degree, Day>::new(0.9856),
        JulianDate::new(2451545.0),
    )
    .expect("valid mean-motion orbit");
    let mean_motion_pos = mean_motion.position_at(jd).expect("valid propagation");
    println!("2. MeanMotionOrbit");
    println!(
        "   heliocentric position = ({:.12}, {:.12}, {:.12}) AU",
        mean_motion_pos.x().value(),
        mean_motion_pos.y().value(),
        mean_motion_pos.z().value(),
    );
    println!(
        "   radius = {:.12} AU\n",
        mean_motion_pos.distance().value()
    );

    // 3. ConicOrbit
    let halley_like = ConicOrbit::try_new(
        AstronomicalUnits::new(0.586),
        0.967,
        siderust::qtty::Degrees::new(162.3),
        siderust::qtty::Degrees::new(58.4),
        siderust::qtty::Degrees::new(111.3),
        siderust::qtty::Degrees::new(38.4),
        JulianDate::new(2451545.0),
    )
    .expect("valid halley-like conic orbit");
    let hyperbolic = ConicOrbit::try_new(
        AstronomicalUnits::new(0.255),
        1.2,
        siderust::qtty::Degrees::new(122.7),
        siderust::qtty::Degrees::new(24.6),
        siderust::qtty::Degrees::new(241.8),
        siderust::qtty::Degrees::new(12.0),
        JulianDate::new(2451545.0),
    )
    .expect("valid hyperbolic conic orbit");
    let halley_pos = halley_like
        .position_at(jd)
        .expect("valid halley propagation");
    let hyperbolic_pos = hyperbolic
        .position_at(jd)
        .expect("valid hyperbolic propagation");
    println!("3. ConicOrbit");
    println!(
        "   halley-like kind = {}",
        if halley_like.kind() == ConicKind::Elliptic {
            "Elliptic"
        } else {
            "Hyperbolic"
        }
    );
    println!(
        "   hyperbolic kind  = {}",
        if hyperbolic.kind() == ConicKind::Elliptic {
            "Elliptic"
        } else {
            "Hyperbolic"
        }
    );
    println!(
        "   elliptic radius  = {:.12} AU | hyperbolic radius = {:.12} AU\n",
        halley_pos.distance().value(),
        hyperbolic_pos.distance().value(),
    );

    // 4. PreparedOrbit
    let prepared = PreparedOrbit::try_from(kepler).expect("valid prepared orbit");
    let prepared_pos = prepared.position_at(jd);
    let delta = ((kepler_pos.x().value() - prepared_pos.x().value()).powi(2)
        + (kepler_pos.y().value() - prepared_pos.y().value()).powi(2)
        + (kepler_pos.z().value() - prepared_pos.z().value()).powi(2))
    .sqrt();
    println!("4. PreparedOrbit");
    println!(
        "   prepared radius = {:.12} AU",
        prepared_pos.distance().value()
    );
    println!("   chord delta vs direct Keplerian = {:.3e}", delta);
}