mlsanityData to 14 Oct 2026

Operator Learning at Machine Precision

arXiv 2511.19980 ↗First posted 25 Nov 2025cs.LG cs.NA math.NAcited by 7

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

  1. 2017Attention Is All You Needcited by 60,361
  2. 2020An Image is Worth 16x16 Words: Transformers for Image Recognition at Scalecited by 30,186
  3. 2023LLaMA: Open and Efficient Foundation Language Modelscited by 14,132
  4. 2020Fourier Neural Operator for Parametric Partial Differential Equationscited by 2,908
  5. 2019Optuna: A Next-generation Hyperparameter Optimization Frameworkcited by 2,501
  6. 2019DeepONet: Learning nonlinear operators for identifying differential equations based on the universal approximation theorem of operatorscited by 2,018
  7. 2018GPyTorch: Blackbox Matrix-Matrix Gaussian Process Inference with GPU Accelerationcited by 609
  8. 2022Scientific Machine Learning through Physics-Informed Neural Networks: Where we are and What's nextcited by 608
  9. 2021Learning the solution operator of parametric partial differential equations with physics-informed DeepOnetscited by 540
  10. 2021Physics-Informed Neural Operator for Learning Partial Differential Equationscited by 527
  11. 2019Algorithm Unrolling: Interpretable, Efficient Deep Learning for Signal and Image Processingcited by 421
  12. 2021Choose a Transformer: Fourier or Galerkincited by 343
31 more
  1. 2023GNOT: A General Neural Operator Transformer for Operator Learningcited by 336
  2. 2022Transformer for Partial Differential Equations' Operator Learningcited by 283
  3. 2024Poseidon: Efficient Foundation Models for PDEscited by 200
  4. 2021Solving and Learning Nonlinear PDEs with Gaussian Processescited by 129
  5. 2022Physics-Informed Deep Neural Operator Networkscited by 113
  6. 2020Estimates on the generalization error of Physics Informed Neural Networks (PINNs) for approximating a class of inverse problems for PDEscited by 102
  7. 2020The Random Feature Model for Input-Output Maps between Banach Spacescited by 90
  8. 2022Learning Operators with Coupled Attentioncited by 90
  9. 2024Operator Learning: Algorithms and Analysiscited by 85
  10. 2023A Mathematical Guide to Operator Learningcited by 76
  11. 2023Neural Inverse Operators for Solving PDE Inverse Problemscited by 65
  12. 2015Operator-valued Kernels for Learning from Functional Response Datacited by 60
  13. 2024Pretraining Codomain Attention Neural Operators for Solving Multiphysics PDEscited by 59
  14. 2023Kernel Methods are Competitive for Operator Learningcited by 50
  15. 2018Kernel Flows: from learning kernels from data into the abysscited by 49
  16. 2021OpenFWI: Large-Scale Multi-Structural Benchmark Datasets for Seismic Full Waveform Inversioncited by 44
  17. 2020Sparse Cholesky factorization by Kullback-Leibler minimizationcited by 41
  18. 2023MgNO: Efficient Parameterization of Linear Operators via Multigridcited by 40
  19. 2024Continuum Attention for Neural Operatorscited by 40
  20. 2024Numerical analysis of physics-informed neural networks and related models in physics-informed machine learningcited by 39
  21. 2022Randomly pivoted Cholesky: Practical approximation of a kernel matrix with few entry evaluationscited by 32
  22. 2023Sparse Cholesky Factorization for Solving Nonlinear PDEs via Gaussian Processescited by 31
  23. 2021Computational Graph Completioncited by 28
  24. 2024Gauss-Newton Natural Gradient Descent for Physics-Informed Computational Fluid Dynamicscited by 28
  25. 2017Compression, inversion, and approximate PCA of dense kernel matrices at near-linear computational complexitycited by 24
  26. 2016Nested Kriging predictions for datasets with large number of observationscited by 18
  27. 2024Operator Learning with Gaussian Processescited by 18
  28. 2020Do ideas have shape? Idea registration as the continuous limit of artificial neural networkscited by 16
  29. 2024RandONet: Shallow-Networks with Random Projections for learning linear and nonlinear operatorscited by 16
  30. 2015Multigrid with rough coefficients and Multiresolution operator decomposition from Hierarchical Information Gamescited by 15
  31. 2025Data-Efficient Kernel Methods for Learning Differential Equations and Their Solution Operators: Algorithms and Error Analysiscited by 7

Built on by

Cited by year (2026 to 14 Oct)

Citing papers per year: 2025: 1; 2026: 6.

7 papers.