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

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

use siderust::coordinates::cartesian::{Displacement, Position};
use siderust::coordinates::centers::Heliocentric;
use siderust::coordinates::frames::EclipticMeanJ2000;
use siderust::qtty::*;

type DispAu = Displacement<EclipticMeanJ2000, Au>;
type PosAu = Position<Heliocentric, EclipticMeanJ2000, Au>;

#[test]
fn displacement_basic_operations() {
    let v1 = DispAu::new(1.0, 0.0, 0.0);
    let v2 = DispAu::new(0.0, 1.0, 0.0);

    let sum = v1 + v2;
    assert_eq!(sum.x().value(), 1.0);
    assert_eq!(sum.y().value(), 1.0);

    let diff = v1 - v2;
    assert_eq!(diff.x().value(), 1.0);
    assert_eq!(diff.y().value(), -1.0);

    let dist = v1.magnitude();
    assert!((dist.value() - 1.0).abs() < 1e-12);
}

#[test]
fn position_affine_operations() {
    let p1 = PosAu::new(1.0, 0.0, 0.0);
    let p2 = PosAu::new(0.0, 1.0, 0.0);

    // Position - Position = Displacement
    let disp = p1 - p2;
    assert_eq!(disp.x().value(), 1.0);
    assert_eq!(disp.y().value(), -1.0);

    let dist = p1.distance();
    assert!((dist.value() - 1.0).abs() < 1e-12);

    let between = p1.distance_to(&p2);
    assert!((between.value() - 2.0_f64.sqrt()).abs() < 1e-12);
}