stem-branch 0.2.0

Native Rust port of stem-branch's solar ephemeris core (full VSOP87D Earth series + JPL DE441-fitted correction + IAU2000B nutation).
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

// @generated by scripts/gen-rust-solar-data.mjs from the TypeScript source.
// Do not edit by hand; re-run the generator to update.
// VSOP87 phase coefficients legitimately sit near math constants (e.g. π),
// which is not a sloppy approximation here.
#![allow(clippy::approx_constant)]

/// IAU2000B nutation-in-longitude terms (77 lunisolar).
/// Columns: `[l, l', F, D, Omega, dpsi_sin, dpsi_sin_T]`, the last two in
/// units of 0.1 microarcseconds. arg = sum(col[0..5] * delaunay), and
/// dpsi += (col[5] + col[6]*T) * sin(arg).
#[rustfmt::skip]
pub static NUT_COEFFS: &[[f64; 7]] = &[
    [0.0, 0.0, 0.0, 0.0, 1.0, -172064161.0, -174666.0],
    [0.0, 0.0, 2.0, -2.0, 2.0, -13170906.0, -1675.0],
    [0.0, 0.0, 2.0, 0.0, 2.0, -2276413.0, -234.0],
    [0.0, 0.0, 0.0, 0.0, 2.0, 2074554.0, 207.0],
    [0.0, 1.0, 0.0, 0.0, 0.0, 1475877.0, -3633.0],
    [0.0, 1.0, 2.0, -2.0, 2.0, -516821.0, 1226.0],
    [1.0, 0.0, 0.0, 0.0, 0.0, 711159.0, 73.0],
    [0.0, 0.0, 2.0, 0.0, 1.0, -387298.0, -367.0],
    [1.0, 0.0, 2.0, 0.0, 2.0, -301461.0, -36.0],
    [0.0, -1.0, 2.0, -2.0, 2.0, 215829.0, -494.0],
    [0.0, 0.0, 2.0, -2.0, 1.0, 128227.0, 137.0],
    [-1.0, 0.0, 2.0, 0.0, 2.0, 123457.0, 11.0],
    [-1.0, 0.0, 0.0, 2.0, 0.0, 156994.0, 10.0],
    [1.0, 0.0, 0.0, 0.0, 1.0, 63110.0, 63.0],
    [-1.0, 0.0, 0.0, 0.0, 1.0, -57976.0, -63.0],
    [-1.0, 0.0, 2.0, 2.0, 2.0, -59641.0, -11.0],
    [1.0, 0.0, 2.0, 0.0, 1.0, -51613.0, -42.0],
    [-2.0, 0.0, 2.0, 0.0, 1.0, 45893.0, 50.0],
    [0.0, 0.0, 0.0, 2.0, 0.0, 63384.0, 11.0],
    [0.0, 0.0, 2.0, 2.0, 2.0, -38571.0, -1.0],
    [0.0, -2.0, 2.0, -2.0, 2.0, 32481.0, 0.0],
    [-2.0, 0.0, 0.0, 2.0, 0.0, -47722.0, 0.0],
    [2.0, 0.0, 2.0, 0.0, 2.0, -31046.0, -1.0],
    [1.0, 0.0, 2.0, -2.0, 2.0, 28593.0, 0.0],
    [-1.0, 0.0, 2.0, 0.0, 1.0, 20441.0, 21.0],
    [2.0, 0.0, 0.0, 0.0, 0.0, 29243.0, 0.0],
    [0.0, 0.0, 2.0, 0.0, 0.0, 25887.0, 0.0],
    [0.0, 1.0, 0.0, 0.0, 1.0, -14053.0, -25.0],
    [-1.0, 0.0, 0.0, 2.0, 1.0, 15164.0, 10.0],
    [0.0, 2.0, 2.0, -2.0, 2.0, -15794.0, 72.0],
    [0.0, 0.0, -2.0, 2.0, 0.0, 21783.0, 0.0],
    [1.0, 0.0, 0.0, -2.0, 1.0, -12873.0, -10.0],
    [0.0, -1.0, 0.0, 0.0, 1.0, -12654.0, 11.0],
    [-1.0, 0.0, 2.0, 2.0, 1.0, -10204.0, 0.0],
    [0.0, 2.0, 0.0, 0.0, 0.0, 16707.0, -85.0],
    [1.0, 0.0, 2.0, 2.0, 2.0, -7691.0, 0.0],
    [-2.0, 0.0, 2.0, 0.0, 0.0, -11024.0, 0.0],
    [0.0, 1.0, 2.0, 0.0, 2.0, 7566.0, -21.0],
    [0.0, 0.0, 2.0, 2.0, 1.0, -6637.0, -11.0],
    [0.0, -1.0, 2.0, 0.0, 2.0, -7141.0, 21.0],
    [0.0, 0.0, 0.0, 2.0, 1.0, -6302.0, -11.0],
    [1.0, 0.0, 2.0, -2.0, 1.0, 5800.0, 10.0],
    [2.0, 0.0, 2.0, -2.0, 2.0, 6443.0, 0.0],
    [-2.0, 0.0, 0.0, 2.0, 1.0, -5774.0, -11.0],
    [2.0, 0.0, 2.0, 0.0, 1.0, -5350.0, 0.0],
    [0.0, -1.0, 2.0, -2.0, 1.0, -4752.0, -11.0],
    [0.0, 0.0, 0.0, -2.0, 1.0, -4940.0, -11.0],
    [-1.0, -1.0, 0.0, 2.0, 0.0, 7350.0, 0.0],
    [2.0, 0.0, 0.0, -2.0, 1.0, 4065.0, 0.0],
    [1.0, 0.0, 0.0, 2.0, 0.0, 6579.0, 0.0],
    [0.0, 1.0, 2.0, -2.0, 1.0, 3579.0, 0.0],
    [1.0, -1.0, 0.0, 0.0, 0.0, 4725.0, 0.0],
    [-2.0, 0.0, 2.0, 0.0, 2.0, -3075.0, 0.0],
    [3.0, 0.0, 2.0, 0.0, 2.0, -2904.0, 0.0],
    [0.0, -1.0, 0.0, 2.0, 0.0, 4348.0, 0.0],
    [1.0, -1.0, 2.0, 0.0, 2.0, -2878.0, 0.0],
    [0.0, 0.0, 0.0, 1.0, 0.0, -4230.0, 0.0],
    [-1.0, -1.0, 2.0, 2.0, 2.0, -2819.0, 0.0],
    [-1.0, 0.0, 2.0, 0.0, 0.0, -4056.0, 0.0],
    [0.0, -1.0, 2.0, 2.0, 2.0, -2647.0, 0.0],
    [-2.0, 0.0, 0.0, 0.0, 1.0, -2294.0, 0.0],
    [1.0, 1.0, 2.0, 0.0, 2.0, 2481.0, 0.0],
    [2.0, 0.0, 0.0, 0.0, 1.0, 2179.0, 0.0],
    [-1.0, 1.0, 0.0, 1.0, 0.0, 3276.0, 0.0],
    [1.0, 1.0, 0.0, 0.0, 0.0, -3389.0, 0.0],
    [1.0, 0.0, 2.0, 0.0, 0.0, 3339.0, 0.0],
    [-1.0, 0.0, 2.0, -2.0, 1.0, -1987.0, 0.0],
    [1.0, 0.0, 0.0, 0.0, 2.0, -1981.0, 0.0],
    [-1.0, 0.0, 0.0, 1.0, 0.0, 4026.0, 0.0],
    [0.0, 0.0, 2.0, 1.0, 2.0, 1660.0, 0.0],
    [-1.0, 0.0, 2.0, 4.0, 2.0, -1521.0, 0.0],
    [-1.0, 1.0, 0.0, 1.0, 1.0, 1314.0, 0.0],
    [0.0, -2.0, 2.0, -2.0, 1.0, -1283.0, 0.0],
    [1.0, 0.0, 2.0, 2.0, 1.0, -1331.0, 0.0],
    [-2.0, 0.0, 2.0, 2.0, 2.0, 1383.0, 0.0],
    [-1.0, 0.0, 0.0, 0.0, 2.0, 1405.0, 0.0],
    [1.0, 1.0, 2.0, -2.0, 2.0, 1290.0, 0.0]
];