lox-frames 0.1.0-alpha.16

Reference frame transformations for the Lox ecosystem
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22

// SPDX-FileCopyrightText: 2024 Helge Eichhorn <git@helgeeichhorn.de>
//
// SPDX-License-Identifier: MPL-2.0

use std::f64::consts::{FRAC_PI_2, TAU};

use lox_core::glam::{DMat3, DVec3};

use crate::rotations::Rotation;

/// Builds the ICRF-to-IAU rotation from rotational elements and their rates.
pub fn icrf_to_iau(
    (right_ascension, declination, rotation_angle): (f64, f64, f64),
    (right_ascension_rate, declination_rate, rotation_rate): (f64, f64, f64),
) -> Rotation {
    let m1 = DMat3::from_rotation_z(-(right_ascension + FRAC_PI_2));
    let m2 = DMat3::from_rotation_x(-(FRAC_PI_2 - declination));
    let m3 = DMat3::from_rotation_z(-(rotation_angle % TAU));
    let m = m3 * m2 * m1;
    let v = DVec3::new(right_ascension_rate, -declination_rate, rotation_rate);
    Rotation::new(m).with_angular_velocity(v)
}