kshana 0.20.0 → 0.23.0
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- package/README.md +65 -1025
- package/kshana-sbom.cdx.json +1135 -0
- package/kshana_bg.wasm +0 -0
- package/package.json +6 -4
package/README.md
CHANGED
|
@@ -1,1065 +1,105 @@
|
|
|
1
|
-
|
|
2
|
-
|
|
3
|
-
|
|
1
|
+
<!-- Surface README for the npm WebAssembly package. Copied into web/pkg/README.md by
|
|
2
|
+
web/build.sh after wasm-pack runs (web/pkg is gitignored/generated). Images/links are
|
|
3
|
+
ABSOLUTE (pinned to /main) because npm does not rewrite relative paths or render Mermaid.
|
|
4
|
+
The canonical, full README lives at README.md on GitHub. To re-pin images to an immutable
|
|
5
|
+
release tag at publish time, replace `/main/` with `/vX.Y.Z/` across this file. -->
|
|
4
6
|
|
|
5
7
|
<p align="center">
|
|
6
|
-
<img src="docs/assets/kshana-wordmark.png" alt="Kshana" width="300">
|
|
8
|
+
<img src="https://raw.githubusercontent.com/AshfordeOU/kshana/main/docs/assets/kshana-wordmark.png" alt="Kshana" width="300">
|
|
7
9
|
</p>
|
|
8
10
|
|
|
9
11
|
<p align="center">
|
|
10
12
|
<strong>क्षण</strong> — Sanskrit for <em>the precise instant</em>, the smallest measure of time.<br>
|
|
11
|
-
Open, reproducible PNT-resilience simulation
|
|
13
|
+
Open, reproducible PNT-resilience simulation, compiled to WebAssembly — the whole engine, in the browser.
|
|
12
14
|
</p>
|
|
13
15
|
|
|
14
16
|
<p align="center">
|
|
15
|
-
<a href="https://
|
|
16
|
-
<a href="
|
|
17
|
-
<a href="
|
|
18
|
-
<a href="https://
|
|
19
|
-
<a href="https://github.com/
|
|
20
|
-
<a href="https://github.com/ashfordeOU/kshana/releases"><img src="https://img.shields.io/badge/release-v0.20.0-c79e63" alt="Release v0.20.0"></a>
|
|
21
|
-
<a href="https://plugins.jetbrains.com/plugin/32181-kshana--pnt-simulator"><img src="https://img.shields.io/badge/JetBrains-Marketplace-c79e63" alt="Kshana on the JetBrains Marketplace"></a>
|
|
22
|
-
<a href="https://glama.ai/mcp/servers/ashfordeOU/kshana"><img src="https://glama.ai/mcp/servers/ashfordeOU/kshana/badges/score.svg" alt="kshana-mcp on Glama — MCP server quality score"></a>
|
|
23
|
-
<a href="LICENSE"><img src="https://img.shields.io/badge/License-AGPL_v3-blue.svg" alt="License: AGPL-3.0-only"></a>
|
|
24
|
-
<a href="LICENSING.md"><img src="https://img.shields.io/badge/commercial_licence-available-2ea043" alt="Commercial licence available from Ashforde OÜ"></a>
|
|
25
|
-
<a href="Cargo.toml"><img src="https://img.shields.io/badge/rust-1.75%2B-orange.svg" alt="Rust 1.75+"></a>
|
|
17
|
+
<a href="https://github.com/AshfordeOU/kshana/blob/main/tests/sgp4_verification.rs"><img src="https://img.shields.io/badge/SGP4-666%2F666%20AIAA%20vectors%20%C2%B7%204.12mm-3fb950" alt="SGP4 validated against all 666 AIAA 2006-6753 vectors, worst 4.12 mm"></a>
|
|
18
|
+
<a href="https://github.com/AshfordeOU/kshana#validation-at-a-glance"><img src="https://img.shields.io/badge/validated-40%20external%20oracles-3fb950" alt="40 of 91 capabilities validated against independent external oracles"></a>
|
|
19
|
+
<a href="https://github.com/AshfordeOU/kshana/releases"><img src="https://img.shields.io/badge/release-v0.23.0-c79e63" alt="Release v0.23.0"></a>
|
|
20
|
+
<a href="https://ashforde.org"><img src="https://img.shields.io/badge/playground-try%20in%20browser-c79e63" alt="Live playground — run in your browser, no install"></a>
|
|
21
|
+
<a href="https://github.com/AshfordeOU/kshana/blob/main/LICENSE"><img src="https://img.shields.io/badge/License-AGPL_v3-blue.svg" alt="License: AGPL-3.0-only"></a>
|
|
26
22
|
<a href="https://doi.org/10.5281/zenodo.20528627"><img src="https://img.shields.io/badge/DOI-10.5281%2Fzenodo.20528627-blue.svg" alt="DOI 10.5281/zenodo.20528627"></a>
|
|
27
23
|
</p>
|
|
28
24
|
|
|
29
|
-
|
|
30
|
-
|
|
31
|
-
|
|
32
|
-
|
|
33
|
-
|
|
34
|
-
|
|
35
|
-
contrast decay, and vibration coupling) that <em>derives</em> the noise coefficient
|
|
36
|
-
rather than looking it up; it is not yet a full quantum-physics simulator (Coriolis and
|
|
37
|
-
light-shift systematics remain coefficient-level — see
|
|
38
|
-
<a href="docs/QUANTUM.md">docs/QUANTUM.md</a> and
|
|
39
|
-
<a href="docs/QUANTUM-MODELS.md">docs/QUANTUM-MODELS.md</a>).
|
|
40
|
-
</p>
|
|
41
|
-
|
|
42
|
-
It quantifies, in hard and reproducible numbers, what quantum clocks, quantum
|
|
43
|
-
inertial sensors, and optical time-transfer buy a navigation system over classical
|
|
44
|
-
PNT — scored against the operational figures of merit that matter for resilient
|
|
45
|
-
navigation. Every result is reproducible from `scenario + seed + engine version`,
|
|
46
|
-
and every sensor parameter is traceable to a published source — consolidated in one
|
|
47
|
-
citable table in [`docs/PROVENANCE.md`](docs/PROVENANCE.md).
|
|
48
|
-
|
|
49
|
-
*Free and open source under the GNU AGPL-3.0 — with a commercial licence available
|
|
50
|
-
from Ashforde OÜ for proprietary/closed integration (see [`LICENSING.md`](LICENSING.md)).
|
|
51
|
-
Professionally developed and maintained by Ashforde OÜ; commercial support,
|
|
52
|
-
integration, and proprietary extensions available.*
|
|
53
|
-
|
|
54
|
-
> **Status: v0.20.0 · a simulation substrate, not yet a product.** A validated,
|
|
55
|
-
> fully reproducible engine spanning the PNT stack — orbit geometry and constellation
|
|
56
|
-
> design, a numerical (Cowell) propagator with a six-perturbation force model, maneuver
|
|
57
|
-
> and trajectory design, time systems, inertial navigation (incl. map-aided and
|
|
58
|
-
> gravity-map-matching alt-PNT), GNSS/INS fusion (loose, tight, UKF, coupled
|
|
59
|
-
> clock+position, 17-state), orbit determination, ARAIM integrity, clocks, advanced
|
|
60
|
-
> time-and-frequency transfer, the GNSS measurement domain, resilience (jamming +
|
|
61
|
-
> multi-layer spoofing), and an open **deep-space / Mars radiometric navigation**
|
|
62
|
-
> engine (light-time + Shapiro, CCSDS-TDM, reduced-dynamic SRIF, one-/two-way fusion);
|
|
63
|
-
> plus first-order **mission-analysis** budgets (launch / re-entry / EO-coverage / pointing /
|
|
64
|
-
> ground-station passes / link), a **space-weather** environment model, an **AI/ML
|
|
65
|
-
> RF-impairment** evaluation testbed, and the versioned **Kshana Interchange Format (KIF)**.
|
|
66
|
-
> Honest by design: every figure of merit is labelled *validated* or *modelled*, and
|
|
67
|
-
> optical-clock figures are space goals on ground hardware (no strontium optical clock has flown).
|
|
68
|
-
>
|
|
69
|
-
> **Validation ladder** (maturity is *not* uniform across domains — and saying so is the point):
|
|
70
|
-
> | Domain | Tier |
|
|
71
|
-
> |---|---|
|
|
72
|
-
> | Earth PNT (orbit, frames, time, clocks, IMU, integrity) | **Real-data validated** — ESA SP3 (Galileo 0.13 m, Swarm-A 0.10 m), NIST SP1065, SOFA/ERFA, heritage vectors |
|
|
73
|
-
> | Deep-space / Mars navigation | **Simulation-validated** — synthetic closed-loop OD + analytic self-consistency; Sun-central dynamics cross-checked vs JPL **DE440** (137 m @ 1-day arc) |
|
|
74
|
-
> | Real-mission deep-space OD | **Roadmap** — pending real DSN/ESTRACK tracking-data validation |
|
|
75
|
-
>
|
|
76
|
-
> Deep-space figures (Mars-LMO OD ≈ 0.2 m; relay-PNT orbiter 0.4 m / rover 5.1 m) are **simulation / covariance figures of merit**, not real-mission results.
|
|
77
|
-
> See **[Capabilities](#capabilities)** for what it does, **[What it is / is not](#what-it-is--is-not)**
|
|
78
|
-
> for scope, and [`docs/CAPABILITY.md`](docs/CAPABILITY.md) / [`docs/VALIDATION.md`](docs/VALIDATION.md)
|
|
79
|
-
> for per-capability maturity. The overclaim closure ledger
|
|
80
|
-
> [`docs/CLAIMS-VS-REALITY.md`](docs/CLAIMS-VS-REALITY.md) tracks every historical overclaim,
|
|
81
|
-
> how it was resolved, and a CI guard (`tests/no_overclaims.rs`) that keeps it resolved.
|
|
82
|
-
|
|
83
|
-
> **Try it in your browser:** the [playground](web/) runs the engine client-side as
|
|
84
|
-
> WebAssembly — pick a scenario, edit the parameters, and see the result, with nothing
|
|
85
|
-
> uploaded. Build it locally with `./web/build.sh` (see [`web/README.md`](web/README.md)),
|
|
86
|
-
> or publish it to GitHub Pages via the `pages` workflow.
|
|
87
|
-
|
|
88
|
-
> **New to this?** In plain terms: GPS-style satellite signals tell things *where they
|
|
89
|
-
> are* and *what time it is*. When those signals are lost (jammed, blocked, or out of
|
|
90
|
-
> view in space), a system has to keep going on its own onboard clock and motion
|
|
91
|
-
> sensors — and they slowly drift. "Quantum" clocks and sensors drift far more slowly.
|
|
92
|
-
> Kshana measures, in honest numbers, **how much longer a quantum-equipped system can
|
|
93
|
-
> coast** before it exceeds its accuracy limits. New readers should start with the
|
|
94
|
-
> [plain-language primer](docs/CONCEPTS.md) and the [glossary](docs/GLOSSARY.md).
|
|
95
|
-
|
|
96
|
-
---
|
|
97
|
-
|
|
98
|
-
## Contents
|
|
99
|
-
|
|
100
|
-
- [Why](#why) · [What it is / is not](#what-it-is--is-not) · [Capabilities](#capabilities) · [Results](#results)
|
|
101
|
-
- [Install & build](#install--build) · [Usage](#usage) ([Python](#python), [WebAssembly](#webassembly))
|
|
102
|
-
- [Scenario format](#scenario-format) · [Output](#output) · [Architecture](#architecture)
|
|
103
|
-
- [Repository layout](#repository-layout) · [Validation & honesty](#validation-reproducibility--honesty)
|
|
104
|
-
- [Documentation](#documentation) · [FAQ](#faq) · [Troubleshooting](#troubleshooting)
|
|
105
|
-
- [Roadmap](#roadmap) · [Contributing](#contributing) · [Citing](#citing) · [Versioning & releases](#versioning--releases) · [License](#license)
|
|
106
|
-
- [Support & professional services](#support--professional-services) · [References](#key-references)
|
|
107
|
-
|
|
108
|
-
## Why
|
|
109
|
-
|
|
110
|
-
Resilient PNT depends on holding position and time when GNSS is denied or jammed.
|
|
111
|
-
Quantum sensors promise far slower drift during those outages. There is no good
|
|
112
|
-
**open** tool to quantify that advantage honestly and reproducibly — so primes,
|
|
113
|
-
agencies, and labs each rebuild private one-offs. Kshana aims to be the neutral,
|
|
114
|
-
citable reference for exactly this question.
|
|
115
|
-
|
|
116
|
-
The engine knows nothing about "quantum" vs "classical": each sensor is an
|
|
117
|
-
**error model** plugged into a common pipeline, so a quantum and a classical
|
|
118
|
-
device are compared *apples-to-apples* on the same scenario, with independent
|
|
119
|
-
noise realizations.
|
|
120
|
-
|
|
121
|
-
## What it is / is not
|
|
122
|
-
|
|
123
|
-
**It is:** a deterministic, dependency-light engine spanning the PNT stack — orbit
|
|
124
|
-
geometry, inertial navigation, GNSS/INS fusion, integrity, clocks, and timing. It
|
|
125
|
-
runs a scenario (often a GNSS outage), evolves calibrated sensor error models
|
|
126
|
-
through the appropriate estimator, and scores the result against the operational
|
|
127
|
-
figures of merit — emitting a reproducible JSON result and an SVG chart, from a
|
|
128
|
-
Rust library, a CLI, a Python extension, an in-browser WebAssembly module, a
|
|
129
|
-
**Model Context Protocol (MCP) server** for AI agents, or a **JetBrains IDE plugin**.
|
|
130
|
-
|
|
131
|
-
**It is not:** flight hardware, a quantum-payload design, a full GNSS signal
|
|
132
|
-
receiver, or a certified avionics product. Quantum-hardware fidelity comes from
|
|
133
|
-
published error models, not from this tool. The granular maturity of each
|
|
134
|
-
capability is documented in [`docs/CAPABILITY.md`](docs/CAPABILITY.md).
|
|
25
|
+
**Kshana** is an open, reproducible **PNT-resilience simulator with quantum-sensor
|
|
26
|
+
performance models** — positioning, navigation, and timing. This package is the Rust
|
|
27
|
+
engine compiled to **WebAssembly**: it runs entirely client-side — pass a scenario TOML
|
|
28
|
+
string in, get a reproducible JSON result and an SVG chart back, with nothing uploaded.
|
|
29
|
+
Every result is reproducible from `scenario + seed + engine version`, and every sensor
|
|
30
|
+
parameter is traceable to a published source.
|
|
135
31
|
|
|
136
|
-
|
|
137
|
-
|
|
138
|
-
|
|
139
|
-
|
|
140
|
-
|
|
141
|
-
a full GNSS *signal-acquisition* receiver (it now solves a single-point **PVT** position
|
|
142
|
-
fix from real RINEX code observations — validated on real IGS data — but does **not**
|
|
143
|
-
acquire or track raw signal); or a full mission-design suite (it has Lambert / porkchop /
|
|
144
|
-
maneuver / orbit-determination building blocks, but is the performance-simulation layer
|
|
145
|
-
*above* GMAT/Orekit, not a replacement). Owning this scope is deliberate. If you need first-principles cold-atom
|
|
146
|
-
interferometer error budgets (e.g. CARIOQA-PMP-grade or X-37B-style validation), see
|
|
147
|
-
the P2 roadmap and [get in touch](#support--professional-services) to collaborate.
|
|
148
|
-
|
|
149
|
-
## Capabilities
|
|
150
|
-
|
|
151
|
-
| Domain | Capability |
|
|
152
|
-
|--------|------------|
|
|
153
|
-
| **Orbit & geometry** | SGP4/SDP4 propagation (validated to 4.12 mm against all 666 AIAA 2006-6753 vectors); real two-line elements (a committed, date-stamped Celestrak `gps-ops` snapshot) or synthetic Walker-delta constellations whose mean elements realise the `i:T/P/F` formula to under 1 km over a 24 h propagation; multi-constellation visibility, **dilution of precision (GDOP/PDOP/HDOP/VDOP/TDOP, validated to 1e-6 against gnss_lib_py 1.0.4, Stanford NAV Lab)**, and GNSS availability; a gradient-free constellation-design optimiser, streets-of-coverage minimum-satellite sizing, a multi-constellation comparison tool, and a Walker **design sweep** that tabulates coverage / PDOP / revisit-time over a planes × satellites grid and reports the Pareto-optimal designs. |
|
|
154
|
-
| **Numerical propagator** | A **Cowell** numerical propagator (`src/propagator.rs`) complementing the analytic SGP4/SDP4 path, with a hierarchical **six-perturbation** force model (`src/forces.rs`): two-body + the full **J2–J6 zonal** field (the exact analytic gradient of its disturbing potential), an optional **EGM2008 tesseral spherical-harmonic geopotential to degree/order 70** (`src/gravity_sh.rs`; real NGA coefficients, Holmes–Featherstone normalized-Legendre recurrence, cross-checked against the closed-form Legendre functions and the analytic ∇V identity), **epoch-driven Sun and Moon third-body** gravity (a built-in low-precision ephemeris, no DE/SPK kernel), **solar-radiation pressure** (cannonball model with a conical umbra+penumbra shadow), **atmospheric drag** (Vallado piecewise-exponential density, co-rotating atmosphere), the **post-Newtonian Schwarzschild relativistic correction**, and the **Lense–Thirring frame-dragging** term (IERS 2010 §10, linear in Earth's angular momentum, ~1–2 orders below Schwarzschild) — driven by a choice of two adaptive integrators (RK4 step-doubling or the **Dormand–Prince RK5(4)** embedded pair). Validated against analytic truth stronger than a cross-tool would give: the unperturbed orbit matches the exact universal-variable Kepler solution to **sub-metre over 24 h**, energy/angular-momentum conserve to ~1e-9, and each perturbation matches a hand-derived closed-form signature. |
|
|
155
|
-
| **Maneuvers & trajectory design** | Impulsive ΔV nodes with 6×6 covariance propagation (ECI / LVLH execution-error frames), finite-burn integration checked against the closed-form **Tsiolkovsky** rocket equation to < 0.01 %, an **Izzo-2015** single-revolution **Lambert** solver, an exact universal-variable **Kepler** propagator, and a **porkchop** (launch × arrival) C3 / arrival-V∞ sweep emitted as a JSON contour grid — the performance-simulation layer above GMAT/Orekit, with every Lambert output round-tripped against two-body truth and the porkchop minimum checked against the analytic Hohmann floor. |
|
|
156
|
-
| **Time systems & reference frames** | IERS leap-second **UTC / TAI / TT / UT1** scales, a Julian-date API, the IAU-2000 **Earth Rotation Angle**, GMST-based **TEME ↔ ECEF** with WGS-84 geodetic frames, IAU 2006 precession (Fukushima–Williams), full **IAU 2000A/2000B nutation**, IERS **polar motion**, and the equinox-free **CIO-based IAU 2006/2000A GCRS↔ITRS** reduction — all validated **bit-for-bit** against the SOFA/ERFA vectors, and **independently cross-checked against ANISE** (the pure-Rust NAIF/SPICE reimplementation): kshana's GCRS→ITRS vs ANISE's ITRF93 from JPL's `earth_latest_high_prec.bpc`, the same IERS Earth-orientation parameters fed to both, agree to **≤ 0.86 m on the ground / ≤ 3.6 m at GNSS orbit** (max 0.028″) across eight epochs 2020–2023. |
|
|
157
|
-
| **Inertial** | Three-axis strapdown INS — quaternion attitude, WGS-84 NED mechanization, coning/sculling compensation, and a deterministic IMU error model (scale-factor, misalignment, g-sensitivity, quantization, drift); a **first-principles cold-atom-interferometer accelerometer** (Mach–Zehnder phase, quantum projection noise, contrast decay, vibration coupling) that *derives* the velocity-random-walk coefficient; and a sequential-importance-resampling **particle filter** for map-aided (terrain-/gravity-referenced) GPS-denied navigation. |
|
|
158
|
-
| **Alt-PNT (GPS-denied)** | A cold-atom **gravimeter measurement model** whose white-noise floor (`σ = ASD/√τ`) is derived from the CAI accelerometer physics; a low-degree, fully-normalised **spherical-harmonic gravity-anomaly field** (validated against the closed-form Legendre functions and a hand-derived single-term anomaly) plus synthetic mascons; and a **gravity-map-matching particle filter** that recovers a GPS-denied track from the anomaly sequence it flies through. It extends to **terrain-referenced navigation** (TERCOM/SITAN against an SRTM `.hgt` DEM, `src/altpnt/terrain.rs`), an **IGRF-14 geomagnetic main field** to degree/order 13 (`src/igrf.rs`, validated against the tilted-dipole closed form and ∇V finite differences), and a **combined gravity + magnetic + terrain** navigator that fuses all three scalar channels through one particle filter (information is additive — no channel makes the fix worse). A **60-minute GPS-denied benchmark** (a ~700 km / one-hour outage where the inertial solution drifts to ~70 km) is recovered to **~145 m (< 500 m)** by a hierarchical coarse-to-fine matcher — the ESA NAVISP *Quantum Wayfarer* target. |
|
|
159
|
-
| **Fusion** | Loosely-coupled 15-state GNSS/INS error-state EKF with closed-loop feedback (the `gnss-ins` pack); a **tightly-coupled** pseudorange update that keeps correcting with fewer than four satellites; a coupled **clock + position** filter; a general **unscented (sigma-point) Kalman** estimator for strongly nonlinear measurements; a tightly-coupled GNSS/INS **UKF navigator** (pseudorange + Doppler) whose force-model orbital coast is validated to **0.77 m RMS** over a 30-minute curving LEO pass that includes a 120-second GNSS outage; and a full **17-state tightly-coupled GNSS/INS UKF** (position, velocity, attitude error, accelerometer and gyro biases, clock bias and drift) whose **quantum-CAI dead-reckoning** coasts a 120-second outage on the cold-atom accelerometer's derived velocity-random-walk. |
|
|
160
|
-
| **Orbit determination** | Recovery of an orbital state `[r, v]` from ground-station range tracking, composing the two-body + J2 force model and RK4 integrator with a **Gauss–Newton batch** corrector (`determine_orbit_batch`, sub-metre / mm·s⁻¹ from noiseless ranges, ~2 m at a 5 m noise floor) and a **sequential** unscented-filter variant (`determine_orbit_sequential`). |
|
|
161
|
-
| **Lunar & cislunar** | An Earth–Moon **circular restricted three-body (CR3BP)** propagator in the rotating frame — conserved Jacobi constant and all five Lagrange points (`src/cr3bp.rs`) — now with a **6×6 state-transition matrix and a single-shooting differential corrector** (`cr3bp_jacobian`, `propagate_state_stm`, `differential_correct_halo`) that produces genuinely periodic **halo / NRHO** orbits: the STM is validated against finite differences, corrected orbits close to machine precision, and seeding the published apolune state reproduces the **L2 southern 9:2 NRHO** (the Gateway orbit) at period ≈ 6.57 d / perilune ≈ 3,250 km, consistent with the published ≈ 6.56 d / ≈ 3,370 km (a CR3BP — circular, Sun-free — solution, **not** validated against a real LANS/Gateway ephemeris; the selenocentric MCI/MCMF transform of the corrected orbit is a follow-on); plus **LunaNet / LNIS** cislunar PNT geometry (MCI↔MCMF reduction, selenographic coordinates) with a **lunar south-pole ARAIM** pass that honestly surfaces the integrity gap: a ~30 m σ_URE drives the protection level well above a 50 m alert limit (`src/lunar.rs`, `scenarios/lunanet-araim.toml`). |
|
|
162
|
-
| **Deep-space & Mars PNT** | An open **radiometric navigation engine**: iterative light-time + **Shapiro** relativistic delay, two-/one-/three-way **Doppler & range** (Moyer two-leg), coherent transponder turnaround ratios, regenerative/PN ranging (CCSDS 414), and **Δ-DOR** plane-of-sky (CCSDS 506), with solar-plasma/tropo/iono media; **CCSDS-TDM (503)** tracking-data-message parse + emit; a **reduced-dynamic Square-Root Information Filter** (RTN empirical accelerations + a 3-state onboard clock + Mars atmospheric drag) that does **Mars-LMO orbit determination to ≈ 0.2 m** in a synthetic closed loop; a joint **one-way + two-way fusion** estimator; a multi-body dynamics core (`Body{μ, re, zonals, gravity, IAU-pole}`, Mars GMM-3 gravity, an IAU body-fixed Mars frame, a pluggable `EphemerisProvider` seam, two-part Julian dates + TT↔TDB); and the **`mars-pnt`** relay-PNT scenario (a MARCONI areostationary relay constellation) with an end-to-end **GSE performance simulator** (geometry → link budget → observables → SRIF → covariance). **Simulation-validated** (covariance / closed-loop figures of merit); the Sun-central Mars dynamics are cross-checked against JPL **DE440** (137 m @ 1-day arc, `xval/anise-mars-od`). Real DSN/ESTRACK tracking-data validation is on the roadmap. |
|
|
163
|
-
| **Integrity** | Snapshot and solution-separation (ARAIM-style) RAIM with horizontal/vertical protection levels (HPL/VPL), fault detection & exclusion, and Stanford integrity diagrams; an explicit integrity-risk-budget (**MHSS**) protection level, including the **dual-/multi-constellation constellation-wide fault mode** (EU ARAIM / DO-316), exercised on a real GPS + Galileo snapshot (`scenarios/araim-gps-galileo.toml`). The protection level applies the one-sided **nominal-bias** projection `b_k = Σ_i|s_i|·b_nom` per fault mode and the **integrity** sigma σ_URA (distinct from the accuracy σ_URE) from the Integrity Support Message — see [`docs/ARAIM_REFERENCE.md`](docs/ARAIM_REFERENCE.md). |
|
|
164
|
-
| **Augmentation (SBAS)** | **SBAS / WAAS protection levels** in the DO-229E weighted-least-squares form (precision-approach and en-route K-factors) and the **L1/L5 dual-frequency ionosphere-free** combination (IS-GPS-705, γ₁₅ ≈ 1.793) that underpins DO-316 — `src/sbas.rs`, with the protection-level geometry cross-checked against a NumPy `inv(GᵀG)` reference. |
|
|
165
|
-
| **Clock & timing** | Two-state Kalman holdover (Joseph-form covariance, NIS/NEES consistency health); Allan-family stability (ADEV / MDEV / TDEV / HDEV) with noise-type-specific confidence intervals and a full **IEEE-1139 five-coefficient power-law fit**; geometric corrections (Sagnac, GNSS common-view); and the operational transfer methods — **TWSTFT** with the BIPM Sagnac closed form, **GNSS common-view**, **PPP** ionosphere-free time transfer, a free-space **optical** link with turbulence scintillation, and an inverse-variance **clock-ensemble (paper) timescale** below the best contributing clock. A **GNSS-denied clock-holdover calculator** (`src/holdover.rs`) exposes the closed-form van-Loan coast-error growth as a *holdover-to-threshold* inversion — how long a clock free-runs before its timing error exceeds budget — across representative classical and quantum-clock classes; **modelled** (cross-checked against the multi-step `clock_state` covariance recursion), and honest that for a very stable clock the holdover to a tight threshold is set by the *assumed* long-tau noise floor, not the cited ADEV. |
|
|
166
|
-
| **GNSS measurement domain** | Forward pseudorange / Doppler synthesis with **Klobuchar** (broadcast) and **IONEX / TEC-grid** (measured) ionosphere — including an IONEX file parser, time interpolation between maps, and the thin-shell slant-obliquity mapping — **Saastamoinen + Niell** troposphere, and snapshot RAIM (HPL/VPL). |
|
|
167
|
-
| **Resilience** | Link-budget **jamming** (J/S → effective C/N₀ → loss of lock, with the anti-jam spectral-separation factor `Q` now **derived from the actual signal and jammer power spectra** via `src/navsignal.rs` — `Q = 1/(R_c·κ)`, cross-checked in CI against the previous representative constant); a stochastic **time-spoof detector** (Neyman–Pearson / χ²₁ energy test with closed-form and Monte-Carlo P_fa/P_md and a Security FoM of 1 − P_md); and a **multi-layer spoof detector** fusing a RAIM-consistency parity test (with the common-mode blind spot modelled honestly), an RF AGC-power monitor, and a signal-quality (SQM early-minus-late) monitor; and a **quantum-inertial dead-reckoning error budget** (`QuantumNavBudget`, `src/inertial/quantum_imu.rs`) composing the cold-atom-interferometer white-noise velocity-random-walk with residual bias (cross-checked against the independent `AccelModel` integrator) and scale-factor error into a position-drift-over-holdover figure — the inertial twin of the clock holdover. |
|
|
168
|
-
| **Nav-signal & code tracking** | The **signal level** between the link budget and the measurement domain (`src/navsignal.rs`): unit-area **power spectral densities** for **BPSK-R(n)** and **sine-BOC(m,n)**; the **spectral-separation coefficient** κ = ∫ G_s·G_i df, which **derives the anti-jam `Q`** the jamming model uses (`Q = 1/(R_c·κ)`) from the actual signal/jammer spectra instead of a representative constant; the **RMS (Gabor) bandwidth** (BOC > BPSK — the ranging-information / Cramér–Rao measure); the **coherent early–late DLL code-tracking thermal-noise jitter** (Kaplan & Hegarty; ~sub-metre for C/A at 45 dB-Hz); and the **multipath error envelope** (coherent EML — narrow-correlator suppression). Validated against closed-form anchors (BPSK self-SSC = 2/(3·R_c), unit-area PSDs, sub-metre C/A jitter). This is signal-**performance** analysis, **not** antenna / RF-payload hardware design (a payload partner's role). |
|
|
169
|
-
| **Interoperability** | **RINEX-3** multi-GNSS broadcast-ephemeris ingestion (GPS, Galileo, QZSS, BeiDou MEO/IGSO via IS-GPS-200; GLONASS via PZ-90 state-vector RK4) usable as a constellation source (RINEX in, PNT geometry out); a **RINEX-3/4** observation parser (pseudorange, carrier phase, Doppler, signal strength) that now **feeds a single-point-positioning solver** (`pvt`) — real code observations in, a real **receiver position** out, validated on IGS data; an **SP3-c/d** precise-ephemeris reader/writer with 9th-order Lagrange interpolation; and **CCSDS OEM 2.0 + OMM** (mean-elements) export for flight-dynamics tools (GMAT, Orekit, STK); and **CCSDS-TDM (503)** tracking-data-message parse + emit for deep-space radiometric tracking. |
|
|
170
|
-
| **Mission analysis (systems engineering)** | First-order mission-design budgets, each a runnable kind: two-body **launch & ascent geometry** (`launch-window` — launch azimuth `sin Az = cos i/cos lat`, minimum inclination, Earth-rotation bonus, dogleg plane-change Δv, daily opportunities; `src/launch.rs`); an **Allen–Eggers ballistic re-entry corridor** (`reentry` — peak deceleration, peak-g velocity/altitude, peak-heating velocity; `src/reentry.rs`); **Earth-observation coverage geometry** (`eo-coverage` — swath / nadir GSD / off-nadir access / revisit via the SMAD space triangle; `src/eo_payload.rs`); a **3-DOF attitude & pointing error budget** (`attitude-budget` — worst-case gravity-gradient torque + RSS pointing budget; `src/attitude_budget.rs`); **ground-station pass prediction** (`passes` — AOS/TCA/LOS, max elevation, access time; `src/passes.rs`); and a **one-way link budget** over the CCSDS 401 / DSN 810-005 link equation (`link-budget` — FSPL, C/N₀, Eb/N₀, margin, closure; `src/linkbudget.rs`). **MODELLED** first-order analytic budgets — the pre-hardware layer below STK/GMAT/Basilisk, not a 6-DoF or radiometric replacement. |
|
|
171
|
-
| **Space environment** | A **space-weather environment model** (`space-weather`, `src/space_weather.rs`): solar (F10.7 / centred-81-day F10.7a) and geomagnetic (Kp, with the definitional Kp↔ap table) activity indices, the **Jacchia-1971** exospheric temperature they drive (validated vs published solar min/mean/max), and the activity-corrected vs static thermospheric neutral density at altitude — the solar-cycle density dependence the static USSA76 atmosphere omits. **MODELLED**: a calibrated first-order scale-height coupling, **not** a data-validated (NRLMSISE) atmosphere. |
|
|
172
|
-
| **AI/ML evaluation & trade** | An **RF-impairment detection evaluation testbed** (`impairment-eval`, `src/impairment_eval.rs`): a labelled, parameter-grounded **synthetic** corpus (nominal / jamming / spoof-time / spoof-position / multipath), a detector-agnostic **ROC/AUC** harness scoring any detector (energy \| agc \| sqm \| parity \| fused) with per-class Pd at a target Pfa, and the in- vs out-of-distribution **optimism gap** (distribution-shift mode). Plus a **quantum-vs-classical PNT trade** (`quantum-trade`, `src/quantum_trade.rs`) quantifying a candidate clock's timing/inertial holdover benefit from a **measured-ADEV** curve vs a classical baseline, with the long-τ floor caveat carried on the artifact and a GNSS-denied resilience-vs-time envelope. The evaluation **metrics** (AUC / confusion / Pd-Pmd) are **validated to an exact match against scikit-learn 1.9.0**, and the trade engine's numerical **kernels** (ADEV NNLS fit, χ² consistency bands, van-Loan clock Q) **against scipy 1.17.1**; the device-benefit numbers built on top stay **MODELLED** operating characteristics — never field/IQ data, no good/bad verdict. |
|
|
173
|
-
| **Frugal engineering & integrity impact** | A **cost-per-coverage ROI** lens (`src/frugal.rs`) — cost per unit of delivered coverage for an architecture trade — and a **detection-miss → integrity-impact** mapping (`src/integrity_impact.rs`) that turns a monitor's missed-detection rate into its integrity-risk contribution. **MODELLED** decision-support budgets, additive. |
|
|
174
|
-
| **Artifact interchange** | The **Kshana Interchange Format (KIF)** (`src/interchange.rs`) — a versioned, self-describing envelope wrapping a scenario result with its kind, schema version, and MODELLED/VALIDATED labels, so a stored artifact stays self-documenting and older envelopes remain forward-compatibly readable. |
|
|
175
|
-
|
|
176
|
-
Each capability is reachable as a Rust API, a runnable scenario `kind`, or both.
|
|
177
|
-
Maturity per capability — *validated*, *runnable*, or *library* — is tracked in
|
|
178
|
-
[`docs/CAPABILITY.md`](docs/CAPABILITY.md). A **machine-checked verification matrix**
|
|
179
|
-
(`src/verification.rs`) renders the requirement → module → test → oracle → status
|
|
180
|
-
cross-reference, with unit-tested honesty invariants that permit a *validated* label
|
|
181
|
-
only where an independent **external** oracle backs it — and that record the
|
|
182
|
-
hardware/PA capabilities Kshana deliberately does **not** provide.
|
|
183
|
-
|
|
184
|
-
## Results
|
|
185
|
-
|
|
186
|
-
Each scenario compares a quantum sensor against its classical counterpart through a
|
|
187
|
-
~1.8 h GNSS outage. Numbers are reproducible (`scenario + seed + version`).
|
|
32
|
+
> ***Validated, not asserted.*** 666/666 AIAA SGP4 vectors to **4.12 mm** · Cowell
|
|
33
|
+
> force model **0.08 m** vs Orekit 12.2 · Galileo **0.61 m** / Swarm-A **0.10 m** vs
|
|
34
|
+
> real ESA precise ephemerides · GCRS→ITRS bit-for-bit vs SOFA/ERFA · ML metrics exact
|
|
35
|
+
> vs scikit-learn · **40 of 91** capabilities validated against independent external
|
|
36
|
+
> oracles; 47 honestly labelled Modelled, 4 partner-owned.
|
|
188
37
|
|
|
189
38
|
<p align="center">
|
|
190
|
-
<img src="docs/assets/
|
|
191
|
-
<br><em>Dead-reckoning position error during a GNSS outage: the quantum sensor (blue)
|
|
192
|
-
stays flat near the spec; the classical sensor (red) diverges to tens of kilometres.
|
|
193
|
-
Generated by Kshana from <code>scenarios/imu-deadreckoning.toml</code>.</em>
|
|
39
|
+
<img src="https://raw.githubusercontent.com/AshfordeOU/kshana/main/docs/assets/diagrams/system-overview.png" alt="Kshana system overview: five front doors (CLI, Python wheel, WebAssembly playground, MCP server, JetBrains plugin) converge on a single api::run_toml dispatch, through the engine, to a reproducible result.json + chart.svg" width="840">
|
|
194
40
|
</p>
|
|
195
41
|
|
|
196
|
-
|
|
197
|
-
|------|----------|---------|-----------|
|
|
198
|
-
| **1 — Clock holdover** | `clock-holdover.toml` (20 ns spec) | optical clock holds the full outage | CSAC breaches the spec mid-outage |
|
|
199
|
-
| **2 — Inertial dead-reckoning** | `imu-deadreckoning.toml` (100 m spec) | cold-atom: **~41 m**, holds full outage | nav-grade: breaches in **~350 s** → tens of km |
|
|
200
|
-
| **3 — Time transfer** (optical inter-satellite link) | `timetransfer.toml` | optical: **~0.3 mm** ranging | RF (TWSTFT): **~150 mm** ranging |
|
|
201
|
-
| **4 — Hybrid fusion** (capstone) | `hybrid-pnt.toml` | full position+timing for the whole outage | **position-limited at ~350 s** |
|
|
202
|
-
|
|
203
|
-
The capstone shows the fusion thesis: optical inter-satellite time-transfer keeps even
|
|
204
|
-
a classical *clock* locked, isolating the *inertial* sensor as the classical suite's
|
|
205
|
-
weak link — i.e. quantum inertial + optical timing together.
|
|
206
|
-
|
|
207
|
-
<p align="center">
|
|
208
|
-
<img src="docs/assets/clock-holdover.svg" alt="Clock holdover: phase error during a GNSS outage — the optical clock stays within the 20 ns spec for the whole outage while the chip-scale clock breaches it mid-outage" width="80%">
|
|
209
|
-
<br><em>Clock holdover through a GNSS outage: the optical clock (blue) stays inside the
|
|
210
|
-
20 ns spec for the full coast; the chip-scale clock (red) breaches it part-way.
|
|
211
|
-
Generated by Kshana from <code>scenarios/clock-holdover.toml</code>.</em>
|
|
212
|
-
</p>
|
|
42
|
+
### Validated against external oracles — every row CI-gated
|
|
213
43
|
|
|
214
|
-
|
|
215
|
-
|
|
216
|
-
|
|
217
|
-
|
|
218
|
-
|
|
219
|
-
|
|
220
|
-
|
|
44
|
+
| | Capability | Result | External oracle |
|
|
45
|
+
|---|---|---|---|
|
|
46
|
+
| ✅ | SGP4/SDP4 propagation | 666/666 vectors, worst **4.12 mm** | AIAA 2006-6753 (Vallado) + independent `sgp4` crate |
|
|
47
|
+
| ✅ | Numerical Cowell force model | **0.08 m** / 24 h, 275 epochs | Orekit 12.2 `DormandPrince853` (CS GROUP) |
|
|
48
|
+
| ✅ | Orbit fit vs precise ephemeris | Galileo **0.61 m** · Swarm-A **0.10 m** | ESA/ESOC SP3 precise orbits |
|
|
49
|
+
| ✅ | GCRS→ITRS frame chain | bit-for-bit vs SOFA; ≤ 0.86 m vs SPICE | ERFA/SOFA + ANISE (pure-Rust SPICE) |
|
|
50
|
+
| ✅ | Allan deviations | reproduce reference deviations | NIST SP 1065 + Stable32 on a real Cs clock |
|
|
51
|
+
| ✅ | GNSS DOP · ML detector metrics | to **1e-6** · to **1e-9** | gnss_lib_py · scikit-learn |
|
|
221
52
|
|
|
222
53
|
<p align="center">
|
|
223
|
-
<img src="docs/assets/
|
|
224
|
-
<br><em>GNSS availability derived from orbital geometry: the visible-satellite count
|
|
225
|
-
(line-of-sight, Earth-occultation, elevation mask) sets the fix state at each step,
|
|
226
|
-
so the clock must coast every gap. Generated by Kshana from
|
|
227
|
-
<code>scenarios/orbit-gnss-challenged.toml</code>.</em>
|
|
54
|
+
<img src="https://raw.githubusercontent.com/AshfordeOU/kshana/main/docs/assets/figures/validation-breakdown.png" alt="Verification status across all 91 capabilities: 40 Validated, 47 Modelled, 4 Partner-owned" width="780">
|
|
228
55
|
</p>
|
|
229
56
|
|
|
230
|
-
|
|
231
|
-
(line 1 + line 2) is propagated with the full **SGP4/SDP4** model — including
|
|
232
|
-
atmospheric drag and the deep-space lunar-solar and 12 h / 24 h resonance terms that
|
|
233
|
-
matter for ~12 h GNSS orbits — validated against the official AIAA 2006-6753 vectors
|
|
234
|
-
to a worst-case ≈ 4 mm. `scenarios/orbit-sgp4-gps.toml` ships a **real Celestrak
|
|
235
|
-
`gps-ops` snapshot** of the operational GPS constellation (2021-07-28, 30 satellites)
|
|
236
|
-
and requires valid TLE checksums — two-line element sets are open data from the US
|
|
237
|
-
Space Force / 18th Space Defense Squadron catalogue, redistributed by Celestrak
|
|
238
|
-
(Dr T. S. Kelso, [celestrak.org](https://celestrak.org)); refresh with
|
|
239
|
-
`scripts/fetch_tles.sh`. A line-2-only block keeps
|
|
240
|
-
the analytic two-body propagation (`scenarios/orbit-real-tle.toml`); the two forms can
|
|
241
|
-
be mixed in one constellation. A constellation can equally be built from a block of
|
|
242
|
-
**RINEX-3 GPS broadcast-ephemeris** records — the format a receiver decodes —
|
|
243
|
-
propagated by the IS-GPS-200 user algorithm and fed through the same geometry
|
|
244
|
-
(`scenarios/orbit-rinex.toml`).
|
|
245
|
-
|
|
246
|
-
## Install & build
|
|
247
|
-
|
|
248
|
-
Requires a Rust toolchain (≥ 1.75; developed on 1.93).
|
|
57
|
+
## Install
|
|
249
58
|
|
|
250
59
|
```bash
|
|
251
|
-
|
|
252
|
-
cd kshana
|
|
253
|
-
cargo build --release
|
|
254
|
-
cargo test # all tests pass
|
|
60
|
+
npm install kshana
|
|
255
61
|
```
|
|
256
62
|
|
|
257
63
|
## Usage
|
|
258
64
|
|
|
259
|
-
|
|
260
|
-
|
|
261
|
-
|
|
262
|
-
```bash
|
|
263
|
-
cargo run -- scenarios/clock-holdover.toml
|
|
264
|
-
cargo run -- scenarios/imu-deadreckoning.toml
|
|
265
|
-
cargo run -- scenarios/timetransfer.toml
|
|
266
|
-
cargo run -- scenarios/hybrid-pnt.toml
|
|
267
|
-
cargo run -- scenarios/orbit-gnss-challenged.toml
|
|
268
|
-
cargo run -- scenarios/orbit-sgp4-gps.toml
|
|
269
|
-
cargo run -- scenarios/orbit-rinex.toml
|
|
270
|
-
cargo run -- scenarios/integrity-raim.toml
|
|
271
|
-
|
|
272
|
-
# Export a propagated constellation to an SP3-c precise-ephemeris file:
|
|
273
|
-
cargo run -- scenarios/orbit-sgp4-gps.toml --export-sp3 gps.sp3
|
|
274
|
-
|
|
275
|
-
# Export the constellation's mean elements to a CCSDS OMM catalogue (one OMM
|
|
276
|
-
# message per TLE-defined satellite, with its real NORAD id / COSPAR designator):
|
|
277
|
-
cargo run -- scenarios/orbit-sgp4-gps.toml --export-omm gps.omm
|
|
278
|
-
|
|
279
|
-
# Export the velocity-carrying state to a CCSDS OEM 2.0 ephemeris (GMAT/Orekit/STK):
|
|
280
|
-
cargo run -- scenarios/orbit-sgp4-gps.toml --export-oem gps.oem
|
|
281
|
-
```
|
|
282
|
-
|
|
283
|
-
**Interoperability role.** Kshana is the *performance-simulation* layer that sits
|
|
284
|
-
alongside the post-processing toolchain, not a replacement for it: feed its **RINEX**
|
|
285
|
-
output into RTKLIB or gLAB for a position solution, and use its **SP3** output as a
|
|
286
|
-
precise-orbit product for tools like Ginan — Kshana answers *what resilience a given
|
|
287
|
-
PNT architecture buys* before you have real signals, in formats those tools already
|
|
288
|
-
ingest (`--export-sp3`, or `export_sp3 = true` in an `orbit` scenario, writes
|
|
289
|
-
`<scenario>.sp3`). The same orbit can be published as standards-track **CCSDS OMM**
|
|
290
|
-
mean elements (`--export-omm`, or `export_omm = true`, writes `<scenario>.omm`) —
|
|
291
|
-
one OMM 502.0 KVN message per TLE-defined satellite, carrying each object's real
|
|
292
|
-
NORAD catalogue number, COSPAR international designator, and epoch, for any
|
|
293
|
-
OMM-aware consumer instead of a bespoke two-line element set.
|
|
294
|
-
|
|
295
|
-
Example output (clock holdover — note the Integrity and Security figures of merit):
|
|
296
|
-
|
|
297
|
-
```
|
|
298
|
-
scenario c827e5d40d25 | quantum holdover 6600s p95 0.0ns integrity 1.000 security 0.997 | classical holdover 2610s p95 19.7ns integrity 1.000 security 0.000
|
|
299
|
-
wrote scenarios/clock-holdover.result.json and scenarios/clock-holdover.chart.svg
|
|
300
|
-
```
|
|
301
|
-
|
|
302
|
-
The optical clock's tight detection floor keeps `security 0.997`; the chip-scale
|
|
303
|
-
clock's own noise over the monitoring window exceeds the 20 ns spec, so it has no
|
|
304
|
-
spoof-detection margin (`security 0.000`). The orbit scenario additionally reports a
|
|
305
|
-
geometry block — fraction of samples with a fix, and best/median PDOP and position
|
|
306
|
-
accuracy — alongside the clock result.
|
|
307
|
-
|
|
308
|
-
> **Read these two numbers carefully.** `security` is an *analytic spoof-detectability
|
|
309
|
-
> bound* derived from each clock's stability — it is meaningful only against a
|
|
310
|
-
> configured spoofing scenario and is **not** a multi-satellite RAIM detector. `integrity`
|
|
311
|
-
> here is the filter's *self-consistency* (fraction of outage samples inside its own k-sigma
|
|
312
|
-
> bound), **not** an aviation HPL/VPL integrity figure. See
|
|
313
|
-
> [`docs/INTEGRITY.md`](docs/INTEGRITY.md).
|
|
314
|
-
>
|
|
315
|
-
> For genuine receiver-autonomous integrity, the **`integrity` scenario kind**
|
|
316
|
-
> (`scenarios/integrity-raim.toml`) runs real snapshot and solution-separation
|
|
317
|
-
> (ARAIM-style) RAIM over the propagated constellation geometry: it computes
|
|
318
|
-
> horizontal/vertical **protection levels (HPL/VPL)** per epoch and reports the
|
|
319
|
-
> fraction of epochs that meet the configured alert limits, with a Stanford
|
|
320
|
-
> integrity diagram for error-vs-PL classification.
|
|
321
|
-
|
|
322
|
-
### Python
|
|
323
|
-
|
|
324
|
-
An optional Python extension (PyO3, abi3) wraps the same engine. Build and install
|
|
325
|
-
it with [maturin](https://www.maturin.rs/):
|
|
326
|
-
|
|
327
|
-
```bash
|
|
328
|
-
pip install maturin
|
|
329
|
-
maturin develop --features python # or: maturin build --features python
|
|
330
|
-
```
|
|
331
|
-
|
|
332
|
-
```python
|
|
333
|
-
import json, kshana
|
|
334
|
-
|
|
335
|
-
result = json.loads(kshana.run(open("scenarios/clock-holdover.toml").read()))
|
|
336
|
-
print(result["quantum"]["fom"]["integrity"])
|
|
337
|
-
|
|
338
|
-
# json, svg, and a one-line summary at once:
|
|
339
|
-
result_json, chart_svg, summary = kshana.run_full(open("scenarios/orbit-gnss-challenged.toml").read())
|
|
340
|
-
print(kshana.version(), summary)
|
|
341
|
-
```
|
|
342
|
-
|
|
343
|
-
Wheels are built for Linux, macOS, and Windows by the `wheels` workflow on each
|
|
344
|
-
release tag.
|
|
345
|
-
|
|
346
|
-
### WebAssembly
|
|
347
|
-
|
|
348
|
-
The engine also runs in the browser via [wasm-pack](https://rustwasm.github.io/wasm-pack/):
|
|
349
|
-
|
|
350
|
-
```bash
|
|
351
|
-
wasm-pack build --target web -- --features wasm
|
|
352
|
-
```
|
|
65
|
+
The package is an ES module with a WebAssembly payload. Initialise it once, then call
|
|
66
|
+
the engine synchronously:
|
|
353
67
|
|
|
354
68
|
```js
|
|
355
|
-
import init, { run, chart_svg, version } from "
|
|
356
|
-
await init();
|
|
357
|
-
const result = JSON.parse(run(tomlText));
|
|
358
|
-
console.log(version(), result.classical.fom.timing_p95_ns);
|
|
359
|
-
```
|
|
360
|
-
|
|
361
|
-
### AI agents (MCP)
|
|
362
|
-
|
|
363
|
-
[](https://glama.ai/mcp/servers/ashfordeOU/kshana)
|
|
364
|
-
|
|
365
|
-
Kshana ships an [MCP](https://modelcontextprotocol.io) server, [`kshana-mcp`](mcp/kshana-mcp/),
|
|
366
|
-
so AI assistants and agents can run the **validated** engine instead of guessing the
|
|
367
|
-
math — usable from **Cursor, JetBrains AI Assistant / Junie, and any MCP-compatible
|
|
368
|
-
assistant or agent**. It exposes `run_scenario`, `list_scenario_kinds`,
|
|
369
|
-
`validate_scenario`, `export_sp3`, and `export_omm` (each a thin wrapper over
|
|
370
|
-
`kshana::api`).
|
|
371
|
-
|
|
372
|
-
```bash
|
|
373
|
-
cargo install kshana-mcp # crates.io
|
|
374
|
-
docker run --rm -i ghcr.io/ashfordeou/kshana-mcp # or OCI, no Rust toolchain
|
|
375
|
-
```
|
|
376
|
-
|
|
377
|
-
Then register `kshana-mcp` in your client's `mcpServers` config — see
|
|
378
|
-
[`mcp/kshana-mcp/README.md`](mcp/kshana-mcp/README.md) for per-client snippets. The
|
|
379
|
-
server is a standalone, workspace-excluded crate (the `rmcp` SDK is edition 2024), so it
|
|
380
|
-
never affects the lean published `kshana` crate or its build.
|
|
381
|
-
|
|
382
|
-
**In a JetBrains IDE** you can also install the
|
|
383
|
-
[**Kshana — PNT simulator**](https://plugins.jetbrains.com/plugin/32181-kshana--pnt-simulator)
|
|
384
|
-
plugin from the JetBrains Marketplace (or *Settings → Plugins → Marketplace → search
|
|
385
|
-
"Kshana"*) to run scenarios from a right-click — see [`ide/jetbrains/`](ide/jetbrains/).
|
|
386
|
-
|
|
387
|
-
## Scenario format
|
|
388
|
-
|
|
389
|
-
Scenarios are declarative TOML. A top-level `kind` selects the pack — **thirty-four** in
|
|
390
|
-
all (`clock` is the default if omitted): `inertial`, `timetransfer`, `hybrid`, `hybrid-ukf`, `fusion`,
|
|
391
|
-
`gnss-ins`, `orbit`, `ephemeris`, `gnss-sim`, `integrity`, `lunar-integrity`, `spoof`,
|
|
392
|
-
`spoof-detect`, `jamming`, `sweep`, `sweep-nd`, `gravity-map`, `terrain-nav`, `terrain-slam`,
|
|
393
|
-
`combined-altpnt`, `pvt`, `mars-pnt`, `impairment-eval` (AI/ML RF-impairment detection
|
|
394
|
-
evaluation testbed — labelled synthetic corpus + detector-agnostic ROC/AUC harness +
|
|
395
|
-
in/out-of-distribution optimism gap), `quantum-trade` (quantum-vs-classical PNT
|
|
396
|
-
trade with measured-ADEV ingestion + GNSS-denied resilience envelope; MODELLED),
|
|
397
|
-
`space-weather` (solar/geomagnetic indices + Jacchia-71 exospheric temperature +
|
|
398
|
-
activity-driven thermospheric density over the static atmosphere; MODELLED),
|
|
399
|
-
`oem-interop` (CCSDS OEM import/round-trip bridge for GMAT/Orekit/STK ephemerides;
|
|
400
|
-
MODELLED), the mission-analysis trio `launch-window` (two-body launch azimuth /
|
|
401
|
-
plane-change / opportunities), `reentry` (Allen-Eggers ballistic re-entry corridor),
|
|
402
|
-
`eo-coverage` (EO swath / GSD / access / revisit geometry), `space-packet` (CCSDS
|
|
403
|
-
133.0 TM/TC Space Packet framing — exact bit layout, round-trip verified), and
|
|
404
|
-
`attitude-budget` (3-DOF gravity-gradient torque + RSS pointing error budget),
|
|
405
|
-
`passes` (ground-station rise/set pass prediction — AOS/TCA/LOS, max elevation,
|
|
406
|
-
access), and `link-budget` (one-way CCSDS/DSN link equation — FSPL / Eb·N₀ /
|
|
407
|
-
margin / closure) — all MODELLED.
|
|
408
|
-
Common fields: `seed`, a `[time]` grid, a `[gnss]` availability timeline (the outage
|
|
409
|
-
driver), and per-sensor blocks with `provenance` strings citing the source of every
|
|
410
|
-
figure. Example (clock):
|
|
411
|
-
|
|
412
|
-
```toml
|
|
413
|
-
seed = 42
|
|
414
|
-
threshold_ns = 20.0
|
|
415
|
-
[time]
|
|
416
|
-
step_s = 10.0
|
|
417
|
-
duration_s = 7200.0
|
|
418
|
-
[gnss]
|
|
419
|
-
windows = [
|
|
420
|
-
{ t0 = 0.0, t1 = 600.0, state = "nominal" }, # 10 min GNSS sync
|
|
421
|
-
{ t0 = 600.0, t1 = 7200.0, state = "denied" }, # ~1.8 h outage
|
|
422
|
-
]
|
|
423
|
-
[clock_quantum]
|
|
424
|
-
id = "optical-sr-lattice"
|
|
425
|
-
provenance = "Strontium optical lattice clock, space-oriented goal sigma_y(1s)=1e-15 (arXiv:1503.08457)"
|
|
426
|
-
y0 = 5.0e-17
|
|
427
|
-
q_wf = 1.0e-30 # white FM: q_wf = sigma_y(1s)^2
|
|
428
|
-
q_rw = 0.0 # random-walk FM
|
|
429
|
-
drift = 0.0 # linear aging (per second)
|
|
430
|
-
[clock_classical]
|
|
431
|
-
id = "csac-sa45s"
|
|
432
|
-
provenance = "Microchip SA65 / SA.45s CSAC datasheet sigma_y(1s)=3e-10"
|
|
433
|
-
y0 = 5.0e-10
|
|
434
|
-
q_wf = 9.0e-20
|
|
435
|
-
q_rw = 0.0
|
|
436
|
-
drift = 0.0
|
|
437
|
-
```
|
|
438
|
-
|
|
439
|
-
Optional fields (off when absent): a clock may add `flicker_floor` (1/f FM Allan
|
|
440
|
-
floor); an inertial sensor may add `gyro_bias` and `q_arw` (gyro bias and angular
|
|
441
|
-
random walk), and `bias_instability` and `q_aa` (the Allan bias-instability floor and
|
|
442
|
-
acceleration random walk) — together a **single-axis (1-DOF) accelerometer error
|
|
443
|
-
budget** (VRW/ARW and bias-instability). This is the error budget the shipped
|
|
444
|
-
`inertial` scenario *pack* runs. Separately, the library now carries a verified
|
|
445
|
-
**3-axis strapdown navigator** (`src/inertial/{attitude,mechanization,imu_errors}.rs`):
|
|
446
|
-
quaternion attitude with coning/sculling compensation, a full NED mechanization
|
|
447
|
-
(Earth-rate and transport-rate terms, WGS-84 Somigliana gravity), and a
|
|
448
|
-
deterministic IMU error model in which **scale-factor, misalignment,
|
|
449
|
-
g-sensitivity, quantization, and rate-ramp are modelled** (IEEE Std 952-1997
|
|
450
|
-
§A.2; Groves 2013 §4.3). That 3-axis path is now **wired into a runnable
|
|
451
|
-
loosely-coupled GNSS/INS pack** (`kind = "gnss-ins"`): a 15-state error-state EKF
|
|
452
|
-
disciplines the strapdown solution against noisy fixes while GNSS is up, then
|
|
453
|
-
coasts through the outage, reporting the fused horizontal error against the
|
|
454
|
-
open-loop free-INS coast. A **tightly-coupled pseudorange** update is also
|
|
455
|
-
available (it forms the innovation in the range domain, so it keeps correcting
|
|
456
|
-
with fewer than four satellites). A
|
|
457
|
-
clock-holdover scenario may add `runs` (> 1) to run a **Monte Carlo ensemble** — each
|
|
458
|
-
figure of merit is then reported as a mean with a 5th–95th-percentile spread and the
|
|
459
|
-
chart shades the error confidence band (see `scenarios/clock-ensemble.toml`).
|
|
69
|
+
import init, { run, run_full, chart_svg, version } from "kshana";
|
|
460
70
|
|
|
461
|
-
|
|
462
|
-
— one for the clock state, one for the position state — disciplined by GNSS and aided by
|
|
463
|
-
optical time transfer, and reports a combined holdover FoM. The two blocks share no
|
|
464
|
-
cross-covariance: this is a stacked pair of error budgets, **not** a true coupled
|
|
465
|
-
clock+position joint filter (cross-block covariance is a roadmap item). See
|
|
466
|
-
`scenarios/fusion-pnt.toml`.
|
|
71
|
+
await init(); // load the wasm
|
|
467
72
|
|
|
468
|
-
|
|
469
|
-
|
|
470
|
-
|
|
471
|
-
two-sided **χ²₁ energy / Neyman–Pearson test** on the clock-aided monitor statistic:
|
|
472
|
-
the threshold is set from a target false-alarm budget `target_pfa`, and the
|
|
473
|
-
**missed-detection probability `P_md`** is reported both closed-form and by
|
|
474
|
-
Monte-Carlo (`mc_runs` trials per hypothesis — the two agree to a few ×1/√N). The
|
|
475
|
-
**Security figure of merit is `1 − P_md`** at the operationally-harmful (spec)
|
|
476
|
-
magnitude, so a quiet clock that catches a spec-sized spoof scores ≈ 1 and a noisy
|
|
477
|
-
one that often misses it scores lower (see `scenarios/spoof-attack.toml`,
|
|
478
|
-
`scenarios/spoof-meaconing.toml`).
|
|
479
|
-
|
|
480
|
-
A `gnss-sim` scenario is a **measurement-domain** simulation: for each visible
|
|
481
|
-
satellite it synthesises the pseudorange `ρ = geometric range + c·δt_rx − c·δt_sv +
|
|
482
|
-
I + T + noise + multipath` and the L1 Doppler, with the **Klobuchar** single-frequency
|
|
483
|
-
ionosphere (`[iono]`, IS-GPS-200 §20.3.3.5.2.5) and the **Saastamoinen** zenith
|
|
484
|
-
troposphere projected by the **Niell (1996)** mapping function (`[tropo]`). The
|
|
485
|
-
residuals feed **snapshot RAIM** for per-epoch HPL/VPL, and every satellite's
|
|
486
|
-
pseudorange, Doppler, C/N₀, and iono/tropo corrections are emitted in the JSON
|
|
487
|
-
`gnss_measurements` array. It is a forward simulator (it generates measurements from
|
|
488
|
-
a known truth), not a receiver/solver — a zero-noise run reproduces geometry plus the
|
|
489
|
-
corrections to sub-millimetre (see `scenarios/gnss-sim-raim.toml`).
|
|
490
|
-
|
|
491
|
-
A `jamming` scenario models RF interference as a **link budget**: a `[jammer]`
|
|
492
|
-
(ECEF position, transmit `power_dbw`, type) raises the jammer-to-signal ratio at a
|
|
493
|
-
`[receiver]` watching a Walker `[constellation]`. From the geometry (free-space
|
|
494
|
-
path loss and the per-direction receive-antenna gain) it computes each satellite's
|
|
495
|
-
`J/S`, the **effective C/N₀** via the standard anti-jam equation (despreading
|
|
496
|
-
processing gain × the spectral-separation factor `Q`; Kaplan & Hegarty §9.4), and
|
|
497
|
-
flags loss of lock below a configurable tracking threshold — reporting an
|
|
498
|
-
`availability_under_jamming` figure of merit. A 10 W broadband jammer at 1 km
|
|
499
|
-
denies the receiver entirely (J/S ≈ 72 dB); the same jammer at 100 km only
|
|
500
|
-
degrades the links (see `scenarios/jamming-demo.toml`).
|
|
501
|
-
|
|
502
|
-
A `sweep` scenario runs a **trade study**: it varies one `parameter` (`threshold_ns`,
|
|
503
|
-
`duration_s`, `quantum_q_wf`, or `classical_q_wf`) from `start` to `stop` over `steps`
|
|
504
|
-
points on a `lin` or `log` `scale`, records a `metric` (e.g. `holdover_s`) for both
|
|
505
|
-
clocks, and charts the two curves. The base scenario goes under `[base]` (see
|
|
506
|
-
`scenarios/sweep-clock-stability.toml`).
|
|
507
|
-
|
|
508
|
-
A `sweep-nd` scenario generalises this to **any pack and any number of axes**: it
|
|
509
|
-
varies dotted TOML keys of a `[base]` scenario (of any `kind`) over the Cartesian
|
|
510
|
-
product of `[[axes]]`, re-runs each grid node, and records `metrics` given as
|
|
511
|
-
dotted JSON paths into the result (e.g. `classical.fom.holdover_s`). It works for
|
|
512
|
-
every pack because it operates at the TOML/result boundary; native runs evaluate
|
|
513
|
-
the grid in parallel (no extra dependency, wasm falls back to sequential) and the
|
|
514
|
-
output is deterministic and row-major (see `scenarios/sweep-nd-inertial.toml`).
|
|
515
|
-
|
|
516
|
-
An `orbit` scenario derives the `[gnss]` timeline from geometry instead of authoring
|
|
517
|
-
it — give a `[user]` orbit, a `[constellation]`, an elevation `mask_deg`, and the two
|
|
518
|
-
clock blocks. It also reports position accuracy from the satellite geometry; the
|
|
519
|
-
optional `sigma_uere_m` (1-sigma user-equivalent range error, default 1 m) scales the
|
|
520
|
-
position dilution of precision into a position sigma. The user orbit may be made
|
|
521
|
-
**eccentric** with `eccentricity` and `argp_deg`, and `j2 = true` adds Earth-oblateness
|
|
522
|
-
secular drift (see `scenarios/orbit-molniya.toml`). The constellation can instead be a
|
|
523
|
-
**real one**: give `[constellation]` a `tle` block of two-line element sets and the
|
|
524
|
-
satellites are parsed from it (see `scenarios/orbit-real-tle.toml`). Add one or more
|
|
525
|
-
`[[constellations]]` blocks for **multi-GNSS** (e.g. GPS + Galileo; see
|
|
526
|
-
`scenarios/orbit-multignss.toml`):
|
|
527
|
-
|
|
528
|
-
```toml
|
|
529
|
-
kind = "orbit"
|
|
530
|
-
seed = 7
|
|
531
|
-
threshold_ns = 5.0
|
|
532
|
-
mask_deg = 10.0
|
|
533
|
-
sigma_uere_m = 1.0 # optional; position sigma = position-DOP * this
|
|
534
|
-
[time]
|
|
535
|
-
step_s = 60.0
|
|
536
|
-
duration_s = 86400.0
|
|
537
|
-
[user] # spacecraft (altitude in km, angles in deg)
|
|
538
|
-
altitude_km = 8000.0
|
|
539
|
-
inclination_deg = 0.0
|
|
540
|
-
[constellation] # Walker-delta GNSS (GPS-like)
|
|
541
|
-
altitude_km = 20180.0
|
|
542
|
-
inclination_deg = 55.0
|
|
543
|
-
planes = 6
|
|
544
|
-
sats_per_plane = 4
|
|
545
|
-
phasing_f = 1.0
|
|
546
|
-
[clock_quantum] # ... as above
|
|
547
|
-
[clock_classical] # ... as above
|
|
548
|
-
```
|
|
549
|
-
|
|
550
|
-
The **GPS-denied alt-PNT** kinds navigate with no GNSS at all, matching a measured field
|
|
551
|
-
sequence against a map through a particle filter. A `gravity-map` scenario flies a track
|
|
552
|
-
through a spherical-harmonic gravity-anomaly field and recovers it from a cold-atom
|
|
553
|
-
gravimeter's reading (`scenarios/gps-denied-gravity-nav.toml`); a `terrain-nav` scenario
|
|
554
|
-
does the same against an SRTM elevation DEM (TERCOM/SITAN, `scenarios/terrain-nav.toml`);
|
|
555
|
-
and a `combined-altpnt` scenario fuses **gravity + IGRF magnetic + terrain** in one filter
|
|
556
|
-
(`scenarios/combined-altpnt.toml`).
|
|
557
|
-
|
|
558
|
-
A `lunar-integrity` scenario evaluates **cislunar** PNT: it runs a lunar south-pole
|
|
559
|
-
ARAIM protection-level pass against a LunaNet/LNIS relay set and honestly reports the
|
|
560
|
-
integrity gap — a ~30 m lunar σ_URE drives the protection level well above a 50 m alert
|
|
561
|
-
limit, so the service is *unavailable* under aviation-style integrity rules
|
|
562
|
-
(`scenarios/lunanet-araim.toml`).
|
|
563
|
-
|
|
564
|
-
See `scenarios/` for one example of every kind.
|
|
565
|
-
|
|
566
|
-
## Output
|
|
567
|
-
|
|
568
|
-
The result artifact is versioned, self-describing JSON: per-step time series, the
|
|
569
|
-
scored figures of merit, the active model specs (with provenance), the seed, a
|
|
570
|
-
**scenario hash** — so any chart can be reproduced from the file — and, for each clock,
|
|
571
|
-
an `adev_curve` (`[{tau_s, adev, n_samples, noise, edf, ci_lo, ci_hi}]`): the overlapping
|
|
572
|
-
Allan deviation across octave-spaced averaging times — the standard way to read a clock's
|
|
573
|
-
stability — now with a **noise-type-specific 95% confidence band** per point (the record's
|
|
574
|
-
power-law type is identified from its modified-Allan slope, and the χ² interval uses the
|
|
575
|
-
matching NIST SP 1065 effective degrees of freedom). The browser playground renders it as a
|
|
576
|
-
log-log "Clock stability (ADEV)" chart. (MDEV, TDEV, and HDEV are available as library
|
|
577
|
-
estimators; the exported result curve is the overlapping ADEV.) Every field, with units and a
|
|
578
|
-
source pointer, is documented in [`docs/SCHEMA.md`](docs/SCHEMA.md).
|
|
579
|
-
|
|
580
|
-
**Every chart is self-describing.** The browser playground, the CLI's `*.chart.svg`
|
|
581
|
-
export, and the HTML scorecard all stamp each chart image with a footer reading
|
|
582
|
-
`Kshana v<version> · scenario <hash> · kshana.dev`. The `scenario <hash>` is the first
|
|
583
|
-
12 hex characters of the run's **scenario hash** — a SHA-256 over the canonical scenario
|
|
584
|
-
definition (seed, thresholds, model parameters, GNSS windows, …); the integrity and lunar
|
|
585
|
-
reports, which carry no hash of their own, fall back to a SHA-256 of the scenario source.
|
|
586
|
-
It is the **same fingerprint** shown in the one-line summary and the result JSON, so a
|
|
587
|
-
saved or pasted chart always carries its version, the exact scenario that produced it (for
|
|
588
|
-
bit-for-bit reproduction), and the source — change any input and the hash changes.
|
|
589
|
-
|
|
590
|
-
The figures of merit follow the standard operational PNT figures of merit:
|
|
591
|
-
|
|
592
|
-
| Figure of merit | How Kshana computes it |
|
|
593
|
-
|-----------------|------------------------|
|
|
594
|
-
| Timing Performance (clock/orbit packs) | clock-phase error RMS + 95th-percentile over the outage, in **nanoseconds** (`timing_rms_ns`) — a timing metric, not position |
|
|
595
|
-
| Positioning Performance (inertial/hybrid packs) | 1-DOF position-error RMS + 95th-percentile over the outage, in **metres** (`pos_rms_m`); single-axis. A single run is flagged `monte_carlo: false`; set `runs = N` for a Monte Carlo ensemble that reports each metric's mean, spread, and bootstrap 95% CI. Still **not** a 2-D CEP/2DRMS or DOP-weighted accuracy (those need the 3-axis model — roadmap) |
|
|
596
|
-
| Autonomy | holdover duration — time in-spec after GNSS loss (grid-quantised: a lower bound) |
|
|
597
|
-
| Resilience | error-growth slope during the outage |
|
|
598
|
-
| Availability | fraction of the run with an in-spec solution |
|
|
599
|
-
| Integrity | filter **self-consistency** — fraction of outage samples whose error stays inside the Kalman filter's own k-sigma bound. **Not** an aviation HPL/VPL/RAIM integrity figure (see [`docs/INTEGRITY.md`](docs/INTEGRITY.md)) |
|
|
600
|
-
| Security | **analytic spoof-*detectability* bound** from clock stability — how small/slow a time-spoof a single-clock consistency monitor could flag. Meaningful only with a configured attack; **not** a multi-satellite RAIM detector |
|
|
601
|
-
|
|
602
|
-
New to these terms? Each is defined in plain language in the [glossary](docs/GLOSSARY.md).
|
|
603
|
-
|
|
604
|
-
## Architecture
|
|
605
|
-
|
|
606
|
-
**One engine, many front doors.** A single Rust core (`kshana`) runs every scenario,
|
|
607
|
-
reached through a CLI, a Python extension, an in-browser WebAssembly module, an **MCP
|
|
608
|
-
server** for AI agents, and a **JetBrains IDE plugin** — all converging on one
|
|
609
|
-
`api::run_toml` dispatch. Inside, the sensor packs plug into a common error-model
|
|
610
|
-
interface; alongside them sit a **reference-frame layer** (IAU 2006/2000A
|
|
611
|
-
precession–nutation and the CIO-based GCRS↔ITRS reduction), an **astrodynamics/numerical
|
|
612
|
-
layer** (analytic SGP4/SDP4 **and** a numerical Cowell propagator with its
|
|
613
|
-
EGM2008/perturbation force model, maneuver design, and orbit determination), an
|
|
614
|
-
**integrity/GNSS layer** (RAIM/ARAIM, SBAS, the measurement domain, jamming, cislunar),
|
|
615
|
-
and a **fusion / alt-PNT layer** (the GNSS/INS estimators and the gravity/terrain/magnetic
|
|
616
|
-
map-matchers).
|
|
617
|
-
|
|
618
|
-
Two standalone, **workspace-excluded** crates sit beside the core — `mcp/kshana-mcp`
|
|
619
|
-
(the MCP server, built on the edition-2024 `rmcp` SDK) and `xval/anise-frames` (the
|
|
620
|
-
ANISE/SPICE frame cross-check, which pulls MPL-2.0 deps) — kept out of the published
|
|
621
|
-
crate's dependency graph, `Cargo.lock`, license gate, and MSRV build by the root
|
|
622
|
-
`Cargo.toml` `exclude` list. The JetBrains plugin (`ide/jetbrains`) is a separate Kotlin
|
|
623
|
-
project. See [`docs/ARCHITECTURE.md`](docs/ARCHITECTURE.md) for the full set of diagrams.
|
|
624
|
-
|
|
625
|
-
```mermaid
|
|
626
|
-
flowchart LR
|
|
627
|
-
SCN["Scenario (.toml)<br/>seed · GNSS timeline · sensor params"] --> ENG
|
|
628
|
-
subgraph ENG["Engine (per step)"]
|
|
629
|
-
direction TB
|
|
630
|
-
M["Error model<br/>step(): evolve noise state"] --> E["Estimator<br/>GNSS-disciplined holdover"]
|
|
631
|
-
E --> F["FoM scoring<br/>vs the 6 figures of merit"]
|
|
632
|
-
end
|
|
633
|
-
ENG --> OUT["result.json + chart.svg<br/>(reproducible: scenario+seed+version)"]
|
|
634
|
-
```
|
|
635
|
-
|
|
636
|
-
```mermaid
|
|
637
|
-
flowchart TD
|
|
638
|
-
cli["CLI · Python · WebAssembly · MCP server · JetBrains plugin"] --> api["api — run_toml: typed dispatch over 34 kinds"]
|
|
639
|
-
subgraph shared["Shared core"]
|
|
640
|
-
types["types · scenario · GNSS timeline"]
|
|
641
|
-
allan["allan — ADEV/MDEV/TDEV/HDEV"]
|
|
642
|
-
end
|
|
643
|
-
subgraph frames["Time & reference frames"]
|
|
644
|
-
ts["timescales · jd2 — UTC/TAI/TT/UT1"]
|
|
645
|
-
pn["precession · nutation — IAU 2006/2000A"]
|
|
646
|
-
cio["cio — GCRS↔ITRS (CIO, SOFA-anchored)"]
|
|
647
|
-
end
|
|
648
|
-
subgraph packs["Sensor packs"]
|
|
649
|
-
p1["clock — models · estimator · kalman · security"]
|
|
650
|
-
p2["inertial — strapdown INS + quantum-CAI"]
|
|
651
|
-
p3["timetransfer — optical/RF/TWSTFT/PPP"]
|
|
652
|
-
p4["hybrid — fused PNT suite"]
|
|
653
|
-
end
|
|
654
|
-
subgraph astro["Astrodynamics & numerical"]
|
|
655
|
-
orbit["orbit · walker — geometry → GNSS + DOP"]
|
|
656
|
-
sgp4["sgp4 · tle — SGP4/SDP4"]
|
|
657
|
-
prop["propagator — Cowell"]
|
|
658
|
-
forces["forces — J2–J6 · 3rd-body · SRP · drag · GR"]
|
|
659
|
-
gsh["gravity_sh — EGM2008 d/o 70"]
|
|
660
|
-
integ["integrator — RK4 · DOPRI"]
|
|
661
|
-
man["maneuver · orbit_determination"]
|
|
662
|
-
cr["cr3bp — Earth–Moon CR3BP + halo/NRHO corrector"]
|
|
663
|
-
end
|
|
664
|
-
subgraph intg["Integrity & GNSS"]
|
|
665
|
-
raim["raim — RAIM/ARAIM · HPL/VPL"]
|
|
666
|
-
sbas["sbas — DO-229E PL · L1/L5"]
|
|
667
|
-
gsim["gnss_sim · ionex — measurement domain"]
|
|
668
|
-
jam["jamming — J/S → C/N₀"]
|
|
669
|
-
nsig["navsignal — PSD · SSC → anti-jam Q · DLL jitter"]
|
|
670
|
-
lun["lunar — cislunar ARAIM"]
|
|
671
|
-
end
|
|
672
|
-
subgraph fnav["Fusion & alt-PNT"]
|
|
673
|
-
fus["fusion — EKF · UKF · 17-state · coupled"]
|
|
674
|
-
grav["gravimeter · mapmatch · particle_filter"]
|
|
675
|
-
terr["altpnt/terrain · igrf — terrain + magnetic"]
|
|
676
|
-
end
|
|
677
|
-
subgraph ds["Deep-space & Mars"]
|
|
678
|
-
dsr["radiometric · ccsds_tdm — light-time · Δ-DOR · TDM"]
|
|
679
|
-
dso["deepspace_od — reduced-dynamic SRIF"]
|
|
680
|
-
dsm["mars_pnt · gse_sim — relay-PNT + GSE sim"]
|
|
681
|
-
end
|
|
682
|
-
subgraph io["Interop formats"]
|
|
683
|
-
iofmt["rinex · sp3 · oem · omm · glonass · ccsds_tdm"]
|
|
684
|
-
end
|
|
685
|
-
api --> packs
|
|
686
|
-
api --> astro
|
|
687
|
-
api --> intg
|
|
688
|
-
api --> fnav
|
|
689
|
-
api --> ds
|
|
690
|
-
packs --> shared
|
|
691
|
-
astro --> frames
|
|
692
|
-
orbit --> sgp4
|
|
693
|
-
prop --> forces
|
|
694
|
-
prop --> gsh
|
|
695
|
-
prop --> integ
|
|
696
|
-
cr --> integ
|
|
697
|
-
nsig --> jam
|
|
698
|
-
fus --> p2
|
|
699
|
-
grav --> p2
|
|
700
|
-
terr --> grav
|
|
701
|
-
orbit --> p1
|
|
702
|
-
orbit --> io
|
|
703
|
-
p4 -. composes .-> p1
|
|
704
|
-
p4 -. composes .-> p2
|
|
705
|
-
p4 -. composes .-> p3
|
|
706
|
-
```
|
|
707
|
-
|
|
708
|
-
**Components & distribution.** The core crate ships through the Rust, Python, and
|
|
709
|
-
JavaScript ecosystems; the MCP server and IDE plugin reach AI agents and JetBrains IDEs.
|
|
710
|
-
Each `vX.Y.Z` tag republishes every channel automatically (see
|
|
711
|
-
[Versioning & releases](#versioning--releases)).
|
|
712
|
-
|
|
713
|
-
```mermaid
|
|
714
|
-
flowchart LR
|
|
715
|
-
subgraph repo["One repository"]
|
|
716
|
-
core["kshana core<br/>library + CLI"]
|
|
717
|
-
mcp["mcp/kshana-mcp<br/>MCP server (excluded crate)"]
|
|
718
|
-
ide["ide/jetbrains<br/>Kotlin IDE plugin"]
|
|
719
|
-
xval["xval/anise-{frames,lunar-od,mars-od}<br/>SPICE/DE440 cross-checks (excluded)"]
|
|
720
|
-
end
|
|
721
|
-
core --> crates["crates.io"]
|
|
722
|
-
core --> pypi["PyPI — wheels"]
|
|
723
|
-
core --> npm["npm — WebAssembly"]
|
|
724
|
-
core --> rel["GitHub Releases<br/>binaries · SBOM · SLSA · validation summary"]
|
|
725
|
-
core --> pages["kshana.dev<br/>GitHub Pages playground"]
|
|
726
|
-
core -. archived .-> zen["Zenodo DOI"]
|
|
727
|
-
mcp --> crates
|
|
728
|
-
mcp --> ghcr["ghcr.io — OCI image"]
|
|
729
|
-
mcp --> reg["official MCP registry"]
|
|
730
|
-
ide --> jb["JetBrains Marketplace"]
|
|
731
|
-
```
|
|
732
|
-
|
|
733
|
-
## Repository layout
|
|
73
|
+
const toml = `kind = "clock_holdover"\n# ... scenario fields ...`;
|
|
74
|
+
const result = JSON.parse(run(toml));
|
|
75
|
+
console.log(version(), result.classical.fom.timing_p95_ns);
|
|
734
76
|
|
|
77
|
+
// JSON result + SVG chart in one call:
|
|
78
|
+
const { json, svg } = run_full(toml);
|
|
735
79
|
```
|
|
736
|
-
kshana/
|
|
737
|
-
├── src/ # the kshana core crate (library + CLI)
|
|
738
|
-
│ ├── api.rs · main.rs · lib.rs # typed dispatch (34 kinds) + CLI + crate root
|
|
739
|
-
│ ├── python.rs · wasm.rs # optional PyO3 / wasm-bindgen bindings
|
|
740
|
-
│ ├── types.rs · scenario.rs · allan.rs # shared core (time grid, GNSS timeline, Allan)
|
|
741
|
-
│ │
|
|
742
|
-
│ ├── models.rs · estimator.rs · kalman.rs # Pack 1 — clock holdover + integrity
|
|
743
|
-
│ ├── security.rs · detection.rs · spoof.rs · spoof_monitors.rs # spoof detection
|
|
744
|
-
│ ├── filter_health.rs · fom.rs · report.rs · chart.rs · run.rs # health · scoring · output
|
|
745
|
-
│ ├── inertial/ # Pack 2 — strapdown INS (attitude · mechanization · imu_errors · quantum_imu)
|
|
746
|
-
│ ├── timetransfer.rs · timetransfer_adv.rs · timegeo.rs # Pack 3 — TWSTFT/CV/PPP/optical, Sagnac
|
|
747
|
-
│ ├── hybrid.rs · ensemble.rs · sweep.rs # Pack 4 — fused PNT, Monte-Carlo, trade sweeps
|
|
748
|
-
│ │
|
|
749
|
-
│ ├── timescales.rs · jd2.rs · ephem.rs # time systems, two-part JD, Sun/Moon ephemeris
|
|
750
|
-
│ ├── precession.rs · nutation.rs · cio.rs # IAU 2006/2000A precession-nutation + CIO GCRS↔ITRS
|
|
751
|
-
│ ├── frames.rs · *_data.rs # TEME↔ECEF + generated nutation/CIO/EGM2008/IGRF tables
|
|
752
|
-
│ │
|
|
753
|
-
│ ├── orbit.rs · sgp4.rs · tle.rs · walker.rs # geometry, SGP4/SDP4, TLE, Walker design
|
|
754
|
-
│ ├── propagator.rs · forces.rs · gravity_sh.rs · integrator.rs # Cowell + perturbations (EGM2008 d/o70, GR) + RK4/DOPRI
|
|
755
|
-
│ ├── maneuver.rs · batch_ls.rs · orbit_determination.rs # burns/Lambert/porkchop, Gauss-Newton, OD
|
|
756
|
-
│ ├── cr3bp.rs · lunar.rs # Earth–Moon CR3BP + halo/NRHO STM corrector, cislunar/LunaNet ARAIM
|
|
757
|
-
│ ├── body.rs · mars_frame.rs · ephem_provider.rs · radiometric.rs · ccsds_tdm.rs # deep-space: multi-body · Mars frame · ephemeris seam · radiometric obs + CCSDS-TDM
|
|
758
|
-
│ ├── deepspace_od.rs · clock_state.rs · mars_atmos.rs · mars_pnt.rs · linkbudget.rs · gse_sim.rs # SRIF OD · onboard clock · Mars drag · relay-PNT · link budget · GSE sim
|
|
759
|
-
│ │
|
|
760
|
-
│ ├── fusion/ # GNSS/INS — EKF · UKF · tightly_coupled(17) · coupled · closed_loop
|
|
761
|
-
│ ├── raim.rs · sbas.rs # RAIM/ARAIM HPL/VPL, SBAS DO-229E PLs + L1/L5 iono-free
|
|
762
|
-
│ ├── gnss_sim.rs · ionex.rs · pvt.rs · jamming.rs # measurement domain · ionosphere maps · single-point positioning · jamming
|
|
763
|
-
│ ├── navsignal.rs # nav-signal PSD (BPSK-R/BOC) · spectral-separation → anti-jam Q · DLL code-tracking jitter · multipath envelope
|
|
764
|
-
│ ├── gravimeter.rs · igrf.rs · mapmatch.rs · particle_filter.rs · altpnt/ # gravity/magnetic/terrain alt-PNT
|
|
765
|
-
│ ├── rinex.rs · rinex_obs.rs · glonass.rs · sp3.rs · oem.rs · omm.rs · permalink.rs # interop formats
|
|
766
|
-
│ ├── launch.rs · reentry.rs · eo_payload.rs · attitude_budget.rs · passes.rs · space_packet.rs # mission-analysis budgets + CCSDS Space Packet
|
|
767
|
-
│ ├── space_weather.rs · holdover.rs # space-weather environment · GNSS-denied clock-holdover calculator
|
|
768
|
-
│ ├── impairment_eval.rs · quantum_trade.rs · frugal.rs · integrity_impact.rs # AI/ML RF-impairment eval · PNT trade · cost-per-coverage ROI · integrity impact
|
|
769
|
-
│ ├── interchange.rs · verification.rs # KIF artifact envelope · machine-checked verification matrix
|
|
770
|
-
│ └── bin/validation_report.rs # builds the release validation-summary HTML
|
|
771
|
-
│
|
|
772
|
-
├── mcp/kshana-mcp/ # standalone, workspace-EXCLUDED crate — the MCP server (+ Dockerfile, server.json)
|
|
773
|
-
├── ide/jetbrains/ # standalone Kotlin/Gradle IntelliJ-Platform plugin
|
|
774
|
-
├── xval/anise-{frames,lunar-od,mars-od}/ # standalone, workspace-EXCLUDED ANISE/SPICE cross-checks (frames · lunar DE440 · Mars DE440)
|
|
775
|
-
│
|
|
776
|
-
├── scenarios/ # one cited .toml per kind + geometry-driven + GPS-denied
|
|
777
|
-
├── scripts/ # reproducibility + repo-hygiene + SBOM guards
|
|
778
|
-
├── docs/ # CONCEPTS, ARCHITECTURE, CAPABILITY, VALIDATION, PROVENANCE, GLOSSARY, …
|
|
779
|
-
├── web/ # the WebAssembly playground + kshana.dev site
|
|
780
|
-
├── tools/ # table generators (EGM2008 · IGRF · nutation · CIO) + fetch_tles.sh
|
|
781
|
-
├── .github/workflows/ # ci · release · publish · wheels · pages · mcp-publish · jetbrains-plugin · frame-xval
|
|
782
|
-
├── pyproject.toml # Python packaging (maturin)
|
|
783
|
-
├── CHANGELOG.md # Keep a Changelog + SemVer
|
|
784
|
-
└── CITATION.cff · ROADMAP.md · CONTRIBUTING.md · SECURITY.md
|
|
785
|
-
```
|
|
786
|
-
|
|
787
|
-
## Documentation
|
|
788
|
-
|
|
789
|
-
| Document | For whom | What's in it |
|
|
790
|
-
|----------|----------|--------------|
|
|
791
|
-
| [Concepts primer](docs/CONCEPTS.md) | everyone, start here | what Kshana does and why, from zero to the physics |
|
|
792
|
-
| [Playground](web/README.md) | everyone | run the engine in your browser (WebAssembly); build & deploy notes |
|
|
793
|
-
| [Glossary](docs/GLOSSARY.md) | everyone | plain-language definitions of every term |
|
|
794
|
-
| [Architecture](docs/ARCHITECTURE.md) | developers / reviewers | module map, engine pipeline, dispatch, and diagrams |
|
|
795
|
-
| [Validation status](docs/VALIDATION.md) | reviewers / citers | what is `validated` vs `not modeled`, with evidence |
|
|
796
|
-
| [Provenance](docs/PROVENANCE.md) | reviewers / citers | every sensor parameter, model, and dataset traced to its published source, in one citable table |
|
|
797
|
-
| [Reproducibility & provenance](docs/REPRODUCIBILITY.md) | reviewers / packagers | determinism guarantees, golden-pinning, SBOM, build provenance |
|
|
798
|
-
| [Wheel platform tags](docs/WHEEL_TAGS.md) | packagers | the abi3 Python wheel matrix — which platform tag `pip install kshana` resolves |
|
|
799
|
-
| [Positioning](docs/POSITIONING.md) | evaluators | where Kshana sits vs RTKLIB/gLAB (complementary), and the zero-install browser tier |
|
|
800
|
-
| [Technical report](paper/kshana-technical-report.md) · [JOSS paper](paper/paper.md) | reviewers / citers / evaluators | the full extended research paper — architecture, per-domain models, validation, case studies, and limitations — plus the concise JOSS submission |
|
|
801
|
-
| [SGP4 validation](docs/SGP4-VALIDATION.md) | reviewers / citers | agreement with the AIAA 2006-6753 reference (666 states, ~4 mm) **and** a head-to-head against the independent `sgp4` crate (agree to sub-micron / 4.12 mm) |
|
|
802
|
-
| [Force-model validation](docs/AGENCY-ORBIT-VALIDATION.md) | reviewers / citers | the full-force engine (`src/precise_od.rs`) fit to agency ephemerides — methodology and validated residuals |
|
|
803
|
-
| [Real TLE guide](docs/REAL_TLE_GUIDE.md) | users | driving scenarios from real Celestrak / Space-Track constellation TLEs (vs the bundled synthetic Walker set) |
|
|
804
|
-
| [Integrity FoM](docs/INTEGRITY.md) | evaluators | what the `integrity` / `security` figures mean — and what they are **not** vs aviation HPL/VPL |
|
|
805
|
-
| [ARAIM reference](docs/ARAIM_REFERENCE.md) | reviewers / integrators | the open MHSS ARAIM protection-level implementation — the `b_k` nominal-bias projection, σ_URA vs σ_URE, and the fault-mode priors |
|
|
806
|
-
| [Quantum models](docs/QUANTUM.md) · [details](docs/QUANTUM-MODELS.md) | reviewers | the cold-atom-interferometer physics layer, and where coefficients are still looked up |
|
|
807
|
-
| [Compliance](docs/COMPLIANCE.md) | evaluators | DO-229E / DO-316 algorithm scope, and what is **not** a conformance claim |
|
|
808
|
-
| [Standards & interoperability](docs/STANDARDS.md) | integrators | the GNSS / flight-dynamics / agency interchange formats Kshana reads and writes (RINEX, SP3, CCSDS OEM/OMM/TDM/Space-Packet, …) |
|
|
809
|
-
| [Result schema](docs/SCHEMA.md) | integrators | every field of the result JSON, with units and a source pointer |
|
|
810
|
-
| [Python API](docs/PYTHON_API.md) | Python users | the PyO3 binding surface — calling the engine, the scenario/result types, and examples |
|
|
811
|
-
| [Claims vs reality](docs/CLAIMS-VS-REALITY.md) | reviewers | the overclaim-closure ledger + the CI guard (`tests/no_overclaims.rs`) that keeps it resolved |
|
|
812
|
-
| [Roadmap](ROADMAP.md) | everyone | the phased roadmap — what has shipped and what is next |
|
|
813
|
-
| [MCP server](mcp/kshana-mcp/README.md) · [JetBrains plugin](ide/jetbrains/README.md) | agents / IDE users | run Kshana from an AI assistant or a JetBrains IDE |
|
|
814
|
-
| [Changelog](CHANGELOG.md) | everyone | released history (Keep a Changelog + SemVer) |
|
|
815
|
-
| [Contributing](CONTRIBUTING.md) | contributors | build, guards, test/citation discipline, DCO |
|
|
816
|
-
| [Governance](GOVERNANCE.md) | contributors / community | how Kshana is governed — who decides, how, and the open/closed boundary |
|
|
817
|
-
| [Code of Conduct](CODE_OF_CONDUCT.md) | community | expected conduct (Contributor Covenant) |
|
|
818
|
-
| [Security policy](SECURITY.md) | reporters | how to report a vulnerability; dual-use note |
|
|
819
|
-
|
|
820
|
-
## Validation, reproducibility & honesty
|
|
821
|
-
|
|
822
|
-
- Every noise term is calibrated to a **published, cited** figure and validated
|
|
823
|
-
against the standard relation (Allan deviation for clocks; Groves' dead-reckoning
|
|
824
|
-
error growth for inertial; the timing→ranging conversion for time transfer). Status
|
|
825
|
-
per term is tracked in [`docs/VALIDATION.md`](docs/VALIDATION.md) as `validated` or
|
|
826
|
-
`not modeled` — nothing is presented as validated that is not.
|
|
827
|
-
- **Reproducible by construction:** `scenario + seed + engine version → identical
|
|
828
|
-
bits`. `scripts/check-reproducible.sh` enforces it; quantum and classical runs use
|
|
829
|
-
independent seeds so their noise is uncorrelated.
|
|
830
|
-
- Maturity is stated honestly: optical-clock and optical-link figures are *targets /
|
|
831
|
-
ground-demonstrator* results, not flown.
|
|
832
|
-
|
|
833
|
-
### Validation at a glance
|
|
834
|
-
|
|
835
|
-
Every row is enforced by a named test in CI; the full evidence (and what is *not*
|
|
836
|
-
modelled) is in [`docs/VALIDATION.md`](docs/VALIDATION.md) and the per-release
|
|
837
|
-
[`kshana-validation-summary.html`](https://github.com/AshfordeOU/kshana/releases)
|
|
838
|
-
artifact (generated by `cargo run --bin validation_report`, SLSA-attested).
|
|
839
|
-
|
|
840
|
-
The **Status** column states the *kind* of evidence, matching the validation ladder above: **VALIDATED** = checked against an independent external oracle (real data, an independent library, or published reference vectors); **MODELLED** = checked against analytic truth or simulation self-consistency (no independent external dataset). VALIDATED describes the *method* of checking, not a pass/fail — an honest miss against real data (the LRO row) is still VALIDATED. `CI` rows are process guards, not figures of merit.
|
|
841
|
-
|
|
842
|
-
| Status | Capability | Agreement | Reference / oracle |
|
|
843
|
-
|--------|------------|-----------|--------------------|
|
|
844
|
-
| **VALIDATED** | SGP4/SDP4 propagation | 666/666 vectors, worst **4.12 mm** | AIAA 2006-6753 (Vallado `tcppver.out`) + head-to-head vs the independent `sgp4` crate |
|
|
845
|
-
| **VALIDATED** | Reference frames — IAU 2000A/B nutation, IAU 2006/2000A CIO chain, ERA | **bit-for-bit** (X,Y to 1e-14, s to 1e-18, ERA to 1e-12) | ERFA/SOFA `eraXys06a` · `eraC2ixys` · `eraEra00` · `eraNut00a/b` |
|
|
846
|
-
| **VALIDATED** | GCRS→ITRS vs an independent SPICE engine | max **0.028″** → ≤ 0.86 m ground, ≤ 3.6 m GNSS orbit | ANISE (pure-Rust NAIF/SPICE), same IERS `finals2000A` EOP, 8 epochs 2020–2023 |
|
|
847
|
-
| **MODELLED** | EGM2008 geopotential (degree/order 70) | acceleration = ∇V to **< 1e-6**; zonal collapse to validated J2 | NGA EGM2008 coefficients + analytic ∇V identity |
|
|
848
|
-
| **VALIDATED** | Allan estimators (ADEV/MDEV/TDEV/HDEV) + confidence bands | reproduce reference deviations; χ² bands match | NIST SP 1065 (Riley), 1000-point Table 31/32 |
|
|
849
|
-
| **MODELLED** | IMU error model — ARW / VRW / bias-instability | recovered to **< 5 %** (bias-instability < 15 %) | Analog Devices ADIS16465 datasheet; NaveGo reference profile |
|
|
850
|
-
| **MODELLED** | Numerical (Cowell) propagator, unperturbed | **sub-metre over 24 h**; energy/momentum conserve ~1e-9 | exact universal-variable Kepler |
|
|
851
|
-
| **MODELLED** | Lambert · Tsiolkovsky · porkchop | round-trip to two-body truth; ΔV **< 0.01 %** | Izzo 2015 · rocket equation · analytic Hohmann floor |
|
|
852
|
-
| **MODELLED** | Orbit determination (Gauss–Newton batch) | sub-m / mm·s⁻¹ noiseless; ~2 m at a 5 m noise floor | two-body + J2 over an RK4 arc |
|
|
853
|
-
| **VALIDATED** | Force-model fit vs Galileo precise ephemeris (full-arc) | **0.61 m** 3-D RMS, 24 h, d/o-70, force-only | ESA/ESOC `ESA0MGNFIN` final orbit (E11), real `finals2000A` EOP |
|
|
854
|
-
| **VALIDATED** | Force-model fit vs Swarm-A precise ephemeris (reduced-dynamic) | **0.10 m** 3-D RMS (empirical-tier bound, not a measure) | ESA `SW_OPER_SP3ACOM_2_` precise orbit |
|
|
855
|
-
| **VALIDATED** | Force-model fit vs LRO lunar (honest miss) | **6.6 m** reduced-dynamic, *above* the 5 m target | JPL Horizons LRO (NAIF −85) + GRAIL `GRGM660PRIM` |
|
|
856
|
-
| **MODELLED** | Deep-space Mars OD (reduced-dynamic SRIF) | **≈ 0.2 m** Mars-LMO (simulation FoM, *not* real-mission) | synthetic closed-loop OD — estimator-machinery validation |
|
|
857
|
-
| **VALIDATED** | Sun-central Mars dynamics vs JPL DE440 | **137 m @ 1-day arc** (grows with arc = unmodelled n-body) | JPL DE440 via ANISE (`xval/anise-mars-od`, kernel-gated) |
|
|
858
|
-
| **VALIDATED** | Single-point positioning vs a surveyed IGS coordinate (real observations) | **5.7 m** 3-D RMS / **1.1 m** horizontal, dual-frequency iono-free code SPP | IGS station ABMF survey + GPS broadcast ephemeris, 2018-05-13 (`tests/pvt_abmf.rs`) |
|
|
859
|
-
| **MODELLED** | Tightly-coupled GNSS/INS UKF | **0.77 m RMS** over a 30-min LEO pass incl. a 120 s outage | force-model coast, hand-derived |
|
|
860
|
-
| **MODELLED** | GPS-denied gravity-map navigation | ~70 km INS drift → **~145 m** recovered | ESA NAVISP *Quantum Wayfarer* target |
|
|
861
|
-
| **MODELLED** | Terrain-referenced navigation (TERCOM/SITAN) | 70 km drift → **< 500 m** (grid-resolution floor ~140 m) | SRTM `.hgt` DEM; hand-injected drift (non-circular check) |
|
|
862
|
-
| **MODELLED** | IGRF-14 main field (degree/order 13) | pole ~80.7°N, dipole ~29.7 µT, physical 22–67 µT band | IAGA `igrf14coeffs.txt` (Schmidt semi-normalised) |
|
|
863
|
-
| **MODELLED** | Nav-signal modulation & code tracking | BPSK self-SSC = **2/(3·R_c)**; unit-area PSDs; **sub-metre** C/A DLL jitter @ 45 dB-Hz | Closed-form SSC/PSD anchors + Kaplan & Hegarty DLL thermal-noise formula |
|
|
864
|
-
| **MODELLED** | CR3BP halo/NRHO differential corrector | STM = finite differences; orbit closes to **machine precision**; L2 9:2 NRHO **≈ 6.57 d / perilune ≈ 3,250 km** | finite-difference STM check + published L2 southern 9:2 NRHO (≈ 6.56 d / ≈ 3,370 km) — CR3BP, not a real Gateway ephemeris |
|
|
865
|
-
| **VALIDATED** | ARAIM dual-constellation integrity | constellation-wide fault mode on real GPS + Galileo | EU ARAIM TR / DO-316; Celestrak `gps-ops` 2021-07-28 |
|
|
866
|
-
| **VALIDATED** | GNSS geometry / DOP (GDOP/PDOP/HDOP/VDOP/TDOP) | match to **1e-6 relative** across 8 geometries (well-conditioned → near-singular) | gnss_lib_py 1.0.4 (Stanford NAV Lab) — independent library (`tests/dop_reference.rs`) |
|
|
867
|
-
| **VALIDATED** | ML detector-evaluation metrics (AUC/ROC/confusion/Pd-Pmd/precision/F1) | **exact counts + < 1e-9** over 5 datasets × 24 thresholds | scikit-learn 1.9.0 (Pedregosa et al., JMLR 2011) — independent library (`tests/eval_metrics_reference.rs`) |
|
|
868
|
-
| **VALIDATED** | Quantum-trade numerical kernels (ADEV NNLS fit · χ² consistency bands · van-Loan clock Q) | NNLS + Q **exact**; χ² **< 5e-4** at operating dof ≥ 48 | scipy 1.17.1 — `optimize.nnls` / `stats.chi2.ppf` / `linalg.expm` (`tests/scipy_reference.rs`) |
|
|
869
|
-
| CI | Cross-platform reproducibility | bit-identical input + shape goldens on 3 OSes | Linux / macOS / Windows CI matrix, SHA-256 goldens |
|
|
870
|
-
| CI | Test coverage | **~96 % line** on `src/`, gated ≥ 85 % | cargo-tarpaulin (LLVM engine) |
|
|
871
|
-
|
|
872
|
-
## FAQ
|
|
873
|
-
|
|
874
|
-
**Do I need to understand quantum physics to use this?**
|
|
875
|
-
No. If you can run a command line you can run Kshana. Start with the
|
|
876
|
-
[plain-language primer](docs/CONCEPTS.md); look terms up in the [glossary](docs/GLOSSARY.md).
|
|
877
|
-
|
|
878
|
-
**Is this a quantum-hardware design or flight software?**
|
|
879
|
-
No. It is a performance *simulator*. Quantum-hardware fidelity comes from published
|
|
880
|
-
error models, not from this tool. See [What it is / is not](#what-it-is--is-not).
|
|
881
|
-
|
|
882
|
-
**Are the quantum results realistic, or marketing?**
|
|
883
|
-
Every parameter is cited to a datasheet or paper, every model is validated against a
|
|
884
|
-
textbook relation, and maturity is labelled honestly in
|
|
885
|
-
[VALIDATION.md](docs/VALIDATION.md) — including that no strontium optical clock has
|
|
886
|
-
flown. The engine is neutral: quantum and classical are the same code with different
|
|
887
|
-
published numbers.
|
|
888
|
-
|
|
889
|
-
**Can I trust two runs to agree?**
|
|
890
|
-
Yes — runs are deterministic: `scenario + seed + engine version → bit-identical output`,
|
|
891
|
-
enforced by `scripts/check-reproducible.sh`.
|
|
892
|
-
|
|
893
|
-
**Can I use it from Python or in a browser?**
|
|
894
|
-
Yes — see [Python](#python) and [WebAssembly](#webassembly). Both call the same engine.
|
|
895
|
-
|
|
896
|
-
**How do I model my own sensor?**
|
|
897
|
-
Write a scenario `.toml` with your sensor's published figures in the `provenance`
|
|
898
|
-
fields. See [Scenario format](#scenario-format) and the examples in `scenarios/`.
|
|
899
|
-
|
|
900
|
-
**Is it free for commercial use?**
|
|
901
|
-
Yes — under the AGPL-3.0, including in commercial settings, as long as you honour the
|
|
902
|
-
AGPL's copyleft (notably: if you modify Kshana and offer it over a network, you must
|
|
903
|
-
offer those users your modified source). If that does not suit you — e.g. you need to
|
|
904
|
-
embed Kshana in a proprietary product or run a closed network service — a commercial
|
|
905
|
-
licence is available from Ashforde OÜ; see [`LICENSING.md`](LICENSING.md) and
|
|
906
|
-
[Support](#support--professional-services).
|
|
907
|
-
|
|
908
|
-
## Troubleshooting
|
|
909
|
-
|
|
910
|
-
**`cargo build` fails on an old toolchain.** Kshana needs Rust ≥ 1.75. Update with
|
|
911
|
-
`rustup update`.
|
|
912
|
-
|
|
913
|
-
**Building the Python extension fails to link on macOS** (`Undefined symbols … _Py…`).
|
|
914
|
-
A Python extension resolves its symbols at load time. `maturin` sets the right linker
|
|
915
|
-
flag automatically — use `maturin develop --features python` rather than a bare
|
|
916
|
-
`cargo build`.
|
|
917
|
-
|
|
918
|
-
**The Python build complains the interpreter is newer than PyO3 knows.** Set
|
|
919
|
-
`PYO3_USE_ABI3_FORWARD_COMPATIBILITY=1` (abi3 wheels are forward-compatible across
|
|
920
|
-
CPython versions).
|
|
921
|
-
|
|
922
|
-
**WebAssembly build can't find the target.** Install it once with
|
|
923
|
-
`rustup target add wasm32-unknown-unknown`, then `wasm-pack build --target web -- --features wasm`.
|
|
924
|
-
|
|
925
|
-
**Where did my output go?** Each run writes `<scenario>.result.json` and
|
|
926
|
-
`<scenario>.chart.svg` next to the input `.toml`. These are git-ignored by design.
|
|
927
|
-
|
|
928
|
-
## Roadmap
|
|
929
|
-
|
|
930
|
-
See [`ROADMAP.md`](ROADMAP.md) for the phased roadmap, [`CHANGELOG.md`](CHANGELOG.md)
|
|
931
|
-
for released history, and [`docs/CAPABILITY.md`](docs/CAPABILITY.md) for the
|
|
932
|
-
per-capability roadmap. The **ITRF-precise frame reduction** is now delivered — the
|
|
933
|
-
full CIO-based IAU 2006/2000A GCRS↔ITRS chain (polar motion + sub-arcsecond nutation),
|
|
934
|
-
validated bit-for-bit against SOFA/ERFA and independently cross-checked against ANISE
|
|
935
|
-
(pure-Rust SPICE) to ≤ 3.6 m at GNSS orbit. Near-term items include tightly-coupled carrier-phase fusion and surfacing the
|
|
936
|
-
loosely-/tightly-coupled GNSS/INS navigator across more packs; the **deep-space / Mars
|
|
937
|
-
radiometric-navigation** engine landed in v0.17.0 (simulation-validated). The
|
|
938
|
-
**quantum physics layer** is a **P2** item: the CAI accelerometer is now simulated from
|
|
939
|
-
first principles (Mach–Zehnder phase, projection noise, contrast decay, vibration
|
|
940
|
-
coupling), while the clock/time-transfer sensors are still driven by published
|
|
941
|
-
Allan/noise-budget coefficients. GMST-based TEME↔ECEF, the IERS
|
|
942
|
-
leap-second time systems (UTC/TAI/TT/UT1), SGP4/SDP4 orbit propagation (v0.7.0,
|
|
943
|
-
validated against the AIAA 2006-6753 vectors), and the runnable `gnss-ins` fusion
|
|
944
|
-
pack have all **shipped**, and the inertial velocity is exposed downstream. An active
|
|
945
|
-
stochastic time-spoof detector (Neyman–Pearson / χ²₁ energy test with Monte-Carlo
|
|
946
|
-
P_fa/P_md and a Security FoM of 1−P_md), a link-budget jamming model (J/S → effective
|
|
947
|
-
C/N₀ → loss of lock), multi-constellation availability, a single-axis (1-DOF)
|
|
948
|
-
IMU error budget, two independent (clock + position) Kalman estimators reported as a
|
|
949
|
-
combined FoM, real constellation geometry from TLEs, an HTML scorecard report,
|
|
950
|
-
geometry-derived GNSS availability
|
|
951
|
-
*and* dilution of precision from Keplerian orbits with eccentricity and J2 drift,
|
|
952
|
-
Monte Carlo confidence bands, trade-study parameter sweeps, an in-browser WebAssembly
|
|
953
|
-
playground, and optional Python (PyO3) and WebAssembly (wasm-bindgen) bindings have
|
|
954
|
-
landed on `main`.
|
|
955
|
-
|
|
956
|
-
## Contributing
|
|
957
|
-
|
|
958
|
-
See [`CONTRIBUTING.md`](CONTRIBUTING.md). In short: tests pass (`cargo test`), the
|
|
959
|
-
two guard scripts pass, Conventional Commits, and a `CHANGELOG.md` `[Unreleased]`
|
|
960
|
-
entry for every user-visible change. Participation is governed by our
|
|
961
|
-
[Code of Conduct](CODE_OF_CONDUCT.md). To report a security issue, see the
|
|
962
|
-
[Security policy](SECURITY.md) — please do not open a public issue for vulnerabilities.
|
|
963
|
-
|
|
964
|
-
## Citing
|
|
965
|
-
|
|
966
|
-
If you use Kshana in academic or technical work, please cite it. Machine-readable
|
|
967
|
-
metadata is in [`CITATION.cff`](CITATION.cff) (GitHub renders a "Cite this repository"
|
|
968
|
-
button from it); cite the version you used (e.g. `v0.20.0`) together with the
|
|
969
|
-
scenario and seed for full reproducibility. Every release is archived on Zenodo with
|
|
970
|
-
a citable DOI — the concept DOI [10.5281/zenodo.20528627](https://doi.org/10.5281/zenodo.20528627)
|
|
971
|
-
always resolves to the latest version.
|
|
972
|
-
|
|
973
|
-
> Baweja, C. (2026). *Kshana — a PNT-resilience simulator with quantum-sensor performance models*. Ashforde OÜ. https://doi.org/10.5281/zenodo.20528627
|
|
974
|
-
|
|
975
|
-
## Versioning & releases
|
|
976
|
-
|
|
977
|
-
Kshana follows [Semantic Versioning](https://semver.org). While pre-1.0 the public
|
|
978
|
-
scenario/result schema may still change; breaking changes are called out explicitly in
|
|
979
|
-
the [`CHANGELOG.md`](CHANGELOG.md). Every result is reproducible from
|
|
980
|
-
`scenario + seed + engine version`.
|
|
981
|
-
|
|
982
|
-
**Every `vX.Y.Z` tag publishes all channels automatically** — one CI pipeline fans out to:
|
|
983
|
-
|
|
984
|
-
| Channel | Install / get | Contents |
|
|
985
|
-
|---------|---------------|----------|
|
|
986
|
-
| [crates.io](https://crates.io/crates/kshana) | `cargo install kshana` · `kshana = "0.18"` | Rust library + CLI |
|
|
987
|
-
| [crates.io](https://crates.io/crates/kshana-mcp) | `cargo install kshana-mcp` | the MCP server |
|
|
988
|
-
| [PyPI](https://pypi.org/project/kshana/) | `pip install kshana` | abi3 wheels (Linux/macOS/Windows) + sdist |
|
|
989
|
-
| [npm](https://www.npmjs.com/package/kshana) | `npm install kshana` | WebAssembly module + JS wrapper |
|
|
990
|
-
| [ghcr.io](https://github.com/AshfordeOU/kshana/pkgs/container/kshana-mcp) | `docker run -i ghcr.io/ashfordeou/kshana-mcp` | multi-arch OCI image — no toolchain needed |
|
|
991
|
-
| official MCP registry | auto-discovered by MCP clients | `io.github.ashfordeOU/kshana-mcp` |
|
|
992
|
-
| [JetBrains Marketplace](https://plugins.jetbrains.com/plugin/32181-kshana--pnt-simulator) | IDE → Plugins → search "Kshana" | the **Kshana — PNT simulator** IDE plugin |
|
|
993
|
-
| [GitHub Releases](https://github.com/AshfordeOU/kshana/releases) | download | `kshana` + `kshana-mcp` binaries, a CycloneDX **SBOM**, **SLSA** build provenance, and an HTML validation summary |
|
|
994
|
-
| [Zenodo](https://doi.org/10.5281/zenodo.20528627) | DOI | a citable archive of every release |
|
|
995
|
-
| [kshana.dev](https://kshana.dev) | open in a browser | the WebAssembly playground (redeployed from `main`) |
|
|
996
|
-
|
|
997
|
-
The MCP server's crate / image / registry version tracks the engine (it bundles the
|
|
998
|
-
library); the JetBrains plugin versions independently (it shells out to your installed
|
|
999
|
-
`kshana` binary).
|
|
1000
|
-
|
|
1001
|
-
## License
|
|
1002
|
-
|
|
1003
|
-
**Dual-licensed.** Use Kshana under **either** the GNU **AGPL-3.0-only** (see
|
|
1004
|
-
[`LICENSE`](LICENSE)) **or** a **commercial licence** from Ashforde OÜ for
|
|
1005
|
-
proprietary/closed integration that the AGPL does not suit. Which one applies, and
|
|
1006
|
-
why it is set up this way, is explained in [`LICENSING.md`](LICENSING.md).
|
|
1007
|
-
|
|
1008
|
-
Contributions are licensed inbound under the AGPL **and** grant Ashforde OÜ the right
|
|
1009
|
-
to include them in the commercially-licensed edition (so the dual-licence keeps
|
|
1010
|
-
working) — see [`CONTRIBUTING.md`](CONTRIBUTING.md). Sign off each commit per the
|
|
1011
|
-
Developer Certificate of Origin with `git commit -s`.
|
|
1012
|
-
|
|
1013
|
-
**Trademark.** "Kshana" and its marks are trademarks of Ashforde OÜ. The licence
|
|
1014
|
-
covers the code, not the name — please rename forks and derivative distributions.
|
|
1015
|
-
|
|
1016
|
-
## Support & professional services
|
|
1017
|
-
|
|
1018
|
-
Kshana is free and open source under the AGPL-3.0 and **professionally developed and
|
|
1019
|
-
maintained by Ashforde OÜ** (Estonia). The open engine is complete and usable on its
|
|
1020
|
-
own. For organisations that need more, Ashforde OÜ offers:
|
|
1021
|
-
|
|
1022
|
-
- **Commercial support & integration** — embedding Kshana in your toolchain, custom
|
|
1023
|
-
scenarios, and priority fixes.
|
|
1024
|
-
- **Custom sensor models** — calibrated to your hardware, including export-sensitive
|
|
1025
|
-
resilience models maintained in a private overlay.
|
|
1026
|
-
- **Kshana Pro** — proprietary model-based systems-engineering and programme tooling
|
|
1027
|
-
that plugs into the open engine to complete the workflow.
|
|
1028
|
-
- **Training & consulting** on quantum/classical PNT performance analysis.
|
|
1029
|
-
|
|
1030
|
-
This is the open-core model: the engine is, and stays, openly licensed; the sustaining
|
|
1031
|
-
business is expertise, support, and the proprietary extensions — not license fees.
|
|
1032
|
-
Contact **contact@ashforde.org**.
|
|
1033
|
-
|
|
1034
|
-
## Key references
|
|
1035
80
|
|
|
1036
|
-
|
|
81
|
+
Beyond `run` / `run_full` / `version`, the module also exports `summary` (the one-line
|
|
82
|
+
result string), `list_kinds` / `error_kind` (introspection), and
|
|
83
|
+
`encode_permalink` / `decode_permalink` — the shareable-URL codec the
|
|
84
|
+
[playground](https://ashforde.org) uses to round-trip a whole scenario through the
|
|
85
|
+
address-bar fragment.
|
|
1037
86
|
|
|
1038
|
-
|
|
1039
|
-
|
|
1040
|
-
|
|
1041
|
-
|
|
1042
|
-
- Pedregosa et al. — *scikit-learn: Machine Learning in Python*, [JMLR 12 (2011)](https://jmlr.org/papers/v12/pedregosa11a.html): the reference ROC/AUC, confusion-matrix and precision/recall/F1 implementations the RF-impairment evaluation testbed is matched to exactly (`tests/eval_metrics_reference.rs`).
|
|
1043
|
-
- Virtanen et al. — *SciPy 1.0*, [Nature Methods 17 (2020)](https://doi.org/10.1038/s41592-019-0686-2): `optimize.nnls`, `stats.chi2` and `linalg.expm` — the reference routines the quantum-trade measured-ADEV NNLS fit, the χ² consistency bands, and the van-Loan clock process-noise covariance are validated against (`tests/scipy_reference.rs`).
|
|
1044
|
-
- Knowles, Kanhere, Neamati & Gao — *gnss_lib_py*, [SoftwareX 27 (2024)](https://doi.org/10.1016/j.softx.2024.101811): used both as open prior art (see *Comparison & open prior art* below) and as the **independent DOP oracle** the GDOP/PDOP/HDOP/VDOP/TDOP computation is matched to 1e-6 (`tests/dop_reference.rs`).
|
|
1045
|
-
- Montenbruck & Gill — *Satellite Orbits: Models, Methods and Applications* ([Springer](https://doi.org/10.1007/978-3-642-58351-3)): the force models behind the force-model fit to agency precise ephemerides.
|
|
1046
|
-
- Howell — *Three-dimensional, periodic, halo orbits*, Celestial Mechanics 32(1) (1984), [doi:10.1007/BF01358403](https://doi.org/10.1007/BF01358403); Zimovan-Spreen, Howell & Davis — *Near rectilinear halo orbits and nearby higher-period dynamical structures*, Astrodynamics 6 (2022), [doi:10.1007/s42064-021-0125-x](https://doi.org/10.1007/s42064-021-0125-x) (the halo/NRHO families the CR3BP differential corrector reproduces).
|
|
87
|
+
Every figure of merit is labelled **validated** or **modelled**; optical-clock figures
|
|
88
|
+
are space goals on ground hardware (no strontium optical clock has flown). Maturity is
|
|
89
|
+
*not* uniform across domains — Earth PNT is real-data validated; deep-space / Mars
|
|
90
|
+
navigation is simulation-validated; real-mission deep-space OD is on the roadmap.
|
|
1047
91
|
|
|
1048
|
-
|
|
92
|
+
## Learn more
|
|
1049
93
|
|
|
1050
|
-
-
|
|
1051
|
-
-
|
|
1052
|
-
-
|
|
1053
|
-
-
|
|
1054
|
-
- Giorgetta et al., *Nature Photonics* 7, 434 (2013) — [arXiv:1211.4902](https://arxiv.org/abs/1211.4902); Deschênes et al., *Phys. Rev. X* 6, 021016 (2016) — [APS](https://journals.aps.org/prx/abstract/10.1103/PhysRevX.6.021016) (optical two-way time-frequency transfer; the optical inter-satellite link models its non-reciprocity budget after these).
|
|
1055
|
-
- Betz — *Binary Offset Carrier Modulations for Radionavigation*, NAVIGATION 48(4) (2001), [doi:10.1002/j.2161-4296.2001.tb00247.x](https://doi.org/10.1002/j.2161-4296.2001.tb00247.x) (the BOC modulation and spectral-separation theory behind `src/navsignal.rs`).
|
|
1056
|
-
- Kaplan & Hegarty (eds.) — *Understanding GPS/GNSS: Principles and Applications* (3rd ed., Artech House, 2017): the anti-jam effective-C/N₀ equation and the early–late DLL code-tracking thermal-noise jitter the nav-signal and jamming models use.
|
|
94
|
+
- **Full README & validation matrix** → <https://github.com/AshfordeOU/kshana>
|
|
95
|
+
- **Live playground** → <https://ashforde.org>
|
|
96
|
+
- **Capabilities** → [docs/CAPABILITY.md](https://github.com/AshfordeOU/kshana/blob/main/docs/CAPABILITY.md)
|
|
97
|
+
- **Validation & provenance** → [docs/VALIDATION.md](https://github.com/AshfordeOU/kshana/blob/main/docs/VALIDATION.md) · [docs/PROVENANCE.md](https://github.com/AshfordeOU/kshana/blob/main/docs/PROVENANCE.md)
|
|
1057
98
|
|
|
1058
|
-
|
|
99
|
+
## Licence
|
|
1059
100
|
|
|
1060
|
-
|
|
1061
|
-
|
|
1062
|
-
|
|
1063
|
-
|
|
1064
|
-
|
|
1065
|
-
- Bertone et al. — *Earth and Space Science* 8(6) (2021), [doi:10.1029/2020EA001454](https://doi.org/10.1029/2020EA001454): GRAIL reduced-dynamic OD, the empirical-acceleration floor the LRO fit reproduces.
|
|
101
|
+
Free and open source under the **GNU AGPL-3.0-only**. A **commercial licence** is
|
|
102
|
+
available from [Ashforde OÜ](https://ashforde.org) for proprietary/closed integration
|
|
103
|
+
— see [LICENSING.md](https://github.com/AshfordeOU/kshana/blob/main/LICENSING.md).
|
|
104
|
+
Professionally developed and maintained by Ashforde OÜ; commercial support, integration,
|
|
105
|
+
and proprietary extensions available.
|