Section C. GeoNum v2.2 Methodology

The substrate precision model underlying the measured benchmarks of Section A and the Phase I plan of Section B. Per-value drift tracking, formal stability bounds, and the cross-domain validation footprint.

C.1 Representation

GeoNum represents a real number as a five-element record:

FieldTypeRole
zoneintegerLucas-tier boundary in log10 space
intensityreal in [0, 1]Continuous log-fractional position within the zone
shadeinteger in [0, 2047]Discrete tessellation; 211 levels match IEEE 754 mantissa width
driftreal, non-negativeAccumulated quantization error in shade units
sign±1Sign bit

The default Theory zone boundary set spans approximately 30 orders of magnitude (Lucas sequence boundaries 1, 2, 3, 4, 7, 11, 18, 29, 47, 76, ...). Domain-specific configurations adjust the boundaries while preserving the arithmetic semantics.

C.2 Arithmetic and drift propagation

Multiplication adds zones, composes intensities in log-space, and propagates drift zone-normalized: dcda · (wa/wc) + db · (wb/wc) + dquant.

Addition uses log-sum-exp on same-sign operands; on different-sign operands it uses log-diff-exp via _logDiffExp with derived cancellation amplification factor 1 / |1 − exp(Δ)|, where Δ is the log-magnitude gap. This makes cancellation observable: as two values of comparable magnitude subtract, drift growth is bounded by the amplification factor, and at exact cancellation the result is flagged as an uncertain zero with the magnitude scale at which the cancellation occurred retained.

The conditional stability theorem bounds drift growth such that EXACT, PRECISE, and VALID trust bands are composable. Any operation between two VALID-or-better operands produces a result whose drift is bounded by the operand drifts and the derived amplification factor.

C.3 Trust bands

TrustDrift bound (shade units)Interpretation
EXACT0Freshly encoded; no accumulated error
PRECISE(0, 1)Sub-shade quantization; less than 1 ULP-equivalent
VALID[1, 204.8)Bounded, composable, safe for solver inputs
DEGRADED[204.8, 1024)Significant accumulated loss
UNRELIABLE[1024, inf)Error spans half a zone width
CANCELLEDuncertain zeroNear-exact cancellation; the result is meaningless at the cancellation scale

The threshold values arise from the substrate s shade count: SHADES / 10 = 204.8 and SHADES / 2 = 1024. The conditional stability theorem applies on the EXACT through VALID range; DEGRADED and UNRELIABLE are outside the theorem s precondition.

C.4 Measured drift on Section A constraint norms

From geonum_drift_on_constraints.json:

ResolutionHam L2 measuredDrift (shades)Trust
1637.82 × 10-10.67PRECISE
2431.52 × 10-10.93PRECISE
3234.96 × 10-20.31PRECISE

All three resolutions return sub-shade drift; the substrate trusts each measured value to better than one ULP-equivalent. The PRECISE tag is the substrate s native precision claim, additive to the standard constraint-violation diagnostic.

C.5 Distinction from f64-only NR codes

HARM, HARM3D, Cactus, the Einstein Toolkit, BAM, AMSS-NCKU, and similar f64-based codes report constraint violation as the primary precision diagnostic. There is no per-value precision tag carried alongside the data; precision loss is inferred after the fact via growing constraint residuals or ad hoc spectral analysis.

GDBS carries four pieces of information that f64-only codes do not provide natively:

ItemWhat it gives
Per-snapshot drift (today, Section A) Every reported constraint value comes with a trust band
Per-RK4-substage drift (Phase I, Section B.5) Every intermediate state of the time integrator carries a tag
Cancellation detection (today) When two large quantities subtract to near zero, the substrate flags CANCELLED and reports the magnitude scale at which the cancellation occurred. The precision floor of the result is explicit, not inferred.
Drift composition theorem (today) Bounds on the drift of any sequence of operations. Relevant for waveform extraction where the precision claim must propagate through many operations from BSSN state to ψ4 to strain to matched filter.

C.6 Cross-domain validation

GeoNum v2.0 supports per-domain zone tessellation: the same drift propagation engine, different boundary configuration tuned to the domain s typical dynamic range. Domains that ship validated benchmarks today:

DomainZone configurationCanonical benchmarkStatus
TheoryLucas (log-spaced)Hawking radiation across 60 orders of magnitude0.27% error, 0.345 shade drift
QuantumTight eV-scaleH2 ground state versus NISTWithin NIST tolerance
FluidsUniform CFD gridTaylor-Green vortex decayRK3-TVD, ensemble L1
PlasmaFrequency-alignedLandau damping rateStability tests
MaterialsLattice-symmetricSilicon lattice constantWithin published tolerance
GeophysicsSpherical harmonicEGM2008 gravitational fieldMulti-extraction radius
BallisticsMach-number scaledHypersonic derivativesCivilian engineering

The Theory and Quantum domains are the load-bearing references for the present validation package; the others demonstrate that the substrate is not bespoke to numerical relativity. The substrate is one engine; the seven configurations differ only in the zone boundary table.

C.7 Auto-tuner: substrate-validated closed-form recognition

The substrate ships cascade_self_tune_mesh_prior (wasm/src/geonum.rs:1004-1053), a search loop over rational combinations of substrate primitives (ζ values, multiple zeta values, π, e, rationals) with built-in false-positive baseline telemetry:

TuneStats {
  candidates_evaluated:        u64,
  matches_in_tolerance:        u64,
  uniform_null_fp_probability: f64,
  ...
}

Used today in three places:

For numerical relativity applications, the auto-tuner becomes the substrate-validated answer to the question: is this measured quantity recognizable as an algebraic constant, or is it transcendental? For phase-drift comparisons in Phase I (Section B.4), the auto-tuner gates whether residual deviations are recognizable closed-form (indicating systematic, not numerical) versus FP-baseline noise.

C.8 Operational implication for the SBIR proposal

The GDBS Phase I claim is not that the platform replaces HARM3D on speed or that it replicates HARM3D BBH simulations in a browser. The claim is precise:

Browser-native BSSN+Z4c evolution with per-value precision tags, validated against industry-standard convergence and stability benchmarks, with a roadmap to LISA-precision waveform extraction grounded in published methods.

The deployed platform at gdbs.getvaultsync.com is the substrate; this validation package is the credibility layer; the Phase I plan delivers the LISA-precision waveform claim.

References

  1. GeoNum v2.2 formal specification: wasm/src/geonum.rs lines 1-1135. Conditional stability theorem at lines 217-250 (Trust enum); arithmetic propagation at 334-428.
  2. cascade_self_tune_mesh_prior source: wasm/src/geonum.rs:1004-1053.
  3. Per-domain zone configurations: core/GdbsWeb.Api/wwwroot/js/geo-num-zones.js.
  4. QED cascade pins (c1-c5): wasm/src/closed_form.rs:1-500.
  5. Theoretical Foundations module validation: core/GdbsWeb.Api/wwwroot/js/hpc-validation.js lines 60-141.