spicekit 0.1.0__tar.gz → 0.2.1__tar.gz

{spicekit-0.1.0 → spicekit-0.2.1}/Cargo.lock

@@ -58,7 +58,7 @@ dependencies = [
  "proc-macro2",
  "quote",
  "regex",
- "rustc-hash",
+ "rustc-hash 1.1.0",
  "shlex",
  "which",
 ]
@@ -81,6 +81,12 @@ version = "3.20.2"
 source = "registry+https://github.com/rust-lang/crates.io-index"
 checksum = "5d20789868f4b01b2f2caec9f5c4e0213b41e3e5702a50157d699ae31ced2fcb"
 
+[[package]]
+name = "bytemuck"
+version = "1.25.0"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "c8efb64bd706a16a1bdde310ae86b351e4d21550d98d056f22f8a7f7a2183fec"
+
 [[package]]
 name = "bytes"
 version = "1.11.1"
@@ -89,9 +95,9 @@ checksum = "1e748733b7cbc798e1434b6ac524f0c1ff2ab456fe201501e6497c8417a4fc33"
 
 [[package]]
 name = "cc"
-version = "1.2.60"
+version = "1.2.61"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "43c5703da9466b66a946814e1adf53ea2c90f10063b86290cc9eb67ce3478a20"
+checksum = "d16d90359e986641506914ba71350897565610e87ce0ad9e6f28569db3dd5c6d"
 dependencies = [
  "find-msvc-tools",
  "shlex",
@@ -619,9 +625,9 @@ dependencies = [
 
 [[package]]
 name = "idna_adapter"
-version = "1.2.1"
+version = "1.2.2"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "3acae9609540aa318d1bc588455225fb2085b9ed0c4f6bd0d9d5bcd86f1a0344"
+checksum = "cb68373c0d6620ef8105e855e7745e18b0d00d3bdb07fb532e434244cdb9a714"
 dependencies = [
  "icu_normalizer",
  "icu_properties",
@@ -672,9 +678,9 @@ checksum = "8f42a60cbdf9a97f5d2305f08a87dc4e09308d1276d28c869c684d7777685682"
 
 [[package]]
 name = "js-sys"
-version = "0.3.95"
+version = "0.3.97"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "2964e92d1d9dc3364cae4d718d93f227e3abb088e747d92e0395bfdedf1c12ca"
+checksum = "a1840c94c045fbcf8ba2812c95db44499f7c64910a912551aaaa541decebcacf"
 dependencies = [
  "cfg-if",
  "futures-util",
@@ -702,9 +708,9 @@ checksum = "09edd9e8b54e49e587e4f6295a7d29c3ea94d469cb40ab8ca70b288248a81db2"
 
 [[package]]
 name = "libc"
-version = "0.2.185"
+version = "0.2.186"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "52ff2c0fe9bc6cb6b14a0592c2ff4fa9ceb83eea9db979b0487cd054946a2b8f"
+checksum = "68ab91017fe16c622486840e4c83c9a37afeff978bd239b5293d61ece587de66"
 
 [[package]]
 name = "libloading"
@@ -878,7 +884,7 @@ dependencies = [
  "num-integer",
  "num-traits",
  "pyo3",
- "rustc-hash",
+ "rustc-hash 1.1.0",
 ]
 
 [[package]]
@@ -1183,6 +1189,12 @@ version = "1.1.0"
 source = "registry+https://github.com/rust-lang/crates.io-index"
 checksum = "08d43f7aa6b08d49f382cde6a7982047c3426db949b1424bc4b7ec9ae12c6ce2"
 
+[[package]]
+name = "rustc-hash"
+version = "2.1.2"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "94300abf3f1ae2e2b8ffb7b58043de3d399c73fa6f4b73826402a5c457614dbe"
+
 [[package]]
 name = "rustix"
 version = "0.38.44"
@@ -1362,17 +1374,19 @@ dependencies = [
 
 [[package]]
 name = "spicekit"
-version = "0.1.0"
+version = "0.2.1"
 dependencies = [
+ "bytemuck",
  "memmap2",
  "rayon",
+ "rustc-hash 2.1.2",
  "tempfile",
  "thiserror",
 ]
 
 [[package]]
 name = "spicekit-bench"
-version = "0.1.0"
+version = "0.2.1"
 dependencies = [
  "cspice-sys",
  "spicekit",
@@ -1381,7 +1395,7 @@ dependencies = [
 
 [[package]]
 name = "spicekit-py"
-version = "0.1.0"
+version = "0.2.1"
 dependencies = [
  "ndarray",
  "numpy",
@@ -1659,9 +1673,9 @@ dependencies = [
 
 [[package]]
 name = "wasm-bindgen"
-version = "0.2.118"
+version = "0.2.120"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "0bf938a0bacb0469e83c1e148908bd7d5a6010354cf4fb73279b7447422e3a89"
+checksum = "df52b6d9b87e0c74c9edfa1eb2d9bf85e5d63515474513aa50fa181b3c4f5db1"
 dependencies = [
  "cfg-if",
  "once_cell",
@@ -1672,9 +1686,9 @@ dependencies = [
 
 [[package]]
 name = "wasm-bindgen-futures"
-version = "0.4.68"
+version = "0.4.70"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "f371d383f2fb139252e0bfac3b81b265689bf45b6874af544ffa4c975ac1ebf8"
+checksum = "af934872acec734c2d80e6617bbb5ff4f12b052dd8e6332b0817bce889516084"
 dependencies = [
  "js-sys",
  "wasm-bindgen",
@@ -1682,9 +1696,9 @@ dependencies = [
 
 [[package]]
 name = "wasm-bindgen-macro"
-version = "0.2.118"
+version = "0.2.120"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "eeff24f84126c0ec2db7a449f0c2ec963c6a49efe0698c4242929da037ca28ed"
+checksum = "78b1041f495fb322e64aca85f5756b2172e35cd459376e67f2a6c9dffcedb103"
 dependencies = [
  "quote",
  "wasm-bindgen-macro-support",
@@ -1692,9 +1706,9 @@ dependencies = [
 
 [[package]]
 name = "wasm-bindgen-macro-support"
-version = "0.2.118"
+version = "0.2.120"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "9d08065faf983b2b80a79fd87d8254c409281cf7de75fc4b773019824196c904"
+checksum = "9dcd0ff20416988a18ac686d4d4d0f6aae9ebf08a389ff5d29012b05af2a1b41"
 dependencies = [
  "bumpalo",
  "proc-macro2",
@@ -1705,9 +1719,9 @@ dependencies = [
 
 [[package]]
 name = "wasm-bindgen-shared"
-version = "0.2.118"
+version = "0.2.120"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "5fd04d9e306f1907bd13c6361b5c6bfc7b3b3c095ed3f8a9246390f8dbdee129"
+checksum = "49757b3c82ebf16c57d69365a142940b384176c24df52a087fb748e2085359ea"
 dependencies = [
  "unicode-ident",
 ]
@@ -1748,9 +1762,9 @@ dependencies = [
 
 [[package]]
 name = "web-sys"
-version = "0.3.95"
+version = "0.3.97"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "4f2dfbb17949fa2088e5d39408c48368947b86f7834484e87b73de55bc14d97d"
+checksum = "2eadbac71025cd7b0834f20d1fe8472e8495821b4e9801eb0a60bd1f19827602"
 dependencies = [
  "js-sys",
  "wasm-bindgen",

{spicekit-0.1.0 → spicekit-0.2.1}/Cargo.toml

@@ -3,7 +3,7 @@ members = ["crates/*"]
 resolver = "2"
 
 [workspace.package]
-version = "0.1.0"
+version = "0.2.1"
 edition = "2021"
 rust-version = "1.85"
 license = "MIT"

spicekit-0.2.1/LICENSE

@@ -0,0 +1,21 @@
+The MIT License (MIT)
+
+Copyright (c) 2026, Asteroid Institute
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

spicekit-0.2.1/LICENSE-NOTICES

@@ -0,0 +1,30 @@
+Third-Party Attributions
+========================
+
+src/naif_builtin_table.rs
+-------------------------
+
+The data in this file (the 692-entry built-in body-name ↔ NAIF ID map)
+was extracted from NASA's CSpice toolkit — specifically from the NPERM
+table in `zzidmap.c`. CSpice is distributed by the Navigation and
+Ancillary Information Facility (NAIF) at the Jet Propulsion Laboratory,
+California Institute of Technology.
+
+The NAIF toolkit is distributed without restriction; NAIF asks only for
+attribution when the toolkit or derived artifacts are redistributed.
+The canonical distribution statement (from the NAIF toolkit's
+distribution-terms documentation) is:
+
+    Copyright (C) [year], California Institute of Technology.
+    U.S. Government sponsorship acknowledged.
+
+spicekit acknowledges NAIF as the source of the body-name ↔ ID data
+embedded in `naif_builtin_table.rs`. Users who rely on that data in
+redistributed work should verify the current NAIF distribution-terms
+document, available with any CSpice toolkit release, for any updated
+requirements.
+
+spicekit is not a port of CSpice. All executable logic in this crate is
+an independent implementation written from the NAIF SPK / PCK / DAF file
+format specifications (which are published separately from the toolkit,
+in the NAIF "required reading" documents).

{spicekit-0.1.0 → spicekit-0.2.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: spicekit
-Version: 0.1.0
+Version: 0.2.1
 Classifier: Development Status :: 4 - Beta
 Classifier: Intended Audience :: Science/Research
 Classifier: License :: OSI Approved :: MIT License
@@ -19,6 +19,8 @@ Requires-Dist: naif-de440 ; extra == 'test'
 Requires-Dist: naif-earth-itrf93 ; extra == 'test'
 Requires-Dist: naif-eop-high-prec ; extra == 'test'
 Provides-Extra: test
+License-File: LICENSE
+License-File: LICENSE-NOTICES
 Summary: Python bindings for spicekit: pure-Rust NAIF kernel readers (DAF/SPK/PCK/text kernels).
 Author: Asteroid Institute
 Author-email: Alec Koumjian <alec@b612foundation.org>
@@ -33,7 +35,7 @@ Project-URL: Repository, https://github.com/B612-Asteroid-Institute/spicekit
 
 [![PyPI](https://img.shields.io/pypi/v/spicekit.svg)](https://pypi.org/project/spicekit/)
 [![CI](https://github.com/B612-Asteroid-Institute/spicekit/actions/workflows/ci.yml/badge.svg)](https://github.com/B612-Asteroid-Institute/spicekit/actions/workflows/ci.yml)
-[![License: MIT](https://img.shields.io/badge/License-MIT-yellow.svg)](../../LICENSE)
+[![License: MIT](https://img.shields.io/badge/License-MIT-yellow.svg)](./LICENSE)
 
 Python bindings for [`spicekit`](../spicekit), a pure-Rust reader for
 NASA/NAIF SPICE kernel formats (DAF, SPK, PCK, text kernels). Independent

{spicekit-0.1.0/crates/spicekit-py → spicekit-0.2.1}/README.md

@@ -2,7 +2,7 @@
 
 [![PyPI](https://img.shields.io/pypi/v/spicekit.svg)](https://pypi.org/project/spicekit/)
 [![CI](https://github.com/B612-Asteroid-Institute/spicekit/actions/workflows/ci.yml/badge.svg)](https://github.com/B612-Asteroid-Institute/spicekit/actions/workflows/ci.yml)
-[![License: MIT](https://img.shields.io/badge/License-MIT-yellow.svg)](../../LICENSE)
+[![License: MIT](https://img.shields.io/badge/License-MIT-yellow.svg)](./LICENSE)
 
 Python bindings for [`spicekit`](../spicekit), a pure-Rust reader for
 NASA/NAIF SPICE kernel formats (DAF, SPK, PCK, text kernels). Independent

{spicekit-0.1.0 → spicekit-0.2.1}/crates/spicekit/Cargo.toml

@@ -15,8 +15,10 @@ categories = ["science", "aerospace", "parser-implementations"]
 path = "src/lib.rs"
 
 [dependencies]
+bytemuck = "1.16"
 memmap2 = "0.9"
 rayon = "1.11"
+rustc-hash = "2.0"
 thiserror = "2.0"
 
 [dev-dependencies]

spicekit-0.2.1/crates/spicekit/LICENSE

@@ -0,0 +1,21 @@
+The MIT License (MIT)
+
+Copyright (c) 2026, Asteroid Institute
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

spicekit-0.2.1/crates/spicekit/LICENSE-NOTICES

@@ -0,0 +1,30 @@
+Third-Party Attributions
+========================
+
+src/naif_builtin_table.rs
+-------------------------
+
+The data in this file (the 692-entry built-in body-name ↔ NAIF ID map)
+was extracted from NASA's CSpice toolkit — specifically from the NPERM
+table in `zzidmap.c`. CSpice is distributed by the Navigation and
+Ancillary Information Facility (NAIF) at the Jet Propulsion Laboratory,
+California Institute of Technology.
+
+The NAIF toolkit is distributed without restriction; NAIF asks only for
+attribution when the toolkit or derived artifacts are redistributed.
+The canonical distribution statement (from the NAIF toolkit's
+distribution-terms documentation) is:
+
+    Copyright (C) [year], California Institute of Technology.
+    U.S. Government sponsorship acknowledged.
+
+spicekit acknowledges NAIF as the source of the body-name ↔ ID data
+embedded in `naif_builtin_table.rs`. Users who rely on that data in
+redistributed work should verify the current NAIF distribution-terms
+document, available with any CSpice toolkit release, for any updated
+requirements.
+
+spicekit is not a port of CSpice. All executable logic in this crate is
+an independent implementation written from the NAIF SPK / PCK / DAF file
+format specifications (which are published separately from the toolkit,
+in the NAIF "required reading" documents).

{spicekit-0.1.0 → spicekit-0.2.1}/crates/spicekit/README.md

@@ -77,11 +77,13 @@ import spicekit
 spk = spicekit.NaifSpk("/path/to/de440.bsp")
 ```
 
-## Parity testing
+## Parity testing and benchmarks
 
 `spicekit-bench` keeps an FFI-based CSpice reference installation and
 asserts that every spicekit numeric path bit-matches CSpice at a
-science-grade tolerance.
+science-grade tolerance. It also runs side-by-side latency
+microbenchmarks on matched inputs; the latest numbers (regenerated by
+CI on every push to `main`) live in [`BENCHMARKS.md`](./BENCHMARKS.md).
 
 See each module's rustdoc for details.
 

{spicekit-0.1.0 → spicekit-0.2.1}/crates/spicekit/src/daf.rs

@@ -1,11 +1,34 @@
 //! NAIF DAF (Double-precision Array File) container reader.
 //!
-//! All adam-core NAIF data — SPK, PCK, and CK — live inside DAF files. This
-//! module parses the file record and summary/name record chain, yielding
-//! typed `Summary` descriptors. Payload interpretation (e.g. SPK Type 2
+//! DAF is the shared on-disk container behind every binary SPICE
+//! kernel format — SPK, PCK, and CK — so this module is the entry
+//! point for the per-kernel-type parsers above it. It parses the
+//! file record and the summary/name record chain, yielding typed
+//! [`Summary`] descriptors; payload interpretation (e.g. SPK Type 2
 //! Chebyshev coefficients) is the caller's responsibility.
 //!
 //! Reference: NAIF "DAF Required Reading" (daf.req).
+//!
+//! # Host-endianness assumption
+//!
+//! [`DafFile::doubles_native`] reinterprets `LTL-IEEE` on-disk bytes as
+//! native `f64` with no byte-swap, which is only correct on
+//! little-endian hosts. We currently fail to compile on big-endian
+//! targets rather than silently produce byte-reversed values. The
+//! supported targets (Linux/macOS/Windows on `x86_64` and `aarch64`)
+//! are all little-endian; adding big-endian support requires either a
+//! byte-swap path here or routing the hot paths back through the
+//! existing endian-agnostic [`DafFile::read_doubles`].
+
+#[cfg(not(target_endian = "little"))]
+compile_error!(
+    "spicekit currently requires a little-endian host (DafFile::doubles_native \
+     reinterprets LTL-IEEE bytes as native f64 with no byte-swap). The supported \
+     targets — Linux/macOS/Windows on x86_64 and aarch64 — are all little-endian. \
+     Adding big-endian support requires either a byte-swap path in doubles_native \
+     or routing SPK/PCK hot paths back through read_doubles (which uses \
+     f64::from_le_bytes and is endian-agnostic)."
+);
 
 use std::path::Path;
 use std::sync::Arc;
@@ -67,8 +90,9 @@ pub struct Summary {
 impl DafFile {
     pub fn open<P: AsRef<Path>>(path: P) -> Result<Self, DafError> {
         let file = std::fs::File::open(path)?;
-        // SAFETY: the mapped file is read-only; adam-core does not modify
-        // SPK/PCK kernels after they've been delivered to the NAIF package.
+        // SAFETY: the mapping is read-only; SPK/PCK kernels are
+        // immutable artifacts — NAIF distributes them as fixed-content
+        // releases and consumers must not modify them in place.
         let mmap = unsafe { Mmap::map(&file)? };
         Self::from_mmap(mmap)
     }
@@ -234,6 +258,36 @@ impl DafFile {
         Ok(out)
     }
 
+    /// Zero-copy view of a double range as a `&[f64]` slice into the
+    /// memory-mapped file.
+    ///
+    /// This is the hot-path accessor used by SPK/PCK segment
+    /// evaluators: it does **no allocation** and **no per-double
+    /// byte-decoding loop** — the bytes are reinterpreted in place.
+    ///
+    /// Correctness rests on four invariants:
+    /// 1. **Host endianness**: the host is little-endian, so the native
+    ///    `f64` byte order matches the on-disk LTL-IEEE byte order.
+    ///    Enforced at compile time by the module-level
+    ///    `compile_error!` above; big-endian builds fail to compile.
+    /// 2. **File endianness**: we rejected non-`LTL-IEEE` files at
+    ///    [`DafFile::open`] (via the private `from_mmap` helper), so
+    ///    on-disk bytes are LE IEEE-754 f64.
+    /// 3. **Alignment**: DAF addresses are 1-indexed counts of 8-byte
+    ///    doubles; `mmap` returns a base pointer that is page-aligned
+    ///    (≥ 4096 B) so every DAF double address lands on an 8-byte
+    ///    boundary.
+    /// 4. **Length**: the byte length is `(end_addr - start_addr + 1) * 8`,
+    ///    exactly a multiple of 8.
+    ///
+    /// `bytemuck::cast_slice` validates (3) and (4) at runtime, so this
+    /// function is fully safe; on a malformed mmap it would panic before
+    /// returning bad data.
+    pub fn doubles_native(&self, start_addr: u32, end_addr: u32) -> Result<&[f64], DafError> {
+        let bytes = self.double_slice(start_addr, end_addr)?;
+        Ok(bytemuck::cast_slice(bytes))
+    }
+
     /// Zero-copy view of a double range as a &[u8] slice.
     pub fn double_slice(&self, start_addr: u32, end_addr: u32) -> Result<&[u8], DafError> {
         if start_addr == 0 || end_addr < start_addr {

spicekit-0.2.1/crates/spicekit/src/lib.rs

@@ -0,0 +1,33 @@
+//! Pure-Rust readers (and writer) for NASA/NAIF SPICE kernel formats.
+//!
+//! Implements the subset of SPICE this crate explicitly supports — DAF
+//! containers, SPK ephemerides, binary PCK rotations, and the
+//! body-name slice of text kernels — with memory-mapped I/O and no
+//! FFI. The DAF parser ([`daf`]) is the shared container reader; the
+//! per-kernel-type parsers ([`spk`], [`pck`], [`text_kernel`]) build
+//! on top of it. [`spk_writer`] emits SPK Type 3 and Type 9 segments.
+//!
+//! This crate does not link CSPICE. Bit-for-bit parity against CSPICE
+//! is asserted in the sibling `spicekit-bench` crate, which links
+//! `cspice-sys` behind a feature flag and exercises every numeric
+//! code path (`spkez`, `pxform`, `sxform`, `bodc2n`, `bodn2c`, plus
+//! the text-kernel precedence semantics) at machine-epsilon
+//! tolerance.
+
+pub mod daf;
+pub mod frame;
+pub mod naif_ids;
+pub mod pck;
+pub mod spk;
+pub mod spk_writer;
+pub mod text_kernel;
+
+pub use daf::{DafError, DafFile, Summary};
+pub use frame::{rotate_state, NaifFrame, OBLIQUITY_J2000_RAD};
+pub use naif_ids::{bodc2n, bodn2c, NaifIdError};
+pub use pck::{PckError, PckFile, PckSegment};
+pub use spk::{SpkError, SpkFile, SpkSegment};
+pub use spk_writer::{SpkWriter, SpkWriterError, Type3Record, Type3Segment, Type9Segment};
+pub use text_kernel::{
+    parse_body_bindings, parse_body_bindings_from_str, BodyBinding, TextKernelError,
+};

{spicekit-0.1.0 → spicekit-0.2.1}/crates/spicekit/src/naif_ids.rs

@@ -12,7 +12,7 @@
 //! per-backend body-name registry, which takes precedence over this
 //! table (last-loaded-wins, matching CSpice).
 //!
-//! Reference: https://naif.jpl.nasa.gov/pub/naif/toolkit_docs/C/req/naif_ids.html
+//! Reference: <https://naif.jpl.nasa.gov/pub/naif/toolkit_docs/C/req/naif_ids.html>
 
 use std::collections::HashMap;
 use std::sync::OnceLock;

{spicekit-0.1.0 → spicekit-0.2.1}/crates/spicekit/src/pck.rs

@@ -1,12 +1,15 @@
 //! Binary PCK (Planetary Constants Kernel) reader.
 //!
-//! Adam-core only needs time-varying planetary body orientations for the
-//! Earth-associated `ITRF93` frame. Binary PCK shares the DAF container
-//! with SPK, but the summary layout differs (NI=5 instead of 6, with
-//! `center` absent: segments describe a single body-fixed frame wrt a
-//! reference inertial frame). The payload is the standard Type 2
-//! Chebyshev record, but with 3 Euler-angle channels rather than 3
-//! position channels. The angles follow the NAIF 3-1-3 convention:
+//! This module reads time-varying planetary body orientations from
+//! binary PCKs. Scope is currently the Earth-associated `ITRF93` frame
+//! and any other body-fixed frames stored as PCK Type 2 Chebyshev
+//! Euler-angle records — the format used by NAIF's standard Earth EOP
+//! kernels. Binary PCK shares the DAF container with SPK, but the
+//! summary layout differs (NI=5 instead of 6, with `center` absent:
+//! segments describe a single body-fixed frame wrt a reference
+//! inertial frame). The payload is the standard Type 2 Chebyshev
+//! record, but with 3 Euler-angle channels rather than 3 position
+//! channels. The angles follow the NAIF 3-1-3 convention:
 //! `(t1, t2, t3)` such that `R_ref→body = Rz(t3) · Rx(t2) · Rz(t1)`.
 //! The `ref` frame is the segment's reference inertial frame — the
 //! standard Earth PCKs use `ECLIPJ2000` (NAIF frame ID 17), so
@@ -20,18 +23,19 @@
 //!
 //! Kernel precedence: three PCKs (`predict`, `historical`, `high_prec`)
 //! may cover overlapping epoch ranges. Last-loaded-wins matches
-//! CSPICE `furnsh` semantics: the caller is expected to load kernels in
-//! precedence order (least- to most-precise) and the file they opened
-//! most recently "wins" for overlapping coverage. This crate offers a
-//! per-file reader; multi-file composition is handled at the Python
-//! layer with a simple list-of-readers dispatcher.
+//! CSPICE `furnsh` semantics: the caller is expected to load kernels
+//! in precedence order (least- to most-precise) and the file they
+//! opened most recently "wins" for overlapping coverage. This crate
+//! offers a per-file reader; multi-file composition (a list-of-
+//! readers dispatcher) is the caller's responsibility — see the
+//! `spicekit-py` Python bindings for one such implementation.
 
 use std::path::Path;
 
 use thiserror::Error;
 
 use crate::daf::{DafError, DafFile};
-use crate::spk::cheby_val_and_deriv;
+use crate::spk::cheby3_val_and_deriv;
 
 #[derive(Debug, Error)]
 pub enum PckError {
@@ -110,7 +114,7 @@ impl PckType2 {
         let idx = raw_idx.clamp(0, self.n_records as isize - 1) as usize;
         let rec_start = self.start_addr + (idx * self.rsize) as u32;
         let rec_end = rec_start + self.rsize as u32 - 1;
-        let rec = self.file.read_doubles(rec_start, rec_end)?;
+        let rec = self.file.doubles_native(rec_start, rec_end)?;
 
         let mid = rec[0];
         let radius = rec[1];
@@ -124,12 +128,19 @@ impl PckType2 {
         let dec_c = &rec[2 + n..2 + 2 * n];
         let w_c = &rec[2 + 2 * n..2 + 3 * n];
 
-        let (ra, dra) = cheby_val_and_deriv(ra_c, s);
-        let (dec, ddec) = cheby_val_and_deriv(dec_c, s);
-        let (w, dw) = cheby_val_and_deriv(w_c, s);
-
+        // Three-channel evaluation with shared T_k(s) / dT_k/ds basis
+        // recurrence — mathematically equivalent to three independent
+        // calls but ~3x cheaper on the recurrence half of the work.
+        let (ang, dang) = cheby3_val_and_deriv(ra_c, dec_c, w_c, s);
         let inv_r = 1.0 / radius;
-        Ok([ra, dec, w, dra * inv_r, ddec * inv_r, dw * inv_r])
+        Ok([
+            ang[0],
+            ang[1],
+            ang[2],
+            dang[0] * inv_r,
+            dang[1] * inv_r,
+            dang[2] * inv_r,
+        ])
     }
 }