Boundary-Depth Sparsity
Implementation & benchmark of arXiv:2607.21355 — A free boundary problem driven by boundary distance in the coincidence set
Source paper: A free boundary problem driven by boundary distance in the coincidence set arXiv:2607.21355 ⓘ · analyzed Aug 30, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Add a nonlocal regularizer to nonnegative CNN feature maps that rewards activation-zero pixels lying deep inside a contiguous zero region. Unlike L1 sparsity, this penalizes isolated holes and favors block-like inactive areas that are more amenable to spatial skipping, channel gating, or structured pruning.
Formulas
Mathematical statement
The paper minimizes the functional J(u)=∫_Ω(1/2|∇u|²−fu)dx−∫_Z F(dist(x,∂Z))dx, where u:Ω→R_+ is the nonnegative field, Z={x∈Ω:u(x)=0} is its zero phase, ∂Z is the topological boundary of that phase, dist(x,A)=inf_{a∈A}|x−a| is Euclidean distance to a set, F:R→R_+ is nondecreasing, and f is a forcing term. The transferable mathematical asset is the depth-weighted zero-phase measure ∫_Z F(dist(x,∂Z))dx: a zero pixel receives little reward near the phase boundary and a larger reward when it lies deep inside a connected zero component. For a neural activation map a_c(i)≥0 on grid sites i∈Ω, use a soft zero indicator q_c(i)≈1[a_c(i)=0], a differentiable distance proxy d_c(i), and add R_depth=−λ∑_{c,i}q_c(i)F(d_c(i)) to the task loss. The adaptation is a smooth surrogate of the paper's topological zero set and distance functional, not a claim that the free-boundary existence theory directly applies to stochastic network training.
Implementation notes
Integrate the regularizer after a ReLU, Softplus, or absolute-value activation in a CNN, before the next convolution. Let a[b,c,h,w]≥0 be the feature tensor and let Ω be the H×W grid for one channel. First compute the soft zero membership q=sigmoid((τ−a)/ε), where τ is a small activation threshold, ε controls softness, and δ=10^-6 prevents division by zero. Approximate the zero-phase boundary by b=q(1−q)+δ; this is large for activations near the threshold and small in confidently active or confidently zero pixels. Compute the soft distance d at every pixel with the log-sum-exp soft minimum shown above, using a separable or truncated Gaussian neighborhood for efficiency. Choose F(d)=d^p with p=1 or 2, or F(d)=log(1+d), and add the negative depth reward to the task loss with coefficient λ. Pseudocode is: a=relu(conv(x)); q=sigmoid((τ−a)/ε); b=q*(1−q)+δ; d=soft_distance_to_boundary(b,rho); L=L_task−λ*mean(q*F(d)); backpropagate L and update weights. Anneal ε from 0.2 to 0.05 and λ from zero during the first few epochs to avoid immediate all-zero collapse. The paper contributes the depth-weighted zero-phase objective; the soft mask, boundary proxy, and distance transform are engineering approximations. Estimate sparsity, connected-component size, mean zero-pixel depth, validation accuracy, and hardware runtime. The first experiment should train ResNet-18 on CIFAR-10 with this term versus identical training with L1 activation regularization, group lasso, and no regularizer. At equal accuracy, success means fewer fragmented nonzeros, larger zero blocks, and measurable speedup after compiling block-sparse or spatial-skipping convolution; a useful target is 20–30% fewer active spatial blocks with less than 0.5 percentage-point accuracy loss.
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.