astronomy 0.1.5

A library for astronomical calculations in Rust
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
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115

use rust_decimal::Decimal;
use rust_decimal::prelude::{FromPrimitive, ToPrimitive};
use std::fmt::{self, Display};
use std::ops::Add;

#[derive(Debug, Clone, Copy, PartialEq, Eq, PartialOrd, Ord)]
pub struct Time(Decimal);

impl Time {
    pub fn from_gps_seconds(seconds: f64) -> Self {
        Time(Decimal::from_f64(seconds).expect("Failed to convert f64 to Decimal for Time."))
    }

    pub fn as_gps_seconds_f64(&self) -> f64 {
        self.0
            .to_f64()
            .expect("Failed to convert Decimal to f64 for Time.")
    }
}

// Implementing the Display trait for Time
// This is useful for printing Time instances.
impl Display for Time {
    fn fmt(&self, f: &mut fmt::Formatter<'_>) -> fmt::Result {
        write!(f, "{}", self.0)
    }
}
impl Add for Time {
    type Output = Self;

    fn add(self, other: Self) -> Self::Output {
        Time(self.0 + other.0)
    }
}

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

    #[test]
    fn test_time_from_gps_seconds() {
        let time = Time::from_gps_seconds(123456.789);
        assert_eq!(time.as_gps_seconds_f64(), 123456.789);
    }

    #[test]
    fn test_time_display() {
        let time = Time::from_gps_seconds(98765.4321);
        assert_eq!(format!("{}", time), "98765.4321");
    }

    #[test]
    fn test_time_equality() {
        let time1 = Time::from_gps_seconds(100.0);
        let time2 = Time::from_gps_seconds(100.0);
        let time3 = Time::from_gps_seconds(200.0);

        assert_eq!(time1, time2);
        assert_ne!(time1, time3);
        assert!(time1 < time3);
    }

    #[test]
    fn test_time_addition() {
        let time1 = Time::from_gps_seconds(50.0);
        let time2 = Time::from_gps_seconds(25.0);
        let result = time1 + time2;

        assert_eq!(result.as_gps_seconds_f64(), 75.0);
    }

    #[test]
    fn test_time_ordering() {
        let time1 = Time::from_gps_seconds(300.0);
        let time2 = Time::from_gps_seconds(400.0);
        let time3 = Time::from_gps_seconds(300.0);

        assert!(time1 < time2);
        assert!(time1 <= time3);
        assert!(time2 > time1);
        assert!(time2 >= time3);
    }

    #[test]
    fn test_time_conversion() {
        let time = Time::from_gps_seconds(123456.789);
        let gps_seconds = time.as_gps_seconds_f64();
        assert_eq!(gps_seconds, 123456.789);
    }

    #[test]
    fn test_time_from_decimal() {
        let decimal_time = Decimal::from_f64(123456.789).expect("Failed to create Decimal");
        let time = Time(decimal_time);
        assert_eq!(time.as_gps_seconds_f64(), 123456.789);
    }

    #[test]
    #[allow(clippy::excessive_precision)]
    fn test_time_creation_and_conversion() {
        let gps_time = 1126259446.0;
        let time_obj = Time::from_gps_seconds(gps_time);
        assert_eq!(time_obj.as_gps_seconds_f64(), gps_time);

        let precise_gps_time = 123456789.123456789;
        let time_obj_precise = Time::from_gps_seconds(precise_gps_time);
        // Due to f64 precision, we might not get the exact value back
        // For Decimal, it shoiuld be exact

        assert_eq!(
            time_obj_precise.0,
            Decimal::from_f64(precise_gps_time).unwrap()
        );
    }
}