Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Multidimensional quadrature evidence

This page connects each public Phase B claim to the artifact that owns it. Passing one row does not promote the other rows, and none establishes universal superiority over another library.

Evidence map

QuestionEvidence ownerWhat it warrantsWhat it does not warrant
Do values and work records match independent truths?truth artifactDeclared analytic and numerical cases at recorded tolerancesAccuracy for unresolved discontinuities or arbitrary integrands
Does first-order replay match the accepted formula?replay artifactDeclared parameter, bound, measure, and quantity derivativesDerivatives of discrete refinement or higher derivatives
Do randomized intervals meet the declared calibration experiment?RQMC artifactReal-scalar fixed-look and bounded sequential coverage in the frozen campaignPer-integral certainty or complex/vector confidence intervals
How were external capabilities compared?comparison artifactOnly the recorded exact, strong-match, node-match, family-match, or capability labelAlgorithmic identity from a broad family label
Did a predeclared optimization trigger fire?performance baselineA dimension-16 analytic memory proxy exceeded the family-median thresholdAn observed process/device-memory bottleneck

Estimator meanings

FamilyQuadError.kindInterpretation
Fixed tensorUNAVAILABLENo truncation-error estimate is supplied
Adaptive tensorSUCCESSIVE_LEVELDifference between accepted global nested levels
Adaptive cubatureEMBEDDED_RULEShared-node high/low Genz-Malik disagreement
Sparse gridSPARSE_GRID_SURPLUSHierarchical frontier or accepted-index evidence
Deterministic SobolUNAVAILABLENo randomization, therefore no confidence interval
Randomized SobolCONFIDENCE_INTERVALFixed-look Student-t or bounded sequential empirical-Bernstein evidence

Replay and quantity boundary

Replay differentiates

I^(θ)=i=1Nacceptedwi(θ)f(xi(θ),θ),\widehat{I}(\theta) = \sum_{i=1}^{N_{\mathrm{accepted}}} w_i(\theta)f(\boldsymbol{x}_i(\theta),\theta),

with the accepted structure frozen. Opt-in quantity normalization may use heterogeneous coordinate units and restores the result unit after the raw JAX kernel. It remains alpha: downstream package adoption and direct Quantity-PyTree quotient-unit Jacobians are not implied.

Performance disposition

The immutable baseline fired only the analytic memory-proxy trigger at dimension 16:

36,864 bytes16,384 bytes=2.25.\frac{36{,}864\ \mathrm{bytes}}{16{,}384\ \mathrm{bytes}} =2.25.

This is expected linear coordinate storage and a small absolute proxy, not an observed peak-memory failure. The optimization addendum therefore authorizes no runtime change until a measured campaign demonstrates a material case.

Reproduce the focused evidence

The following checks are intentionally narrower than the full release gate:

uv run --no-sync python scripts/generate_quad_multidim_evidence.py --check
uv run --no-sync python scripts/generate_quad_rqmc_evidence.py --check
uv run --no-sync python scripts/benchmark_quad_multidim.py --suite baseline --check

The exhaustive test, documentation, clean-wheel, and observed-memory campaign remain separate release gates.