siderust 0.8.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

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

use siderust::qtty::*;

#[test]
fn hour_angles_from_hms() {
    let ha = HourAngles::from_hms(5, 30, 0.0);
    assert!((ha.value() - 5.5).abs() < 1e-12);

    let neg = HourAngles::from_hms(-5, 15, 30.0);
    assert!((neg.value() + 5.258333333333333).abs() < 1e-12);
}

#[test]
fn degrees_from_dms() {
    let deg = Degrees::from_dms(-33, 52, 0.0);
    assert!((deg.value() + 33.86666666666667).abs() < 1e-12);

    let with_sign = Degrees::from_dms_sign(-1, 33, 52, 0.0);
    assert!((with_sign.value() + 33.86666666666667).abs() < 1e-12);
}

#[test]
fn meters_kilometers_conversion() {
    let m = Meters::new(5000.0);
    let km = m.to::<Km>();
    assert!((km.value() - 5.0).abs() < 1e-12);

    let m_back = km.to::<Meter>();
    assert!((m_back.value() - 5000.0).abs() < 1e-12);
}

#[test]
fn astronomical_unit_lightyear_roundtrip() {
    let au = AstronomicalUnits::new(2.0);
    let ly: LightYears = au.into();
    let au_back: AstronomicalUnits = ly.into();
    assert!((au_back.value() - au.value()).abs() < 1e-10);
}

#[test]
fn quantity_arithmetic_and_simplify() {
    let distance = Meters::new(120.0);
    let extra = Meters::new(5.0);
    let sum = distance + extra;
    assert_eq!(sum.value(), 125.0);

    let diff = sum - extra;
    assert!((diff.value() - distance.value()).abs() < 1e-12);

    let time = Seconds::new(10.0);
    let velocity = distance / time; // Quantity<Per<Meter, Second>>
    let dist_back = velocity * time;
    assert!((dist_back.value() - distance.value()).abs() < 1e-12);

    let time_back = distance / velocity;
    assert!((time_back.value() - time.value()).abs() < 1e-12);

    let unitless_ratio = distance / distance;
    assert!((unitless_ratio.asin() - std::f64::consts::FRAC_PI_2).abs() < 1e-12);
}

#[test]
fn unitless_from_same_unit_division_and_display() {
    let a = Meters::new(84.0);
    let b = Meters::new(2.0);
    let ratio: f64 = a / b;
    assert!((ratio - 42.0).abs() < 1e-12);
    assert_eq!(ratio, 42.0);
}

#[test]
fn erase_unit_raw_and_display() {
    let meters = Meters::new(42.0);
    let raw: f64 = meters.erase_unit_raw();
    assert_eq!(raw, 42.0);
}