geodate 0.5.0

Geodate computes geocentric expressions of points in time using a natural lunisolar calendar with metric time based on decimal fractions of the mean solar day.
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
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76

use core::f64::consts::PI;
use core::ops::Div;
#[cfg(not(feature = "std"))]
use num_traits::Float;

// FIXME: Use builtin rust to_radians()
pub fn rad(num: f64) -> f64 {
    num * PI / 180.0
}

// FIXME: Use builtin rust to_degrees()
pub fn deg(num: f64) -> f64 {
    num * 180.0 / PI
}

pub fn cos_deg(num: f64) -> f64 {
    rad(num).cos()
}

pub fn sin_deg(num: f64) -> f64 {
    rad(num).sin()
}

pub fn tan_deg(num: f64) -> f64 {
    rad(num).tan()
}

pub fn acos_deg(num: f64) -> f64 {
    deg(num.acos())
}

pub fn asin_deg(num: f64) -> f64 {
    deg(num.asin())
}

pub fn atan2_deg(x: f64, y: f64) -> f64 {
    deg(x.atan2(y))
}

pub fn dec_deg(d: f64, m: f64, s: f64) -> f64 {
    d + m / 60.0 + s / 3600.0
}

pub fn modulo(x: f64, y: f64) -> f64 {
    (y + x % y) % y
}

pub fn interpolate(y1: f64, y2: f64, y3: f64, n: f64) -> f64 {
    // From "Astronomical Algorithms" by Jean Meeus
    // Formula 3.3
    let a = y2 - y1; // First difference
    let b = y3 - y2; // First difference
    let c = b - a;   // Second difference

    y2 + (n.div(2.0)) * (a + b + n * c)
}

#[cfg(test)]
mod tests {
    use super::*;

    #[test]
    fn cos_deg_test() {
        assert_eq!(-1.0, cos_deg(180.0));
    }

    #[test]
    fn acos_deg_test() {
        assert_eq!(180.0, acos_deg(-1.0));
    }

    #[test]
    fn dec_deg_test() {
        assert_eq!(1.3958333333333333, dec_deg(1.0, 23.0, 45.0));
    }
}