Abstract:Operator learning, the approximation of mappings between infinite-dimensional function spaces using ideas from machine learning, has gained increasing research attention in recent years. Approximate operators, learned from data, hold promise to serve as efficient surrogate models for problems in computational science and engineering, complementing traditional numerical methods. However, despite their empirical success, our understanding of the underpinning mathematical theory is in large part still incomplete. In this paper, we study the approximation of Lipschitz operators in expectation with respect to Gaussian measures. We prove higher Gaussian Sobolev regularity of Lipschitz operators and establish lower and upper bounds on the Hermite polynomial approximation error. We further consider the reconstruction of Lipschitz operators from $m$ arbitrary (adaptive) linear samples. A key finding is the tight characterization of the smallest achievable error for all possible (adaptive) sampling and reconstruction maps in terms of $m$. It is shown that Hermite polynomial approximation is an optimal recovery strategy, but we have the following curse of sample complexity: No method to approximate Lipschitz operators based on finitely many samples can achieve algebraic convergence rates in $m$. On the positive side, we prove that a sufficiently fast spectral decay of the covariance operator of the Gaussian measure guarantees convergence rates which are arbitrarily close to any algebraic rate in the large data limit $m \to \infty$. Finally, we focus on the recovery of Lipschitz operators from finitely many point samples. We consider Christoffel sampling and weighted least-squares approximation, and present an algorithm which provably achieves near-optimal sample complexity.

Higher Order Lipschitz Greedy Recombination Interpolation Method (HOLGRIM)

Greedy Recombination Interpolation Method (GRIM)

A Directional Interpolation Pattern for Multilevel Green's Function Interpolation Method

A generalization of Floater--Hormann interpolants

IG-LSPIA: Least Squares Progressive Iterative Approximation for Isogeometric Collocation Method

Learning Lipschitz Operators with respect to Gaussian Measures with Near-Optimal Sample Complexity

Interpolation via weighted $l_1$ minimization

Probabilistic Richardson Extrapolation

Splitting the Conditional Gradient Algorithm

Interpolating Greedy and Reluctant Algorithms

Algorithms for computing the optimal Lipschitz constant of interpolants with Lipschitz derivative

Greedy Adaptive Local Recovery of Functions in Sobolev Spaces

On the Optimality of Target-Data-Dependent Kernel Greedy Interpolation in Sobolev Reproducing Kernel Hilbert Spaces

Proximal-Like Incremental Aggregated Gradient Method with Linear Convergence under Bregman Distance Growth Conditions

Convergence of ease-controlled Random Reshuffling gradient Algorithms under Lipschitz smoothness

On the convergence of generalized kernel-based interpolation by greedy data selection algorithms

Construction techniques for highly accurate quasi-interpolation operators

Interpolatory H-infinity Model Reduction

Restarted Nonconvex Accelerated Gradient Descent: No More Polylogarithmic Factor in the O(<SUP>-7/4</SUP>) Complexity

Iteratively reweighted $\ell_1$ algorithms with extrapolation

Numerical Gradient Schemes for Heat Equations Based on the Collocation Polynomial and Hermite Interpolation