## A direct approach for function approximation on data defined manifolds

## H. N. Mhaskar ∗

## Abstract

In much of the literature on function approximation by deep networks, the function is assumed to be defined on some known domain, such as a cube or a sphere. In practice, the data might not be dense on these domains, and therefore, the approximation theory results are observed to be too conservative. In manifold learning, one assumes instead that the data is sampled from an unknown manifold; i.e., the manifold is defined by the data itself. Function approximation on this unknown manifold is then a two stage procedure: first, one approximates the Laplace-Beltrami operator (and its eigen-decomposition) on this manifold using a graph Laplacian, and next, approximates the target function using the eigen-functions. Alternatively, one estimates first some atlas on the manifold and then uses local approximation techniques based on the local coordinate charts.

In this paper, we propose a more direct approach to function approximation on unknown , data defined manifolds without computing the eigen-decomposition of some operator or an atlas for the manifold, and without any kind of training in the classical sense. Our constructions are universal; i.e., do not require the knowledge of any prior on the target function other than continuity on the manifold. We estimate the degree of approximation. For smooth functions, the estimates do not suffer from the so-called saturation phenomenon. We demonstrate via a property called good propagation of errors how the results can be lifted for function approximation using deep networks where each channel evaluates a Gaussian network on a possibly unknown manifold.

Keywords: Manifold learning, deep networks, Gaussian networks, weighted polynomial approximation.

## 1 Introduction

One of the main problems of machine learning is the following. Given data { ( y j , f ( y j ) + glyph[epsilon1] j ) } M j =1 , where f is an unknown function, y j 's are sampled randomly from a probability distribution µ ∗ defined on a subset of R Q for some typically high dimension Q , and glyph[epsilon1] j 's are realizations of a mean zero random variable, find an approximation P from a class V n to f [21, 15, 16], where { V n } is a nested sequence of subsets of L 2 ( µ ∗ ). In practice, this approximation is found by empirical risk minimization, assuming some prior on f , such as that it belongs to some reproducing kernel Hilbert space with a known kernel, or that it has a certain number of derivatives, or that it satisfies some conditions on its Fourier transform. To set up the minimization problem, one needs to know in advance the complexity of the model P , typically, the number of parameters desired to be estimated. In theory, the usual way of estimating this number is to estimate the so called approximation error, inf P ∈ V n E µ ∗ (( f -P ) 2 ). Necessarily, this results in a fundamental gap in the theory, namely, that the minimizer of the empirical risk may have no connection with the minimizer of the approximation error.

Since the fundamental problem is one of function approximation, it is natural to wonder if appropriate tools in approximation theory can be developed in order to close this gap. One of the difficulties in doing so is that most of the results in classical approximation theory assume that the approximation takes place on a known domain, such as the cube, or Euclidean space, or sphere or similar known manifold. In turn, this requires that the data should be dense on this domain; i.e., the domain should be the (exact) support of µ ∗ . The problem is that µ ∗ being unknown, it is not possible to ensure this requirement.

During this century, manifold learning has sought to ameliorate the situation, with many practical applications. An early introduction to this topic is in the special issue [10] of Applied and Computational Harmonic Analysis, edited by Chui and Donoho. In this theory, one assumes that the support of µ ∗ is an unknown smooth compact connected manifold; for simplicity, even that µ ∗ is the Riemannian volume measure for the manifold, normalized to be a probability measure. Following, e.g., [4, 6, 5, 23, 38], one constructs first a 'graph Laplacian' from the data, and finds its eigen-decomposition. It is proved in the above mentioned papers that as the size of the data tends to infinity, the graph Laplacian converges to the Laplace-Beltrami operator on the manifold and the eigen-values

∗ Institute of Mathematical Sciences, Claremont Graduate University, Claremont, CA 91711. email: hrushikesh.mhaskar@cgu.edu. The research is supported in part by NSF grant DMS 2012355.

(respectively, eigen-vectors) converge to the corresponding quantities on the manifold. A great deal of work is devoted to studying the geometry of this unknown manifold (e.g., [22, 24]), based on the so called heat kernel. The theory of function approximation on such manifolds is also well developed (e.g., [26, 31, 32, 19, 20]).

All this work depends upon a two stage procedure - finding the eigen-decomposition of the graph Laplacian and then using approximation in terms of the eigen-vectors/eigen-functions. Once more, this leads to errors not just from the approximation of the target function but also from the approximation of the eigen-decomposition of the Laplace-Beltrami operator itself. In recent years, there are some efforts to explore alternative approaches using deep networks (e.g., [14, 11, 9, 37]). These papers also take a two-step approach: developing an atlas on the manifold first, and then using some local approximation schemes based on the local coordinate charts.

↦

Our objective in this paper is to develop a single-shot method to solve the problem, knowing only the dimension of the manifold. In particular, we aim not to find any eigen-decomposition nor to learn any atlas on the manifold, but to give a direct construction that starts with the data and constructs an approximation without involving any optimization/training and with guaranteed approximation error estimated in a probabilistic sense. Our approximation can be implemented as a Gaussian network; i.e., a function of the form x → ∑ k a k exp( - | λ x -y k | 2 2 ,Q ), where | · | 2 ,Q denotes the glyph[lscript] 2 norm on R Q . The size of the data set required depends only on the dimension of the manifold and the smoothness of the target function measured in a technical manner as explained in this paper. We will extend our results to approximation by deep Gaussian networks.

## 2 Technical introduction and outline

In this section, let us assume that the data y j is sampled from some unknown manifold, uniformly with respect to the Riemannian volume element of that manifold. One of the fundamental results in manifold learning is the following theorem of Belkin and Niyogi [6].

Theorem 2.1 Let X be a smooth, compact, q -dimensional sub-manifold of R Q , µ ∗ be its Riemannian volume measure, normalized by µ ∗ ( X ) = 1 , and ∆ denote the Laplace-Beltrami operator on X . Then for a smooth function f : X → R ,

$$\Gamma ^ { \ e w }, \quad \lim _ { t \to 0 } \frac { 1 } { t ( 4 \pi t ) ^ { q / 2 } } \int _ { \mathfrak { X } } \exp \left ( - \frac { | \mathfrak { x } - \mathfrak { y } | _ { 2, Q } ^ { 2 } } { t } \right ) ( f ( \mathfrak { y } ) - f ( \mathfrak { x } ) ) d \mu ^ { * } ( \mathfrak { y } ) = \Delta ( f ) ( \mathfrak { x } )$$

uniformly for x ∈ X , where | · | 2 ,Q denotes the glyph[lscript] 2 norm on R Q . Equivalently, uniformly for x ∈ X , we have

$$\left | \frac { 1 } { ( 4 \pi t ) ^ { q / 2 } } \int _ { \mathbb { X } } \exp \left ( - \frac { | \mathbf x - \mathbf y | _ { 2, Q } ^ { 2 } } { t } \right ) ( f ( \mathbf y ) - f ( \mathbf x ) ) d \mu ^ { * } ( \mathbf y ) - t \Delta ( f ) ( \mathbf x ) \right | = o ( t )$$

as t → 0+ .

From an approximation theory point of view, the theorem is more of a saturation theorem for approximating f on X , analogous to the Voronowskaja estimates for Bernstein polynomials ([25, Section 1.6.1], See Appendix A). Thus, (2.2) states that the rate of approximation of f cannot be better than O ( ), even if t f is infinitely differentiable, unless f is in the null space of the Laplace-Beltrami operator. This is to be expected because the Gaussian kernel involved is a positive operator. In particular, this phenomenon holds even if X is a Euclidean space rather than a manifold. Moreover, the curvature of the manifold contributes to the saturation as well. The Gaussian kernel has many advantages, invariance under translations and rotations is one of the them. This plays a major role in the proof of Theorem 2.1. Nevertheless, it is natural to ask whether another kernel can be found that leads directly to the approximation of the target function f on the manifold from the data without knowing the manifold itself and without having to go through an expensive eigen-decomposition. The curvature of the manifold will still affect the rate of convergence, but when applied to an affine space rather than a manifold, such a construction should lead to approximation without any saturation, without knowing what the affine space is (Remark 3.3).

The main objective of this paper is to demonstrate such a construction using certain localized kernels based on Hermite polynomials (Theorem 3.1). This theorem gives an analogue of Theorem 2.1 to obtain function approximation on an unknown manifold based only on noise-corrupted samples on the manifold, and give estimates on the degree of approximation. In the case when the approximation is done on an affine space rather than a manifold, our construction is free of any saturation, and does not need to know what the affine space is (Theorem 7.1).

↦

To recapture the advantage of the Gaussian kernel, we will study approximation by Gaussian networks. A (shallow) Gaussian network with n neurons has the form x → ∑ n k =1 a k exp( - | λ x -y k | 2 2 ,Q ). A deep Gaussian

network is constructed following a DAG structure, where each node (referred to as 'channel' in the literature on deep learning) evaluates a Gaussian network. Using the close connection between Hermite polynomials and Gaussian networks (cf. [27, 29, 13]), we can translate the result about approximation on the manifold into a result on approximation by shallow Gaussian networks, where the input is assumed to lie on an unknown low dimensional manifold of the nominally high dimensional ambient space (Theorem 5.1). In turn, using a property called 'good propagation of errors' (Theorem 5.2), we will 'lift' this theorem to estimate the degree of approximation by deep Gaussian networks, where each channel evaluates a Gaussian network on a similarly manifold-based data (Theorem 5.3). The networks themselves are constructed from certain pre-fabricated networks in the ambient space to approximate the Hermite functions with a correspondingly high number of neurons. However, we will give an explicit formula for such networks (Proposition 6.6), so that there is no training required here . The amount of information used in the final synthesis of the network will depend only on the dimension of the manifold on which the input lives. We consider this to be a step in bringing approximation theory of deep networks closer to the practice, so that the results are proved in the setting of approximation on unknown manifolds analogous to diffusion geometry rather than on known domains.

The statement of the main results in this paper mentioned above require a good deal of background information on the theory of weighted polynomial approximation, which we defer to Section 6. We will state the main results about approximation on a manifold in Section 3, and illustrate them using a simple numerical example in Section 4. We explain our ideas about shallow and deep networks in Section 5. To develop the details required in the constructions and proofs, we start by summarizing the relevant facts from the theory of weighted polynomial approximation in Section 6. Of particular interest is the approximation of a weighted polynomial using pre-fabricated Gaussian networks whose weights and centers do not depend upon the polynomial, as described in Section 6.3. Our main theorem in the context of approximation on unknown affine spaces is stated and proved in Section 7. The proofs of the results in Section 3 and 5 are given in Sections 8 and 9 respectively.

## 3 Approximation on manifolds

In this section, we state our main results on approximation on manifolds. The details and motivations for these constructions will be clearer after reading Sections 6 and 7. The notation on the manifolds is described in Section 3.1, the results themselves are discussed in Section 3.2.

## 3.1 Definitions

Let Q ≥ q ≥ 1 be integers, X be a q dimensional, compact, connected, sub-manifold of R Q (without boundary), with geodesic distance ρ and volume measure µ ∗ , normalized so that µ ∗ ( X ) = 1. We will identify the tangent space at x ∈ X with an affine space T x ( X ) in R Q passing through x . For any x ∈ X , we need to consider in this section three kinds of balls.

$$B _ { Q } ( \mathbf x, r ) \coloneqq \{ y \in \mathbb { R } ^ { Q } \, \colon | \mathbf x - y | _ { 2, Q } \leq r \}, \ B _ { \mathbb { T } } ( \mathbf x, r ) \coloneqq \mathbb { T } _ { \mathbf x } ( \mathbb { X } ) \cap B _ { Q } ( \mathbf x, r ), \ \mathbb { B } ( \mathbf x, r ) \coloneqq \{ y \in \mathbb { X } \, \colon \rho ( \mathbf x, y ) \leq r \}. \quad ( 3. 1 )$$

With this convention, the exponential map E x at x ∈ X (based on the definition in [18, Proposition 2.9]) is a diffeomorphism of an open ball centered at x in T x ( X ) onto its image in X such that ρ ( x , E x ( u )) = | u - | x 2 ,Q . Since X is compact, there exists ι ∗ &gt; 0 such that for every x ∈ X , E x is defined on B T ( x , ι ∗ ), and ρ ( x , E x ( u )) = | u - | x 2 ,Q for all u ∈ B T ( x , ι ∗ ).

We now define the smoothness class W γ ( X ). If f, g : X → R , the function fg : X → R is defined as usual by ( fg )( x ) = f ( x ) g ( x ) for x ∈ X . The space C ( X ) is the space of all continuous real-valued functions on X , equipped with the supremum norm ‖ ◦ ‖ X . The space C ∞ ( X ) is the subspace of C ( X ) comprising all infinitely differentiable functions on X . Let f ∈ C ( X ), γ &gt; 0. We say that f ∈ W γ ( X ) if for every x ∈ X , and φ ∈ C ∞ ( X ), supported on B ( x , ι ∗ / 2), the function F x ,φ : T x ( X ) → R defined by F x ,φ ( u ) := f ( E x ( u )) φ ( E x ( u )) is in W γ ( T x ) in the sense described in Section 7 (See (6.44), (7.3), and (7.4)). We define

$$\| f \| _ { W _ { \gamma } ( \mathbb { X } ) } \coloneqq \sup _ { x \in \mathbb { X }, \| \phi \| _ { x \leq 1 } } \| F _ { \mathfrak { x }, \phi } \| _ { W _ { \gamma } ( \mathbb { T } _ { \mathfrak { x } } ( \mathbb { X } ) ) }.$$

If γ is an integer and f is γ times differentiable on X then f ∈ W γ ( X ). The space W γ ( X ) can contain functions which are not differentiable. For example, we say that f ∈ Lip( X ) if glyph[negationslash]

$$\| f \| _ { \text{Lip} ( \mathbb { X } ) } \coloneqq \sup _ { \mathbf x, \mathbf y \in \mathbb { X }, \mathbf x \neq \mathbf y } \frac { | f ( \mathbf x ) - f ( \mathbf y ) | } { \rho ( \mathbf x, \mathbf y ) } < \infty.$$

We have Lip( X ) ⊂ W 1 ( X ).

Next, we define the approximation operators. The orthonormalized Hermite polynomial h k of degree k is defined recursively by

$$h _ { k } ( x ) \ \coloneqq \ \sqrt { \frac { 2 } { k } } x h _ { k - 1 } ( x ) - \sqrt { \frac { k - 1 } { k } } h _ { k - 2 } ( x ), \quad k = 2, 3, \cdots, \\ h _ { 0 } ( x ) \coloneqq \pi ^ { - 1 / 4 }, \ h _ { 1 } ( t ) \coloneqq \sqrt { 2 } \pi ^ { - 1 / 4 } x.$$

We write

$$\psi _ { k } ( t ) \coloneqq h _ { k } ( t ) \exp ( - t ^ { 2 } / 2 ), \quad t \in \mathbb { R }, \ k \in \mathbb { Z } _ { + }.$$

The functions { ψ k } ∞ k =0 are an orthonormal set with respect to the Lebesgue measure (cf. (6.1)). In the sequel, we fix an infinitely differentiable function H : [0 , ∞ → ) [0 , 1], such that H t ( ) = 1 if 0 ≤ t ≤ 1 / 2, and H t ( ) = 0 if t ≥ 1. We define for x ∈ R , m ∈ Z + :

$$\mathcal { P } _ { m, q } ( x ) \coloneqq \begin{cases} \pi ^ { - 1 / 4 } ( - 1 ) ^ { m } \frac { \sqrt { ( 2 m )! } } { 2 ^ { m } m! } \psi _ { 2 m } ( x ), & \text{if $q=1$}, \\ \frac { 1 } { \pi ^ { ( 2 q - 1 ) / 4 } \Gamma ( ( q - 1 ) / 2 ) } \sum _ { \ell = 0 } ^ { m } ( - 1 ) ^ { \ell } \frac { \Gamma ( ( q - 1 ) / 2 + m - \ell ) } { ( m - \ell )! } \frac { \sqrt { ( 2 \ell )! } } { 2 ^ { \ell } \ell! } \psi _ { 2 \ell } ( x ), & \text{if $q\geq2$}, \end{cases}$$

and the kernel Φ ˜ n,q for x ∈ R , n ∈ Z + by

$$\widetilde { \Phi } _ { n, q } ( x ) \coloneqq \sum _ { m = 0 } ^ { \lfloor n ^ { 2 } / 2 \rfloor } H \left ( \frac { \sqrt { 2 m } } { n } \right ) \mathcal { P } _ { m, q } ( x ).$$

## Constant convention:

In the sequel, c, c 1 , · · · will denote generic positive constants depending upon the dimension and other fixed quantities in the discussion, such as the norm. Their values may be different at different occurrences, even within a single formula. The notation A ∼ B means c 1 A ≤ B ≤ c 2 B . glyph[squaresolid]

## 3.2 Approximation theorems

The traditional machine learning paradigm is to consider data of the form { ( y j , f ( y j ) + glyph[epsilon1] j ) } , where y j 's are drawn randomly with respect to µ ∗ and glyph[epsilon1] j 's are random, mean 0 samples from an unknown distribution. More generally, we assume here a noisy data of the form ( y , glyph[epsilon1] ), with a joint probability distribution τ and assume further that the marginal distribution of y with respect to τ has the form dν ∗ = f dµ 0 ∗ for some f 0 ∈ C ( X ). In place of f ( y ), we consider a noisy variant F ( y , glyph[epsilon1] ), and denote

$$f ( \mathbf y ) \coloneqq \mathbb { E } _ { \tau } ( \mathcal { F } ( \mathbf y, \epsilon ) | \mathbf y ).$$

Remark 3.1 In practice, the data may not lie on a manifold, but it is reasonable to assume that it lies on a tubular neighborhood of the manifold. Our notation accommodates this - if z is a point in a neighborhood of X , we may view it as a perturbation of a point y ∈ X , so that the noisy value of the target function is F ( y , glyph[epsilon1] ), where glyph[epsilon1] encapsulate the noise in both the y variable and the value of the target function. An example is given in Example 4.1. glyph[squaresolid]

Our approximation process is simple: given by

$$\widehat { F } _ { n, \alpha } ( Y ; \mathbf x ) \coloneqq \frac { n ^ { q ( 1 - \alpha ) } } { M } \sum _ { j = 1 } ^ { M } \mathcal { F } ( \mathbf y _ { j }, \epsilon _ { j } ) \tilde { \Phi } _ { n, q } ( n ^ { 1 - \alpha } | \mathbf x - \mathbf y _ { j } | _ { 2, Q } ), \quad \mathbf x \in \mathbb { R } ^ { Q },$$

where 0 &lt; α ≤ 1.

Our main theorem is the following.

Theorem 3.1 Let γ &gt; 0 , τ be a probability distribution on X × Ω for some sample space Ω such the marginal distribution of τ restricted to X is absolutely continuous with respect to µ ∗ with density f 0 ∈ W γ ( X ) . We assume that

$$\sup _ { \mathbf x \in \mathbb { X }, r > 0 } \frac { \mu ^ { * } ( \mathbb { B } ( \mathbf x, r ) ) } { r ^ { q } } \leq c.$$

Let F : X × Ω → R be a bounded function, f defined by (3.7) be in W γ ( X ) , the probability density f 0 ∈ W γ ( X ) . Let M ≥ 1 , Y = ( { y 1 , glyph[epsilon1] 1 ) , · · · , ( y M , glyph[epsilon1] M ) } be a set of random samples chosen i.i.d. from τ . If

$$0 < \alpha < \frac { 4 } { 2 + \gamma }, \quad \alpha \leq 1,$$

then for every n ≥ 1 , 0 &lt; δ &lt; 1 and M ≥ n q (2 -α )+2 αγ √ log( n/δ ) , we have with τ -probability ≥ 1 -δ :

$$\left \| \widehat { F } _ { n, \alpha } ( Y ; \circ ) - f f _ { 0 } \right \| _ { \mathbb { X } } \leq c _ { 1 } \frac { \sqrt { \| f _ { 0 } \| \mathbb { X } } \| \mathcal { F } \| _ { \mathbb { X } \times \Omega } + \| f f _ { 0 } \| _ { W _ { \gamma } ( \mathbb { X } ) } } { n ^ { \alpha \gamma } }.$$

We record two corollaries of Theorem 3.1 as separate theorems. The first is the approximation of f itself, assuming that f 0 ≡ 1.

Theorem 3.2 With the set-up as in Theorem 3.1, let f 0 ≡ 1 (i.e., the marginal distribution of y with respect to τ is µ ∗ ). Then we have with τ -probability ≥ 1 -δ :

$$\left \| \widehat { F } _ { n, \alpha } ( Y ; \circ ) - f \right \| _ { \mathbb { X } } \leq c _ { 1 } \frac { \| \mathcal { F } \| _ { \mathbb { X } \times \Omega } + \| f \| _ { W _ { \gamma } ( \mathbb { X } ) } } { n ^ { \alpha \gamma } }.$$

The second is a consequence analogous to Theorem 2.1.

Theorem 3.3 With the set-up as in Theorem 3.1, we have with τ -probability ≥ 1 -δ :

$$\left \| \frac { n ^ { q ( 1 - \alpha ) } } { M } \sum _ { j = 1 } ^ { M } \left ( \mathcal { F } ( y _ { j }, \epsilon _ { j } ) - f ( \circ ) \right ) \tilde { \Phi } _ { n, q } ( n ^ { 1 - \alpha } | \circ - \mathbf y _ { j } | _ { 2, Q } ) \right \| _ { \mathbf x } \leq c _ { 1 } \frac { \sqrt { \| f _ { 0 } \| _ { \mathbf x } \| \mathcal { F } \| _ { \mathbf x \times \Omega } + \| f f _ { 0 } \| _ { \mathbf w _ { \mathbf y } ( \mathbf x ) } } } { n ^ { \alpha \gamma } }.$$

Remark 3.2 To compare the estimate (3.13) with (2.2), which is applicable with γ = 2, we are tempted to take any α ∈ (0 1), set , t = n -2(1 -α ) , and obtain the upper bound t A with A = α/ (1 -α ). Clearly, this bound tends to 0 arbitrarily fast with t . However, the estimate (2.2) uses a fixed kernel, while the estimate (3.13) uses a kernel depending upon t . glyph[squaresolid]

Remark 3.3 Although the curvature of the manifold forces us to put limitations on the rate of convergence in (3.12), this is not a saturation phenomenon. Thus, it is not ruled out that the rate can be much better than that given in (3.12) for non-trivial functions. glyph[squaresolid]

Remark 3.4 If γ &lt; 2, we may choose α = 1 without knowing the value of γ . The formula (3.8) itself does not require any prior knowledge of the smoothness of f . glyph[squaresolid]

## 4 Numerical example

We illustrate the theory using the following simple example. We let X ⊂ R 3 to be the helix defined by

$$\mathbf x ( t ) \coloneqq ( \cos ( \pi t ), \sin ( \pi t ), \pi t ), \quad 0 \leq t \leq 2 \pi.$$

This does not satisfy the conditions of the theorems in Section 3, and we will see an 'end point effect' in the errors, but we find it easy to work with this example because of the ease in computing the various quantities like the volume measure (arc-length) : dµ ∗ = ( √ 8 π 2 ) -1 dt . The target function f is given by

$$f ( { \mathbf x } ( t ) ) \coloneqq \cos ( x _ { 1 } ( t ) - x _ { 2 } ( t ) + x _ { 3 } ( t ) / 2 ) = \cos ( \cos ( \pi t ) - \sin ( \pi t ) - \pi t / 2 ), \quad 0 \leq t \leq 2 \pi. \quad \ \ ( 4. 2 )$$

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00156v4-figure-86e2d7ae-d43c-4875-8bd2-60096ceef211.jpg)

Figure 1: In all figures, black continuous line is the approximation, red dashed line is the target function (4.2). Left: Reconstruction without noise using 256 random training points, 2048 equidistant test points, Middle: Estimate in one trial, 256 random training points (blue dots) according to (4.3), 2048 equidistant test points, Right: Average of the estimates in 100 trials, 256 random training points plus noise each, 2048 equidistant test points

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00156v4-figure-530c714e-347c-4237-92f5-1bd0b888d273.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00156v4-figure-c6f7cbf2-9981-4d1a-8b02-1b95cce2f4cd.jpg)

## Example 4.1 We consider data of the form

$$\mathcal { F } ( \mathbf y, \epsilon ) \coloneqq f ( \mathbf y + \epsilon ) \exp ( 1. 1 2 5 ),$$

where glyph[epsilon1] is a random normal variable with mean 0 and standard deviation 1 5. . The factor exp(1 125) ensures that . the expected value of F is f . This example illustrates a multiplicative noise as well as additive noise. We may also consider this to be an example where every point y on the helix is perturbed by a normal noise with mean 0 and standard deviation 1, although we cannot deal directly with the perturbed points in the calculation of ̂ F n,α . We took M = 256, n = 64, α = 1. The results are reported in Figure 1 on one trial, as well as the average of ̂ F n,α over 100 trials. glyph[squaresolid]

Example 4.2 We consider data of the form

$$\mathcal { F } ( \mathbf y, \epsilon ) \coloneqq f ( \mathbf y ) + \epsilon,$$

where glyph[epsilon1] is a random normal variable with mean 0 and standard deviation 0 3. We take . M = 1024, and M samples of y distributed uniformly according to µ ∗ . We take n = 64, α = 1, and compute the quantity ̂ F 64 1 , ( Y, x ) for x = x ( ), where t t ranges over 2048 equidistant points on [0 , 2 π ]. The results are shown in Figure 2. glyph[squaresolid]

Figure 2: With f and F as in (4.2) and (4.4) respectively, M = 1024, n = 64, α = 1. On the x axis are 2048 equidistant samples on [0 , 2 π ]. Left: The function f in red, the sampled values F in blue dots for one trial, Middle: The function f in red, the reconstruction ̂ F 64 1 , for one trial in black, Right: A cumulative histogram of errors over 50 trials, the point ( p, y ) signifies that the error is 0 3 . y at p % of the test data.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00156v4-figure-057f74d6-58f3-4657-bac9-4639aa95f226.jpg)

## 5 Gaussian networks

In this section, we describe the consequences of Theorem 3.1 for Gaussian networks. In the case of shallow networks, we can give an explicit construction and error bounds in Section 5.1. In the case of deep networks (Section 5.2), we give only an existence theorem, explaining when the theorem can be described more constructively.

## 5.1 Shallow networks

Since P m,q and hence Φ ˜ n,q are even polynomials of degree &lt; n 2 , ˜ Φ n,q ( n 1 -α | ◦ | 2 ,Q ) ∈ Π . We will see in Remark 6.2 Q n that Φ ˜ n,q ( n 1 -α | ◦ | 2 ,Q ) = Φ n,q,Q ( 0 , n 1 -α ( ◦ )) for a polynomial kernel Φ n,q,Q on R Q . We may then define a prefabricated Gaussian network using (6.41)

$$\mathbb { G } _ { n, q, Q } ^ { * } \coloneqq \mathfrak { G } _ { Q } ( \Phi _ { n, q, Q } ( 0, ( n ^ { 1 - \alpha } ( \circ ) _ { 2, Q } ) ).$$

Using Corollary 6.2, we then deduce easily the following theorem about Gaussian networks. We note again that there is no training involved here. Even though the number of non-linearities in the network in the following theorem is O ( Mn 2 Q ), this potentially large number of non-linearities is not as much of a problem as it would be if we were to use an optimization procedure to train the network.

Theorem 5.1 Let (3.9) be satisfied, γ &gt; 0 , τ be a probability distribution on X × Ω for some sample space Ω such the marginal distribution of τ restricted to X is ν ∗ with dν ∗ = f dµ 0 ∗ for some f 0 ∈ W γ ( X ) . Let F : X × Ω → R be a bounded function, and f defined by (3.7) be in W γ ( X ) . Let 0 &lt; δ &lt; 1 , α satisfy (3.10). Let M ≥ 1 , Y = ( { y 1 , glyph[epsilon1] 1 ) , · · · , ( y M , glyph[epsilon1] M ) } be a set of random samples chosen i.i.d. from τ . If

$$M \geq n ^ { q ( 2 - \alpha ) + 2 \alpha \gamma } \sqrt { \log ( n / \delta ) }$$

we have with τ -probability ≥ 1 -δ :

$$\left \| \frac { 1 } { M } \sum _ { j = 1 } ^ { M } \left ( \mathcal { F } ( \mathbf y _ { j }, \epsilon _ { j } ) - f ( \circ ) \right ) \mathbb { G } _ { n, q, Q } ^ { * } ( \circ - \mathbf y _ { j } ) \right \| _ { \mathbb { X } } \leq c _ { 1 } \frac { \sqrt { \| f _ { 0 } \| _ { \mathbb { X } } \| } \mathcal { F } \| _ { \mathbb { X } \times \Omega } + \| f \| _ { W _ { r } ( \mathbb { X } ) } } { n ^ { \circ \gamma } }.$$

In particular, let

$$\dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots,\dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots, \dots }$$

If f 0 ≡ 1 , we have with τ -probability ≥ 1 -δ :

$$\| \mathbb { G } _ { n, q, Q } ( Y ; \mathcal { F } ) - f \| _ { \mathbb { X } } \leq c _ { 1 } \frac { \| \mathcal { F } \| _ { \mathbb { X } \times \Omega } + \| f \| _ { W _ { \gamma } ( \mathbb { X } ) } } { n ^ { \alpha \gamma } }.$$

## 5.2 Deep networks

The following discussion about the terminology about the deep networks is based on (almost taken from) the discussion in [35, 36], and elaborates upon the same. In particular, Figure 3 is taken from the arxiv version of [35]. A commonly used definition of a deep network is the following. Let φ : R → R be an activation function; applied to a vector x = ( x , 1 · · · , x q ), φ ( x ) = ( φ x ( 1 ) , · · · , φ ( x q )). Let L ≥ 2 be an integer, for glyph[lscript] = 0 , · · · , L , let q glyph[lscript] ≥ 1 be an integer ( q 0 = q ), T glyph[lscript] : R q glyph[lscript] → R q glyph[lscript] +1 be an affine transform, where q L +1 = 1. A deep network with L -1 hidden layers is defined as the compositional function

↦

$$\mathbf x \mapsto T _ { L } ( \phi ( T _ { L - 1 } ( \phi ( T _ { L - 2 } \cdots \phi ( T _ { 0 } ( \mathbf x ) ) \cdots ). }$$

This definition has several shortcomings. First, it does not distinguish between a function and the network architecture. As demonstrated in [36], a function may have more than one compositional representation, so that the affine transforms and L are not determined uniquely by the function itself. Second, this notion does not capture the connection between the nature of the target function and its approximation. Third, the affine transforms T glyph[lscript] define a special directed acyclic graph (DAG). It is cumbersome to describe notions of weight sharing, convolutions, sparsity, skipping of layers, etc. in terms of these transforms. Therefore, we have proposed in [35] to separate the architecture from the function itself, and describe a deep network more generally as a directed acyclic graph (DAG) architecture.

Let G be a DAG, with the set of nodes V ∪ S , where S is the set of source nodes, and V that of non-source nodes. For each node v ∈ V ∪ S , we denote its in-degree by d v ( ). Associated with each v ∈ V ∪ S is a compact, connected, Riemmanian submanifold X v of R d v ( ) with dimension q v , metric ρ v and volume element µ ∗ v . We assume further that (3.9) is satisfied with q v in place of q . Each of the in-edges to each node in V ∪ S represents an input real variable. If v ∈ V , u ∈ V ∪ S , u is called the child of v if there is an edge from u to v . The notion of the level

of a node is defined as follows. The level of a source node is 0. The level of v ∈ V is the length of the longest path from the nodes in S to v .

Each node v is supposed to evaluate a function f v on its input variables, supplied via the in-edges for v . The value of this function is propagated along the out-edges of v . Each of the source nodes obtains an input from some smooth manifold as described in Section 3. Other nodes can also obtain such an input, but by introducing dummy nodes, it is convenient to assume that only the source nodes obtain an input from the manifold.

Intuitively, we wish to say that the DAG structure implies a compositional structure for the functions involved; for example, if u , 1 · · · , u d v ( ) are children of v , then the function evaluated at v is f v ( f u 1 , · · · , f u d v ( ) ). To make this meaningful, we have to assume some 'pooling' operation on the input variables to make sure that the output of the vector valued function ( f u 1 , · · · , f u d v ( ) ) belongs to X v . Thus, for example, if the domain of f v is the cube [ -1 1] , d v ( ) , some clipping operation is required; if the domain is the torus in d v ( ) dimensions then some standard substitutions need to be made (e.g., [36]). We do not know how to specify the pooling operation in the general case of an unknown manifold, but assume that this pooling operation π v : R d v ( ) → X v has the following property: For any two sets of functions { f v ∈ C ( X v ) } v ∈ V , { g v ∈ C ( X v ) } v ∈ V ,

$$\rho _ { v } \left ( \pi _ { u _ { 1 } } ( x _ { u _ { 1 } } ), \cdots, f _ { u _ { d ( v ) } } ( x _ { u _ { d ( v ) } } ) ), \pi _ { v } ( g _ { u _ { 1 } } ( x _ { u _ { 1 } } ), \cdots, g _ { u _ { d ( v ) } } ( x _ { u _ { d ( v ) } } ) ) \right ) \leq c ( v ) \sum _ { k = 1 } ^ { d ( v ) } \| f _ { u _ { k } } - g _ { u _ { k } } \| _ { x _ { u _ { k } } }, \quad v \in V. \ \ ( 5. 7 )$$

A G -function is defined to be a set of functions { f v } v ∈ V ∪ S such that each f v ∈ C ( X v ), and if v ∈ V , u , 1 · · · , u d v ( ) are children of v , then the function evaluated at v is f v ( π v ( f u 1 , · · · , f u d v ( ) )). The individual functions f v will be called constituent functions .

For example, the DAG G in Figure 3 ([35]) represents the compositional function

$$f ^ { * } ( x _ { 1 }, \cdots, x _ { 9 } ) \ = \ h _ { 1 9 } ( h _ { 1 7 } ( h _ { 1 3 } ( h _ { 1 0 } ( x _ { 1 }, x _ { 2 }, x _ { 3 }, h _ { 1 6 } ( h _ { 1 2 } ( x _ { 6 }, x _ { 7 }, x _ { 8 }, x _ { 9 } ) ) ), h _ { 1 1 } ( x _ { 4 }, x _ { 5 } ) ), \\ h _ { 1 4 } ( h _ { 1 0 }, h _ { 1 1 } ), h _ { 1 6 } ), h _ { 1 8 } ( h _ { 1 5 } ( h _ { 1 1 }, h _ { 1 2 } ), h _ { 1 6 } ) ).$$

The G -function is { h 10 , · · · , h 19 = f ∗ } .

Figure 3: An example of a G -function ( f ∗ given in (5.8)). The vertices of the DAG G are the channels of the network. The input to the various channels is indicated by the in-edges of the nodes, and the output of the sink node h 19 indicates the output value of the G -function, f ∗ in this example.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00156v4-figure-211c6ca1-aa7f-4af5-a1df-8ba3c51b4548.jpg)

We assume that there is only one sink node, v ∗ (or v ∗ ( G )) whose output is denoted by f v ∗ (the target function ). Technically, there are two functions involved here: one is the final output as a function of all the inputs to all source nodes, the other is the final output as a function of the inputs to the node v ∗ . We will use the symbol f v ∗ to denote

both with comments on which meaning is intended when we feel that it may not be clear from the context. A similar convention is followed with respect to each of the constituent functions as well. For example, in the DAG of Figure 3, the function h 15 can be thought of both as a function of two variables, namely the outputs of h 11 and h 12 as well as a function of six variables x , 4 · · · , x 9 . In particular, if each constituent function is a neural network, h 15 is a shallow network receiving two inputs.

We define the notion of the variables 'seen' by a node. If u ∈ S , then these are the variables input to u . Let v ∈ V , and u , 1 · · · , u d v ( ) be the children of v . If x 1 , · · · , x d v ( ) are the inputs seen by u , 1 · · · , u d v ( ) , then the inputs seen by v are ( x 1 , · · · , x d v ( ) ), where the order is respected. For example, consider the function

$$f ^ { * } ( x _ { 1 }, x _ { 2 }, x _ { 3 }, x _ { 4 } ) = f ( f _ { 1 } ( x _ { 1 }, x _ { 2 } ), f _ { 2 } ( x _ { 4 }, x _ { 2 } ), f _ { 3 } ( x _ { 3 }, x _ { 1 } ) ).$$

The inputs seen by the leaves f 1 , f 2 , f 3 are ( x , x 1 2 ), ( x , x 4 2 ), ( x , x 3 1 ) respectively (not ( x , x 1 2 ), ( x , x 2 4 ), ( x , x 1 3 )). The inputs seen by f ∗ are ( x , x 1 2 , x 3 , x 4 ).

The following theorem enables us to 'lift' a theorem about shallow networks to that about deep networks.

Theorem 5.2 Let G be a DAG as described above, { f v } v ∈ V ∪ S , { g v } v ∈ V ∪ S be G -functions, and

$$\| f _ { v } - g _ { v } \| _ { \mathbf x _ { v } } \leq \varepsilon, \quad v \in V \cup \mathbf S.$$

Further assume that for each v ∈ V , f v ∈ Lip( X v ) , with L = max v ∈ V ‖ f v ‖ Lip ( X v ) . Then for the target function, thought of as a compositional function of all the input variables x to all the nodes in S , we have

$$| f _ { v ^ { * } } ( { \mathbf x } ) - g _ { v ^ { * } } ( { \mathbf x } ) | \leq c ( L, \mathcal { G } ) \varepsilon.$$

Theorem 5.2 allows us to lift Theorem 5.1 to deep networks. In general, we do not know the constituent functions. Also, for any given function and a DAG structure, it may not be possible to devise an algorithm to find the constituent functions uniquely. For example, (cos 2 x ) 2 and (1 / 4)(1 + cos(2 x )) 2 both have the structures g 1 ( g 2 ( x )) or f 1 ( f 2 ( x )), both representing the same DAG but with different constituent functions. Thus, even if we may assume that the noise occurs only in the approximation of the target function at the sink node and not in the constituent functions, it seems to be an extremely difficult problem to determine theoretically for any target function what the optimal DAG structure and the input/output for the constituent functions ought to be. Therefore, we have to state our theorem for deep networks only as an existence theorem, in the non-noisy case, not to complicate the notations too much. We assume also that at each node v , the input data is distributed according to the volume measure of X v .

Theorem 5.3 Let G be a DAG as described above, { f v } be a G -function, and we assume that each of the constituent functions f v ∈ W γ ( X v ) ∩ Lip( X v ) for some γ &gt; 0 , α satisfy (3.10) . Let n ≥ 1 . Then there exists a G -function { g v } such that each g v is a Gaussian network constructed using O ( n q v (2 -α )+2 αγ log n ) samples of its inputs, such that for any x seen by v ∗ ,

$$| f _ { v ^ { * } } ( \mathbf x ) - g _ { v ^ { * } } ( \mathbf x ) | \leq c ( L, \mathcal { G } ) n ^ { - \alpha \gamma }.$$

## 6 Background on weighted polynomials

## 6.1 Weighted polynomials

A good preliminary source of many identities regarding Hermite polynomials is the book [39] of Szeg¨ or the o Bateman manuscript [3].

We denote the class of all univariate algebraic polynomials of degree &lt; n by P n . The orthonormalized Hermite polynomial h k of degree k is defined recursively by (3.3). With ψ k ( x ) = h k ( x ) exp( -x / 2 2), one has the orthogonality relation for k, j ∈ Z + , glyph[negationslash]

$$\int _ { \mathbb { R } } \psi _ { k } ( x ) \psi _ { j } ( x ) d x = \begin{cases} \ 1, & \text{if $k=j$,} \\ 0, & \text{if $k\neq j$.} \end{cases}$$

Using (3.3), it is easy to deduce by induction that

$$\psi _ { \ell } ( 0 ) = \begin{cases} \ \pi ^ { - 1 / 4 } ( - 1 ) ^ { \ell / 2 } \frac { \sqrt { \ell }! } { 2 ^ { \ell / 2 } ( \ell / 2 )! }, & \text{if $\ell$ is even,} \\ 0, & \text{if $\ell$ is odd,} \end{cases}$$

The Hermite polynomial h m has m real and simple zeros x k,m . Writing

$$\lambda _ { k, m } \coloneqq \left ( \sum _ { j = 0 } ^ { m - 1 } h _ { j } ( x _ { k, m } ) ^ { 2 } \right ) ^ { - 1 },$$

it is well known (cf. [39, Section 3.4]) that

$$\sum _ { k = 1 } ^ { m } \lambda _ { k, m } P ( x _ { k, m } ) = \int _ { \mathbb { R } } P ( x ) \exp ( - x ^ { 2 } ) d x, \quad P \in \mathbb { P } _ { 2 m }.$$

It is also known (cf. [27, Theorem 8.2.7], applied with p = 2, b = 0) that

$$\sum _ { k = 1 } ^ { m } \lambda _ { k, m } \exp ( x _ { k, m } ^ { 2 } ) \leq c m ^ { 1 / 2 }.$$

The Mehler formula [1, Formula (6.1.13)] states that

$$& \sum _ { j = 0 } ^ { \infty } \psi _ { j } ( y ) \psi _ { j } ( z ) w ^ { j } = \frac { 1 } { \sqrt { \pi ( 1 - w ^ { 2 } ) } } \exp \left ( \frac { 2 y z w - ( y ^ { 2 } + z ^ { 2 } ) w ^ { 2 } } { 1 - w ^ { 2 } } \right ) \exp ( - ( y ^ { 2 } + z ^ { 2 } ) / 2 ), \quad y, z \in \mathbb { R }, \, w \in \mathbb { C }, \, | w | < 1. \\ & \quad \text{Nav+} \quad \text{no individual and return the non-interactive of Normalize in the minimal indicator of the} \, \text{$\mathbb{A}$ will} \, \text{no}.$$

Next, we introduce and review the properties of Hermite polynomials in the multivariate setting. We will need to use spaces with many different dimensions. Therefore, in this section, we will use the symbol d to denote a generic dimension, which will be replaced later by q , Q q , v , etc.

If d ≥ 2 is an integer, we define Hermite polynomials on R d using tensor products. We adopt the notation x = ( x , 1 · · · , x d ). The orthonormalized Hermite function is defined by

$$\psi _ { \mathbf k } ( \mathbf x ) \coloneqq \prod _ { j = 1 } ^ { d } \psi _ { k _ { j } } ( x _ { j } ).$$

In general, when univariate notation is used in multivariate context, it is to be understood in the tensor product sense as above; e.g., k ! = ∏ d j =1 ( k j !), x k = ∏ d j =1 x k j j , etc. The notation | · | p,d will denote the glyph[lscript] p norm on R d .

For any set A ⊂ R d and f : A → R , we denote by C A ( ) the space of all uniformly continuous and bounded functions on A , with the norm ‖ f ‖ A = sup x ∈ A | f ( x ) . | The space C 0 ( A ) is the subspace of all f ∈ C A ( ) vanishing at infinity.

We will often use (without mentioning it explicitly) the fact deduced from the univariate bounds proved in [2] that

$$| \psi _ { \mathbf k } ( \mathbf x ) | \leq c.$$

We will denote by Π d n the span of { ψ k : √ | k | 1 ,d &lt; n } and by P d n the space of all algebraic polynomials of total degree &lt; n . Thus, if P ∈ Π , then d n P ( x ) = R ( x ) exp( -| x | 2 2 ,d / 2) for some R ∈ P d n 2 . The following proposition lists a few important properties of these spaces (cf. [27, 30, 33]).

## Proposition 6.1 Let n &gt; 0 , P ∈ Π d n .

- (a) (Infinite-finite range inequality) For any δ &gt; 0 , there exists c = ( ) c δ such that

$$\| P \| _ { \mathbb { R } ^ { d } \vee [ - \sqrt { 2 } n ( 1 + \delta ), \sqrt { 2 } n ( 1 + \delta ) ] ^ { d } } \leq c _ { 1 } e ^ { - c n ^ { 2 } } \| P \| _ { [ - \sqrt { 2 } n ( 1 + \delta ), \sqrt { 2 } n ( 1 + \delta ) ] ^ { d } }$$

- (b) (MRS identity) We have

$$\| P \| _ { \mathbb { R } ^ { d } } = \| P \| _ { [ - \sqrt { 2 } n, \sqrt { 2 } n ] ^ { d } }.$$

- (c) (Bernstein inequality) There is a positive constant B depending only on d such that

$$\| | \nabla P | _ { 2, d } \| _ { \mathbb { R } ^ { d } } \leq B n \| P \| _ { \mathbb { R } ^ { d } }.$$

Let m ≥ 1. For a multi-integer j , 1 ≤ j ≤ m , we write x j ,m,d := ( x j 1 ,m , · · · , x j d ,m ), and λ j ,m,d := ∏ d glyph[lscript] =1 λ j glyph[lscript] ,m . We observe further that if P , P 1 2 ∈ Π d m , then P 1 ( x ) P 2 ( x ) = R ( x ) exp( -| x | 2 ) for some R ∈ P d 2 m 2 . Therefore, (6.4) and (6.5) lead to the following fact, which we formulate as a proposition.

Proposition 6.2 For m ≥ 1 , we have

$$\sum _ { 1 \leq j \leq m ^ { 2 } } \lambda _ { \mathbb { j }, m ^ { 2 }, d } \exp ( | \mathbf x _ { \mathbb { j }, m ^ { 2 }, d } | _ { 2, d } ^ { 2 } ) P _ { 1 } ( \mathbf x _ { \mathbb { j }, m ^ { 2 }, d } ) P _ { 2 } ( \mathbf x _ { \mathbb { j }, m ^ { 2 }, d } ) = \int _ { \mathbb { R } ^ { d } } P _ { 1 } ( \mathbf x ) P _ { 2 } ( \mathbf x ) d \mathbf x, \quad P _ { 1 }, P _ { 2 } \in \Pi _ { m } ^ { d },$$

and

$$\sum _ { 1 \leq j \leq m ^ { 2 } } \lambda _ { j, m ^ { 2 }, d } \exp ( | \mathbf x _ { j, m ^ { 2 }, d } | _ { 2, d } ^ { 2 } ) \leq c m ^ { d / 2 }.$$

## 6.2 Applications of Mehler identity

The Mehler identity for multivariate Hermite polynomials is expressed conveniently by writing

$$\text{Proj} _ { m, d } ( x, y ) \coloneqq \sum _ { | k | _ { 1, d } = m } \psi _ { k } ( x ) \psi _ { k } ( y ).$$

Using the univariate Mehler identity (6.6), it is then easy to deduce that for w ∈ C , | w &lt; | 1,

$$\text{Using the univariate Mehler identity } ( 6. 6 ), \text{it is then easy to deduce that for } w \in \mathbb { C }, | w | < 1, \\ \sum _ { \kappa \in \mathbb { Z } ^ { 4 } } \psi _ { k } ( x ) \psi _ { k } ( y ) w ^ { | k | _ { 1 }, d } = \sum _ { m = 0 } ^ { \infty } w ^ { m } \text{Proj} _ { m, d } ( x, y ) \\ = \frac { 1 } { ( \pi ( 1 - w ^ { 2 } ) ) ^ { d / 2 } } \exp \left ( \frac { 4 w x \cdot y - ( 1 + w ^ { 2 } ) ( | x | _ { 2, d } ^ { 2 } + | y | _ { 2, d } ^ { 2 } ) } { 2 ( 1 - w ^ { 2 } ) } \right ) \\ = \frac { 1 } { ( \pi ( 1 - w ^ { 2 } ) ) ^ { d / 2 } } \exp \left ( \frac { - 1 + w | x - y | _ { 2, d } ^ { 2 } } { 1 - w } \frac { 1 - w | x + y | _ { 2, d } ^ { 2 } } { 1 + w } \right ) \\ = \frac { 1 } { ( \pi ( 1 - w ^ { 2 } ) ) ^ { d / 2 } } \exp \left ( - \frac { 1 + w ^ { 2 } } { 2 ( 1 - w ^ { 2 } ) } \left | x - \frac { 2 w } { 1 + w ^ { 2 } } y ^ { | } _ { 2, d } \right ) \exp \left ( - \frac { 1 - w ^ { 2 } } { 2 ( 1 + w ^ { 2 } ) } | y | _ { 2, d } ^ { 2 } \right ). \\ \text{We note an identity } ( 6. 8 ) \text{ which follows immediately from } ( 6. 1 5 ) \text{ by setting } x = y = 0. \text{ For integer } d \text{ (not}$$

We note an identity (6.18) which follows immediately from (6.15) by setting x = y = 0 . For integer d (not necessarily positive), we define the sequence D d r ; by

$$\cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cd. \quad D _ { d ; r } \colon = \begin{cases} \begin{array} { c } \pi ^ { - d / 2 } ( - 1 ) ^ { r / 2 } \frac { \Gamma ( 1 - d / 2 ) } { \Gamma ( 1 - d / 2 - r / 2 ) ( r / 2 )! } \end{array} & \text{if $r$ is even, $d\leq 0$,} \\ \begin{array} { c } \pi ^ { - d / 2 } \frac { \Gamma ( d / 2 + r / 2 ) } { \Gamma ( d / 2 ) ( r / 2 )! }, & \text{if $r$ is even, $d\geq 1$,} \\ 0, & \text{if $r$ is odd.} \end{cases} &$$

This sequence is chosen so as to satisfy

$$\pi ^ { - d / 2 } ( 1 - w ^ { 2 } ) ^ { - d / 2 } = \sum _ { r = 0 } ^ { \infty } D _ { d ; r } w ^ { r }, \quad | w | < 1.$$

Using the Mehler identity (6.15), we deduce that for any integer d ≥ 1

$$\sum _ { r = 0 } ^ { \infty } w ^ { 2 r } \sum _ { | k | _ { 1, d } = 2 r } | \psi _ { \mathbf k } ( 0 ) | ^ { 2 } = ( \pi ( 1 - w ^ { 2 } ) ) ^ { - d / 2 } = \pi ^ { - d / 2 } \sum _ { r = 0 } ^ { \infty } \frac { \Gamma ( d / 2 + r ) } { \Gamma ( d / 2 ) r! } w ^ { 2 r } = \sum _ { \ell = 0 } ^ { \infty } D _ { d ; \ell } w ^ { \ell }.$$

In this section, we point out the invariance and localization properties of certain kernels using the Mehler identity.

## 6.2.1 Rotation invariance

An interesting consequence of the Mehler identity is that the projection Proj m,d is invariant under rotations. For d ≥ 2 and any x y , ∈ R d , we may therefore use an appropriate rotation to write

$$\text{Proj} _ { m, d } ( x, y ) = \sum _ { j = 0 } ^ { m } \text{Proj} _ { j, 2 } ( ( | x | _ { 2, d }, 0 ), ( | y | _ { 2, d } \cos \theta, | y | _ { 2, d } \sin \theta ) ) \sum _ { | k | _ { 1, d - 2 } \leq m - j } | \psi _ { k } ( 0 ) | ^ { 2 },$$

where cos θ = x · y / ( | x y || | ), with obvious modifications if y = 0 . Hence, we obtain from (6.19) and (6.18) (used with d -2 in place of d ),

$$\text{Proj} _ { m, d } ( x, y ) = \sum _ { j = 0 } ^ { m } \text{Proj} _ { j, 2 } ( ( | x | _ { 2, d }, 0 ), ( | y | _ { 2, d } \cos \theta, | y | _ { 2, d } \sin \theta ) ) D _ { d - 2 ; m - j }.$$

In the case when d = 1, (6.19) takes the form

$$\text{proj} _ { m, 1 } ( x, y ) = \psi _ { m } ( | x | ) \psi _ { m } ( | y | \cos \theta ),$$

where cos θ = xy/ ( | x y || | ) = sgn( xy ), (sgn(0) = 0).

Let Q ≥ q ≥ 1 be integers. We can extend the definition of Proj m,q to x y , ∈ R Q by

$$\text{Proj} _ { m, q, Q } ( x, y ) \coloneqq \underset { j = 0 } { \overset { m } { \sum } } \text{Proj} _ { j, 2 } ( ( | x | _ { 2, Q }, 0 ), ( | y | _ { 2, Q } \cos \theta, | y | _ { 2, Q } \sin \theta ) ) D _ { q - 2 ; m - j }, \quad \text{if $q \geq 2$}, \\ \psi _ { m } ( | x | _ { 2, Q } ) \psi _ { m } ( | y | _ { 2, Q } \cos \theta ), \quad \text{if $q = 1$}.$$

The relationship between Proj m,q,Q and Proj m,Q , both defined on R Q is given by the following proposition.

Proposition 6.3 Let Q &gt; q ≥ 2 be integers. Let m ≥ 0 , and x y , ∈ R Q .

- (a) We have

$$\Pr o j _ { m, Q } ( { \mathbf x }, { \mathbf y } ) = \pi ^ { ( q - Q ) / 2 } \sum _ { \ell = 0 } ^ { \lfloor m / 2 \rfloor } \left ( \begin{smallmatrix} ( Q - q ) / 2 + \ell - 1 \\ \ell \end{smallmatrix} \right ) \Pr o j _ { m - 2 \ell, q, Q } ( { \mathbf x }, { \mathbf y } ).$$

- (b) We have

$$\Pr o j _ { m, q, Q } ( x, y ) = \pi ^ { ( Q - q ) / 2 } \sum _ { \ell = 0 } ^ { \lfloor m / 2 \rfloor } ( - 1 ) ^ { \ell } \Big { ( } \begin{pmatrix} ( Q - q ) / 2 \\ \ell \end{pmatrix} \Pr o j _ { m - 2 \ell, Q } ( x, y ).$$

Hence, Proj m,q,Q ( x y , ) is a weighted polynomial in Π Q m as a function of x and y .

(c) If x is a scalar multiple of y , then (6.23) and (6.24) both hold also when q = 1 .

Proof. In this proof, let x ′ = ( | x | 2 ,Q , 0 , 0 , · · · , 0 ), ︸ ︷︷ ︸ q -2 times y ′ = ( | y | 2 ,Q cos θ, | y | 2 ,Q sin θ, 0 , · · · , 0 ). ︸ ︷︷ ︸ q -2 times In view of (6.19), we observe that

$$\Pr o j _ { m, q, Q } ( x, y ) = \Pr o j _ { m, q } ( x ^ { \prime }, y ^ { \prime } ).$$

Further, | x -y | 2 ,Q = | x ′ -y ′ | 2 ,q , | x + y | 2 ,Q = | x ′ + y ′ | 2 ,q . Therefore, the Mehler identity (6.15) shows that

$$\Gamma \widehat { u } \widehat { u }, | A - y | 2, Q - | A - y | 2, q, | A + y | 2, Q - | A + y | 2, q } \cdots \Gamma \widehat { u } \widehat { u } \widehat { u } \\ \sum _ { m = 0 } ^ { \infty } w ^ { m } \Pr o j _ { m, Q } ( x, y ) & = \frac { 1 } { ( \pi ( 1 - w ^ { 2 } ) ) ^ { Q / 2 } } \exp \left ( - \frac { 1 + w } { 1 - w } \frac { | x - y | ^ { 2 } _ { 2, Q } } { 4 } - \frac { 1 - w } { 1 + w } \frac { | x + y | ^ { 2 } _ { 2, Q } } { 4 } \right ) \\ & = \frac { 1 } { ( \pi ( 1 - w ^ { 2 } ) ) ( Q - q ) / 2 } \sum _ { m = 0 } ^ { \infty } w ^ { m } \Pr o j _ { m, q } ( x ^ { \prime }, y ^ { \prime } ) \\ & = \frac { 1 } { ( \pi ( 1 - w ^ { 2 } ) ) ( Q - q ) / 2 } \sum _ { m = 0 } ^ { \infty } w ^ { m } \Pr o j _ { m, q, Q } ( x, y ). \\ \intertext { W a n o w o o l l t h o M a C l a n i n o w o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i m o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i N o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w O i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n O w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n u n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o i n o w o ier r e e,$$

We now recall the McClaurin expansion for (1 -w 2 ) -( Q q / -) 2 (cf. (6.17)), multiply the two power series using the Cauchy-Leibnitz formula, and compare the coefficients to arrive at (6.23). Part (b) is proved similarly by observing that

$$\sum _ { m = 0 } ^ { \infty } w ^ { m } \Pr o j _ { m, q, Q } ( { \mathbf x }, { \mathbf y } ) = \pi ^ { ( Q - q ) / 2 } ( 1 - w ^ { 2 } ) ^ { ( Q - q ) / 2 } \sum _ { m = 0 } ^ { \infty } w ^ { m } \Pr o j _ { m, Q } ( { \mathbf x }, { \mathbf y } ).$$

If x is a scalar multiple of y , then sin θ = 0, so that x ′ = ( | x | 2 ,Q , 0 , · · · , 0 ), ︸ ︷︷ ︸ q -1 times y ′ = ( | y | 2 ,Q cos θ, 0 , · · · , 0 ). ︸ ︷︷ ︸ q -1 times Part

- (c) is then proved using the same calculations as above. glyph[squaresolid]

↦

Remark 6.1 Clearly, for every x y , ∈ R Q , Proj m,Q,Q ( x y , ) = Proj m,Q ( x y , ), Proj m,q,Q ( x y , ) = Proj m,q,Q ( -x , -y ) and the kernel ( x y , ) → Proj m,q,Q ( 0 x , -y ) is both rotation invariant and translation invariant. Using (6.19), (6.2), and (6.18) (used with d = q -1) show that for all q ≥ 1 and m ∈ Z + , Proj 2 m -1 ,q,Q ( 0 x , ) = 0, and

$$\Pr j _ { 2 m, q, Q } ( 0, x ) = \Pr j _ { 2 m, q } ( 0, x ^ { \prime } ) = \sum _ { j = 0 } ^ { 2 m } \psi _ { j } ( | x | _ { 2, Q } ) \psi _ { j } ( 0 ) \sum _ { | k | _ { 1, q - 1 } \leq 2 m - j } | \psi _ { k } ( 0 ) | ^ { 2 } = \mathcal { P } _ { m } ( | x | _ { 2, Q } ).$$

glyph[squaresolid]

## 6.2.2 Localized kernels

In this section, we recall the localization properties of certain kernels. In the sequel, H : [0 , ∞ → ) [0 , 1] is a fixed, infinitely differentiable function, with H t ( ) = 1 if 0 ≤ t ≤ 1 / 2, H t ( ) = 0 if t ≥ 1. All constants may depend upon H as well. We define

$$\Phi _ { n, d } ( H ; \mathbf x, \mathbf y ) \coloneqq \Phi _ { n, d } ( \mathbf x, \mathbf y ) \coloneqq \sum _ { \mathbf k \in \mathbb { Z } _ { + } ^ { d } } H \left ( \frac { \sqrt { | \mathbf k | _ { 1, d } } } { n } \right ) \psi _ { \mathbf k } ( \mathbf x ) \psi _ { \mathbf k } ( \mathbf y ) = \sum _ { m = 0 } ^ { n ^ { 2 } } H \left ( \frac { \sqrt { m } } { n } \right ) \Pr o j _ { m, d } ( \mathbf x, \mathbf y ), \quad \mathbf x, \mathbf y \in \mathbb { R } ^ { d }.$$

(6.29)

Using Mehler identity and the Tauberian theorem in [34, Theorem 4.3], we proved in [12, Lemma 4.1] the following proposition.

Proposition 6.4 For n ≥ 1 , x y , ∈ R d , we have

$$| \Phi _ { n, d } ( { \mathbf x }, { \mathbf y } ) | \leq \frac { c n ^ { d } } { \max ( 1, ( n | { \mathbf x } - { \mathbf y } | _ { 2, d } ) ^ { S } } ) }.$$

In particular, and for 1 ≤ p &lt; ∞ ,

$$\sup _ { x \in \mathbb { R } ^ { d } } \int _ { \mathbb { R } ^ { d } } | \Phi _ { n, d } ( x, y ) | ^ { p } d y \leq c n ^ { d ( p - 1 ) }.$$

We extend the definition of Φ n,d as follows. Let Q ≥ q ≥ 1 be integers. We define

$$\Phi _ { n, q, Q } ( \mathbf x, \mathbf y ) \coloneqq & \sum _ { m = 0 } ^ { \infty } H \left ( \frac { \sqrt { m } } { n } \right ) \Pr o j _ { m, q, Q } ( \mathbf x, \mathbf y ), \quad \mathbf x, \mathbf y \in \mathbb { R } ^ { Q }.$$

Remark 6.2 In view of Remark 6.1, the kernel ˜ Φ n,q defined in (3.6) satisfies ˜ Φ n,q ( | x | 2 ,Q ) = Φ n,q,Q ( 0 x , ). In particular, Φ ˜ n,q ( | ◦ | 2 ,Q ) ∈ Π . Q n glyph[squaresolid]

Proposition 6.5 Let S &gt; Q ≥ q ≥ 2 be integers. The kernel Φ n,q,Q ( x y , ) ∈ Π Q n as a function of x and y . For x y , ∈ R Q , n = 1 2 , , · · · ,

$$| \Phi _ { n, q, Q } ( \mathbf x, y ) | \leq \frac { c n ^ { q } } { \max ( 1, ( n | \mathbf x - y | _ { 2, Q } ) ^ { S } ) }.$$

In particular,

$$| \Phi _ { n, q, Q } ( \mathbf x, y ) | \leq c n ^ { q }.$$

If x is a scalar multiple of y , then

$$| \Phi _ { n, 1, Q } ( { \mathbf x }, { \mathbf y } ) | \leq \frac { c n } { \max ( 1, ( n | { \mathbf x } - { \mathbf y } | _ { 2, Q } ) ^ { S } } ) }, \quad | \Phi _ { n, 1, Q } ( { \mathbf x }, { \mathbf y } ) | \leq c n.$$

Proof. Let x ′ , y ′ be as in the proof of Proposition 6.3. Since Φ n,q,Q ( x y , ) = Φ n,q ( x ′ , y ′ ), this proposition follows directly from Proposition 6.4. glyph[squaresolid]

Corollary 6.1 The kernel ˜ Φ n,q defined in (3.6) satisfies each of the following properties.

$$| \widetilde { \Phi } _ { n, q } ( | \mathbf x | _ { 2, Q } ) ) | \leq \frac { c n ^ { q } } { \max ( 1, ( n | \mathbf x | _ { 2, Q } ) ^ { S } ) }, \quad \mathbf x \in \mathbb { R } ^ { q },$$

$$| \widetilde { \Phi } _ { n, q } ( | x | _ { 2, Q } ) ) | \leq c n ^ { q }, \ \ | \widetilde { \Phi } _ { n, q } ( | x | _ { 2, Q } ) ) - \widetilde { \Phi } _ { n, q } ( | y | _ { 2, Q } ) ) | \leq c n ^ { q + 1 } \, | | x | _ { 2, Q } - | y | _ { 2, Q } | \,, \quad x, y \in \mathbb { R } ^ { Q }.$$

$$| \Phi _ { n, d } ( \mathbf x, \mathbf y ) | \leq c n ^ { d },$$

Proof. The estimate (6.37) and the first estimate in (6.38) follows from Proposition 6.5 and the fact that ˜ Φ n,q ( | x | 2 ,Q ) = Φ n,q,Q ( 0 x , ). The second estimate in (6.38) follows from the Bernstein inequality (6.11) applied with d = 1 to the univariate polynomial Φ ˜ n,q . glyph[squaresolid]

## 6.3 From Hermite polynomials to Gaussian networks

We discuss in this section the close connection between Hermite polynomials and Gaussian networks.

Proposition 6.6 Let m ≥ 1 , k ∈ Z d + , and for | k | 1 ,d &lt; m 2 , x ∈ R d ,

$$\mathfrak { G } _ { k, m, d } ( x ) \coloneqq \left ( \frac { 3 } { 2 \pi } \right ) ^ { d / 2 } 3 ^ { | k | _ { 1, d } / 2 } \sum _ { 1 \leq 1 \leq 2 m ^ { 2 } } \lambda _ { 1, 2 m ^ { 2 }, d } \exp ( 3 | x _ { 1, 2 m ^ { 2 }, d | _ { 2, d } / 4 } ) \psi _ { k } ( x _ { 1, 2 m ^ { 2 }, d } ) \exp \left ( \left | x - \frac { \sqrt { 3 } } { 2 } x _ { 1, 2 m ^ { 2 }, d } \right | _ { 2, d } \right ).$$

Then

$$\max _ { | k | _ { 1, d } < m } \| \psi _ { k } - \mathfrak { G } _ { k, m, d } \| _ { \mathbb { R } ^ { d } } \leq c m ^ { d - 2 } 3 ^ { - m ^ { 2 } / 2 }.$$

Clearly, the number of neurons in the network G k ,m,d is O ( m 2 d ) .

Proof. This proof is the same as that in [13, Lemma 4.2] and [29, Lemma 4.1]. Using the last expression in (6.15) with w = 1 / √ 3, we obtain

$$\sum _ { k \in \mathbb { Z } _ { + } ^ { d } } \psi _ { k } ( x ) \psi _ { k } ( \mathbf u ) 3 ^ { - | k | _ { 1, d } / 2 } = \left ( \frac { 3 } { 2 \pi } \right ) ^ { d / 2 } \exp \left ( - \left | \mathbf x - \frac { \sqrt { 3 } } { 2 } \mathbf u \right | _ { 2, d } ^ { 2 } \right ) \exp ( - | \mathbf u | ^ { 2 } / 4 ).$$

In this proof, we denote by ν ∗ m,d the measure that associates the mass λ j , 2 m ,d 2 exp( | x j , 2 m ,d 2 | 2 2 ,d ) with the point x j , 2 m ,d 2 for 1 ≤ j 1 , · · · , j d ≤ 2 m 2 . Therefore, using Proposition 6.2 with m √ 2 in place of m , we obtain

$$\ddot { \psi } _ { 1 }, e m ^ { *, a } \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \colon \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \ddots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots.$$

The first term on the right hand side above is G k ,m,d . The second term is estimated using (6.8) and (6.13) (applied with m √ 2 in place of m ) exactly as in the proof of [13, Lemma 4.2]. We omit the details. glyph[squaresolid]

The following corollary is easy to deduce (cf. [29, Proposition 4.1]). If P = ∑ | k | 1 ,d &lt;m 2 b k ψ k ∈ Π d m , we define

$$\mathfrak { G } _ { d } ( P ) \coloneqq \sum _ { | k | _ { 1, d } < m ^ { 2 } } b _ { k } \mathfrak { G } _ { k, m, d }.$$

Corollary 6.2 Let m ≥ 1 , P ∈ Π d m . Then

$$\| P - \mathfrak { G } _ { d } ( P ) \| _ { \mathbb { R } ^ { d } } & \leq c _ { 1 } m ^ { c _ { 3 } - m ^ { 2 } / 2 } \| P \| _ { \mathbb { R } ^ { d } }.$$

We note that the centers and the number of neurons in the network ∑ | k | 1 ,d &lt;m 2 b k G k ,m,d are independent of P . In particular, the number of neurons is O ( m 2 d ) .

## 6.4 Function approximation

In this section, we describe some results on approximation of functions on R d . If f ∈ C 0 ( R d ), we define its degree of approximation by

$$E _ { n } ( \mathbb { R } ^ { d } ; f ) \coloneqq \min _ { P \in \Pi _ { n } ^ { d } } \| f - P \| _ { \mathbb { R } ^ { d } }.$$

For γ &gt; 0, the smoothness class W γ ( R d ) comprises f ∈ C 0 ( R d ) for which

$$\| f \| _ { W _ { \gamma } ( \mathbb { R } ^ { d } ) } \coloneqq \| f \| _ { \mathbb { R } ^ { d } } + \sup _ { n \geq 0 } 2 ^ { n \gamma } E _ { 2 ^ { n } } ( \mathbb { R } ^ { d } ; f ) < \infty.$$

We need some results from [27, 28], reformulated in the form stated in Theorem 6.1 below. To state this theorem, we need some notation first. First, for δ ∈ (0 , 1], x ∈ R d , 1 ≤ k ≤ d , we write

$$Q _ { k, \delta } ^ { \prime } ( \mathbf x ) \coloneqq \min ( \delta ^ { - 1 }, | x _ { k } | ).$$

For t &gt; 0 and integer j ≥ 0, the forward difference of a function f : R d → R is defined by

$$\Delta _ { k, t } ^ { j } ( f ) ( { \mathbf x } ) \coloneqq \sum _ { \ell = 0 } ^ { j } ( - 1 ) ^ { j - \ell } \binom { j } { \ell } f ( x _ { 1 }, \cdots, x _ { k - 1 }, x _ { k } + \ell t, x _ { k + 1 }, \cdots, x _ { d } ).$$

and for integers r ≥ 1

$$\coloneqq \, ^ { * } \, \\ \omega _ { r } ( f, \delta ) \coloneqq \sum _ { k = 1 } ^ { d } \sum _ { j = 0 } ^ { r } \delta ^ { r - j } \sup _ { | t | \leq \delta } \| ( Q ^ { \prime } _ { k, \delta } ) ^ { r - j } \Delta ^ { j } _ { k, t } ( f ) \| _ { \mathbb { R } ^ { d } }.$$

↦

Remark 6.3 If λ &gt; 0, f λ ( x ) = f ( x /λ ), then ∆ j k,t ( f λ )( x ) = ∆ j k,t/λ ( f )( x /λ ), and Q ′ k,δ ( x ) = λQ ′ k,λδ ( x /λ ). Using the fact that δ → δQ k,δ ( x ) is non-decreasing for every x , it is not difficult to deduce that

$$\omega _ { r } ( f _ { \lambda }, \delta ) \quad & = \quad \sum _ { k = 1 } ^ { d } \sum _ { j = 0 } ^ { r } \delta ^ { r - j } \sup _ { | l | \leq \delta } \| ( Q ^ { \prime } _ { k, \delta } ) ^ { r - j } \Delta ^ { j } _ { k, t } ( f _ { \lambda } ) \| _ { \mathbb { R } ^ { d } } = \sum _ { k = 1 } ^ { d } \sum _ { j = 0 } ^ { r } ( \lambda \delta ) ^ { r - j } \sup _ { | \omega | \leq \delta / \lambda } \| ( Q ^ { \prime } _ { k, \lambda \delta } ) ^ { r - j } \Delta ^ { j } _ { k, t } ( f ) \| _ { \mathbb { R } ^ { d } } \\ & \leq \quad \omega _ { r } ( f, \max ( \lambda, 1 / \lambda ) \delta ).$$

Theorem 6.1 Let f ∈ C 0 ( R d ) , r ≥ 1 , 0 &lt; γ &lt; r . Then (a) For n ≥ 1 ,

$$E _ { n } ( \mathbb { R } ^ { d } ; f ) \leq c \omega _ { r } ( f, 1 / n ).$$

(b)

The function

f

∈

W

γ

(

R

d

)

if and only if

ω

r

(

f, δ

) =

O

(

δ

γ

)

for

0

&lt; δ

≤

1

.

In fact,

$$\| f \| _ { W _ { \gamma } ( \mathbb { R } ^ { d } ) } \sim \| f \| _ { \mathbb { R } ^ { d } } + \sup _ { 0 < \delta \leq 1 } \delta ^ { - \gamma } \omega _ { r } ( f, \delta ).$$

Proof. The theorem is already contained in the results in [28], but we need to reconcile notation and explain why. In [28, Formulas (42),(43)] we have defined a univariate K -functional and a pre-modulus of smoothness for g ( x ) = exp( | x | 2 2 ,d / 2) f ( x ) applied to the k -th component of x , k = 1 , · · · , d . The K -functional obtained in this way is denoted in [28, Formula (21)] by K r,k . Likewise, the quantity denoted by ω r in [28] is the k -th summand of the right hand side of (6.46). Our definition of Q ′ k,δ is slightly different from that in [28] (where it is defined to be min( δ -1 , (1 + x 2 k ) 1 / 2 )). However, our Q ′ k,δ as defined in (6.45) satisfies Q ′ k,δ ∼ min( δ -1 , (1 + x 2 k ) 1 / 2 ), Therefore, [28, Theorem 5.1, Proposition 4.5] lead to the statement of this theorem. glyph[squaresolid]

Remark 6.4 If γ = r + β , where r ≥ 0 is an integer and 0 &lt; β ≤ 1, f ∈ C r 0 ( R d ) and satisfies

$$\sup _ { | \mathfrak { u } | _ { 2, d } \leq \delta } \| f ^ { ( r ) } ( \circ + \mathfrak { u } ) - f ^ { ( r ) } \| _ { \mathbb { R } ^ { d } } + \delta \left \| \min ( \delta ^ { - 1 }, | \circ | _ { 2, d } ) f ^ { ( r ) } \right \| _ { \mathbb { R } ^ { d } } \leq c ( f ) \delta ^ { \beta },$$

for every derivative f ( r ) of order r , then ω r ( f, δ ) = O ( δ γ ) for 0 &lt; δ ≤ 1, and f ∈ W γ ( R d ). If f ∈ C r 0 ( R d ) is compactly supported, and every derivative f ( r ) of order r satisfies

$$\sup _ { | \mathfrak { u } | _ { 2, d } \leq \delta } \| f ^ { ( r ) } ( \circ + \mathfrak { u } ) - f ^ { ( r ) } \| _ { \mathbb { R } ^ { d } } \leq c ( f ) \delta ^ { \beta },$$

then f ∈ W γ ( R d ). In particular, if f is compactly supported and satisfies a Lipschitz condition, then f ∈ W 1 ( R d ), and therefore, also f ∈ W γ ( R d ) for every γ ∈ (0 , 1). glyph[squaresolid]

We define

$$\ e \ e n u e \\ \ e \ e n u e \\ \sigma _ { n } ( \mathbb { R } ^ { d } ; f ) ( \mathbf x ) \coloneqq \int _ { \mathbb { R } ^ { g } } \Phi _ { n, d } ( \mathbf x, \mathbf y ) f ( \mathbf y ) d \mathbf y, \quad f \in C _ { 0 } ( \mathbb { R } ^ { d } ), \ n > 0, \ \mathbf x \in \mathbb { R } ^ { d }.$$

The following proposition is routine to prove using Proposition 6.4:

Proposition 6.7 (a) If n &gt; 0 and P ∈ Π d n/ √ 2 , then σ n ( R d ; P ) = P .

$$( \mathbf b ) \ I f \, f \in C _ { 0 } ( \mathbb { R } ^ { d } ), \, n > 0, \, t h e n$$

$$\| \sigma _ { n } ( \mathbb { R } ^ { d } ; f ) \| _ { \mathbb { R } ^ { d } } \leq c \| f \| _ { \mathbb { R } ^ { d } }, \ \ E _ { n } ( \mathbb { R } ^ { d } ; f ) \leq \| f - \sigma _ { n } ( \mathbb { R } ^ { d } ; f ) \| _ { \mathbb { R } ^ { d } } \leq c E _ { n / \sqrt { 2 } } ( \mathbb { R } ^ { d } ; f ).$$

## 7 Approximation on affine spaces

In the sequel, we fix integers Q ≥ q ≥ 1.

Let Y be a q -dimensional affine subspace of R Q , passing through a point x 0 ∈ R Q . Then there exists a rotation operator R on R Q depending only on Y such that any point x ∈ Y can be expressed in the form (with 0 Q q -= (0 , · · · , 0) ∈ R Q q -)

$$x = \colon x _ { 0 } + \mathcal { R } ( \mathbf u, 0 _ { Q - q } ) \quad \mathbf u \coloneqq \mathbf u ( x ) \coloneqq ( u _ { 1 } ( x ), \cdots, u _ { q } ( x ) ).$$

With an abuse of notation, we will write this as x = x 0 + R ( u ). In this section only, the function F : R q → R is defined by

$$F ( \mathbf u ) \coloneqq f \left ( \mathbf x _ { 0 } + \mathcal { R } ( \mathbf u ) \right ),$$

we define

$$E _ { n } ( \mathbb { Y } ; f ) \coloneqq E _ { n } ( \mathbb { R } ^ { q } ; F ).$$

Similarly, if γ &gt; 0, then f ∈ W γ ( Y ) if F ∈ W γ ( R q ); i.e., f ∈ W γ ( Y ) if f ∈ C 0 ( Y ) and

$$\| f \| _ { W _ { \gamma } ( Y ) } \coloneqq \| F \| _ { W _ { \gamma } ( \mathbb { R } ^ { q } ) } < \infty.$$

↦

In terms of the points x = x 0 + R ( u 0 , Q q -) ∈ Y , the class of approximants of functions on Y have the form x → P ( x ) exp( -| x -x 0 | 2 / 2), where P ∈ P Q n 2 . If we are interested only in approximation on Y , we may decide to use some standard point, such as the best approximation to 0 ∈ R Q from Y . This section is meant to be preparatory to Section 8 where the results in this section will be used with Y replaced by the tangent space T x 0 ( X ) to a manifold X . With this goal in mind, our definition is more natural. We note that if f is supported on a compact neighborhood of x 0 , then F is supported on a compact neighborhood of 0 ∈ R q . Therefore, for such functions, we may use Theorem 6.1 (and Remark 6.4) with F and get the estimates where the constants do not depend upon x 0 , although the space of approximants does.

Our goal in this section is to study the analogue of Proposition 6.7 in the context of approximation on Y . We denote the volume measure of Y by µ Y , and for f ∈ C 0 ( Y ), λ &gt; 0, x = x 0 + R ( u ),

$$\sigma _ { n, \lambda } ( \mathbb { Y } ; f ) ( \mathbf x ) \coloneqq \sigma _ { n, \lambda } ( \mathbf x _ { 0 }, \mathbb { Y } ; f ) ( \mathbf x ) \coloneqq \lambda ^ { q } \int _ { \mathbb { Y } } \Phi _ { n, q, Q } ( \lambda ( \mathbf x - \mathbf x _ { 0 } ), \lambda ( \mathbf y - \mathbf x _ { 0 } ) ) f ( \mathbf y ) d \mu _ { \mathbb { Y } } ( \mathbf y ).$$

Theorem 7.1 Let Q ≥ q ≥ 1 be integers, Y be a q -dimensional affine subspace of R Q , passing through x 0 ∈ R Q , f ∈ C 0 ( Y ) , λ &gt; 0 . Then

$$\| \sigma _ { n, \lambda } ( \mathbb { Y } ; f ) - f \| _ { \mathbb { Y } } \leq c E _ { n / \sqrt { 2 } } ( \mathbb { Y } ; f ( \mathbf x _ { 0 } + \mathcal { R } ( ( \circ - \mathbf x _ { 0 } ) / \lambda ) ).$$

In particular, if γ &gt; 0 , f ∈ W γ ( Y ) , λ ≥ 1 , then

$$\| \sigma _ { n, \lambda } ( \mathbb { Y } ; f ) - f \| _ { \mathbb { Y } } \leq c \| f \| _ { W _ { \gamma } ( \mathbb { Y } ) } ( \lambda / n ) ^ { \gamma }.$$

Here, all the constants are independent of λ .

Proof. Since the kernel Φ n,q,Q is invariant under rotations, it is easy to verify that for x = x 0 + R u ∈ Y ,

$$\sigma _ { n, \lambda } ( \mathbb { Y } ; f ) ( \mathbf x ) = \int _ { \mathbb { R } ^ { q } } \Phi _ { n, q } ( \lambda \mathbf u, \mathbf v ) F ( \mathbf v / \lambda ) d \mathbf v.$$

Hence, (7.6) follows from Proposition 6.7. The estimate (7.7) follows from Remark 6.3. glyph[squaresolid]

## 8 Proofs of the theorems in Section 3

For any x ∈ X , we need to consider in this section three kinds of balls, defined in (3.1):

$$B _ { Q } ( { \mathbf x }, { r } ) \coloneqq \{ { \mathbf y } \in \mathbb { R } ^ { Q } \, \colon | { \mathbf x } - { \mathbf y } | _ { 2, Q } \leq { r } \}, \ B _ { \mathbb { T } } ( { \mathbf x }, { r } ) \coloneqq \mathbb { T } _ { \mathbf x } ( { \mathbf X } ) \cap B _ { Q } ( { \mathbf x }, { r } ), \ \mathbb { B } ( { \mathbf x }, { r } ) \coloneqq \{ { \mathbf y } \in { \mathbb { X } } \, \colon \rho ( { \mathbf x }, { \mathbf y } ) \leq { r } \}.$$

Clearly, if r ≤ ι ∗ , then B ( x , r ) = E x ( B T ( x , r )).

The following proposition is not difficult to prove using definitions and Taylor expansions (cf. [6]). In this section, we will simplify the notation to write d u in place of dµ T x ( X ) ( u ).

Proposition 8.1 There exists a constant C ∗ &gt; 0 depending only on X such that each of the following statements holds for every x ∈ X .

(a)

We have

$$| | x - \mathcal { E } _ { x } ( \mathfrak { u } ) | _ { 2, Q } - \rho ( x, \mathcal { E } _ { x } ( \mathfrak { u } ) ) | = | | x - \mathcal { E } _ { x } ( \mathfrak { u } ) | _ { 2, Q } - | x - \mathfrak { u } | _ { 2, Q } | \leq C ^ { * } \rho ( x, \mathcal { E } _ { x } ( \mathfrak { u } ) ) ^ { 3 }, \quad \mathcal { E } _ { x } ( \mathfrak { u } ) \in \mathbb { B } ( x, \iota ^ { * } ). \quad ( 8. 1 )$$

(b) If δ ≤ ι ∗ then

$$| \mathcal { E } _ { \mathbf x } ( \mathbf u ) - \mathbf u | _ { 2, Q } \leq C ^ { * } \delta ^ { 2 }, \quad \mathcal { E } _ { \mathbf x } ( \mathbf u ) \in \mathbb { B } ( \mathbf x, \delta ),$$

(c) If δ ≤ ι ∗ then

$$\int _ { \mathbb { B } ( \mathbf x, \delta ) } | d \mu ^ { * } ( \mathcal { E } _ { \mathbf x } ( \mathbf u ) ) - d \mathbf u | \leq C ^ { * } \delta ^ { q + 2 }.$$

Proof. In this proof only, let r be any geodesic passing through x , parametrized by the arclength s from x , and g be the metric tensor of X . Then, using the fact that | r ′ ( s ) | 2 ,Q = 1, and r ′ ( s ) · r ′′ ( s ) = 0, it is easy to deduce using Taylor expansions that for | s | ≤ ι ∗ ,

$$| | \mathbf r ( s ) - \mathbf x | _ { 2, Q } ^ { 2 } - s ^ { 2 } | \leq c s ^ { 4 } ; \ \text{i.e., } 1 - \frac { | \mathbf r ( s ) - \mathbf x | _ { 2, Q } ^ { 2 } } { s ^ { 2 } } \leq c s ^ { 2 }.$$

Since 1 - | r ( s ) -x | 2 ,Q /s ≤ 1 - | r ( s ) -x | 2 2 ,Q /s 2 , this proves (8.1). The estimate (8.2) follows from the fact that r ( s ) = E x ( x + s r ′ (0)) and a simple estimate using Taylor theorem. The estimate (8.3) follows from the well known fact that in exponential coordinates √ det( g ) = 1 + O ( δ 2 ) in B ( x , δ ) if δ ≤ ι ∗ . glyph[squaresolid]

Corollary 8.1 There exists C ∗ 1 &gt; 0 depending only on X such that for every x y , ∈ X ,

$$| \mathbf x - \mathbf y | _ { 2, Q } \leq \rho ( \mathbf x, \mathbf y ) \leq C _ { 1 } ^ { * } | \mathbf x - \mathbf y | _ { 2, Q }.$$

In particular, for r &gt; 0 , and (3.9) is equivalent to

$$\sup _ { x \in \mathbb { X }, r > 0 } \frac { \mu ^ { * } ( \mathbb { B } _ { Q } ( \mathbf x, r ) ) } { r ^ { q } } \leq c.$$

Proof. In this proof only, let a = min((2 C ∗ ) -1 / 2 , ι ∗ / 2). Then for ρ ( x y , ) ≤ a , (8.1) shows that

$$0 \leq 1 - \frac { | { \mathbf x } - { \mathbf y } | _ { 2, Q } } { \rho ( { \mathbf x }, { \mathbf y } ) } \leq C ^ { * } \rho ( { \mathbf x }, { \mathbf y } ) ^ { 2 } \leq 1 / 2.$$

Therefore,

$$| \mathbf x - \mathbf y | _ { 2, Q } \leq \rho ( \mathbf x, \mathbf y ) \leq 2 | \mathbf x - \mathbf y | _ { 2, Q }, \quad \text{if } \rho ( \mathbf x, \mathbf y ) \leq a.$$

↦

In this proof only, let A = { ( x y , ) ∈ X × X : ρ ( x y , ) ≥ a } . Then A is a compact set and the function ( x y , ) → | x - | y 2 ,Q /ρ ( x y , ), being continuous on A , attains its (necessarily positive) minimum. Thus, there exists c such that

$$| { \mathbf x } - { \mathbf y } | _ { 2, Q } \leq \rho ( { \mathbf x }, { \mathbf y } ) \leq c | { \mathbf x } - { \mathbf y } | _ { 2, Q }, \quad \text{if } \rho ( { \mathbf x }, { \mathbf y } ) \geq a.$$

Together with (8.7), this leads to (8.4), and hence to (8.5). glyph[squaresolid]

To motivate the construction of the operator for approximation, our idea is to transfer the target function locally at each point to the tangent space at that point. Therefore, we use the operator defined as in Section 7. In the present situation, at any point x at which the approximation is desired, the affine space passes through the point x itself, which plays the dual role of x 0 in Section 7. While there is only one parameter t in Theorem 2.1, our construction allows us to have two parameters to control localization: the parameter n controlling the degree of the polynomials involved and an additional parameter to control scaling. Recalling that Φ n,q,Q ( x y , ) = Φ n,q,Q ( -x , -y ) we can define our operator as a convolution as follows.

$$\sigma _ { n, \lambda } ( \mathbb { X } ; f ) ( \mathbf x ) \coloneqq \lambda ^ { q } \int _ { \mathbb { X } } \Phi _ { n, q, Q } ( 0, \lambda ( \mathbf y - \mathbf x ) ) f ( \mathbf y ) d \mu ^ { * } ( \mathbf y ) = \lambda ^ { q } \int _ { \mathbb { X } } \widetilde { \Phi } _ { n, q, Q } ( \lambda | \mathbf x - \mathbf y | _ { 2, Q } ) f ( \mathbf y ) d \mu ^ { * } ( \mathbf y ). \quad \quad ( 8. 8 )$$

Our first theorem is the analogue of Theorem 7.1 when X is a manifold instead of an affine space.

$$\mathbb { B } ( \mathbf x, r ) \subseteq \mathbb { B } _ { Q } ( \mathbf x, r ) \subseteq \mathbb { B } ( \mathbf x, C _ { 1 } ^ { * } r ),$$

Theorem 8.1 Let γ &gt; 0 , f ∈ W γ ( X ) , 0 &lt; α ≤ 1 , α &lt; 4 / γ ( +2) . Then for n ≥ 1 , λ = n 1 -α ,

$$\| f - \sigma _ { n, \lambda } ( \mathbb { X } ; f ) \| _ { \mathbb { X } } \leq c n ^ { - \alpha \gamma } \| f \| _ { W _ { \gamma } ( \mathbb { X } ) }.$$

It is convenient to summarize some details of the proof of this theorem in the form of the following lemma.

Lemma 8.1 Let x ∈ X , g ∈ C ( X ) be supported on B ( x , ι ∗ / 8) , G ( u ) = g ( E x ( u )) , γ &gt; 0 , 0 &lt; α ≤ 1 , α &lt; 4 / γ ( +2) . Then for n ≥ 1 , λ = n 1 -α ,

$$\left | \lambda ^ { q } \int _ { \chi } \widetilde { \Phi } _ { n, q } ( \lambda | x - \mathbf y | _ { 2, Q } ) g ( \mathbf y ) d \mu ^ { * } ( \mathbf y ) - \lambda ^ { q } \int _ { \mathbb { T } _ { x } ( \mathbf x ) } \widetilde { \Phi } _ { n, q } ( \lambda | x - \mathbf u | _ { 2, Q } ) G ( \mathbf u ) d \mathbf u \right | \leq c n ^ { - \alpha \gamma } \| g \| _ { \mathbf x },$$

where G is extended outside B T ( x , ι ∗ / 8) as a zero function.

Proof.

Without loss of generality, we assume that ‖ g ‖ X = 1. First, we summarize our choices of various parameters. In this proof only, let

$$\delta = n ^ { - ( ( 2 - \alpha ) ( q + 1 ) + \alpha \gamma ) / ( q + 3 ) },$$

so that for sufficiently large n ,

$$\delta < \min ( 1, \iota ^ { * } / 6 ), \ \ n ^ { q + 1 } \lambda ^ { q + 1 } \delta ^ { q + 3 } = n ^ { - \alpha \gamma }, \ \ n \lambda \delta = n ^ { ( 4 - \alpha \gamma - 2 \alpha ) / ( q + 3 ) } \uparrow \infty.$$

We choose

$$S \geq \frac { ( q ( 2 - \alpha ) + \alpha \gamma + 1 ) ( q + 3 ) } { 4 - \alpha \gamma - 2 \alpha }, \ \ ( \Rightarrow n ^ { q } \lambda ^ { q } ( n \lambda \delta ) ^ { - S } \leq n ^ { - \alpha \gamma - 1 } ).$$

We now assume further that n is large enough so that with C ∗ as in Proposition 8.1, C δ ∗ 2 ≤ δ/ 2.

Next, we summarize the implications of our choices on the distances on the manifold, tangent space, and the ambient space.

If y ∈ B ( x , ι ∗ / 8) ∩ B Q ( x , δ ), u ∈ T x ( X ), y = E x ( u ), then (8.2) shows that

$$| { \mathbf x } - { \mathbf u } | _ { 2, Q } \leq | { \mathbf x } - { \mathbf y } | _ { 2, Q } + | { \mathcal { E } } _ { \mathbf x } ( { \mathbf u } ) - { \mathbf u } | _ { 2, Q } \leq \delta + C ^ { * } \delta ^ { 2 } \leq ( 3 / 2 ) \delta, \quad \rho ( { \mathbf x }, { \mathbf y } ) \leq 3 \delta < \iota ^ { * } / 2. \quad \ ( 8. 1 3 )$$

Thus,

$$E _ { \delta } \coloneqq \mathcal { E } _ { \mathbf x } ^ { - 1 } ( \mathbb { B } ( \mathbf x, \iota ^ { * } / 8 ) \cap B _ { Q } ( \mathbf x, \delta ) ) \subseteq B _ { \mathbb { T } } ( \mathbf x, 3 \delta / 2 ).$$

If u ∈ B T ( x , ι ∗ / 8) then E x ( u ) is well defined. If u ∈ B T ( x , ι ∗ / 8) \ E δ , then (8.2), (8.1) show that

$$| x - \mathbf u | _ { 2, Q } \geq | x - \mathcal { E } _ { \mathbf x } ( \mathbf u ) | _ { 2, Q } - | \mathcal { E } _ { \mathbf x } ( \mathbf u ) - \mathbf u | _ { 2, Q } \geq \delta - C ^ { * } \delta ^ { 2 } \geq \delta / 2.$$

With this preparation, we are now ready to start with the main estimates. Since g is supported on B ( x , ι ∗ / 8), we find that (cf. (8.12), (6.37))

$$\int _ { X \subset B _ { Q } ( \mathbf x, \delta ) } | \widetilde { \Phi } _ { n, q } ( \lambda | \mathbf x - \mathbf y | _ { 2, Q } ) g ( \mathbf y ) | d \mu ^ { * } ( \mathbf y ) & = \int _ { \mathbb { B } ( \mathbf x, t ^ { * } / 8 ) \mathcal { M } _ { Q } ( \mathbf x, \delta ) } | \widetilde { \Phi } _ { n, q } ( \lambda | \mathbf x - \mathbf y | _ { 2, Q } ) g ( \mathbf y ) | d \mu ^ { * } ( \mathbf y ) \\ & \leq \ c n ^ { q } ( n \lambda \delta ) ^ { - S } \leq c n ^ { - \alpha \gamma - 1 } \lambda ^ { - q }.$$

Using (6.38) and (8.1), we deduce that for y = E x ( u ) ∈ B ( x , ι ∗ / 8) ∩ B Q ( x , δ ),

$$| \widetilde { \Phi } _ { n, q } ( \lambda | { \mathbf x } - \mathcal { E } _ { \mathbf x } ( { \mathbf u } ) | _ { 2, Q } ) - \widetilde { \Phi } _ { n, q } ( \lambda | { \mathbf x } - { \mathbf u } | _ { 2, Q } ) | \leq c n ^ { q + 1 } \lambda | { \mathbf x } - \mathcal { E } _ { \mathbf x } ( { \mathbf u } ) | _ { 2, Q } - | { \mathbf x } - { \mathbf u } | _ { 2, Q } | \leq n ^ { q + 1 } \lambda \delta ^ { 3 }.$$

The estimates (8.13) and (3.9) lead further to

$$\left | \int _ { \mathbb { B } ( x, \mu ^ { * } / 8 ) \cap B _ { Q } ( x, \delta ) } d \mu ^ { * } ( \mathbf y ) - \int _ { E _ { \delta } } d \mathbf u \right | \leq \left | \int _ { E _ { \delta } } | d \mu ^ { * } ( \mathcal { E } _ { \mathbf x } ( \mathbf u ) ) - d \mathbf u | \right | \leq c \delta ^ { q + 2 }.$$

In view of (8.11), (8.17) and (8.18), we deduce that

$$& \text{In view of (8.11), (8.17) and (8.18), we deduce that} \\ & \left | \int _ { B ( x, + ^ { \ast } / 8 ) \cap B _ { Q } ( x, \delta ) } \widetilde { \Phi } _ { n, q } ( \lambda | x - y | _ { 2, Q } ) g ( y ) d \mu ^ { \ast } ( y ) - \int _ { E _ { \delta } } \widetilde { \Phi } _ { n, q } ( \lambda | x - \mathfrak { u } | _ { 2, Q } ) G ( \mathfrak { u } ) d \mathfrak { u } \right | \\ & \quad = \ \left | \int _ { E _ { \delta } } \widetilde { \Phi } _ { n, q } ( \lambda | x - \mathcal { E } _ { \mathfrak { x } } ( \mathfrak { u } ) | _ { 2, Q } ) G ( \mathfrak { u } ) d \mu ^ { \ast } ( \mathcal { E } _ { \mathfrak { x } } ( \mathfrak { u } ) ) - \int _ { E _ { \delta } } \widetilde { \Phi } _ { n, q } ( \lambda | x - \mathfrak { u } | _ { 2, Q } ) G ( \mathfrak { u } ) d \mathfrak { u } \right | \\ & \quad \leq \ \left | \int _ { E _ { \delta } } \left ( \widetilde { \Phi } _ { n, q } ( \lambda | x - \mathcal { E } _ { \mathfrak { x } } ( \mathfrak { u } ) | _ { 2, Q } ) - \widetilde { \Phi } _ { n, q } ( \lambda | x - \mathfrak { u } | _ { 2, Q } ) \right ) G ( \mathfrak { u } ) d \mathfrak { u } \right | + c \mathfrak { m } ^ { q } \delta ^ { q + 2 } \\ & \quad \leq \ c n ^ { q + 1 } \lambda \delta ^ { q + 3 } = c n ^ { - \alpha \gamma } \lambda ^ { - q }.$$

The localization estimate (6.37) shows (cf. (8.12)) that

$$\left | \int _ { \mathcal { T } _ { \mathfrak { x } } ( \mathbb { X } ) \, \mathcal { B } _ { \mathfrak { T } } ( \mathfrak { x }, t ^ { * } / 8 ) } \widetilde { \Phi } _ { n, q } ( \lambda | \mathfrak { x } - \mathbf u | _ { 2, Q } ) G ( \mathbf u ) d \mathbf u \right | \leq c n ^ { q } ( n \lambda ) ^ { - S } \leq c n ^ { - \alpha \gamma - 1 } \lambda ^ { - q }.$$

Invoking the localization estimate (6.37) and (8.11), (8.15) again, we deduce that

$$\left | \int _ { B _ { r } ( x, i ^ { * } / 8 ) \L E _ { \delta } } \widetilde { \Phi } _ { n, q } ( \lambda | x - \mathbf u | _ { 2, Q } ) G ( \mathbf u ) d \mathbf u \right | \leq c n ^ { q } ( n \lambda \delta ) ^ { - S } \leq c n ^ { - \alpha \gamma - 1 } \lambda ^ { - q }.$$

The estimates (8.16), (8.19), (8.20) and (8.21) lead to (8.10). glyph[squaresolid]

We are now in a position to prove Theorem 8.1.

## Proof of Theorem 8.1.

Let x ∈ X . Let φ ∈ C ∞ ( X ) be chosen so that φ ( y ) = 1 if y ∈ B ( x , ι ∗ / 16), φ ( y ) = 0 if y ∈ X \ B ( x , ι ∗ / 8), and 0 ≤ φ ( y ) ≤ 1 for y ∈ X . Then the function fφ is supported on B ( x , ι ∗ / 8), and hence, the function F : T x ( X ) → R defined by F ( u ) := f ( E x ( u )) φ ( E x ( u )) is in W γ ( T x ( X )). Clearly,

$$\| F \| _ { \mathbb { T } _ { \times } ( \mathbb { X } ) } \leq \| f \| _ { \mathbb { X } }, \quad \| F \| _ { W _ { \gamma } ( \mathbb { T } _ { \times } ( \mathbb { X } ) ) } \leq \| f \| _ { W _ { \gamma } ( \mathbb { X } ) }.$$

We choose S &gt; q +( αγ +1) (2 / -α ), and write a = ι ∗ / (16 C ∗ 1 ), where C ∗ 1 is the constant defined in Corollary 8.1. Then, the inclusion (8.5) and the localization property (6.37) show that

$$\int _ { X } \widetilde { \Phi } _ { n, q } ( \lambda | x - y | _ { 2, Q } ) & \left ( 1 - \phi ( y ) \right ) f ( y ) d \mu ^ { * } ( y ) \right | = \left | \int _ { X \mathbb { V } ( x, \omega ^ { \ast } / 16 ) } \widetilde { \Phi } _ { n, q } ( \lambda | x - y | _ { 2, Q } ) ( 1 - \phi ( y ) ) f ( y ) d \mu ^ { * } ( y ) \right | \\ & \leq \int _ { X \mathbb { V } \mathbb { B } _ { Q } ( x, \omega ) } \left | \widetilde { \Phi } _ { n, q } ( \lambda | x - y | _ { 2, Q } ) ( 1 - \phi ( y ) ) f ( y ) \right | d \mu ^ { * } ( y ) \\ & \leq c n ^ { q - S } n ^ { - ( 1 - \alpha ) S } \| f \| _ { X } \leq c n ^ { - \alpha \gamma - 1 } \lambda ^ { - q } \| f \| _ { X }.$$

In view of Lemma 8.1,

$$\left | \int _ { \mathbb { X } } \widetilde { \Phi } _ { n, q } ( \lambda | x - \mathbf y | _ { 2, Q } ) \phi ( \mathbf y ) f ( \mathbf y ) d \mu ^ { * } ( \mathbf y ) - \int _ { \mathbb { T } _ { x } ( \mathbb { X } ) } \widetilde { \Phi } _ { n, q } ( \lambda | x - \mathbf u | _ { 2, Q } ) F ( \mathbf u ) d \mathbf u \right | \leq c n ^ { - \alpha \gamma } \lambda ^ { - q } \| f \| _ { \mathbb { X } },$$

so that

$$\text{so that} \\ \left | \int _ { \mathbb { X } } \widetilde { \Phi } _ { n, q } ( \lambda | x - y | _ { 2, Q } ) f ( y ) d \mu ^ { * } ( y ) - \int _ { \mathbb { T } _ { x } ( \mathbb { X } ) } \widetilde { \Phi } _ { n, q } ( \lambda | x - \mathbf u | _ { 2, Q } ) F ( \mathbf u ) d \mathbf u \right | \leq c n ^ { - \alpha \gamma } \lambda ^ { - q } \| f \| _ { \mathbb { X } }.$$

Since F ( x ) = f ( x ), (7.7) in Theorem 7.1 now shows that

$$\left | \lambda ^ { q } \int _ { \chi } \widetilde { \Phi } _ { n, q } ( \lambda | x - \mathbf y | _ { 2, Q } ) f ( \mathbf y ) d \mu ^ { * } ( \mathbf y ) - f ( x ) \right | \leq c ( n / \lambda ) ^ { - \gamma } \| F \| _ { W _ { \mathbf y } ( \mathbf T _ { x } ( \mathbf x ) ) } \leq c n ^ { - \alpha \gamma } \| f \| _ { W _ { \mathbf y } ( \mathbf x ) }.$$

This proves (8.9). glyph[squaresolid]

Our next objective in this section is to obtain the following discretization of Theorem 8.1 based on noise-corrupted random samples of f as in Theorem 3.1.

The proof of Theorem 3.1 is included in that of the following theorem, together with Theorem 8.1 applied with ff 0 in place of f .

Theorem 8.2 We assume the set up as in Theorem 3.1. Then for every n ≥ 1 and M ≥ n q (2 -α )+2 αγ log( n/δ ) we have with λ = n 1 -α ,

$$\text{Prob} _ { \tau } \left ( \left \| \widehat { F } _ { n, \alpha } ( Y ; \circ ) - \sigma _ { n, \lambda } ( \mathbb { X } ; f f _ { 0 } ) \right \| _ { \mathbb { R } ^ { Q } } \geq c \sqrt { \| f _ { 0 } \| _ { \mathbb { X } } } \| \mathcal { F } \| _ { \mathbb { X } \times \Omega } n ^ { - \alpha \gamma } \right ) \leq \delta.$$

The proof of Theorem 8.2 requires some preparation. We start with the following concentration inequality [7, Section 2.7].

Proposition 8.2 ( Bernstein concentration inequality ) Let Z , 1 · · · , Z M be independent real valued random variables such that for each j = 1 , · · · , M , | Z j | ≤ R , and E ( Z 2 j ) ≤ V . Then for any t &gt; 0 ,

$$\text{Prob} \left ( \left | \frac { 1 } { M } \sum _ { j = 1 } ^ { M } ( Z _ { j } - \mathbb { E } ( Z _ { j } ) ) \right | \geq t \right ) \leq 2 \exp \left ( - \frac { M t ^ { 2 } } { 2 ( V + R t / 3 ) } \right ).$$

In order to apply Proposition 8.2, we need to estimate the second moment of F ( y , glyph[epsilon1] )Φ ˜ n,q,Q ( λ | x -y | 2 ,Q ) = F ( y , glyph[epsilon1] )Φ n,q,Q ( 0 , λ ( x -y )) for every x ∈ R Q . This is done in the following lemma.

## Lemma 8.2 We have

$$\lambda ^ { 2 q } \sup _ { x \in \mathbb { R } ^ { Q } } \int _ { \mathbb { X } \times \Omega } | \mathcal { F } ( \mathbf y, \epsilon ) \Phi _ { n, q, Q } ( 0, \lambda ( x - y ) ) | ^ { 2 } d \tau ( \mathbf y, \epsilon ) \leq c ( n \lambda ) ^ { q } \| \mathcal { F } \| _ { \mathbb { X } \times \Omega } ^ { 2 } \| f _ { 0 } \| _ { \mathbb { X } }.$$

Proof. Let x ∈ R Q . We need only to estimate

$$\int _ { \mathbb { X } \times \Omega } | \mathcal { F } ( \mathbf y, \epsilon ) \Phi _ { n, q, Q } ( \mathbf 0, \lambda ( x - \mathbf y ) ) | ^ { 2 } \, d \tau ( \mathbf y, \epsilon ) \leq \| \mathcal { F } \| _ { \mathbb { X } \times \Omega } ^ { 2 } \| f _ { 0 } \| _ { \mathbb { X } } \int _ { \mathbb { X } } \Phi _ { n, q, Q } ( \mathbf 0, \lambda ( x - \mathbf y ) ) ^ { 2 } d \mu ^ { * } ( \mathbf y ).$$

Using Proposition 6.3 and (8.6), and keeping in mind that λ ≥ 1, we deduce that

$$& \text{Using Proposition 6.3 and (8.6), and keeping in mind that } \lambda \geq 1, \, \text{we deduce that} \\ & \int _ { \mathfrak { G } _ { n, q } ( 0, \lambda ( x - y ) ) ^ { 2 } d \mu ^ { * } ( y ) } = \int _ { \mathfrak { X } \cap \mathfrak { B } _ { Q } ( x, 1 / ( n \lambda ) ) } \Phi _ { n, q, Q } ( 0, \lambda ( x - y ) ) ^ { 2 } d \mu ^ { * } ( y ) \\ & \quad + \sum _ { k = 0 } ^ { \infty } \int _ { \mathfrak { X } \cap ( \mathfrak { B } _ { Q } ( x, 2 ^ { k + 1 } / ( n \lambda ) ) \mathfrak { B } _ { Q } ( x, 2 ^ { k } / ( n \lambda ) ) ) } \Phi _ { n, q, Q } ( 0, \lambda ( x - y ) ) ^ { 2 } d \mu ^ { * } ( y ) \\ & \quad \leq c n ^ { 2 q } \left \{ \mu ^ { * } ( \mathfrak { B } _ { Q } ( x, 1 / ( n \lambda ) ) ) + \sum _ { k = 0 } ^ { \infty } 2 ^ { - k ^ { 2 } k ^ { 2 } } \mu ^ { * } ( \mathfrak { B } _ { Q } ( x, 2 ^ { k + 1 } / ( n \lambda ) ) \, \mathfrak { B } _ { Q } ( x, 2 ^ { k } / ( n \lambda ) ) ) \right \} \\ & \quad \leq c n ^ { q } \lambda ^ { - q } \left \{ 1 + \sum _ { k = 0 } ^ { \infty } 2 ^ { - k ( 2 ^ { k } - q ) } \right \} \leq c n ^ { q } \lambda ^ { - q }. \\ & \quad \text{$$

glyph[squaresolid]

The proof of Theorem 8.2 requires an estimation of a quantity of the form

$$\sup _ { y _ { 1 }, \cdots, \mathcal { Y } _ { M } \in \mathbb { X } } \left \| \frac { \lambda ^ { q } } { M } \sum _ { j = 1 } ^ { M } \mathcal { F } ( \mathbf y _ { j }, \epsilon _ { j } ) \Phi _ { n, q, Q } ( 0, \lambda ( \circ - \mathbf y _ { j } ) ) - \sigma _ { n, \lambda } ( \mathbb { X } ; f f _ { 0 } ) \right \| _ { \mathbb { R } ^ { Q } }$$

in terms of the maximum of the function involved at finitely many points. The following lemma accomplishes this by considering the difference between two measures on X : one that associates the mass (1 /M ) F ( y j , glyph[epsilon1] j ) with each y j , and other given by f ( y ) dν ∗ ( y ) = f ( y ) f 0 ( y ) dµ ∗ ( y ). We will denote the total variation of a measure ν by | | | ν | | | TV . The total variation of the difference between the two measures mentioned above is clearly ≤ ‖F‖ 2 X × Ω .

Lemma 8.3 Let S &gt; Q +2 , λ be as in Theorem 8.1. There exists c ∗ = c ∗ ( S ) &gt; 0 and a finite set D ∗ ⊂ R Q with |D | ∼ ∗ n c ∗ such that for any measure ν on X ,

$$\left \| \lambda ^ { \omega } \int _ { \chi } \Phi _ { n, q, Q } ( 0, \lambda ( \circ - \bar { \mathbf y } ) ) d \nu ( \bar { \mathbf y } ) \right \| _ { \mathbb { R } ^ { Q } } \leq \max _ { \bar { \mathbf x } \in \mathcal { D } ^ { * } } \left | \lambda ^ { \omega } \int _ { \bar { \mathbf x } } \Phi _ { n, q, Q } ( 0, \lambda ( \bar { \mathbf x } - \bar { \mathbf y } ) ) d \nu ( \bar { \mathbf y } ) \right | + c n ^ { - S } \| \nu \| _ { T V }.$$

Proof. We assume n to be large enough so that X ⊂ -[ √ 2 n, √ 2 n ] Q . Then Proposition 6.5 (used with 2 S in place of S ) shows that

$$\sup _ { x \in \mathbb { R } ^ { Q } \bigcup [ - 2 n, 2 n ] ^ { Q } } \left | \int _ { \mathbb { X } } \Phi _ { n, q, Q } ( 0, \lambda ( x - y ) ) d \nu ( y ) \right | \leq c n ^ { Q - 2 S } \| \nu \| _ { T V } \leq c n ^ { - S } \| \nu \| _ { T V }.$$

Therefore,

$$\left \| \int _ { \mathbb { Z } } \Phi _ { n, q, Q } ( 0, \lambda ( \circ - y ) ) d \nu ( y ) \right \| _ { \mathbb { R } ^ { Q } } \leq \sup _ { \mathbb { x } \in [ - 2 n, 2 n ] ^ { Q } } \left | \int _ { \mathbb { X } } \Phi _ { n, q, Q } ( 0, \lambda ( \mathbb { x } - y ) ) d \nu ( y ) \right | + c n ^ { - S } \| \nu \| _ { T V }.$$

Next, we observe that for any y ∈ R Q

$$\nabla _ { x } \left ( \Phi _ { n, q, Q } ( 0, \lambda ( x - y ) ) \right ) = \lambda \left ( \nabla _ { x } \Phi _ { n, q, Q } ( 0, \circ ) \right ) ( \lambda ( x - y ) ).$$

Therefore, for any x ∈ R Q ,

$$\left | \nabla _ { \mathbf x } \left ( \int _ { \mathbf X } \Phi _ { n, q, Q } ( 0, \lambda ( \mathbf x - \mathbf y ) ) d \nu ( \mathbf y ) \right ) \right | \leq \lambda \int _ { \mathbf X } \left | \left ( \nabla _ { \mathbf x } \Phi _ { n, q, Q } ( 0, \circ ) \right ) \right ) \left ( \lambda ( \mathbf x - \mathbf y ) \right ) \right | d | \nu | ( \mathbf y ).$$

Using the Bernstein inequality Proposition 6.1(c), we conclude that

$$\sup _ { x \in \mathbb { R } ^ { Q } } \left | \nabla _ { x } \left ( \int _ { \mathbb { X } } \Phi _ { n, q, Q } ( 0, \lambda ( x - y ) ) d \nu ( y ) \right ) \right | \leq c n ^ { q + 1 } \lambda \| \nu \| _ { T V } = c n ^ { q + 2 - \alpha } \| \nu \| _ { T V }.$$

and hence, for any x z , ∈ R Q ,

$$\left | \lambda ^ { q } \int _ { X } \Phi _ { n, q, Q } ( 0, \lambda ( x - y ) ) d \nu ( \mathbf y ) - \lambda ^ { q } \int _ { X } \Phi _ { n, q, Q } ( 0, \lambda ( z - y ) ) d \nu ( \mathbf y ) \right | \leq c n ^ { ( q + 1 ) ( 2 - \alpha ) } \| \nu \| _ { T V } | x - z | _ { \infty, Q }.$$

We now let D ∗ be a finite subset of [ -2 n, 2 n ] Q such that

$$\max _ { x \in [ - 2 n, 2 n ] ^ { Q } } \min _ { z \in \mathcal { D } ^ { * } } | x - z | _ { \infty, Q } \leq n ^ { - ( q + 1 ) ( 2 - \alpha ) - S },$$

and observe that |D | ∼ ∗ n Q q (( +1)(2 -α )+ S ) . The estimate (8.30) is easy to deduce using (8.32), (8.33), and (8.34). glyph[squaresolid]

With this preparation, we now prove Theorem 8.2, and hence, Theorem 3.1.

Proof of Theorem 8.2 (and Theorem 3.1).

Let x ∈ R Q . We consider the random variables

$$Z _ { j } ( x ) & = \lambda ^ { q } \mathcal { F } ( y _ { j }, \epsilon _ { j } ) \Phi _ { n, q, Q } ( 0, \lambda ( x - y _ { j } ) ).$$

It is easy to verify using Fubini's theorem that if F is integrable with respect to τ then for any x ∈ R Q ,

$$\mathbb { E } _ { \tau } ( \lambda ^ { q } \mathcal { F } ( y, \epsilon ) \Phi _ { n } ( 0, \lambda ( x - y ) ) = \sigma _ { n, \lambda } ( \mathbb { X } ; f f _ { 0 } ) ( \mathbf x ).$$

The estimate (6.35) implies that | Z j | ≤ c nλ ( ) q ‖F‖ X × Ω . Further, Lemma 8.2 yields E τ ( Z 2 j ) ≤ c nλ ( ) q ‖F‖ 2 X × Ω ‖ f 0 ‖ X . Therefore, we deduce using Proposition 8.2 that for any t ∈ (0 1), ,

$$\text{Prob} _ { \tau } \left ( \left | \frac { 1 } { M } \sum _ { j = 1 } ^ { M } Z _ { j } ( \mathbf x ) - \sigma _ { n, \lambda } ( \mathbf X ; f f _ { 0 } ) ( \mathbf x ) \right | \geq t \| \mathcal { F } \| _ { \mathbf X \times \Omega } \| f _ { 0 } \| _ { \mathbf x / 2 } \right ) \leq 2 \exp \left ( - c \frac { M \| f _ { 0 } \| _ { \mathbf x } t ^ { 2 } } { ( n \lambda ) ^ { q } } \right ).$$

In view of Lemma 8.3, we have for S ≥ Q +2+ αγ ,

$$\text{Prob} _ { \tau } \left ( \left \| \frac { 1 } { M } \sum$$

We recall that nλ = n 2 -α and choose

$$t = c _ { 3 } \sqrt { \frac { n ^ { q ( 2 - \alpha ) } } { M \$$

for a suitable constant to make the right hand side of (8.38) to be ≤ δ , to obtain

$$P r o b _ { r } \left ( \left \| \frac { 1 } { M } \sum _ { j =$$

Wenow observe that since 1 = ∫ X f dµ 0 ∗ , and µ ∗ ( X ) = 1, ‖ f 0 ‖ X ≥ 1. Therefore, choosing M ≥ n q (2 -α )+2 αγ √ log( n/δ ), we arrive at (8.26). glyph[squaresolid]

Theorem 3.2 is obtained immediately from Theorem 3.1 by setting f 0 ≡ 1. To obtain Theorem 3.3, we use Theorem 3.1 once as stated and again with F ( Y ; ◦ ) ≡ 1 to get an approximation to f 0 .

## 9 Proof of the theorems in Section 5

Proof of Theorem 5.1

Theorem 5.1 follows easily from Theorem 8.2 and Corollary 6.2.

Proof of Theorem 5.2 .

Let v ∈ V , and u , 1 · · · , u d v ( ) be the children of v , and x 1 , · · · , x d v ( ) be the inputs seen by these in that order. Let x be the corresponding input seen by v . Then using the Lipschitz condition on f v and the property (5.7), we obtain

$$\text{ obtain} \quad \begin{array} {$$

We now use induction on the level of v . Thus, if v ∗ ∈ S , then the 'shallow network' estimate implied in Theorem 5.1 is already the one which we want. Suppose the theorem is proved for the DAGs for which the sink node is at level glyph[lscript] ≥ 0. If v ∈ V , so that its level glyph[lscript] ≥ 1, then its children are at level glyph[lscript] -1 ≥ 0. For each of the children, say u , we consider the subgraph G u of G comprising only those nodes and edges that culminate in u as the sink node. We then apply the theorem to each of these subgraphs, and then use (9.1) to conclude that the statement is true for the subgraph G v of G comprising only those nodes and edges that culminate in v as the sink node. glyph[squaresolid]

Remark 9.1 Suppose we consider a shallow Gaussian network acting on a 2 s dimensional manifold of R Q . The number of samples required to obtain an accuracy of n -αγ predicted by Theorem 5.1 is O ( n 2 (2 s -α )+2 αγ log n ). On the other hand, suppose the target function has a compositional structure according to a binary tree, but in addition, for any v ∈ V with children u , u 1 2 , the image of ( f u 1 , f u 2 ) forms a curve in R 2 . Then the number of samples required to get the same accuracy with the corresponding network is only O ( n 2 -α +2 αγ log n ) at each level. In fact, it seems likely that this is the number of samples in the orignal submanifold of R Q itself, since the input variables external to the machine are given only at the source nodes. glyph[squaresolid]

## 10 Conclusions

We have given a direct solution to the problem of function approximation if the data is sampled from a compact, smooth, connected Riemannian manifold, without knowing the manifold itself , except for its dimension. Our construction avoids the evaluation of an eigen-decomposition of a matrix or otherwise the need to compute the local charts on the manifold. Also, the construction avoids any optimization/training in the classical paradigm.

Our construction is universal; i.e., can be used for any target function without any assumption on its prior. The approximation error is estimated in the probabilistic sense, and of course, depends upon the smoothness of the target function. In the case when the data is taken from an affine space, our approximation error does not suffer from any saturation, but can be as small as the smoothness of the target function allows. In the general case, the curvature of the manifold imposes some limitations on how well we can estimate the degree of approximation, but there is no saturation in the sense that if the degree of approximation is better for a function, then it must be 'trivial' in some sense.

We have extended our results to the case of deep Gaussian networks. However, in this context, they are not completely constructive unless the constituent functions in the DAG defining the deep network are known.

## A Saturation phenomenon

The notation in this section is not the same as that in the rest of the paper, except that ‖ · ‖ A will denote the supremum norm on a set A . A detailed discussion of saturation phenomena in approximation theory can be found in [8]. Intuitively, an approximation process on a metric space A is a sequence of operators U n : C A ( ) → C A ( ) such that U n ( f ) → f uniformly on A . The process is saturated with the rate { δ n } if ‖ U n ( f ) -f ‖ A = o δ ( n ) as n →∞ implies that f is trivial in some sense (classically U n ( f ) = f ) and there exists a non-trivial function f for which ‖ U n ( f ) - ‖ f A = O ( δ n ). We are unable to find in the literature a precise definition that covers the many applications where this phenomenon holds. As remarked earlier, Theorem 2.1 is one example. We give two other examples.

Example A.1 For f ∈ C ([ -1 1]), the Bernstein polynomial is defined by ,

$$B _ { n } ( f ) ( x ) \coloneqq \sum _ { k = 0 } ^ { n } \bin$$

The Voronowskaja theorem ([25, Section 1.6.1]) states that if f ∈ C 2 ([ -1 1]) then uniformly in , x ∈ -[ 1 1], ,

$$\lim _ { n \to \infty } \left | n \left ( B _ { n } ( f ) ( x ) - f ( x ) \right ) - f ^ { \prime \prime } ( x ) \frac { x ( 1 - x ) } { 2 } \right | = 0.$$

Thus, f ∈ C 2 ([ -1 1]), , ‖ B n ( f ) - ‖ f [ -1 1] , = O (1 /n ) and if ‖ B n ( f ) - ‖ f [ -1 1] , = (1 o /n ) then f ′′ ( x ) = 0 for x ∈ -( 1 1), , so that f is a linear function. glyph[squaresolid]

Example A.2 A function S : [ -1 1] , → R is called piecewise constant with n break-points if there are points t 0 = -1 &lt; t 1 &lt; · · · &lt; t n +1 = 1 such that S is a constant on each ( t j , t j +1 ), j = 0 , · · · , n . We denote the class of all piecewise constants with n break-points by S n , and define for f ∈ C ([ -1 1]), ,

$$\sigma _ { n } ( f ) \coloneqq \inf _ { S \in \mathcal { S } _ { n } } \| f - S \| _ { [ - 1, 1 ] }.$$

We note that the break-points of the approximating function may depend upon the target function f . It is known ([17, Chapter 12, Theorem 4.3, Corollary 4.4]) that if f has a bounded total variation on [ -1 1] then , σ n ( f ) = O (1 /n ). Moreover, if f ∈ C ([ -1 1]) and , σ n ( f ) = o (1 /n ) then f is a constant. glyph[squaresolid]

## List of Symbols

B Q ( x , r ), B T ( x , r ), B ( x , r ) Defined in (3.1)

$$\Delta _ { k, \ell } ^ { j }, Q _ { k, \delta } ^ { \prime }, \omega _ { r } \, \text{Section 6.4}$$

ι ∗ Inradius of X

λ Scaling factor, typically, n 1 -α

λ k,m , λ k ,m Quadrature weights, Section 6.1

G ∗ n,q,Q , G n,q,Q Special Gaussian network (5.1), (5.4)

P d n , Π d n Polynomial spaces Section 6.1

E x Exponential map at x ∈ X , E x : T x ( X ) → X

G DAG for deep networks, Section 5.2

P m,q , ˜ Φ n,q Univariate polynomials defined in (3.5), (3.6)

G k ,m,d , G Q Basic Gaussian networks (6.39), (6.41)

Proj m,d , Proj m,q,Q Projection kernels (6.14), (6.22)

Lip( X ) Lipschitz functions on X

- µ ∗ Volume measure on X
- Φ n,d , Φ n,q,Q Localized kernels (6.29), (6.33)
- π v Pooling operation Section 5.2
- ρ Metric on X
- σ n , σ n,λ Approximation operators (6.51), (7.5), (8.8)
- τ Probability distribution for the data
- T x ( X ) Tangent space to X at x
- X Manifold
- Y Affine space
- d Generic dimension, Section 6
- d v ( ) Ambient dimension at vertex v Section 5.2
- E n ( A f ; ) Degree of approximation of f on A
- f 0 Density of the marginal distribution
- H Low pass filter Section 3.1
- h k , ψ k , ψ k Orthonormalized Hermite polynomial, Hermite function, tensor product Hermite function
- n , α Parameters in approximation
- Q Dimension of the ambient space
- q Dimension of affine space or manifold
- q v Dimension in Section 5.2
- S Large integer controlling localization
- V , S Non-source, source vertices Section 5.2
- v v Constituent function at v Section 5.2
- W A γ ( ) Smoothness class on A
- f , F , F n,α
- ̂ Target function, observations, and estimator x k,m , x k ,m Quadrature nodes Section 6.1

## References

- [1] G. E. Andrews, R. Askey, and R. Roy. Special functions , volume 71. Cambridge university press, 1999.
- [2] R. Askey and S. Wainger. Mean convergence of expansions in Laguerre and Hermite series. American Journal of Mathematics , 87(3):695-708, 1965.
- [3] H. Bateman, A. Erd´ elyi, W. Magnus, F. Oberhettinger, and F. G. Tricomi. Higher transcendental functions , volume 2. McGraw-Hill New York, 1955.
- [4] M. Belkin and P. Niyogi. Laplacian eigenmaps for dimensionality reduction and data representation. Neural computation , 15(6):1373-1396, 2003.
- [5] M. Belkin and P. Niyogi. Semi-supervised learning on Riemannian manifolds. Machine learning , 56(1-3):209239, 2004.
- [6] M. Belkin and P. Niyogi. Towards a theoretical foundation for Laplacian-based manifold methods. Journal of Computer and System Sciences , 74(8):1289-1308, 2008.
- [7] S. Boucheron, G. Lugosi, and P. Massart. Concentration inequalities: A nonasymptotic theory of independence . Oxford university press, 2013.
- [8] P. L. Butzer and R. J. Nessel. Fourier analysis and approximation , volume 40. Academic Press, 2011.
- [9] M. Chen, H. Jiang, W. Liao, and T. Zhao. Efficient approximation of deep ReLU networks for functions on low dimensional manifolds. arXiv preprint arXiv:1908.01842 , 2019.
- [10] C. K. Chui and D. L. Donoho. Special issue: Diffusion maps and wavelets. Appl. and Comput. Harm. Anal. , 21(1), 2006.
- [11] C. K. Chui and H. N. Mhaskar. Deep nets for local manifold learning. Frontiers in Applied Mathematics and Statistics , 4:12, 2018.
- [12] C. K. Chui and H. N. Mhaskar. A Fourier-invariant method for locating point-masses and computing their attributes. Appl. Comput. Harmon. Anal. , 45:436-452, 2018.
- [13] C. K. Chui and H. N. Mhaskar. A unified method for super-resolution recovery and real exponential-sum separation. Appl. Comput. Harmon. Anal. , 46(2):431-451, March 2019.

| [14]   | A. Cloninger, R. R. Coifman, N. Downing, and H. M. Krumholz. Bigeometric organization of deep nets. arXiv preprint arXiv:1507.00220 , 2015.                                    |
|--------|--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [15]   | F. Cucker and S. Smale. On the mathematical foundations of learning. Bulletin of the American Mathematical Society , 39:1-49, 2002.                                            |
| [16]   | F. Cucker and D. X. Zhou. Learning theory: an approximation theory viewpoint , volume 24. Cambridge University Press, 2007.                                                    |
| [17]   | R. A. DeVore and G. G. Lorentz. Constructive approximation , volume 303. Springer Science & Business Media, 1993.                                                              |
| [18]   | M. P. do Carmo Valero. Riemannian geometry . Birkh¨user, a 1992.                                                                                                               |
| [19]   | M. Ehler, F. Filbir, and H. N. Mhaskar. Locally learning biomedical data using diffusion frames. Journal of Computational Biology , 19(11):1251-1264, 2012.                    |
| [20]   | F. Filbir and H. N. Mhaskar. Marcinkiewicz-Zygmund measures on manifolds. Journal of Complexity 27(6):568-596, 2011.                                                           |
| [21]   | F. Girosi and T. Poggio. Networks and the best approximation property. Biological cybernetics , 63(3):169-176, 1990.                                                           |
| [22]   | P. W. Jones, M. Maggioni, and R. Schul. Universal local parametrizations via heat kernels and eigenfunctions of the Laplacian. Ann. Acad. Sci. Fenn. Math. , 35:131-174, 2010. |
| [23]   | S. S. Lafon. Diffusion maps and geometric harmonics . PhD thesis, Yale University, Yale, 2004.                                                                                 |
| [24]   | W. Liao and M. Maggioni. Adaptive geometric multiscale approximations for intrinsically low-dimensional data. arXiv preprint arXiv:1611.01179 , 2016.                          |
| [25]   | G. G. Lorentz. Bernstein polynomials . American Mathematical Soc., 2013.                                                                                                       |
| [26]   | M. Maggioni and H. N. Mhaskar. Diffusion polynomial frames on metric measure spaces. Applied and Com- putational Harmonic Analysis , 24(3):329-353, 2008.                      |
| [27]   | H. N. Mhaskar. Introduction to the theory of weighted polynomial approximation , volume 56. World Scientific Singapore, 1996.                                                  |
| [28]   | H. N. Mhaskar. On the degree of approximation in multivariate weighted approximation. In Advanced Problems in Constructive Approximation , pages 129-141. Springer, 2003.      |
| [29]   | H. N. Mhaskar. When is approximation by Gaussian networks necessarily a linear process? Neural Networks 17(7):989-1001, 2004.                                                  |
| [30]   | H. N. Mhaskar. A Markov-Bernstein inequality for Gaussian networks. In Trends and applications in construc- tive approximation , pages 165-180. Springer, 2005.                |
| [31]   | H. N. Mhaskar. Eignets for function approximation on manifolds. Applied and Computational Harmonic Analysis , 29(1):63-87, 2010.                                               |
| [32]   | H. N. Mhaskar. A generalized diffusion frame for parsimonious representation of functions on data defined manifolds. Neural Networks , 24(4):345-359, 2011.                    |
| [33]   | H. N. Mhaskar. Local approximation using Hermite functions. In Progress in Approximation Theory and Applicable Complex Analysis , pages 341-362. Springer, 2017.               |
| [34]   | H. N. Mhaskar. A unified framework for harmonic analysis of functions on directed graphs and changing data. Appl. Comput. Harm. Anal. , 44(3):611-644, 2018.                   |
| [35]   | H. N. Mhaskar and T. Poggio. Deep vs. shallow networks: An approximation theory perspective. Analysis and Applications , 14(06):829-848, 2016.                                 |
| [36]   | H. N. Mhaskar and T. Poggio. An analysis of training and generalization errors in shallow and deep networks. Neural Networks , 121:229-241, 2020.                              |

- [37] J. Schmidt-Hieber. Deep ReLU network approximation of functions on a manifold. arXiv preprint arXiv:1908.00695 , 2019.
- [38] A. Singer. From graph to manifold Laplacian: The convergence rate. Applied and Computational Harmonic Analysis , 21(1):128-134, 2006.
- [39] G. Szeg¨. o Orthogonal polynomials. In Colloquium publications/American mathematical society , volume 23. Providence, 1975.
Hrushikesh Mhaskar
A direct approach for function approximation on data defined manifolds
In much of the literature on function approximation by deep networks, the
function is assumed to be defined on some known domain, such as a cube or a
sphere. In practice, the data might not be dense on these domains, and
therefore, the approximation theory results are observed to be too
conservative. In manifold learning, one assumes instead that the data is
sampled from an unknown manifold; i.e., the manifold is defined by the data
itself. Function approximation on this unknown manifold is then a two stage
procedure: first, one approximates the Laplace-Beltrami operator (and its
eigen-decomposition) on this manifold using a graph Laplacian, and next,
approximates the target function using the eigen-functions. Alternatively, one
estimates first some atlas on the manifold and then uses local approximation
techniques based on the local coordinate charts.
  In this paper, we propose a more direct approach to function approximation on
\emph{unknown}, data defined manifolds without computing the
eigen-decomposition of some operator or an atlas for the manifold, and without
any kind of training in the classical sense. Our constructions are universal;
i.e., do not require the knowledge of any prior on the target function other
than continuity on the manifold. We estimate the degree of approximation. For
smooth functions, the estimates do not suffer from the so-called saturation
phenomenon. We demonstrate via a property called good propagation of errors how
the results can be lifted for function approximation using deep networks where
each channel evaluates a Gaussian network on a possibly unknown manifold.