Every idea extracted from recent arXiv mathematics papers — verified and unverified. Click an idea to open its full card; badges show the empirical verdict.
Approximate a long hidden-state or key/value trajectory by a small sum of exponential modes, storing only the mode nodes and vector weights instead of every timestep. Recover the modes from a Hankel matrix through an annihilating polynomial, then reconstruct the trajectory or a compact recurrent state during inference.
Apply a Birkhoff-normal-form-inspired monitor to momentum optimization and recurrent-state updates, where oscillatory modes are identified from recent parameter or hidden-state trajectories. When two dominant frequencies approach a low-order ratio such as 2:1 or 3:1, increase damping before nonlinear mode coupling produces large oscillations; away from resonance, retain the faster low-damping update.
Track where the loss Hessian's eigenvectors are most sensitive to the current minibatch perturbation, rather than using only eigenvalues or a global learning-rate estimate. Apply extra damping only to spectral bands with high geometric response, allowing flat and well-separated curvature modes to retain a larger step size.
Maintain an SPD matrix preconditioner with the paper's deterministic Ornstein–Uhlenbeck covariance recursion rather than estimating an inverse through Newton–Schulz or an explicit matrix inverse. Apply this preconditioner to gradients from a small layer block, using damping and a conservative step size to preserve positive definiteness. The method is most plausible for low-rank, per-layer, or blockwise curvature matrices where dense matrix storage is affordable.
Use a rapidly cycling preconditioner or learning-rate vector during optimization, but construct a static averaged optimizer with the same mean update. When the parameter dynamics are locally contractive, the averaged optimizer should track the periodic optimizer while requiring less schedule bookkeeping and potentially fewer expensive state updates.
Replace the purely diagonal preconditioner in AdamW or SGD with a blockwise, single-secant BFGS inverse-curvature metric. Use spectral damping and clipping relative to the diagonal RMS metric so the learned metric cannot become arbitrarily ill-conditioned, mirroring the paper's uniform comparison between its conjugate-free scaling and the primal barrier Hessian.
Replace the scalar learning rate of SGD or Adam's outer update by a blockwise Barzilai--Borwein estimate computed from consecutive parameters and gradients. Use gradient smoothing, denominator checks, and clipping so that the curvature estimate remains usable with stochastic neural-network gradients.
Replace real diagonal state-space channels with complex damped oscillators whose hidden states encode both amplitude and phase. Train with parallel causal convolution and deploy with the equivalent one-step recurrence, allowing the same layer to support efficient batched training and low-memory streaming inference.
Replace dense token-to-token attention in selected layers with communication through a small number of multi-token hyperedges. Each hyperedge aggregates its incident token states and broadcasts the resulting message back to those tokens, allowing higher-order interactions while reducing the number of pairwise links. Reconstruct hyperedges periodically from cumulative token displacement so stable tokens retain useful groups while rapidly changing tokens are regrouped.
Replace a learned dense token-mixing matrix or residual-state transition with a sparse diffusive mixer whose Laplacian has a deliberately small largest Jordan block. Balance the two chain lengths around the central coupling/core, because the paper proves that this minimizes the worst defective transient among the tridiagonal family. Use a scalar residual step size to move the non-consensus spectrum inside the unit disk while preserving the sparse structure.
Implement f(A)b inside a neural network with a short Lanczos recurrence instead of an eigendecomposition or dense matrix-function operation. Use an SPD operator A such as a regularized graph Laplacian or feature covariance matrix, and choose the number of iterations by monitoring successive approximations. For Stieltjes functions, Lanczos is guaranteed to be close to the best vector in the same Krylov subspace.
Replace fixed-period federated averaging or distributed all-reduce with a Bernoulli communication decision whose probability is selected from estimated network connectivity and optimization conditioning. Local workers continue making corrected updates between communication events, while the contraction theorem exposes when communication is worth its cost.
Represent periodic input-output behavior using a compact real vector of Fourier coefficients and learn an invertible neural map from input coefficients to output coefficients. Inference then obtains the input representation for a desired periodic output by a single inverse pass instead of iterative optimization through a nonlinear forward model, while the Fourier representation reduces sequence dimensionality when high-rate signals are spectrally sparse.
Use the Bregman objective's exact residual-dependent curvature to build a positive-semidefinite Gauss-Newton preconditioner for a neural network's scalar regression head. Negative curvature weights are clipped or damped before solving the update, preserving the original gradient while preventing residual patterns from producing unstable parameter steps.
Train a finite-element surrogate by minimizing the assembled discrete potential energy rather than a loss against solved displacement labels. The objective uses only the sparse stiffness matrix and load vector, while its exact energy-gap identity makes it equivalent to supervised regression in the stiffness norm.
Replace repeated full-dimensional matrix-exponential or ODE solves in a conditioned continuous-time state-space layer with contour quadrature evaluated in a projection basis. The same reduced basis and contour nodes can serve many conditioning vectors, while shifted reduced resolvents provide a stable and differentiable approximation over a prescribed time window.
Add a periodic coarse optimization phase to SGD or Adam that operates on a compressed parameterization and returns a prolongated correction to the full network. Retain nonsmooth constraints or regularizers explicitly through a primal-dual update instead of relying on penalty smoothing. Accept the correction only when it improves a cheap fine-batch merit test, making the method useful even when the coarse objective is only approximately coherent with the fine objective.
Use enumerated weighing matrices as sparse orthogonal channel-mixing operators inside MLPs or residual blocks. Their ternary entries reduce multiplication to signed additions, while exact orthogonality prevents amplification or attenuation of feature norms; a trainable fixed-support version can recover expressivity without giving up computational sparsity.
Replace exhaustive optimization of N binary gates by the geometrically admissible sign patterns induced by projections onto a common direction. For two-dimensional gate vectors, enumerate angular cells exactly; for higher-dimensional vectors, sample directions and evaluate only the induced configurations.
Replace an explicit gradient step by an implicit correction using the trajectory derivative \(Dg(\theta)g(\theta)=H(\theta)g(\theta)\), where \(g=\nabla f\) and \(H=\nabla^2 f\). The update should strongly damp high-curvature or stiff modes while preserving fourth-order matching of the local linearized dynamics. Start with a self-contained fourth-order L-stable rational prototype, then compare it with the paper's exact two-stage coefficients after recovering those coefficients from the full…
Replace ordinary reverse-mode differentiation through a long sequence of stiff relaxation updates with a projected adjoint that separates slow conserved features from rapidly relaxing residual features. The neural layer can use large outer time steps even when its internal relaxation time is very small, while reconstructing only the microscopic gradient component required by the preceding layer.
Replace a delayed optimizer state or gradient by a causal lower-triangular history transformation that predicts the current descent direction from recently stored states and inputs. Use Fredholm terms to incorporate the recent history and Volterra terms to preserve causal invertibility, then apply the optimizer update in transformed coordinates. This targets oscillation and divergence caused by concurrent delays in distributed or asynchronous training.
Use a strictly upper-triangular block operator to represent interactions between ordered layers, experts, or token groups, and compute its inverse exactly with a finite Neumann series. Because nilpotency truncates the series after a known number of blocks, the module avoids an iterative solver while retaining controlled long-range interactions.
Replace heuristic magnitude pruning in a layerwise convex reconstruction problem with safe screening based on a perspective relaxation of the cardinality constraint. A channel is removed only when a lower bound for every solution containing that channel exceeds the loss of a feasible incumbent; conversely, a channel is forced to remain when every solution excluding it is provably worse.