sidereon 3.0.4__tar.gz → 3.0.5__tar.gz

{sidereon-3.0.4 → sidereon-3.0.5}/.github/workflows/ci.yml

@@ -169,6 +169,8 @@ jobs:
         run: uv run --locked --extra test --group dev mypy
       - name: Build the extension (maturin develop)
         run: uv run --locked --extra test --group dev maturin develop --locked --features exact-cache-test-failpoints
+      - name: Release consumer contract
+        run: uv run --locked --extra test --group dev python scripts/smoke-release.py
       - name: stub accuracy
         env:
           SIDEREON_CORE_FIXTURES: ${{ github.workspace }}/.core-checkout/crates/sidereon-core/tests/fixtures

{sidereon-3.0.4 → sidereon-3.0.5}/.github/workflows/release.yml

@@ -153,10 +153,10 @@ jobs:
             done
             for source in mod.rs ocean.rs pole.rs; do
               unzip -p "$wheel" \
-                "third_party_source/sidereon-core-3.0.4/tides/$source" \
+                "third_party_source/sidereon-core-3.0.5/tides/$source" \
                 >"$RUNNER_TEMP/$source"
               cmp "$RUNNER_TEMP/$source" \
-                "third_party_source/sidereon-core-3.0.4/tides/$source"
+                "third_party_source/sidereon-core-3.0.5/tides/$source"
             done
           done
           notice="$RUNNER_TEMP/sdist.notices"
@@ -173,10 +173,10 @@ jobs:
           done
           for source in mod.rs ocean.rs pole.rs; do
             tar -xOf "${sdists[0]}" --wildcards \
-              "*/third_party_source/sidereon-core-3.0.4/tides/$source" \
+              "*/third_party_source/sidereon-core-3.0.5/tides/$source" \
               >"$RUNNER_TEMP/sdist-$source"
             cmp "$RUNNER_TEMP/sdist-$source" \
-              "third_party_source/sidereon-core-3.0.4/tides/$source"
+              "third_party_source/sidereon-core-3.0.5/tides/$source"
           done
       - name: Smoke-test installed x86_64 Linux wheel
         shell: bash

{sidereon-3.0.4 → sidereon-3.0.5}/CHANGELOG.md

@@ -4,6 +4,20 @@ All notable changes to the Sidereon Python interface are documented here.
 
 ## [Unreleased]
 
+## [3.0.5] - 2026-10-30
+
+### Added
+
+- Add opt-in SP3 column rounding writers that report maximum absolute SI changes; strict writers remain unchanged.
+- Preserve strictly writable SP3 columns for grid-aligned merge and alignment products.
+
+### Fixed
+
+- Normalize extreme-scale quaternions robustly while preserving ordinary quaternion rounding.
+- Reject zero probability values where probability inputs must be strictly between zero and one.
+- Describe the seven raw LNAV broadcast angular fields as semicircles or
+  semicircles/second, matching their unchanged transmitted-integer scaling.
+
 ## [3.0.4] - 2026-10-18
 
 ### Fixed

{sidereon-3.0.4 → sidereon-3.0.5}/Cargo.lock

@@ -611,9 +611,9 @@ checksum = "f8fadd59c855ef2080decdef8ff161eb6661b86933c9d82e5ba29dc602a55aba"
 
 [[package]]
 name = "sidereon"
-version = "3.0.4"
+version = "3.0.5"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "c042d74f42d0e2a9f8b2645109fd06f407e42ecab964fb91c492b18a2c034096"
+checksum = "fa3b2d18f88640b1bc0a86ed17bbacd9298eeea5b202174ec8da5ce7ae9a46e2"
 dependencies = [
  "crc32fast",
  "flate2",
@@ -622,9 +622,9 @@ dependencies = [
 
 [[package]]
 name = "sidereon-core"
-version = "3.0.4"
+version = "3.0.5"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "3634af272103625414a0e63bacbf2b76e181c33f8652690823acdfa0e25c938e"
+checksum = "3ba2c37c8b12bc20827e8e90971d635895e8a45fc715ad81c38d45cc3a4d3591"
 dependencies = [
  "approx",
  "cc",
@@ -646,7 +646,7 @@ dependencies = [
 
 [[package]]
 name = "sidereon-py"
-version = "3.0.4"
+version = "3.0.5"
 dependencies = [
  "nalgebra",
  "numpy",

{sidereon-3.0.4 → sidereon-3.0.5}/Cargo.toml

@@ -1,6 +1,6 @@
 [package]
 name = "sidereon-py"
-version = "3.0.4"
+version = "3.0.5"
 authors = []
 edition = "2021"
 description = "Python (PyO3) bindings over the sidereon ergonomic API: SP3 loading and SPP/RTK/PPP positioning solves"
@@ -25,8 +25,8 @@ exact-cache-test-failpoints = ["sidereon-core/exact-cache-test-failpoints"]
 # The thin ergonomic engine surface this binding marshals to. No modeling logic
 # lives here: every solve delegates to a `sidereon` entry point, which in turn
 # delegates to `sidereon-core`.
-sidereon = "=3.0.4"
-sidereon-core = { version = "=3.0.4", features = ["mmap"] }
+sidereon = { version = "=3.0.5" }
+sidereon-core = { version = "=3.0.5", features = ["mmap"] }
 # The standalone general-purpose trust-region least-squares solver: a data-driven
 # `least_squares` whose residual and Jacobian are evaluated entirely in Rust, so
 # no per-iteration Python callback is paid.

{sidereon-3.0.4 → sidereon-3.0.5}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: sidereon
-Version: 3.0.4
+Version: 3.0.5
 Classifier: Programming Language :: Rust
 Classifier: Programming Language :: Python :: 3
 Classifier: Topic :: Scientific/Engineering
@@ -280,7 +280,10 @@ The Python package mirrors the full breadth of the engine.
   defines, including exact validity seconds (see `docs/antex.md`).
 - **Data acquisition:** the `sidereon.data` module downloads and caches GNSS
   products (SP3 and IONEX from IGS/MGEX analysis centers, including merged
-  multi-center SP3) and DTED terrain tiles. Exact-product requests can select
+  multi-center SP3) and DTED terrain tiles. Merges from column-grid inputs
+  retain SP3-writable records; higher-precision inputs can still need the
+  opt-in rounding writer. Agreement metrics use combined values before final
+  column quantization. Exact-product requests can select
   direct archives, NASA CDDIS/Earthdata, local files, or caller-provided bytes
   without changing product identity. `ExactSp3Request` validates the declared
   cadence, regular epoch grid, exact half-open or inclusive coverage, producer,

{sidereon-3.0.4 → sidereon-3.0.5}/README.md

@@ -261,7 +261,10 @@ The Python package mirrors the full breadth of the engine.
   defines, including exact validity seconds (see `docs/antex.md`).
 - **Data acquisition:** the `sidereon.data` module downloads and caches GNSS
   products (SP3 and IONEX from IGS/MGEX analysis centers, including merged
-  multi-center SP3) and DTED terrain tiles. Exact-product requests can select
+  multi-center SP3) and DTED terrain tiles. Merges from column-grid inputs
+  retain SP3-writable records; higher-precision inputs can still need the
+  opt-in rounding writer. Agreement metrics use combined values before final
+  column quantization. Exact-product requests can select
   direct archives, NASA CDDIS/Earthdata, local files, or caller-provided bytes
   without changing product identity. `ExactSp3Request` validates the declared
   cadence, regular epoch grid, exact half-open or inclusive coverage, producer,

{sidereon-3.0.4 → sidereon-3.0.5}/THIRD-PARTY-NOTICES.md

@@ -330,9 +330,9 @@ endorsed by the IERS Conventions Center. The Rust source describes the origin,
 the renamed routines, and the differences from the Fortran implementation.
 The complete non-test tide source distributed in the compiled extension is
 included under
-`third_party_source/sidereon-core-3.0.4/tides/` in wheels and source
-distributions. These files match sidereon-core 3.0.4 at commit 5cfcdf929b57a5998e304df786735b171ba8afaf.
-Source: [sidereon-core tides at the pinned commit](https://github.com/neilberkman/sidereon/tree/5cfcdf929b57a5998e304df786735b171ba8afaf/crates/sidereon-core/src/tides).
+`third_party_source/sidereon-core-3.0.5/tides/` in wheels and source
+distributions. These files match sidereon-core 3.0.5 at commit da59d0eba8e273cc37821eac5535662e5d991126.
+Source: [sidereon-core tides at the pinned commit](https://github.com/neilberkman/sidereon/tree/da59d0eba8e273cc37821eac5535662e5d991126/crates/sidereon-core/src/tides).
 The authoritative original is available from:
 
 <https://iers-conventions.obspm.fr/content/chapter7/software/dehanttideinel/DEHANTTIDEINEL.F>

{sidereon-3.0.4 → sidereon-3.0.5}/docs/public-gnss-distribution-sources.md

@@ -278,11 +278,11 @@ except sidereon.Sp3WriteError as error:
 assert data.verify_merge_report(json.loads(persisted_report_json))
 ```
 
-The file helpers write SP3 text only when every value restates exactly. The
-default mean merge of two centers holds positions finer than the millimetre
-record column wherever the centers differ, and a clock moved onto the reference
-clock datum can be finer than its column; `Sp3WriteError` names the value
-instead of the helper rounding it.
+The merge path preserves grid-in/grid-out for record values derived from
+on-grid inputs, including means between column values. The file helpers write
+those merged products directly. Other products can still hold higher-precision
+values than the record columns; the strict writer names those with
+`Sp3WriteError`, while the explicit rounding writer reports its SI changes.
 
 `report.to_dict()` writes merged-SP3 report schema 3. A merge keeps a
 satellite whose orbit is missing but whose clock reached a consensus as a

{sidereon-3.0.4 → sidereon-3.0.5}/docs/sp3.md

@@ -25,6 +25,7 @@ parsing, rounding or repair of its own.
 | `clock_record(satellite, epoch_index)` | `Sp3ClockRecord`: a clock beside a missing orbit |
 | `clock_records_at(epoch_index)` | dict from satellite token to `Sp3ClockRecord` at one epoch, ascending |
 | `to_sp3_string()` | the SP3 text, or `Sp3WriteError` |
+| `to_sp3_string_with_rounding()` | `(text, Sp3RoundingReport)`, or `Sp3WriteError` |
 
 `state` and `clock_record` raise `IndexError` for an epoch index past the end
 and `KeyError` when the satellite has no such record at that epoch.
@@ -86,11 +87,20 @@ back unchanged; it can still be refused on a header base finer than its field,
 on header text wider than its columns, or on velocity records under a `P`
 header.
 
-A mean or median merge usually holds positions finer than the millimetre
-`F14.6` column, and a clock shifted onto the reference clock datum can hold a
-clock finer than its column, so writing such a merge raises
-`RecordValueNotRepresentable` naming the value. `data.write_sp3` and
-`data.fetch_merged_sp3_file` write nothing in that case.
+`to_sp3_string_with_rounding()` is an explicit alternative for products whose
+record values need the precision of their SP3 columns. It returns
+`(text, Sp3RoundingReport)`. The immutable report gives maximum absolute changes
+per field in SI units (`position_*_m`, `velocity_*_m_s`, `clock_s`, and
+`clock_rate_s_s`). It applies no rounding to header fields or epochs; refusals
+there still raise `Sp3WriteError`. `data.write_sp3_with_rounding()` writes the
+text atomically and returns `(path, report)`.
+
+The merge path preserves grid-in/grid-out for record values derived from
+on-grid inputs, including means between column values. Agreement and dispersion
+metrics are computed from combined values before the final column quantization;
+their arithmetic is unchanged. Other products can still hold record values
+finer than their columns; the strict writer refuses those, and the explicit
+rounding writer can state them while reporting the maximum SI change.
 
 ### Sp3WriteError
 

{sidereon-3.0.4 → sidereon-3.0.5}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "maturin"
 
 [project]
 name = "sidereon"
-version = "3.0.4"
+version = "3.0.5"
 description = "Python (PyO3) bindings over the sidereon engine: SP3 loading and SPP/RTK/PPP positioning solves"
 readme = "README.md"
 license = { text = "MIT" }
@@ -82,7 +82,7 @@ include = [
   "third_party_licenses/ERFA-BSD-3-Clause.txt",
   "third_party_licenses/IERS-CONVENTIONS-SOFTWARE-LICENSE.txt",
   "third_party_licenses/SciPy-BSD-3-Clause.txt",
-  "third_party_source/sidereon-core-3.0.4/tides/mod.rs",
-  "third_party_source/sidereon-core-3.0.4/tides/ocean.rs",
-  "third_party_source/sidereon-core-3.0.4/tides/pole.rs",
+  "third_party_source/sidereon-core-3.0.5/tides/mod.rs",
+  "third_party_source/sidereon-core-3.0.5/tides/ocean.rs",
+  "third_party_source/sidereon-core-3.0.5/tides/pole.rs",
 ]

{sidereon-3.0.4 → sidereon-3.0.5}/python/sidereon/__init__.pyi

@@ -509,6 +509,8 @@ class Sp3:
     def to_sp3_string(self) -> str:
         """Serialize to SP3 text; raises `Sp3WriteError` rather than change a
         value to fit its columns."""
+    def to_sp3_string_with_rounding(self) -> tuple[str, Sp3RoundingReport]:
+        """Write rounded record values and return their maximum SI deltas."""
     def precise_ephemeris_samples(self) -> list[PreciseEphemerisSample]: ...
     def precise_ephemeris_accuracy_samples(
         self,
@@ -518,6 +520,26 @@ class Sp3:
     ) -> bytes: ...
     def __repr__(self) -> str: ...
 
+class Sp3RoundingReport:
+    """Immutable maximum absolute SI changes made by the SP3 rounding writer."""
+    @property
+    def position_x_m(self) -> float: ...
+    @property
+    def position_y_m(self) -> float: ...
+    @property
+    def position_z_m(self) -> float: ...
+    @property
+    def velocity_x_m_s(self) -> float: ...
+    @property
+    def velocity_y_m_s(self) -> float: ...
+    @property
+    def velocity_z_m_s(self) -> float: ...
+    @property
+    def clock_s(self) -> float: ...
+    @property
+    def clock_rate_s_s(self) -> float: ...
+    def __repr__(self) -> str: ...
+
 class ExactSp3ValidationError(ValueError):
     kind: str | None
     detail: dict[str, object] | None
@@ -13177,9 +13199,13 @@ class LnavDecoded:
     @property
     def crs(self) -> float: ...
     @property
-    def delta_n(self) -> float: ...
+    def delta_n(self) -> float:
+        """Mean motion difference, semicircles/second."""
+        ...
     @property
-    def m0(self) -> float: ...
+    def m0(self) -> float:
+        """Mean anomaly at reference time, semicircles."""
+        ...
     @property
     def cuc(self) -> float: ...
     @property
@@ -13197,19 +13223,29 @@ class LnavDecoded:
     @property
     def cic(self) -> float: ...
     @property
-    def omega0(self) -> float: ...
+    def omega0(self) -> float:
+        """Longitude of ascending node at reference time, semicircles."""
+        ...
     @property
     def cis(self) -> float: ...
     @property
-    def i0(self) -> float: ...
+    def i0(self) -> float:
+        """Inclination at reference time, semicircles."""
+        ...
     @property
     def crc(self) -> float: ...
     @property
-    def omega(self) -> float: ...
+    def omega(self) -> float:
+        """Argument of perigee, semicircles."""
+        ...
     @property
-    def omega_dot(self) -> float: ...
+    def omega_dot(self) -> float:
+        """Rate of right ascension, semicircles/second."""
+        ...
     @property
-    def idot(self) -> float: ...
+    def idot(self) -> float:
+        """Rate of inclination, semicircles/second."""
+        ...
     def __repr__(self) -> str: ...
 
 def lnav_tow(bits: Sequence[int] | bytes) -> int | None: ...

{sidereon-3.0.4 → sidereon-3.0.5}/python/sidereon/data.py

@@ -284,6 +284,7 @@ __all__ = [
     "sp3_merge_input_identity",
     "verify_merge_report",
     "write_sp3",
+    "write_sp3_with_rounding",
     "DistributionSource",
     "Distribution",
     "ProductIdentity",
@@ -5872,9 +5873,10 @@ def write_sp3(sp3: "sidereon.Sp3", path: str, *, gzip: bool = False) -> str:
 
     Pass ``gzip=True`` to gzip-compress the output (pair it with a ``.gz``
     extension on ``path``). Returns the written path. Raises
-    :class:`sidereon.Sp3WriteError` before anything is written when the
-    product holds a value SP3 text cannot state, such as a mean-merged position
-    finer than a millimetre; ``detail`` names the value.
+    :class:`sidereon.Sp3WriteError` before anything is written when a product
+    holds record values finer than SP3 columns; ``detail`` names the value.
+    Merges from column-grid inputs keep records on that grid, including means.
+    Higher-precision inputs can still require :func:`write_sp3_with_rounding`.
     """
     data = sp3.to_sp3_string().encode("ascii")
     if gzip:
@@ -5890,6 +5892,29 @@ def write_sp3(sp3: "sidereon.Sp3", path: str, *, gzip: bool = False) -> str:
     return path
 
 
+def write_sp3_with_rounding(
+    sp3: "sidereon.Sp3", path: str, *, gzip: bool = False
+) -> tuple[str, "sidereon.Sp3RoundingReport"]:
+    """Write rounded SP3 record values atomically and return ``(path, report)``.
+
+    Header fields and epochs retain the exact checks of the strict writer. The
+    immutable report gives maximum absolute changes in SI units for each field.
+    """
+    text, report = sp3.to_sp3_string_with_rounding()
+    data = text.encode("ascii")
+    if gzip:
+        data = _gzip.compress(data)
+    directory = _os.path.dirname(_os.path.abspath(path))
+    _ensure_dir(directory)
+    tmp = _write_temp(directory, data)
+    try:
+        _os.replace(tmp, path)
+    except OSError as exc:
+        _silent_remove(tmp)
+        raise CacheNotWritable(f"cannot write {path}: {exc}") from exc
+    return path, report
+
+
 def fetch_merged_sp3_file(
     target: Union[_dt.date, _dt.datetime],
     centers: Sequence[str],
@@ -5911,10 +5936,11 @@ def fetch_merged_sp3_file(
     path by default. Pass ``return_report=True`` to receive ``(path, report)``
     and retain the exact contributor provenance and stable input identity.
     Nothing is written if the fetch/merge step raises, or if the merged product
-    cannot be written as SP3 text (:class:`sidereon.Sp3WriteError`). A mean or
-    median of positions that differ below the column resolution, or a clock
-    shifted onto the reference clock datum, can be finer than its SP3 column;
-    ``detail`` on the exception names the value.
+    cannot be written as SP3 text (:class:`sidereon.Sp3WriteError`). When
+    contributors are on the SP3 column grid, merge results stay on that grid,
+    including means, medians and clock-datum adjustments. Higher-precision
+    contributor values can still require :func:`write_sp3_with_rounding`;
+    ``detail`` on a strict-writer exception names the value.
     """
     merged, report = fetch_merged_sp3(
         target,

{sidereon-3.0.4 → sidereon-3.0.5}/python/sidereon/data.pyi

@@ -494,3 +494,6 @@ def fetch_merged_sp3_file(
     **fetch_opts: object,
 ) -> Union[str, tuple[str, MergeReport]]: ...
 def write_sp3(sp3: sidereon.Sp3, path: str, *, gzip: bool = ...) -> str: ...
+def write_sp3_with_rounding(
+    sp3: sidereon.Sp3, path: str, *, gzip: bool = ...
+) -> tuple[str, sidereon.Sp3RoundingReport]: ...

{sidereon-3.0.4 → sidereon-3.0.5}/scripts/check-release.py

@@ -13,7 +13,7 @@ import tomllib
 ROOT = Path(__file__).resolve().parents[1]
 CORE_GIT_URL = "https://github.com/neilberkman/sidereon"
 FULL_GIT_REVISION = re.compile(r"^[0-9a-f]{40}$")
-CORE_SOURCE_REVISION = "5cfcdf929b57a5998e304df786735b171ba8afaf"
+CORE_SOURCE_REVISION = "da59d0eba8e273cc37821eac5535662e5d991126"
 
 
 def registry_version(name, value, package_version):
@@ -307,7 +307,7 @@ def main() -> None:
         "ocean.rs": "25946677944425671a92717860ac2d70f255de5403eeb1fbf98b361716821d5c",
         "pole.rs": "b4cc4c16bdd8ce1d8f04073602ab47dfb85a002b946ab192e8d4d2d600f0a1f8",
     }
-    tide_root = ROOT / "third_party_source" / "sidereon-core-3.0.4" / "tides"
+    tide_root = ROOT / "third_party_source" / "sidereon-core-3.0.5" / "tides"
     for filename, expected_digest in tide_sources.items():
         source = tide_root / filename
         if not source.is_file():
@@ -316,7 +316,7 @@ def main() -> None:
         if digest != expected_digest:
             raise SystemExit(
                 f"IERS-derived source disclosure {source} has digest {digest}, "
-                f"expected {expected_digest} from sidereon-core 3.0.4 at "
+                f"expected {expected_digest} from sidereon-core 3.0.5 at "
                 f"{CORE_SOURCE_REVISION}"
             )
 

{sidereon-3.0.4 → sidereon-3.0.5}/scripts/core_source_identity.rs

@@ -8,7 +8,7 @@ use std::path::{Path, PathBuf};
 
 const REGISTRY_SOURCE: &str = "registry+https://github.com/rust-lang/crates.io-index";
 const GIT_SOURCE_PREFIX: &str = "git+https://github.com/neilberkman/sidereon?rev=";
-const RELEASE_VERSION: &str = "3.0.4";
+const RELEASE_VERSION: &str = "3.0.5";
 
 #[derive(Debug)]
 struct LockedPackage {

sidereon-3.0.5/scripts/smoke-release.py

@@ -0,0 +1,154 @@
+#!/usr/bin/env python3
+"""Smoke-test public APIs from the installed release artifact."""
+
+from __future__ import annotations
+
+import math
+import os
+import struct
+
+import sidereon
+
+
+def _sp3(text: str):
+    return sidereon.load_sp3(text.encode("ascii"))
+
+
+def _check_sp3_grid_merge() -> None:
+    first_text = """\
+#cP2020  6 25  0  0  0.00000000       1 ORBIT IGS14 FIT  TST
+## 2111 432000.00000000   900.00000000 59025 0.0000000000000
++    1   G01  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0
+++         0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0
+%c G  cc GPS ccc cccc cccc cccc cccc ccccc ccccc ccccc ccccc
+%c cc cc ccc ccc cccc cccc cccc cccc ccccc ccccc ccccc ccccc
+%f  1.2500000  1.025000000  0.00000000000  0.000000000000000
+%f  0.0000000  0.000000000  0.00000000000  0.000000000000000
+%i    0    0    0    0      0      0      0      0         0
+%i    0    0    0    0      0      0      0      0         0
+/* RELEASE CONSUMER SP3-C FIXTURE
+*  2020  6 25  0  0  0.00000000
+PG01  15000.000000 -20000.000000   5000.000000    100.000000
+EOF
+"""
+    second_text = first_text.replace("15000.000000", "15000.000001")
+    merged, _report = sidereon.merge_sp3(
+        [_sp3(first_text), _sp3(second_text)],
+        sidereon.Sp3MergeOptions(min_agree=1, clock_min_common=1),
+    )
+
+    strict_text = merged.to_sp3_string()
+    reread = sidereon.load_sp3(strict_text.encode("ascii"))
+    first_x = _sp3(first_text).state("G01", 0).position_m[0]
+    second_x = _sp3(second_text).state("G01", 0).position_m[0]
+    x_m = reread.state("G01", 0).position_m[0]
+    if x_m != merged.state("G01", 0).position_m[0] or x_m not in (first_x, second_x):
+        raise SystemExit(
+            "strict SP3 merge round-trip did not preserve a parser-scaled grid value"
+        )
+
+    rounded_text, rounding = merged.to_sp3_string_with_rounding()
+    if rounded_text != strict_text or rounding.position_x_m != 0.0:
+        raise SystemExit(
+            "SP3 rounding report disagrees with the writable merged product"
+        )
+
+
+def _clock_sp3(clocks_us: list[float]):
+    satellites = ["G01", "G02", "G03"]
+    satellite_header = "".join(satellites) + "  0" * 14
+    records = "\n".join(
+        f"P{sat}{15000.0:14.6f}{-20000.0:14.6f}{5000.0:14.6f}{clock:14.6f}"
+        for sat, clock in zip(satellites, clocks_us)
+    )
+    text = f"""\
+#cP2020  6 25  0  0  0.00000000       1 ORBIT IGS14 FIT  TST
+## 2111 432000.00000000   900.00000000 59025 0.0000000000000
++    3   {satellite_header}
+++         0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0
+%c G  cc GPS ccc cccc cccc cccc cccc ccccc ccccc ccccc ccccc
+%c cc cc ccc ccc cccc cccc cccc cccc ccccc ccccc ccccc ccccc
+%f  1.2500000  1.025000000  0.00000000000  0.000000000000000
+%f  0.0000000  0.000000000  0.00000000000  0.000000000000000
+%i    0    0    0    0      0      0      0      0         0
+%i    0    0    0    0      0      0      0      0         0
+/* RELEASE CONSUMER SP3 CLOCK FIXTURE
+*  2020  6 25  0  0  0.00000000
+{records}
+EOF
+"""
+    return _sp3(text)
+
+
+def _check_sp3_half_grid_clock_merge() -> None:
+    first = _clock_sp3([-598.410235, 138.983111, 0.0])
+    second = _clock_sp3([-598.410234, 138.983110, 0.0])
+    first.to_sp3_string()
+    second.to_sp3_string()
+    merged, _report = sidereon.merge_sp3(
+        [first, second], sidereon.Sp3MergeOptions(min_agree=1, clock_min_common=3)
+    )
+
+    strict_text = merged.to_sp3_string()
+    reread = sidereon.load_sp3(strict_text.encode("ascii"))
+    for satellite in ("G01", "G02", "G03"):
+        left = first.state(satellite, 0).clock_s
+        right = second.state(satellite, 0).clock_s
+        clock = reread.state(satellite, 0).clock_s
+        if clock != merged.state(satellite, 0).clock_s or clock not in (left, right):
+            raise SystemExit(f"strict SP3 merge round-trip changed {satellite} clock")
+
+    rounded_text, report = merged.to_sp3_string_with_rounding()
+    if (
+        rounded_text != strict_text
+        or report.clock_s != 0.0
+        or report.clock_rate_s_s != 0.0
+    ):
+        raise SystemExit(
+            "SP3 half-grid clock merge did not return a zero rounding report"
+        )
+
+
+def main() -> None:
+    expected_version = os.environ.get("SIDEREON_EXPECTED_VERSION")
+    if expected_version is not None and sidereon.__version__ != expected_version:
+        raise SystemExit(
+            f"installed sidereon {sidereon.__version__}, expected {expected_version}"
+        )
+
+    header = struct.pack(">ddddii", -90.0, -180.0, 0.25, 0.25, 721, 1440)
+    data = header + struct.pack(">f", 2.5) * (721 * 1440)
+    grid = sidereon.GeoidGrid.from_proj_egm96_gtx(data)
+
+    for arithmetic in (
+        sidereon.ProjVgridshiftArithmetic.SEPARATE_MULTIPLY_ADD,
+        sidereon.ProjVgridshiftArithmetic.FUSED_MULTIPLY_ADD,
+    ):
+        got = grid.undulation_proj_rad(
+            math.radians(40.0), math.radians(-105.0), arithmetic
+        )
+        if got != 2.5:
+            raise SystemExit(f"PROJ EGM96 GTX smoke returned {got}, expected 2.5")
+
+    try:
+        grid.undulation_proj_rad(
+            math.nan,
+            0.0,
+            sidereon.ProjVgridshiftArithmetic.FUSED_MULTIPLY_ADD,
+        )
+    except sidereon.ProjVgridshiftNonFiniteCoordinateError as error:
+        if not isinstance(error, (sidereon.ProjVgridshiftError, ValueError)):
+            raise SystemExit(
+                "typed PROJ coordinate exception hierarchy is broken"
+            ) from None
+    else:
+        raise SystemExit("non-finite PROJ coordinate did not raise")
+
+    _check_sp3_grid_merge()
+    _check_sp3_half_grid_clock_merge()
+
+    print(f"sidereon {sidereon.__version__} release smoke passed")
+
+
+if __name__ == "__main__":
+    main()