Organizations: Department of Mathematics, Ludwig-Maximilians-Universität München (LMU Munich), Akademiestr. 7, 80799 München · Department of Mathematics, City University of Hong Kong, Hong Kong, China · School of Mathematics and Statistics, The University of Sydney, Sydney, New South Wales 2006, Australia · Munich Center for Machine Learning (MCML), Germany · University of Tromsø, Norway · German Aerospace Center (DLR), Germany
Deep neural networks have emerged as powerful tools for learning operators defined over infinite-dimensional function spaces. However, existing theories frequently encounter difficulties related to dimensionality and limited interpretability. This work investigates how sparsity can help address these challenges in functional learning, a central ingredient in operator learning. We propose a framework that employs convolutional architectures to extract sparse features from a finite number of samples, together with deep fully connected networks to effectively approximate nonlinear functionals. Using universal discretization methods, we show that sparse approximators enable stable recovery from discrete samples. In addition, both the deterministic and the random sampling schemes are sufficient for our analysis. These findings lead to improved approximation rates and reduced sample sizes in various function spaces, including those with fast frequency decay and mixed smoothness. They also provide new theoretical insights into how sparsity can alleviate the curse of dimensionality in functional learning.
Figures & tables
Figure 1: Functional learning pipeline proposed in Theorem 7 through functional neural networks.
Input functions
Encoder
Rate
Wpα(Sd−1) ( Yang et al., 2026 )
Spherical harmonics
O~((logK)−βα/d−1)
Cα([−1,1]d) ( Song et al., 2023a )
Legendre polynomials
O~((logK)−βα/d)
Analytic functions ( Song et al., 2023b )
O~(exp(−c(logM)d+11))
B2,∞α(Rd),H0α([0,1]d) ( Shi et al., 2025 )
Eigensystem induced by Mercer kernel
O~((logK)−βα/d)
Gaussian RKHS ( Shi et al., 2025 )
O~(exp(−c(logM)d+11))
A1α(DN) ( ours )
CNN
O~((logK)−β(α−3/2))
Table 1: Comparison of approximation rates of Hölder continuous functionals across various input function spaces and encoding schemes. We employ O~ notation in place of standard O notation to emphasize the leading term, suppressing the loglog factor.