Distribution-Preserving Fragmentation Augmentation
Implementation & benchmark of arXiv:2607.04807 — Thresholds, fragmentation and symmetrization in parabolic equations
Source paper: Thresholds, fragmentation and symmetrization in parabolic equations arXiv:2607.04807 ⓘ · analyzed Aug 30, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Augment spatial training examples by replacing a compact active region with several separated components while preserving its exact value histogram, total active area, and amplitude. The augmentation probes the nonlinear interaction between diffusion-like receptive fields and threshold activations, which the paper shows can make fragmented and compact inputs evolve in opposite directions despite identical distributions.
Formulas
Mathematical statement
The paper studies the reaction-diffusion dynamics ∂_t u = ∂_xx u + f(u) for bounded nonnegative u, with a bistable nonlinearity satisfying f(0) = f(θ) = f(1) = 0, f'(0) < 0, f(s) < 0 for 0 < s < θ, f(s) > 0 for θ < s < 1, integral from 0 to 1 of f(s) ds greater than 0, and f(s) < 0 for s > 1. Here θ is the activation threshold. Theorem 1.6 states that when 0 < R_1 < L_infinity < R_2, there exist an amplitude α > θ and an integer n ≥ 2 such that the compact datum u_0^blk(x) = α times the indicator of (-R_1/2, R_1/2) extinguishes, while the fragmented datum u_0^frag(x) = α times the indicator of E_n(R_1,R_2) propagates to 1; both have the same distribution function and active measure R_1. Define E_n(R_1,R_2) as n equal intervals, each of length R_1/n, uniformly spaced inside an envelope of length R_2, so its measure is R_1. The transferable property is that threshold dynamics depend on spatial connectivity and gap size, not only on pointwise statistics.
Implementation notes
Integrate this at the input or intermediate-feature augmentation stage of a CNN, U-Net, neural cellular automaton, or spatial graph network. For a tensor X with shape C by H by W, select one channel or a learned scalar mask M. Threshold it into an active set, preserve the active values, and replace the active set by n separated components with the same total pixel count. Choose n from {2, 4, 8}, place equal-area components inside an envelope expanded by a factor gamma, and copy the original active amplitude α. Apply the transformation with probability p while keeping the original label. Pseudocode: mask = (X[c] > tau); A = mask.sum(); n = sample({2,4,8}); E = place_equal_components(area=A, envelope=gamma * bounding_box(mask), n=n); X_aug[c] = alpha * E + X[c] * (1-mask). The preserved mathematical quantities are the value distribution, active measure R_1, fragmentation count n, and envelope scale R_2. Estimate tau from feature quantiles and choose gamma and p on a validation split. First test on MNIST or CIFAR-10 with a small CNN and synthetic objects whose connectedness varies. Compare ordinary augmentation, histogram-preserving pixel shuffling, and structured fragmentation at equal compute. Measure clean accuracy, accuracy under object splitting, calibration, and performance as n increases. Success means improved accuracy on fragmented layouts without degrading compact-layout accuracy, together with lower prediction variance across spatial arrangements having the same activation histogram.
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.