Training ideas

Curricula, sampling schedules, pseudo-labeling and multi-task mixing, with the standard training recipe as the control.

△ Mechanism confirmed, baseline not beaten 2026

Matrix-Free Krylov Backpropagation Through Solver Layers

Turn an iterative optimization or equilibrium computation inside a neural network into a differentiable layer whose backward pass solves the implicit adjoint system with conjugate gradients or GMRES using only automatic-differentiation matrix-vector products. This avoids storing unrolled iterations and avoids explicit Hessian or Jacobian construction, enabling longer solver horizons and lower-memory implicit architectures.

Useful8/10
Difficulty6/10
Novelty5/10
Paper: PANDA: A Matrix-Free Differentiable NMPC Solver via Proximal Averaged Quasi-Newton with Adaptive Linesearch Algorithm arXiv:2608.16280
✗ Failed on benchmark 2026

Localized Petrov–Galerkin Neural Residuals

Replace the pointwise strong-form PINN loss with a vector of localized weak residuals generated by fixed compactly supported polynomial test functions. Use a neural network or KAN as the trial function, integrate by parts once, and evaluate each test residual with Gauss–Legendre quadrature; this lowers the required derivative order and prevents a few high-curvature collocation points from dominating training.

Useful8/10
Difficulty5/10
Novelty5/10
Paper: PG-KINN: A Physics-Informed Petrov-Galerkin Kolmogorov-Arnold Network for Solving Forward and Inverse PDEs arXiv:2607.20378
✗ Failed on benchmark 2026

Second-Order Brownian Jet Residual

Replace pointwise high-order PINN residuals with a stochastic one-step residual evaluated on Brownian transitions. A single scalar network produces the value, gradient, and Hessian by automatic differentiation, and the quadratic centered increment supplies a stochastic probe of the Hessian. Add a terminal gradient penalty so the learned full jet is constrained at the terminal boundary, not only the scalar value.

Useful8/10
Difficulty5/10
Novelty7/10
Paper: A Deep Second-Order Stochastic Residual Method for Fully Nonlinear Parabolic PDEs arXiv:2607.16730
✗ Mechanism failed 2026

Fenchel-Gap Certified Neural PDE Training

Train a primal state network and a dual flux network jointly, using the convex primal-dual gap as the main loss and as an a posteriori certificate of state error. Unlike a strong residual, the certificate is based on monotonicity and convex duality, so it can remain informative even when differentiating rapidly oscillatory coefficients would amplify noise by $1/\varepsilon$.

Useful8/10
Difficulty6/10
Novelty7/10
Paper: Non-Asymptotic Variational Learning for Monotone Nonlinear Multiscale Elliptic Equations: Scale-Robust Primal-Dual Bounds and Strong-Form Statistical Ill-Conditioning arXiv:2607.15702
△ Mechanism confirmed, baseline not beaten 2026

Derivative-Free Very-Weak Neural PDE Solver

Train a neural trial function for an elliptic PDE using a very-weak residual in which all derivatives act on fixed smooth test functions rather than on the neural network. This eliminates second-order reverse-mode or forward-mode automatic differentiation and allows low-regularity activations while retaining a least-squares objective over many test functions.

Useful8/10
Difficulty4/10
Novelty6/10
Paper: Neural Very Weak Formulations enabling Hardware-Oriented deep PDE solvers arXiv:2607.14498
△ Mechanism confirmed, baseline not beaten 2026

Residual-Christoffel Collocation for Random-Feature PDE Networks

Replace uniform collocation for a fixed random-feature neural PDE solver with sampling from the leverage-score density of the operator-applied features. Whiten the retained residual feature space before solving for output coefficients, so the sampled least-squares matrix has an identity-like expected Gram rather than inheriting severe anisotropy from the differential operator. The same construction can be used for a linearized neural network by treating Jacobian features as the trial functions.

Useful8/10
Difficulty5/10
Novelty7/10
Paper: Residual-Christoffel Sampling for Random Feature Collocation of Linear PDEs arXiv:2607.13382
△ Mechanism confirmed, baseline not beaten 2026

Walk-on-Spheres stochastic target layer

Train a neural network to represent an elliptic solution using Walk-on-Spheres rollouts as stochastic targets instead of evaluating a mesh-based PDE residual. For each input point, recursively jump to a random point on the largest interior sphere, accumulate source contributions, evaluate boundary data at termination, and regress the network output to the resulting Monte Carlo estimate.

Useful7/10
Difficulty5/10
Novelty6/10
Paper: Walk-on-Spheres Monte Carlo and deep neural network approximations of elliptic PDEs with drift and killing arXiv:2608.09494
△ Mechanism confirmed, baseline not beaten 2026

Quasi-uniform residual least squares

Train a frozen-feature or linearized neural network by residual least squares on deterministic quasi-uniform points rather than independently sampled collocation points. The paper's norm-equivalence result predicts that, once the number of residual points is proportional to the number of active features, the empirical residual controls the continuous residual and avoids random undersampling of localized errors.

Useful7/10
Difficulty4/10
Novelty6/10
Paper: Optimal Neural Network Approximation via Empirical Least Squares with Deterministic Samples arXiv:2608.06687
Unverified 2026

Actuator-Aware Envelope Scheduler

Use adaptive performance specifications to prevent a neural controller or policy from demanding output changes that exceed bounded actuator amplitude or action-rate limits. The target error envelope tightens when the policy has control authority and relaxes when saturation or rate clipping persists, instead of allowing the controller to destabilize while chasing an infeasible target. This converts actuator clipping into an explicit slow state that can be used by reinforcement-learning policies…

Useful6/10
Difficulty4/10
Novelty7/10
Paper: Output Feedback Adaptive Performance Control arXiv:2608.15758
Unverified 2026

Scaled-boundary quadrature loss

Replace random interior collocation points in domain-integrated neural losses with deterministic points generated by scaling oriented boundary patches toward an interior center. The resulting quadrature exactly integrates low-degree polynomial integrands on affine polytopes and gives a mesh-free, geometry-aware loss for PINNs and neural fields.

Useful6/10
Difficulty4/10
Novelty7/10
Paper: Scaled boundary cubature scheme in higher dimensions: integration over polytopes and curved regions arXiv:2608.15563
Unverified 2026

Pivot-safe discrete LU initialization

Initialize an invertible neural linear layer from a bounded discrete random matrix only after checking that every leading principal submatrix is nonsingular and that its LU growth factor is below a prescribed threshold. This replaces blind random initialization with a cheap resampling rule designed to prevent zero pivots and excessive finite-precision amplification in reversible or flow-based networks.

Useful6/10
Difficulty4/10
Novelty7/10
Paper: LU Factorization of Discrete Random Matrices arXiv:2608.08998
Unverified 2026

Gap-Graded Neural Collocation

Replace uniform PINN or neural-operator collocation by a graded point distribution concentrated in narrow regions between nearly touching interfaces. Use the paper's distance-dependent mesh scale to determine point spacing, and switch to a gap-dependent minimum scale when the separation becomes too small for the global mesh.

Useful5/10
Difficulty4/10
Novelty6/10
Paper: High-order finite element method for perfect conductivity and linear elasticity with nearly touching inclusions arXiv:2607.22128