rfa 0.5.9

A port ERFA to 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

mod common;
#[cfg(test)]

mod tests{
    use super::common::*;
    use erfa_sys::*;
    use rfa::all::{ursa_epv00, ursa_moon98, ursa_plan94};
    #[test]
    pub fn epv00_test(){
        let dj1 = ERFA_DJ00; let dj2 = 365.0;
        for i in 0..20  {
            let dj_new = dj2 + (i*100) as f64;
            let mut erfa_pvh = [[0.0; 3];2]; let mut erfa_pvb = [[0.0; 3];2];

            unsafe{eraEpv00(dj1, dj_new, erfa_pvh.as_mut_ptr(), erfa_pvb.as_mut_ptr())};

            let mut rfa_pvh = [[0.0; 3];2]; let mut rfa_pvb = [[0.0; 3];2];
            ursa_epv00(dj1, dj_new, &mut rfa_pvh, &mut rfa_pvb);

            for i in 0..rfa_pvh.len() {
                for j in 0..rfa_pvh[i].len() {
                    assert!(rel_equal(rfa_pvh[i][j], erfa_pvh[i][j], 1E-14) );
                    assert!(rel_equal(rfa_pvb[i][j], erfa_pvb[i][j], 1E-14) );

                }
            }
        }
    }

    #[test]
    pub fn moon98_test(){
        let dj1 = ERFA_DJ00; let dj2 = 365.0;
        for i in 0..20  {
            let dj_new = dj2 + (i*100) as f64;
            let mut erfa_pv = [[0.0; 3];2]; 

            unsafe{eraMoon98(dj1, dj_new, erfa_pv.as_mut_ptr())};

            let mut rfa_pv = [[0.0; 3];2]; 
            ursa_moon98(dj1, dj_new, &mut rfa_pv);

            for i in 0..rfa_pv.len() {
                for j in 0..rfa_pv[i].len() {
                    assert!(rel_equal(rfa_pv[i][j], erfa_pv[i][j], 1E-13) );
                }
            }
        }
    }

    #[test]
    pub fn plan94_test(){
        let dj1 = ERFA_DJ00; let dj2 = 365.0;
        for i in 0..20  {
            let dj_new = dj2 + (i*100) as f64;
            let mut erfa_pv = [[0.0; 3];2]; 

            for np in 1..9{
            unsafe{eraPlan94(dj1, dj_new, np, erfa_pv.as_mut_ptr())};

            let mut rfa_pv = [[0.0; 3];2]; 
            ursa_plan94(dj1, dj_new, np, &mut rfa_pv);

            for i in 0..rfa_pv.len() {
                for j in 0..rfa_pv[i].len() {
                    assert!(rel_equal(rfa_pv[i][j], erfa_pv[i][j], 1E-13) );
                }
            }
            }
        }
    }
}