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

use super::{NutationModel, NutationResult};
use crate::{TimeError, TimeResult, TT};
use celestial_core::nutation::NutationIAU2006A as CoreCalculator;

pub fn calculate(tt: &TT) -> TimeResult<NutationResult> {
    let _ = super::tt_to_centuries(tt)?;

    let jd = tt.to_julian_date();

    let calculator = CoreCalculator::new();
    let core_result = calculator.compute(jd.jd1(), jd.jd2()).map_err(|_| {
        TimeError::CalculationError("IAU 2006A nutation calculation failed".to_string())
    })?;

    Ok(NutationResult::new(core_result, NutationModel::IAU2006A))
}