siderust 0.9.1

High-precision astronomy and satellite mechanics in Rust.
Documentation
// SPDX-License-Identifier: AGPL-3.0-or-later
// Copyright (C) 2026 Vallés Puig, Ramon

#![allow(missing_docs)]

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);
}

#[test]
fn position_display_includes_center_frame_and_components() {
    let v = siderust::coordinates::cartesian::position::ICRS::<AstronomicalUnit>::new(
        1.0 * AU,
        2.0 * AU,
        3.0 * AU,
    );

    let rendered = format!("{v}");

    assert!(rendered.contains("Center: Barycentric"));
    assert!(rendered.contains("Frame: ICRS"));
    assert!(rendered.contains("X: 1"));
}