Convex-overlap boundary spectral loss
Implementation & benchmark of arXiv:2608.19806 — Boundary-Weighted Fourier Inequalities for Convex Domains
Source paper: Boundary-Weighted Fourier Inequalities for Convex Domains arXiv:2608.19806 ⓘ · analyzed Sep 1, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Add a Fourier-domain residual loss whose per-frequency weight is determined by the geometric overlap of a convex bandwidth domain with its reflection about that frequency. Frequencies close to the boundary receive larger weight through \(\omega_\Omega^{-d}\), forcing the network to model fragile spectral components instead of optimizing only the high-energy interior. Use clipping or a bounded transform of the singular weight so that a few boundary bins cannot dominate training.
Formulas
Mathematical statement
The paper studies the inequality \(\int_\Omega |\hat f(x)|^p/\omega_\Omega(x)^d\,dx\le C\|f\|_{L^q}^p\), where \(\Omega\subset\mathbb R^n\) is an open convex set containing no affine line, \(f\in L^q(\mathbb R^n)\) has Fourier support in \(\overline\Omega\), \(\hat f\) is its Fourier transform, \(1\le p,q<\infty\), \(d\in\mathbb R\), \(C\) is an inequality constant, and \(m\) denotes Lebesgue measure. The geometric multiplier is \(\omega_\Omega(x)=m(\Omega\cap(2x-\Omega))\) for \(x\in\Omega\), where \(2x-\Omega=\{2x-y:y\in\Omega\}\). The overlap is large for interior frequencies and shrinks toward the boundary, so \(\omega_\Omega(x)^{-d}\) emphasizes boundary frequencies when \(d>0\). For the discrete rectangular domain \(\Omega=\prod_{j=1}^n(-1,1)\), the exact continuous overlap is \(\omega_\Omega(x)=2^n\prod_{j=1}^n(1-|x_j|)\). The neural adaptation replaces \(\hat f\) by the prediction residual spectrum \(\widehat r\), estimates the integral by a normalized FFT sum, and uses a clipped weight \(w_k=[\max(\omega_\Omega(\xi_k),\varepsilon)]^{-d}\), with frequency \(\xi_k\), cutoff \(\varepsilon>0\), and exponent \(d>0\).
Implementation notes
Integrate this at the output-loss stage of a convolutional image-to-image model, such as a small U-Net or FNO-like CNN. Let the prediction be \(\hat y\), target be \(y\), and residual be \(r=y-\hat y\); compute a channelwise 2D FFT of \(r\). Normalize FFT frequencies to \(\xi_k\in[-1,1]^2\), construct the rectangular domain \(\Omega=(-1,1)^2\), and evaluate \(\omega_\Omega(\xi_k)=4(1-|\xi_{k,1}|)(1-|\xi_{k,2}|)\). Set \(w_k=\max(\omega_\Omega(\xi_k),\varepsilon)^{-d}\), initially using \(d=0.25\) or \(0.5\) and \(\varepsilon\) equal to the fifth percentile of nonzero overlap values. The training loss is \(\mathcal L=\mathcal L_{\mathrm{pixel}}+\lambda\mathcal L_{\mathrm{BO}}\), where \(\mathcal L_{\mathrm{pixel}}\) is ordinary L1 or MSE and \(\mathcal L_{\mathrm{BO}}\) is the normalized weighted FFT residual in the displayed formula. Stop gradients through \(w_k\). To prevent corner bins from dominating, cap \(w_k\) at the 99th percentile and warm up \(\lambda\) linearly over the first 5 epochs. The paper supplies the overlap geometry and inverse-power form; no theorem constant \(C\) is estimated. The first cheap experiment is CIFAR-10 or DIV2K 2x super-resolution with a 4-layer U-Net, comparing pixel MSE, unweighted FFT MSE, radial \(|\xi|^2\) weighting, and this boundary-overlap loss at matched training FLOPs. Track PSNR/SSIM, validation loss, spectral residual separately in interior and boundary bands, and gradient-norm variance. Success means lower boundary-band error and improved validation PSNR without exploding gradients or increasing pixel-space error.
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.