Operator Learning at Machine Precision
Not on the map. Why?
Abstract (arXiv)
Neural operator learning methods have garnered significant attention in scientific computing for their ability to approximate infinite-dimensional operators. However, increasing their complexity often fails to substantially improve their accuracy, leaving them on par with much simpler approaches such as kernel methods and more traditional reduced-order models. In this article, we set out to address this shortcoming and introduce CHONKNORIS (Cholesky Newton--Kantorovich Neural Operator Residual Iterative System), an operator learning paradigm that can achieve machine precision. CHONKNORIS draws on numerical analysis: many nonlinear forward and inverse PDE problems are solvable by Newton-type methods. Rather than regressing the solution operator itself, our method regresses the Cholesky factors of the elliptic operator associated with Tikhonov-regularized Newton--Kantorovich updates. The resulting unrolled iteration yields a neural architecture whose machine-precision behavior follows from achieving a contractive map, requiring far lower accuracy than end-to-end approximation of the solution operator. We benchmark CHONKNORIS on a range of nonlinear forward and inverse problems, including a nonlinear elliptic equation, Burgers' equation, a nonlinear Darcy flow problem, the Calder'on problem, an inverse wave scattering problem, and a problem from seismic imaging. We also present theoretical guarantees for the convergence of CHONKNORIS in terms of the accuracy of the emulated Cholesky factors. Additionally, we introduce a foundation model variant, FONKNORIS (Foundation Newton--Kantorovich Neural Operator Residual Iterative System), which aggregates multiple pre-trained CHONKNORIS experts for diverse PDEs to emulate the solution map of a novel nonlinear PDE. Our FONKNORIS model is able to accurately solve unseen nonlinear PDEs such as the Klein--Gordon and Sine--Gordon equations.
Key facts. Cited by 7 papers in the map (data through 14 Oct 2026). Builds on 43 of its 75 references that are in the map. First posted on arXiv 25 Nov 2025.
Builds on
43 of 75 references.
References · 43
- 2017Attention Is All You Needcited by 60,361
- 2020An Image is Worth 16x16 Words: Transformers for Image Recognition at Scalecited by 30,186
- 2023LLaMA: Open and Efficient Foundation Language Modelscited by 14,132
- 2020Fourier Neural Operator for Parametric Partial Differential Equationscited by 2,908
- 2019Optuna: A Next-generation Hyperparameter Optimization Frameworkcited by 2,501
- 2019DeepONet: Learning nonlinear operators for identifying differential equations based on the universal approximation theorem of operatorscited by 2,018
- 2018GPyTorch: Blackbox Matrix-Matrix Gaussian Process Inference with GPU Accelerationcited by 609
- 2022Scientific Machine Learning through Physics-Informed Neural Networks: Where we are and What's nextcited by 608
- 2021Learning the solution operator of parametric partial differential equations with physics-informed DeepOnetscited by 540
- 2021Physics-Informed Neural Operator for Learning Partial Differential Equationscited by 527
- 2019Algorithm Unrolling: Interpretable, Efficient Deep Learning for Signal and Image Processingcited by 421
- 2021Choose a Transformer: Fourier or Galerkincited by 343
31 more
- 2023GNOT: A General Neural Operator Transformer for Operator Learningcited by 336
- 2022Transformer for Partial Differential Equations' Operator Learningcited by 283
- 2024Poseidon: Efficient Foundation Models for PDEscited by 200
- 2021Solving and Learning Nonlinear PDEs with Gaussian Processescited by 129
- 2022Physics-Informed Deep Neural Operator Networkscited by 113
- 2020Estimates on the generalization error of Physics Informed Neural Networks (PINNs) for approximating a class of inverse problems for PDEscited by 102
- 2020The Random Feature Model for Input-Output Maps between Banach Spacescited by 90
- 2022Learning Operators with Coupled Attentioncited by 90
- 2024Operator Learning: Algorithms and Analysiscited by 85
- 2023A Mathematical Guide to Operator Learningcited by 76
- 2023Neural Inverse Operators for Solving PDE Inverse Problemscited by 65
- 2015Operator-valued Kernels for Learning from Functional Response Datacited by 60
- 2024Pretraining Codomain Attention Neural Operators for Solving Multiphysics PDEscited by 59
- 2023Kernel Methods are Competitive for Operator Learningcited by 50
- 2018Kernel Flows: from learning kernels from data into the abysscited by 49
- 2021OpenFWI: Large-Scale Multi-Structural Benchmark Datasets for Seismic Full Waveform Inversioncited by 44
- 2020Sparse Cholesky factorization by Kullback-Leibler minimizationcited by 41
- 2023MgNO: Efficient Parameterization of Linear Operators via Multigridcited by 40
- 2024Continuum Attention for Neural Operatorscited by 40
- 2024Numerical analysis of physics-informed neural networks and related models in physics-informed machine learningcited by 39
- 2022Randomly pivoted Cholesky: Practical approximation of a kernel matrix with few entry evaluationscited by 32
- 2023Sparse Cholesky Factorization for Solving Nonlinear PDEs via Gaussian Processescited by 31
- 2021Computational Graph Completioncited by 28
- 2024Gauss-Newton Natural Gradient Descent for Physics-Informed Computational Fluid Dynamicscited by 28
- 2017Compression, inversion, and approximate PCA of dense kernel matrices at near-linear computational complexitycited by 24
- 2016Nested Kriging predictions for datasets with large number of observationscited by 18
- 2024Operator Learning with Gaussian Processescited by 18
- 2020Do ideas have shape? Idea registration as the continuous limit of artificial neural networkscited by 16
- 2024RandONet: Shallow-Networks with Random Projections for learning linear and nonlinear operatorscited by 16
- 2015Multigrid with rough coefficients and Multiresolution operator decomposition from Hierarchical Information Gamescited by 15
- 2025Data-Efficient Kernel Methods for Learning Differential Equations and Their Solution Operators: Algorithms and Error Analysiscited by 7
Built on by
Citing papers per year: 2025: 1; 2026: 6.
7 papers.
- 2025Convolutional-neural-operator-based transfer learning for solving PDEscited by 2
- 2026Spectrally Safe Neural Operator Warm-Starts for Large-Scale Newton Solverscited by 1
- 2026Deflation-PINNs: Learning Multiple Solutions for PDEs and Landau-de Gennescited by 0
- 2026Learning Contractive Integral Operators with Fredholm Integral Neural Operatorscited by 0
- 2026Kernel Methods for Learning Operators with Multiple Inputs and Outputscited by 0
- 2026A Systematic Analysis of Automatic Differentiation versus Discretization-based Constraints for Physics-Informed PDE Solverscited by 0
- 2026Transferability of Learned States in Neural PDE Solverscited by 0
- 2026Transferability of Learned States in Neural PDE Solvers7 Oct 2026
- 2026A Systematic Analysis of Automatic Differentiation versus Discretization-based Constraints for Physics-Informed PDE Solvers7 Sep 2026
- 2026Kernel Methods for Learning Operators with Multiple Inputs and Outputs12 Aug 2026
- 2026Spectrally Safe Neural Operator Warm-Starts for Large-Scale Newton Solvers20 Jun 2026
- 2026Learning Contractive Integral Operators with Fredholm Integral Neural Operators3 Apr 2026
- 2026Deflation-PINNs: Learning Multiple Solutions for PDEs and Landau-de Gennes30 Mar 2026
- 2025Convolutional-neural-operator-based transfer learning for solving PDEs19 Dec 2025