Saturation-certified checkpoint bracket
Implementation & benchmark of arXiv:2608.16473 — Reference-free logged energy-oracle recovery for neural approximations of symmetric coercive variational problems: conforming Riesz reconstruction and archive-level selection
Source paper: Reference-free logged energy-oracle recovery for neural approximations of symmetric coercive variational problems: conforming Riesz reconstruction and archive-level selection arXiv:2608.16473 ⓘ · analyzed Sep 1, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Use two nested Riesz reconstruction spaces to estimate unresolved residual energy for every checkpoint. Under a measurable saturation assumption, convert the coarse and enriched monitors into lower and upper error bounds, and certify a unique checkpoint whenever its upper bound lies below every competitor's lower bound.
Formulas
Mathematical statement
For candidate \(u_k\), let \(M_m(k)=\|z_{k,m}\|_a\) be the conforming Riesz monitor in \(W_m\), and let \(E(k)=\|u_k-u^*\|_a\) be the unknown energy error. The lower-bound property is \(M_m(k)\leq E(k)\). Assume hierarchical saturation: \(\|z_k-z_{k,m+1}\|_a\leq q\|z_k-z_{k,m}\|_a\) for \(0<q<1\). Let \(\Delta_m(k)=M_{m+1}(k)-M_m(k)\) estimate the unresolved energy. Use \(L_k=M_{m+1}(k)\) as a lower estimate and \(U_k=M_{m+1}(k)+q\Delta_m(k)/(1-q)\) as a practical upper estimate. If \(U_{k_*}<L_j\) for every competitor \(j\neq k_*\), then \(k_*\) is certified as the unique minimum-energy checkpoint, conditional on saturation.
Implementation notes
Add a certification stage after archive scoring for a PINN or neural finite-element solver. Construct a coarse conforming auxiliary space \(W_m\) and a nested refined space \(W_{m+1}\), then assemble and factorize both stiffness matrices. For every checkpoint, evaluate residual vectors \(b_m^{(k)}\) and \(b_{m+1}^{(k)}\), solve the two Riesz systems, and compute \(M_m(k)\) and \(M_{m+1}(k)\). Set \(L_k=M_{m+1}(k)\), \(\Delta_k=M_{m+1}(k)-M_m(k)\), and \(U_k=M_{m+1}(k)+q\Delta_k/(1-q)\). Select the candidate with smallest upper estimate, but issue a certificate only when its upper estimate is below every other candidate's lower estimate. Estimate \(q\) on a calibration subset with a third refinement using \(q_k=\Delta_{m+1}(k)/\Delta_m(k)\); use a high quantile plus safety margin and clip it below one. Log saturation violations rather than silently assuming the condition holds. First test on a 2D Poisson or elasticity PINN with 100 checkpoints and a manufactured exact solution. Compare single-level selection, bracket selection, and the exact oracle. Measure energy-error regret, certificate frequency, saturation violations, and linear-solve overhead. The desired signal is fewer incorrect rankings at modest auxiliary resolution and increasing certificate frequency under refinement.
Verification
This idea has not been verified yet.
Verification happens in two stages: Stage 1 — a mechanism check on a toy system confirms the claimed mathematical phenomenon reproduces; Stage 2 — a benchmark implements the idea on a real (small) neural network task and compares it against a tuned baseline over 8 paired seeds with a permutation test.
Artifacts
Artifacts unavailable.