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

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

use siderust::coordinates::centers::Geodetic;
use siderust::coordinates::frames::ECEF;
use siderust::coordinates::spherical::{direction, position};
use siderust::qtty::*;

const EPS: f64 = 1e-6;

#[test]
fn ecef_normalization_and_altitude() {
    // Input angles outside their canonical ranges
    // new(lon, lat) normalizes lat via wrap_quarter_fold to [-90, 90]
    // and lon via normalize to [0, 360)
    // Note: IAU convention - lon first, lat second
    let dir = direction::EcefDir::new(190.0 * DEG, 95.0 * DEG);
    // 95° lat wraps over the pole via wrap_quarter_fold: lat → 85°,
    // longitude flips by 180° (190° + 180° = 370° → 10°).
    assert!(
        (dir.lat().value() - 85.0).abs() < EPS,
        "lat mismatch: {}",
        dir.lat().value()
    );
    assert!(
        (dir.lon().value() - 10.0).abs() < EPS,
        "lon mismatch: {}",
        dir.lon().value()
    );

    // Geodetic::<ECEF>::new() normalises longitude to [-180, 180).
    // 190° wraps to -170°.
    let coord = Geodetic::<ECEF>::new(190.0 * DEG, 85.0 * DEG, 10_000.0 * M);
    assert!(
        (coord.lat.value() - 85.0).abs() < EPS,
        "geodetic lat mismatch: {}",
        coord.lat.value()
    );
    assert!(
        (coord.lon.value() - (-170.0)).abs() < EPS,
        "geodetic lon mismatch: {}",
        coord.lon.value()
    );
}

#[test]
fn ecliptic_normalization() {
    // Direction::new(lon, lat) normalizes both
    let dir = direction::EclipticMeanJ2000::new(120.0 * DEG, -45.0 * DEG);
    // After normalization: lat = -45° (in [-90, 90]), lon = 120° (in [0, 360))
    assert!(
        (dir.lon().value() - 120.0).abs() < EPS,
        "lon mismatch: {}",
        dir.lon().value()
    );
    assert!(
        (dir.lat().value() - (-45.0)).abs() < EPS,
        "lat mismatch: {}",
        dir.lat().value()
    );

    // Position::new(lon, lat, distance) also normalizes
    let pos =
        position::EclipticMeanJ2000::<AstronomicalUnit>::new(120.0 * DEG, -45.0 * DEG, 2.0 * AU);
    assert!(
        (pos.lon().value() - 120.0).abs() < EPS,
        "pos lon mismatch: {}",
        pos.lon().value()
    );
    assert!(
        (pos.lat().value() - (-45.0)).abs() < EPS,
        "pos lat mismatch: {}",
        pos.lat().value()
    );
    assert!((pos.distance - 2.0 * AU).abs() < EPS * AU);
}

#[test]
fn horizontal_normalization() {
    // Direction is now frame-only (no site parameter)
    // Note: new(alt, az) - IAU Alt-Az convention (altitude first)
    let dir = direction::Horizontal::new(120.0 * DEG, -30.0 * DEG);
    // Altitude 120° wraps over the zenith via wrap_quarter_fold: alt → 60°,
    // azimuth flips by 180° (-30° + 180° = 150°).
    assert!((dir.alt().value() - 60.0).abs() < EPS, "alt={}", dir.alt());
    assert!((dir.az().value() - 150.0).abs() < EPS, "az={}", dir.az());

    // Positions use new_with_site for Topocentric center
    let site = Geodetic::<ECEF>::default();
    // Note: new_with_site takes (site, alt, az, dist) - IAU Alt-Az convention (altitude first)
    let pos = siderust::coordinates::centers::Topocentric::horizontal(
        site,
        120.0 * DEG, // alt - wraps to 60
        -30.0 * DEG, // az - flips by 180° to 150
        2.0 * AU,
    );
    assert!(
        (pos.alt().value() - 60.0).abs() < EPS,
        "pos alt={}",
        pos.alt()
    );
    assert!(
        (pos.az().value() - 150.0).abs() < EPS,
        "pos az={}",
        pos.az()
    );
    assert!((pos.distance - 2.0 * AU).abs() < EPS * AU);
}