siderust 0.9.0

High-precision astronomy and satellite mechanics 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

// SPDX-License-Identifier: AGPL-3.0-or-later
// Copyright (C) 2026 Vallés Puig, Ramon

//! Embedded Sun-Earth Lagrange Chebyshev kernels.
//!
//! The five SCK binary files (`l1.sck`..`l5.sck`) are embedded in
//! [`siderust_archive`] under `src/lagrange/vsop87/` and activated by the
//! `lagrange` feature of that crate (pulled in automatically when this crate
//! is built with `lagrange-centers`).
//!
//! The bytes are lazily parsed on first access so the parsing cost is paid only
//! when the Lagrange API is first called.

use std::sync::LazyLock;

use crate::archive::lagrange::vsop87 as sck_data;
use crate::formats::sck::parse_sck;

pub(crate) const NCOEFF: usize = 8;

static L1: LazyLock<Vec<f64>> = LazyLock::new(|| {
    parse_sck(sck_data::L1_BYTES)
        .expect("embedded L1 SCK must be valid")
        .records
});
static L2: LazyLock<Vec<f64>> = LazyLock::new(|| {
    parse_sck(sck_data::L2_BYTES)
        .expect("embedded L2 SCK must be valid")
        .records
});
static L3: LazyLock<Vec<f64>> = LazyLock::new(|| {
    parse_sck(sck_data::L3_BYTES)
        .expect("embedded L3 SCK must be valid")
        .records
});
static L4: LazyLock<Vec<f64>> = LazyLock::new(|| {
    parse_sck(sck_data::L4_BYTES)
        .expect("embedded L4 SCK must be valid")
        .records
});
static L5: LazyLock<Vec<f64>> = LazyLock::new(|| {
    parse_sck(sck_data::L5_BYTES)
        .expect("embedded L5 SCK must be valid")
        .records
});

/// SHA-256 of the original concatenated f64 payload from which the SCK files were derived.
pub(crate) const CHECKSUM: &str =
    "69c8eddcaeedc6f2e906a2ac5a4a2c2c7d97f491c08ce2b5ac6f4038e1966846";

#[must_use]
pub(crate) fn records_l1() -> &'static [f64] {
    (*L1).as_slice()
}

#[must_use]
pub(crate) fn records_l2() -> &'static [f64] {
    (*L2).as_slice()
}

#[must_use]
pub(crate) fn records_l3() -> &'static [f64] {
    (*L3).as_slice()
}

#[must_use]
pub(crate) fn records_l4() -> &'static [f64] {
    (*L4).as_slice()
}

#[must_use]
pub(crate) fn records_l5() -> &'static [f64] {
    (*L5).as_slice()
}