sofars 0.6.2

Pure Rust implementation of the IAU SOFA library
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
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139

use criterion::{Criterion, criterion_group, criterion_main};
use sofars::erst;
use std::hint::black_box;

fn bench_erst(c: &mut Criterion) {
    let mut group = c.benchmark_group("erst");

    // Common test parameters
    let dj1 = 2400000.5;
    let dj2 = 53736.0;
    let epsa = 0.4090789763356509900;
    let dpsi = -0.9630909107115582393e-5;

    group.bench_function("ee00", |b| {
        b.iter(|| {
            erst::ee00(
                black_box(dj1),
                black_box(dj2),
                black_box(epsa),
                black_box(dpsi),
            )
        })
    });

    group.bench_function("ee00a", |b| {
        b.iter(|| erst::ee00a(black_box(dj1), black_box(dj2)))
    });

    group.bench_function("ee00b", |b| {
        b.iter(|| erst::ee00b(black_box(dj1), black_box(dj2)))
    });

    group.bench_function("ee06a", |b| {
        b.iter(|| erst::ee06a(black_box(dj1), black_box(dj2)))
    });

    group.bench_function("eect00", |b| {
        b.iter(|| erst::eect00(black_box(dj1), black_box(dj2)))
    });

    group.bench_function("eqeq94", |b| {
        b.iter(|| erst::eqeq94(black_box(dj1), black_box(41234.0)))
    });

    group.bench_function("era00", |b| {
        b.iter(|| erst::era00(black_box(dj1), black_box(54388.0)))
    });

    group.bench_function("gmst00", |b| {
        b.iter(|| {
            erst::gmst00(
                black_box(dj1),
                black_box(dj2),
                black_box(dj1),
                black_box(dj2),
            )
        })
    });

    group.bench_function("gmst06", |b| {
        b.iter(|| {
            erst::gmst06(
                black_box(dj1),
                black_box(dj2),
                black_box(dj1),
                black_box(dj2),
            )
        })
    });

    group.bench_function("gmst82", |b| {
        b.iter(|| erst::gmst82(black_box(dj1), black_box(dj2)))
    });

    group.bench_function("gst00a", |b| {
        b.iter(|| {
            erst::gst00a(
                black_box(dj1),
                black_box(dj2),
                black_box(dj1),
                black_box(dj2),
            )
        })
    });

    group.bench_function("gst00b", |b| {
        b.iter(|| erst::gst00b(black_box(dj1), black_box(dj2)))
    });

    let rnpb = [
        [
            0.9999989440476103608,
            -0.1332881761240011518e-2,
            -0.5790767434730085097e-3,
        ],
        [
            0.1332858254308954453e-2,
            0.9999991109044505944,
            -0.4097782710401555759e-4,
        ],
        [
            0.5791308472168153320e-3,
            0.4020595661593994396e-4,
            0.9999998314954572365,
        ],
    ];

    group.bench_function("gst06", |b| {
        b.iter(|| {
            erst::gst06(
                black_box(dj1),
                black_box(dj2),
                black_box(dj1),
                black_box(dj2),
                black_box(&rnpb),
            )
        })
    });

    group.bench_function("gst06a", |b| {
        b.iter(|| {
            erst::gst06a(
                black_box(dj1),
                black_box(dj2),
                black_box(dj1),
                black_box(dj2),
            )
        })
    });

    group.bench_function("gst94", |b| {
        b.iter(|| erst::gst94(black_box(dj1), black_box(dj2)))
    });

    group.finish();
}

criterion_group!(benches, bench_erst);
criterion_main!(benches);