celestial-time 0.1.1-alpha.2

Pure Rust astronomical time scales
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

pub mod constants;
pub mod iau2000;
pub mod iau2006;

use crate::{TimeResult, TT};
use celestial_core::matrix::RotationMatrix3;

#[derive(Debug, Clone, PartialEq)]
pub struct PrecessionResult {
    pub bias_matrix: RotationMatrix3,
    pub precession_matrix: RotationMatrix3,
    pub bias_precession_matrix: RotationMatrix3,
    pub model: PrecessionModel,
}

#[derive(Debug, Clone, Copy, PartialEq, Eq)]
pub enum PrecessionModel {
    IAU1976,
    IAU2000,
    IAU2006,
}

pub trait PrecessionCalculator {
    fn precession_iau2000(&self) -> TimeResult<PrecessionResult>;

    fn precession_iau2006(&self) -> TimeResult<PrecessionResult>;

    fn precession(&self) -> TimeResult<PrecessionResult> {
        self.precession_iau2006()
    }

    fn bias_precession_matrix(&self) -> TimeResult<RotationMatrix3> {
        Ok(self.precession_iau2006()?.bias_precession_matrix)
    }
}

impl PrecessionCalculator for TT {
    fn precession_iau2000(&self) -> TimeResult<PrecessionResult> {
        iau2000::calculate(self)
    }

    fn precession_iau2006(&self) -> TimeResult<PrecessionResult> {
        iau2006::calculate(self)
    }
}