Anti-aliased Sloan refinement
Implementation & benchmark of arXiv:2608.20098 — Superconvergence and aliasing saturation in Sloan iteration for spherical integral equations
Source paper: Superconvergence and aliasing saturation in Sloan iteration for spherical integral equations arXiv:2608.20098 ⓘ · analyzed Sep 1, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Convert a neural operator block into a shared-weight iterative fixed-point refinement scheme that exploits repeated smoothing while avoiding repeated low-resolution projections. Compute all refinement steps at an overresolved latent bandwidth and apply the target-bandwidth projection only at the end, reducing the opportunity for unresolved frequencies to alias into retained channels.
Formulas
Mathematical statement
The paper's central Sloan identity is \(u-\mathcal{S}_{q}(v)=\mathcal{T}^{q}(u-v)\), with \(\mathcal{S}_{0}(v)=v\) and \(\mathcal{S}_{q+1}(v)=f+\mathcal{T}\mathcal{S}_{q}(v)\). Here \(u\) solves \((I-\mathcal{T})u=f\), \(f\) is the forcing term, \(\mathcal{T}\) is a compact smoothing operator, \(v\) is the initial approximation, and \(q\) is the number of iterations. Thus the initial error is filtered by \(\mathcal{T}^{q}\). Let \(P_B\) be projection onto the retained spherical-harmonic or graph-frequency bandwidth \(B\), and let \(K_\theta\) be a learned analogue of \(\mathcal{T}\). Repeated low-resolution refinement applies \((P_BK_\theta P_B)^q\), whereas overresolved refinement approximates \(P_BK_\theta^q\), with projection delayed until the final step. The difference matters because the intermediate high-frequency component \((I-P_B)K_\theta h\) can be folded back into low frequencies by quadrature or sampling. Increasing the overresolved bandwidth \(B'>B\) reduces this aliasing term, so the depth at which accuracy saturates should increase with \(B'\).
Implementation notes
(1) Integration point: replace a spectral convolution, graph message-passing layer, or spherical neural-operator block that maps a feature field z to an output y. Let K_theta be the learned smoothing operator, f_theta(z) be a residual or forcing branch, B be the normal representation bandwidth, and B'>B be an overresolved bandwidth. (2) Pseudocode: h = upsample_or_zero_pad(z, B_prime); for i in 1..q: h = f_theta(upsample_or_zero_pad(z, B_prime), h) + K_theta(h); h = normalization(h); y = project_to_bandwidth(h, B); return y. Share K_theta across iterations to match Sloan refinement; compare against the baseline that projects h to B after every iteration. The baseline therefore implements h = P_B(f + K_theta(P_B(h))), while the proposed block keeps h at B' until the final projection. (3) Compute from the paper's mathematics the two competing mechanisms: repeated smoothing corresponds to K_theta^q, while aliasing is diagnosed by the discrepancy D_q = norm(P_B K_theta^q h0 - (P_B K_theta P_B)^q h0)_2. Estimate D_q with one high-resolution reference run, using zero-padding or a larger spherical quadrature grid to approximate B'. Estimate the effective smoothing spectrum empirically by applying K_theta to random Fourier or spherical-harmonic probes and measuring per-frequency gain. (4) First experiment: train a small spherical spectral network on synthetic band-limited Poisson or Helmholtz data, using B=16, B'=16,32,64 and q=1,2,4,8. Match parameter count and compare test relative L2 error, training loss versus FLOPs, activation memory, and D_q. The expected signal is that q>1 improves error at B'=2B or 4B while the B'=B baseline saturates; increasing B' should shift the saturation point to larger q. A useful negative result is no benefit when K_theta is not smoothing, which tests whether the mechanism is genuinely tied to spectral decay rather than extra depth.
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.