Math: Numerical analysis

Machine-learning ideas tagged Numerical analysis in the Math taxonomy of the Math2NN corpus.

✗ Failed on benchmark 2026

Rankine–Hugoniot Front Tokens

Augment a 1D neural operator or transformer with explicit tokens for detected discontinuities. Advance each front analytically using the local Rankine–Hugoniot speed and train the network only to reconstruct smooth regions and the residual caused by source terms and grid resolution.

Useful7/10
Difficulty6/10
Novelty6/10
Paper: Physics-informed token transformer methodology for nonlinear balance laws. I. Schwarzschild--Burgers fluid flows arXiv:2607.23143
△ Mechanism confirmed, baseline not beaten 2026

Residual-Pivoted Kernel Attention

Replace full PSD self-attention with a pivoted Cholesky/Nyström approximation whose landmarks are sampled from the unexplained diagonal mass. Tokens with large residual self-similarity are more likely to become landmarks, so the rank budget is spent on difficult regions rather than uniformly selected tokens.

Useful7/10
Difficulty5/10
Novelty5/10
Paper: A new analysis of the randomly pivoted Cholesky algorithm arXiv:2608.20633
✗ Mechanism failed 2026

Closure-Decorrelation Memory Scheduler

Choose the neural operator's input-history length from the measured correlation time of the unresolved closure signal produced by coarse-graining. This avoids under-memory, which causes systematic closure error, and over-memory, which increases attention cost and can destabilize training. The same diagnostic can drive adaptive memory truncation across physical regimes.

Useful7/10
Difficulty4/10
Novelty7/10
Paper: Flux-form spatiotemporal neural operators for coarse-grained dynamics of multiscale PDEs arXiv:2608.18148
△ Mechanism confirmed, baseline not beaten 2026

Singular-Value-Robust Projector-Splitting LoRA

Train a fixed-rank neural weight update Y=USV^T with a projector-splitting Runge–Kutta step instead of independently applying Adam or gradient descent to U, S, and V. The update evolves the full low-rank matrix using the neural gradient but performs QR-based factor updates, avoiding S^{-1} and remaining stable when adapter singular values collapse or cross zero. Use a common-base midpoint construction so every internal stage starts from the same U,V basis and remains rank r.

Useful7/10
Difficulty5/10
Novelty6/10
Paper: Robust Projector-Splitting Runge-Kutta Integrators of Orders Two and Three arXiv:2608.17157
✗ Mechanism failed 2026

Coarse-to-fine active-support transport attention

Replace dense cross-attention weights with a balanced transport plan whose nonzero query-key edges are maintained by a multiscale active-set procedure. Solve the coarse token-group problem first, lift its support to the fine token grid, add only edges indicated by local cost or marginal residuals, and warm-start the fine problem from the lifted plan. This should provide a principled sparse attention pattern rather than fixing a global top-k pattern before seeing the transport solution.

Useful7/10
Difficulty7/10
Novelty6/10
Paper: A Multiscale Primal-Dual Interior-Point Relaxation Method for Large-Scale Optimal Transport Problems arXiv:2608.12060
△ Mechanism confirmed, baseline not beaten 2026

Dyadic Hankel Boundary Attention

Replace dense attention between tokens on opposite sides of a one-dimensional boundary or segment split with a dyadic low-rank approximation of a Cauchy/Hankel distance kernel. Each distance-scale block uses O(log(1/\varepsilon)) features, and the number of active scales grows only logarithmically with context length after discarding a narrow boundary layer. This is especially suitable for a relative-position attention branch or state-space-like long-range branch, rather than arbitrary…

Useful7/10
Difficulty5/10
Novelty6/10
Paper: An independent proof of the plunge-region conjecture for time-frequency localization operators in dimension one arXiv:2607.23016
✗ Failed on benchmark 2026

Rank-Revealing Representative Tokens

Compress a transformer KV cache by selecting actual past tokens whose key or hidden-state columns form a stable basis for all cached tokens. Instead of retaining tokens with the largest attention scores or leverage scores independently, compute rank-revealing pivoting of the leading right-singular-vector matrix and retain its pivot columns, then evaluate attention using the representatives plus an optional low-cost residual correction.

Useful7/10
Difficulty5/10
Novelty6/10
Paper: Computing Strong Rank-Revealing Factorizations for Matrices with Orthonormal Rows arXiv:2607.13532
△ Mechanism confirmed, baseline not beaten 2026

Residual-energy cross attention

Replace dense query-key attention with an adaptive cross approximation constructed from selected query and key pivot tokens. At each rank, choose the pivot pair that removes large estimated residual energy, update the residual by a rank-1 cross correction, and stop when the residual estimate reaches a target tolerance. The resulting factorization computes approximate attention using a small number of landmark interactions while adapting to the actual token distribution.

Useful7/10
Difficulty6/10
Novelty5/10
Paper: Continuous Cross Approximation of Matrices Arising Out of Kernel Functions arXiv:2607.12540
△ Mechanism confirmed, baseline not beaten 2026

Branch-Free Double-Word FMA Accumulator

Replace ordinary low-precision multiply-add accumulation in selected neural-network reductions with a two-word floating-point accumulator updated by the paper's branch-free DW-FMA network. The high word retains the main sum and the low word stores the rounding residual, improving cancellation behavior without the control-flow divergence of conditional compensated summation.

Useful7/10
Difficulty5/10
Novelty6/10
Paper: Performance evaluation of branch-free fused multiply-add algorithms for multi-component-type multiple-precision floating-point arithmetic arXiv:2607.11391
✗ Failed on benchmark 2026

Hard-Constrained Bernstein Memory Head

Add a causal memory branch whose lag-response function is represented by a Bernstein polynomial with coefficients constrained to produce a nonnegative, decreasing, convex kernel. The branch aggregates past hidden states using this kernel, giving the model a learnable long-memory profile while preventing oscillatory, negative, or increasing historical influence.

Useful7/10
Difficulty4/10
Novelty7/10
Paper: Neural Discovery of Memory and Nonlocal Kernels in Integro-Differential Equations with Constrained Kolmogorov--Arnold Networks arXiv:2607.11110
△ Mechanism confirmed, baseline not beaten 2026

Tolerance-controlled adaptive low-rank layers

Replace selected dense neural-network operators by low-rank factors whose rank is selected by a randomized residual test at a user-specified tolerance. Construct candidate bases in large blocks for efficient matrix operations, then prune the block to the smallest rank that passes the residual criterion instead of treating the block size as the final rank.

Useful7/10
Difficulty5/10
Novelty6/10
Paper: Adaptive, Matrix-Free Low-Rank Approximation arXiv:2607.06758
✗ Failed on benchmark 2026

Conservative Sparse Mortar Cross-Attention

Replace dense coarse-to-fine cross-attention at multiresolution interfaces with a sparse, nonnegative overlap operator whose weighted feature average is exactly conserved between the two resolutions. Use this operator as a low-order path and blend it with an unrestricted neural cross-attention path through a convex limiter that keeps features inside a box or simplex domain. The construction is especially suitable for adaptive token grids, hierarchical graph neural networks, neural operators…

Useful7/10
Difficulty5/10
Novelty8/10
Paper: Invariant-domain-preserving limiting with Adaptive Mesh Refinement for Legendre-Gauss-Lobatto Discontinuous Galerkin Spectral Element Methods arXiv:2607.06045
△ Mechanism confirmed, baseline not beaten 2026

Star-Delta Hub Elimination

Remove a latent relay or hub token from an attention or graph layer and replace its two-hop influence by direct effective edges between retained tokens. The correction is a normalized rank-one update, so it can preserve hub-mediated communication while reducing the number of stored and processed states.

Useful7/10
Difficulty5/10
Novelty7/10
Paper: The Invariant Measure of Multiscale Markov Chains via Fast Arborescence Factorization arXiv:2606.31596
✓✓ Beats tuned baseline 2026

Spectral-gap local mixing

Replace a dense graph-attention or token-mixing matrix by a resolvent-like interaction operator and truncate it to graph neighborhoods whose radius is selected from an estimated spectral gap. Unlike fixed-window sparse attention, the sparsity level is tied to a measurable stability parameter and has an explicit exponential tail criterion.

Useful7/10
Difficulty6/10
Novelty6/10
Paper: Waveguiding in systems of high contrast resonators: Theory and fast computations arXiv:2608.26906
✗ Mechanism failed 2026

Samplet-compressed kernel interaction layer

Replace a dense coordinate-kernel interaction among N points by an orthogonal samplet transform with a sparse detail-detail matrix and a small polynomial branch. Detail basis vectors have vanishing moments, so smooth low-frequency behavior is represented by a few polynomial coefficients while localized residual interactions become sparse in the transformed domain.

Useful7/10
Difficulty6/10
Novelty7/10
Paper: Samplet compression for conditionally positive definite kernels and universal Kriging arXiv:2608.24283
Unverified 2026

Critical Heavy-Tail Propagation Layer

Insert a sparse heavy-tailed nonlocal mixing operator into a residual sequence, graph, or spatial network so that information can traverse distant positions without stacking many local layers. Use the critical tail exponent s = 1/2, whose truncated first moment grows logarithmically and predicts an effective propagation distance proportional to depth times log depth rather than merely depth.

Useful6/10
Difficulty5/10
Novelty7/10
Paper: Propagation rates in integro-differential equations of ignition type arXiv:2608.28450
Unverified 2026

Concatenated-Neighbor Low-Rank Message Passing

Replace a sparse graph layer's separate edge transformations with one joint low-rank factorization of all transformations entering each target node. For target node i, concatenate the neighbor matrices horizontally, project all neighbor features into a shared low-dimensional receiving basis, and reconstruct one output; retain the self transformation exactly. This can reduce edge-parameter storage and message-passing FLOPs when the incoming block row has rapidly decaying singular values.

Useful6/10
Difficulty5/10
Novelty6/10
Paper: Accelerated S-NFC for Million-Chaff RCS Computation Using Low-Rank Compression of Concatenated Block Rows arXiv:2608.27936
Unverified 2026

Indefinite Sketched Nyström Interaction

Approximate a dense symmetric interaction matrix in a neural layer by \(\widehat A=C\widehat M C^{\top}\), but compute the small core \(\widehat M\) from a two-sided sketched least-squares fit rather than from the landmark principal submatrix. This preserves signed or indefinite directions and avoids exploding outputs caused by an almost-singular \(A(I,I)\).

Useful6/10
Difficulty5/10
Novelty6/10
Paper: Nyström method for symmetric indefinite matrices arXiv:2608.20531
Unverified 2026

Hierarchical Polynomial Kernel Attention

Construct a sparse attention variant whose positional interaction kernel is piecewise polynomial rather than an arbitrary softmax score. Store key-value content in exact monomial moments within a spatial or learned-coordinate box hierarchy, then evaluate all queries from those moments without materializing the query-by-key matrix.

Useful6/10
Difficulty6/10
Novelty6/10
Paper: Exact hierarchical algorithms for accelerating particle--mesh coupling in sparse-grid particle-in-cell methods arXiv:2608.19702
Unverified 2026

Local-Tail Adaptive Attention

Replace a global attention truncation rule with a per-query local-tail budget. For each query token, retain nearby or high-priority keys until the estimated discarded interaction strength is below a target epsilon; this uses the paper's central distinction between local tail mass and the extensive norm of the discarded operator. The resulting attention pattern can allocate long-range computation only to tokens whose local tail is large.

Useful6/10
Difficulty5/10
Novelty5/10
Paper: Local observable errors from truncating interaction tails in gapped quantum lattice systems arXiv:2608.15576
Unverified 2026

Quadratic-budget Toeplitz long-range layer

Replace a dense translation-invariant interaction matrix with a positive-definite Toeplitz kernel K_n(e^f) whose log-spectrum is parameterized by a small number of Fourier coefficients with 1/|k| decay. Use the paper's explicit quadratic term as a spectral-volume budget, allowing long-range structure while discouraging uncontrolled determinant growth and ill-conditioning. Subtracting this term from a log-determinant regularizer leaves a residual intended to capture higher-order deviations from…

Useful6/10
Difficulty6/10
Novelty7/10
Paper: On Toeplitz determinants with slow Fourier decay arXiv:2608.13182
Unverified 2026

Power-of-Two Latency Attention

Make Laplacian latency attention hardware-friendly by quantizing the exponential decay and row normalization to powers of two. Exponentials and divisions can then be implemented with lookup tables and integer shifts, enabling a low-bit attention kernel for neuromorphic and edge accelerators.

Useful6/10
Difficulty4/10
Novelty7/10
Paper: Lapis: Laplacian Spiking Attention via First-Spike Timing and Membrane Leakage arXiv:2608.11865
Unverified 2026

Random-Batch Interaction Layer

Replace an all-pairs interaction or attention-like message field with an unbiased estimator obtained from a uniformly random partition into batches. Keep the partition fixed for \(\tau\) integration steps, then redraw it; the paper's error law predicts that too-long reuse causes a quadratic-in-\(\tau\) error accumulation, while finite-population error decreases as \(N^{-1}\).

Useful6/10
Difficulty4/10
Novelty5/10
Paper: Mean field error estimate of the random batch method for vortex blob dynamics for the 2D Navier--Stokes Equation arXiv:2608.06533
Unverified 2026

Outlier-Spectral Landmark Attention

Construct a symmetric token affinity matrix and approximate only its spectrally outlying token-mixing modes using a small set of sampled landmark columns. The resulting low-rank operator replaces an \(O(n^2)\) dense mixer by two skinny matrix multiplications, while the paper's residual guarantee predicts that large-magnitude global interaction modes are preserved.

Useful6/10
Difficulty5/10
Novelty5/10
Paper: Sublinear Time Eigenvector Approximation via Column Sampling arXiv:2608.00840