Distributed Learning via Filtered Hyperinterpolation on Manifolds

Learning mappings of data on manifolds is an important topic in contemporary machine learning, with applications in astrophysics, geophysics, statistical physics, medical diagnosis, biochemistry, and 3D object analysis. This paper studies the problem of learning real-valued functions on manifolds through filtered hyperinterpolation of input–output data pairs where the inputs may be sampled deterministically or at random and the outputs may be clean or noisy. Motivated by the problem of handling large data sets, it presents a parallel data processing approach which distributes the data-fitting task among multiple servers and synthesizes the fitted sub-models into a global estimator. We prove quantitative relations between the approximation quality of the learned function over the entire manifold, the type of target function, the number of servers, and the number and type of available samples. We obtain the approximation rates of convergence for distributed and non-distributed approaches. For the non-distributed case, the approximation order is optimal.


Introduction
Learning functions over manifolds has become an increasingly important topic in machine learning. The performance of many machine learning algorithms depends strongly on the geometry of the data. In real-world applications, one often has huge data sets with noisy samples. In this paper, we propose distributed filtered hyperinterpolation on manifolds, which combines filtered hyperinterpolation and distributed learning [22,23]. Filtered hyperinterpolation [31,38] provides a constructive approach to modeling mappings between inputs and outputs in a way that can reduce the influence of noise. The distributed strategy assigns the learning task of the input-output mapping to multiple local servers, enabling parallel computing for massive data sets. Each server handles a small fraction of all data by filtered hyperinterpolation. It then synthesizes the local estimators as a global estimator. We show the precise quantitative relation between the approximation error of the distributed filtered hyperinterpolation, the number of the local servers, and the amount of data. The approximation error (over the entire manifold) converges to zero provided the available amount of data increases sufficiently fast with the number of servers.
Filtered hyperinterpolation was introduced by Sloan and Womersley [31] on the two-sphere S 2 , which is a form of filtered polynomial approximation method motivated by hyperinterpolation [29]. Hyperinterpolation uses a Fourier expansion where the integral for the Fourier coefficients is approximated by numerical integration with a quadrature rule. The filtered hyperinterpolation adopts a similar strategy as hyperinterpolation but uses a filter to modify the Fourier expansion. The filter is a restriction on the eigenvalues of the basis functions. Effectively this restricts the capacity of the approximation class and yields a reproducing property for polynomials of a certain degree specified by the filter. It has some similarities to kernel methods. Filtering improves the approximation accuracy of plain hyperinterpolation for noiseless data that is sampled deterministically [18]. With appropriate choice of filter, the filtered hyperinterpolation achieves the best approximation by polynomials of a given degree depending on the amount of data (see Sect. 3.1). As shown in the left part of Fig. 1, one aims at finding the closest approximation of f * within the polynomial space Π n on the manifold M, which, nevertheless, is difficult to achieve. The filtered hyperinterpolation is an approximator V D,n constructed from data D = (x i , y i ) N i=1 which lies in a slightly larger polynomial space Π 2n and whose distance to f * is very close to the distance between f * and Π n .
Motivated by the problem of handling massive amounts of data, we propose a distributed computational strategy based on filtered hyperinterpolation. As shown in the right part of Fig. 1, we can split estimation task of filtered hyperinterpolation into multiple servers j = 1, . . . , m, each of which computes a filtered hyperinterpolation V D j ,n , for a small subset D j of all the training data. It consists of creating a filtered expansion in terms of eigenfunctions of the manifold to best-fit the corresponding fraction of the training data set. The "best-fit" means that the local servers can achieve best approximation for noisy data y i = f * (x i ) + i , i = 1, . . . , N , for any continuous function f * : M → R on the manifold and independent bounded noise i . The central processor then takes a weighted average of the filtered hyperinterpolations obtained in the local servers to synthesize as a global estimator V (m) D,n . We call the global estimator the distributed filtered hyperinterpolation.
The remaining of the paper is organized as follows. In Sect. 2, we introduce the main mathematical settings and notation. Then, we proceed with the study of non-distributed and distributed filtered hyperinterpolation on manifolds, for which we derive upper bounds on the error. Our bounds depend on (1) the dimension d of the manifold and the smoothness r of the Sobolev space that contains the target function, (2) the degree n of the approximating polynomials, which is tied to the number N of available data points, (3) the smoothness of the filter, (4) the presence of noise in the output data points. Here, we base the analysis on properties of the quadrature formulas, which we couple with the arrangement of the input data points (deterministic or random). For the deterministic case, we require the quadrature rule has polynomial exactness of degree 3n − 1; for the random case, the condition that the volume measure on the manifold controls the distribution of the sampling points.
In Sect. 3, we study non-distributed filtered interpolation on manifolds. We obtain an error bound O N −r /d for the noiseless setting on general manifolds (see Theorem 6). This result generalizes the same bound that was previously obtained on the sphere [35]. Since the bound on the sphere is optimal, the new bound is also optimal. We further study learning with noisy output data. The error bound for the noisy case is O N −2r /(2r +d) . Due to the impact of the noise, it does not entirely reduce to the error bound of the noiseless case. To the best of our knowledge, this is the first error upper bound for noisy learning on general Riemannian manifolds. The optimality of this bound remains open at this point.
In Sect. 4, we study distributed learning. We obtain similar rates of convergence as in the non-distributed setting, provided the number of servers satisfies a certain upper bound in terms of the total amount of data. As it turns out, the distributed estimator has the same convergence rate as the non-distributed estimator for a class of functions with given smoothness. Compared with the clean data case, the distributed filtered hyperinterpolation with noisy data has slightly lower convergence order than the non-distributed. See Theorems 13 and 14.
Section 5 illustrates definitions, methods, and convergence results on a concrete numerical example. Section 6 summarizes and compares the convergence rates of the different methods and settings (see Table 1). It also presents a concise description of the implementation (see Algorithm 1). All the proofs are deferred to Appendix A. The proofs utilize the wavelet decomposition of filtered hyperinterpolation, Marcinkiewicz-Zygmund inequality, Nikolskiî-type inequality on a manifold, bounds of best approximation on a manifold, and concentration inequality, estimates of covering number and bounds of sampling operators from learning theory. We also show a table of the notations used throughout the article in Appendix B for readers' reference. In the left part, V D,n is a filtered hyperinterpolation function constructed from a data set D from the target function f * . We show that the distance between f * and V D,n is approximately equal to the distance between f * and f * Π n , which is the optimal approximation in the space Π n . In the right part, V (m) D,n is a weighted average of the individual filtered hyperinterpolations V D j ,n obtained from multiple data sets sampled from the target function f * . Here again, the distance between f * and V (m) D,n is approximately equal to the distance between f * and its optimal approximation f * Π n in Π n

Preliminaries on Approximation on Manifolds
In this section, we discuss L p and Sobolev spaces of functions on manifolds, assumptions on the manifolds and embedding theorems to the space of continuous functions.
We start with a brief description of L p spaces and norms. Let M be a compact and smooth Riemannian manifold of dimension d ≥ 1 with smooth or empty boundary and Riemannian measure μ normalized to have the total volume μ(M) = 1. For 1 ≤ p < ∞, let L p (M) = L p (M, μ) be the complex-valued L p -function space with respect to the measure μ on M, endowed with the L p norm We will write f (x) L p (M),x = f L p (M) to indicate the variable for integration when necessary. For p = 2, L 2 (M) is a Hilbert space with inner product where g is the complex conjugate to g.

Diffusion Polynomial Space
Diffusion polynomials are a generalization of regular polynomials. We will use them to construct approximations of real-valued functions on manifolds. Let N := {1, 2, . . . } be the set of positive integers and let N 0 = N ∪ {0}. Let Δ be the Laplace-Beltrami operator on M, which has a sequence of eigenvalues {λ } ∈N and a corresponding sequence of orthonormal eigenfunctions {φ ∈ L 2 (M) | Δφ = −λ 2 φ , ∈ N}. We let λ 0 := 0 and φ 0 := 1. For n ∈ N 0 , the span Π n := span{φ |λ ≤ n} is called the diffusion polynomial space of degree n on M, and an element of Π n is called a diffusion polynomial of degree n. In the following, we will refer to diffusion polynomials simply as polynomials.
Let ρ(x, y) be the geodesic distance of points x and y induced by the Riemannian metric on M. For x ∈ M and β > α > 0, let B (x, α) := {y ∈ M|ρ(x, y) ≤ α} be the ball with center x and radius α, and let B (x, α, β) . We make the following assumptions for the measure μ and the eigenfunctions of Δ on M. The first is a standard assumption about the regularity of the measure on the manifold.

Assumption 1 (Volume of ball)
There exist positive constants c, c , c ≤ c, depending only upon the measure μ and the dimension d such that for all β > α > 0 and x ∈ M, Here c, c are the constants depending only on d.
The second is an assumption stating that the space of polynomials is closed under multiplication.
Assumption 2 (Product of eigenfunctions) For , ∈ N 0 , the product of eigenfunctions φ ∈ Π n , φ ∈ Π n for the Laplace-Beltrami operator Δ on M is a polynomial of degree n + n , i.e., φ φ ∈ Π n+n . Assumption 2 implies that the product P P of two polynomials P ∈ Π n and P ∈ Π n of degrees n ∈ N 0 and n ∈ N 0 , respectively, is a polynomial of degree n + n . Assumptions 1 and 2 are satisfied by typical manifolds, such as hypercubes [0, 1] d , unit spheres and balls in real or complex Euclidean coordinate spaces [10,19], flat tori T d , d ≥ 1, and Grassmannians [4,5], simplices in R d [41,48], with Lebesgue measures induced by the corresponding Riemannian metric, and also graphs equipped with an atomic measure on the graph nodes with the graph Laplacian defined as the difference of the identity matrix and the adjacency matrix [42].

Generalized Sobolev Spaces
We give a brief introduction to the Sobolev spaces on a Riemannian manifold M. The Fourier coefficients for f in L 1 (M) are For r > 0, the generalized Sobolev space W r p (M) may be defined as the set of all .
We let W 0 p (M) := L p (M). In the context of numerical analysis, we need to use the following Lemma 1, which is an embedding theorem of Sobolev space into the space of continuous functions on a manifold, see, e.g., [1,Section 2.7]. It guarantees that any function in the Sobolev space has a representation by a continuous function so that the numerical integration is valid and the quadrature rule can be applied.

Filtered Approximation on Manifolds
This section defines the filtered polynomial approximation on a compact Riemannian manifold M in terms of the eigenfunctions of the Laplace-Beltrami operator Δ on M. Given a target function f * ∈ L p (M) with 1 ≤ p ≤ ∞, the filtered polynomial approximation converges to functions in L p (M) as the degree n tends to infinity. Filter A real-valued continuous compactly supported function on R + is called a filter. Without loss of generality, we will only consider filters with support a subset of [0, 2]. In this paper, we focus on the following function H on R + as the filter.
and H ∈ C κ (R + ) for some κ ∈ N. Then, H is called a filter.
Definition 2 (Filtered kernel) A filtered kernel of degree n for n ∈ N on M with filter H is defined by Here λ and φ are eigenvalues and eigenfunctions of the Laplace-Beltrami operator on M.
For a kernel G : M × M → R and f ∈ L 1 (M), the convolution of f with G is defined as Definition 3 (Filtered approximation) We can define a filtered approximation V n on L 1 (M) as an integral operator with the filtered kernel K n,H (·, ·): for f ∈ L 1 (M) and x ∈ M, Note that for n = 0 this is just the integral of f . By (1) and (3), The following lemma, as given by Maggioni and Mhaskar [24,Theorem 4.1], shows that a filtered kernel is highly localized when the filter is sufficiently smooth.
where the constant c depends only on d, H and κ and ρ(x, y) is the geodesic distance between x and y.
By (4), if y is not close to x, and |K n (x, y)| decays to zero with rate n κ−d . It means given x, the kernel |K n (x, ·)| is concentrated on a small neighborhood of x, although it is supported on the whole manifold. This localization is essential to the boundedness of the filtered approximation operator. Lemma where the constant c depends only on d, p, H and κ.

Remark 1
For the case that M is a sphere, the above lemma for p = 1 was proved by Wang et al. [38] (see also [27] for κ ≥ d + 1); the case p > 1 can be obtained from the case p = 1 with the fact that K n ∈ Π d 2n and the Nikolskiî inequality for spherical polynomials [25].
Using the interpolation theorem with (3) gives This with Lemma 3 implies the following boundedness of the filtered approximation on L p (M). Polynomial space and best approximation Recall Π n := span{φ |λ ≤ n} is the (diffusion) polynomial space of degree n on manifold M.
The following theorem proves the convergence error for the filtered approximation of f ∈ L p (M).

Theorem 4 Let d ≥ 1, 1 ≤ p ≤ ∞ and M be a compact Riemannian manifold of dimension d. Let V n be the filtered approximation with filter H given by Definition
where the constant c depends only on d, H and κ.

Remark 2
For L p ([0, 1]), the case of filtered approximation with an appropriate filter reduces to a classic result of de la Vallée-Poussin approximation [20]. Stein [32] proved in a general context the convergence of de La Vallée-Poussin approximation to the target function. The sphere case of Theorem 4 was proved by Rustamov [28], Le Gia and Mhaskar [21], and Sloan [30].
The following lemma gives the convergence error of the best approximation for f ∈ W s p (M), see [24].
where the constant c depends only on d, p and r .
Theorem 4 and Lemma 4 imply the following convergence order for the filtered approximation of a smooth function on a compact Riemannian manifold.
where the constant c depends only on d, p, r , H and κ.
In Sects. 3 and 4, we will study the non-distributed and distributed filtered hyperinterpolation which use single and multiple servers to find a global estimator, respectively. For both non-distributed and distributed learning by filtered hyperinterpolation, we need to take account of the data type (noise or noiseless) and the quadrature point type (deterministic or random). There are in total 8 cases for which we have to treat separately.

Non-Distributed Filtered Hyperinterpolation on Manifolds
In this section, we study the non-distributed version of filtered hyperinterpolation (NDFH) on a manifold. We consider the cases when the data is either clean or noisy, and the input samples are either deterministic or random. It turns out that the NDFH for clean data achieves the optimal convergence order of the approximation error, while noise on the data would reduce the convergence order.
Filtered hyperinterpolation is a special type of regression, and the primary tool that we will use. Within this approach, as introduced in Definition 3, a target function f * is approximated by the filtered polynomial approximation Here, H is a filter for the eigenvalues λ to eigenfunctions φ , and ( f * ) are the Fourier coefficients. The Fourier coefficients cannot be computed in practice, because they would require to integrate the unknown target function. Instead, they are estimated from samples. This estimation is conducted via a quadrature formula, We rewrite (5) as After rearranging, our approximation takes the form which is a weighted sum of kernels centered at the data locations x i . In practice, the estimator of (6) is scaled by the observed values y i instead of f * (x i ).
In the following, we define the (non-distributed) filtered hyperinterpolation (approximation) on a compact Riemannian manifold M for a data set D. Besides the traditional deterministic quadrature rule, we also consider the filtered hyperinterpolation with random quadrature rule where the quadrature points are distributed with some probability measure on the manifold. We first introduce some notion about data and quadrature rule.
for noises i . Deterministic and random sampling In this paper, we consider two types of input samples depending on whether they are randomly sampled: the deterministic sampling and random sampling. The data D has random sampling if x i are randomly chosen with respect some probability measure on M. In contrast, D has deterministic sampling if the x i are fixed. Noisy and noiseless data We also distinguish data types by its data values y i . We say D is noiseless data or clean data if y i is equal to the function value of the associated (ideal) target function value f * (x i ) (that is, the noises i ≡ 0). We say D is noisy data if the noises i in (7) are nonzero.

Quadrature rule A set
is said to be a quadrature rule for numerical integration on M. We say Q D is a positive quadrature rule if all weights w i > 0, i = 1, . . . , N . In this paper, we only consider positive quadrature rules. For general Riemannian manifolds, the construction of quadrature formulas depends on the existence of designs on manifolds, which has been proved in [15]. See also [3]. On a flat torus, the weights are equal for a regular grid in [−π, π] d . On the real spheres, numerically constructive spherical designs with double precision (cubature with equal weights) can be found in [45]. However, in general, the weights do not have to be equal. For example, for a hypercube [0, 1] d , non-equal quadrature rule examples can be found in [6]; for spheres or spherical caps, quadrature rules with non-equal area weights examples can be found in [19]; for complex spheres [40]; Grassmannians [5], where the weights are not necessarily equal.
a positive quadrature rule on M and H be a filter in Definition 2 on R + . For n ∈ N, the non-distributed filtered hyperinterpolation (NDFH) for data D and quadrature rule Q D is be the noiseless data for the ideal target function f * and data D, then (8) becomes We call V D * ,n ( f * ) non-distributed filtered hyperinterpolation (NDFH) for clean data set or for quadrature rule Q D , for the function f * .

Remark 3
Non-distributed filtered hyperinterpolation on the sphere was studied by Sloan and Womersley [31].

Non-Distributed Filtered Hyperinterpolation for Clean Data
We first assume that we have a quadrature rule that has polynomial exactness. That is, the weighted sum by the quadrature rule can recover the integral for polynomials on manifolds. The non-distributed filtered hyperinterpolation with polynomial-exact quadrature rule can reach the same optimal convergence order as the filtered approximation in Sect. 2.3, and the convergence rate is optimal.
Let ∈ N 0 . A positive quadrature rule That the quadrature is exact for polynomials is a strong assumption, as the optimalorder number of points is O N d in typical examples of manifolds, see, e.g., [6,19].
The following lemma shows that the filtered hyperinterpolation V D,n with filter H given by Definition 1 reproduces polynomials of degree up to n if the associated quadrature rule Q D is exact for degree 3n − 1.

Lemma 5 Let n ∈ N 0 and M be a d-dimensional compact Riemannian manifold. Let
be a positive quadrature rule on M exact for polynomials of degree up to 3n − 1 and let V D * ,n be a non-distributed filtered hyperinterpolation on M for quadrature rule Q D with filter H given by Definition 1. Then, Theorem 6 (NDFH with clean data and deterministic samples) Let d ≥ 1, 1 ≤ p ≤ ∞ and n ≥ 1. Let M be a compact Riemannian manifold of dimension d. Let H be a filter given by Definition 1 with κ ≥ d + 1 and Q D be a positive quadrature rule exact for polynomials of degree up to 3n − 1. Then, for f * ∈ W r p (M) with r > d/ p, the NDFH for the quadrature rule Q D has the error upper bounded by where the constant c depends only on d, p, r , H and κ.
From the perspective of information-based complexity it is interesting to observe that if the target function f * is in the Sobolev space W s p (M), s > 0, the convergence rate is optimal in the sense of optimal recovery. This is due to the fact that on a real unit sphere when one uses optimal-order number of points N = O n d , the order n −s = N −s/d in (10) is optimal, as proved by Wang and Sloan [35] and Wang and Wang [36]. Theorem 6 can be viewed as the non-distributed filtered hyperinterpolation for clean data, where the estimator uses the whole data set in one machine.
We now introduce the (non-distributed) filtered hyperinterpolation for clean data with random sampling. We say a data set D has random sampling (with distribution ν) if the sampling points x i of D are independent and identically distributed (i.i.d.) random points with distribution ν on M. To construct the filtered hyperinterpolation for random sampling, we need the following lemma, which shows that there exist N quadrature weights given N i.i.d. random points x i such that the resulting quadrature rule is exact for polynomials for degree n with high probability. For 1 ≤ p ≤ ∞, let L p,ν (M) be L p space on manifold M with respect to probability measure ν.
for a positive absolute constant c. Then, for integer n satisfying N /n 2d > c for sufficiently large constant c, there exists a quadrature rule where C is a constant depending only on dimension d.

Definition 5 (Quadrature rule for random samples) We modify the weights of the set
in Lemma 6, as follows.

If the event
the quadrature rule for random samples on the manifold M for measure ν.
The following theorem gives the approximation error of the non-distributed filtered hyperinterpolation with clean data and random sampling for sufficiently smooth functions. Here, we want to obtain an estimated value of the expected error and take the expectation over the distribution of the data P(X )P(Y |X ).  (11). Given some τ , 0 < τ ≤ d, for cn d+τ ≤ |D * | ≤ c n 2d with two positive constants c, c , the filtered hyperinterpolation V D * ,n for clean data set D * with target function f * ∈ W r 2 (M), as given by (9) with the quadrature rule for random samples , where c 5 and c depend only on μ(M), c 1 , d, r , and the filter H and its smoothness κ, and C from Lemma 6.
Theorem 7 shows that the filtered hyperinterpolation with random sampling for clean data can achieve the same optimal convergence rate as the filtered hyperinterpolation with deterministic sampling. We give the proof of Theorem 7 in Sect. A.1.

Non-Distributed Filtered Hyperinterpolation for Noisy Data
In the following, we describe non-distributed filtered hyperinterpolation with deterministic or random sampling for noisy data. The data y i are the values of a function f * on M plus noise. Here, we assume that the noise has mean zero and is bounded. To be precise, we let The D satisfying (12) is then called noisy data set associated with f * . For real data, the f * is an unknown mapping from input to output. We study the performance of the non-distributed filtered hyperinterpolation for a noisy data set D whose data are stored in a sufficiently big machine.
We first consider the case where the locations of sampling points are fixed, which we call filtered hyperinterpolation with deterministic sampling. The kernel K n provides a smoothing method for the function f * using data D. As we shall see below, the approximation error of this filtered hyperinterpolation has the convergence rate depending on the smoothness of function f * . The following assumes that there exists a quadrature rule with N nodes and N "almost equal" weights which are exact for polynomials of degree approximately N 1/d .
Remark 4 For the sphere of any dimension, Assumption 8 always holds [26]. In order to construct the quadrature rule for general Riemannian manifolds, one needs to find weights that make the worst case error vanish. This corresponds to solving a particular equation The following theorem shows that the filtered hyperinterpolation V D,n can approximate f * well, provided that the support of the filtered kernel is appropriately tuned and Assumption 8 holds.
Theorem 9 (NDFH for noisy data and deterministic samples) Let d ≥ 2 and r > d/2. The sampling point set of the data set D satisfies Assumption 8. Then, for c 3 (13), the filtered hyperinterpolation V D,n for noisy data set D with target function f * ∈ W r 2 (M) satisfies where is a constant only depending upon r , d, M and the Sobolev norm of the target function f * , where M is the upper bound for the noise as given by (12), and c 1 , c 2 , c 3 , c 5 are constants depending only on d, r , the filter H and its smoothness r , where c 1 is the constant c in Lemma 3, and c 2 , c 3 are from Assumption 8.
Here, in contrast to Theorem 6, y contains noise. The expectation in (14) is with respect to the noise on y. The variance of the noise enters in C 1 .

Remark 5
Here the condition r > d/2 is the embedding condition such that any function in W r 2 (M) has a representation of a continuous function on M, which makes quadrature rule of filtered hyperinterpolation feasible for numerical computation.
Theorem 9 illustrates that if the scattered data Λ D has polynomial exactness, and the support of the filter η is appropriately chosen, then the filtered hyperinterpolation for noisy data set D can approximate sufficiently smooth target function f * on the manifold in high precision in probabilistic sense. By Györfi et al. [17], the rate |D| −2r /(2r +d) in (14) cannot be essentially improved in the scenario of (12). Theorem 9 thus provides a feasibility analysis of the filtered hyperinterpolation for manifold-structured data with random noise. Now, we introduce the (non-distributed) filtered hyperinterpolation for noisy data with random sampling.
where the x i are i.i.d. random points with distribution ν on M. The following theorem gives the approximation error of the non-distributed filtered hyperinterpolation for sufficiently smooth functions. Here, we want to get an estimated value of the expected error and take the expectation over the distribution P(X )P(Y |X ) of the data. For two sequences {a n } ∞ n=1 , {b n } ∞ n=1 , a n b n means that there exist constants c , c such that c b n ≤ a n ≤ cb n .
Theorem 10 (NDFH for noisy data and random samples) (11). For n |D| 1/(2r +d) , the filtered hyperinterpolation V D,n given by (8) with the quadrature rule for random samples , has the approximation error where M is the upper bound for the noise as given by (12), and c , c 5 , c 5 are constants depending only on μ(M), d, r, and filter H and its smoothness κ, and C from Lemma 6. Theorems 9 and 10 show that the filtered hyperinterpolation approximations with deterministic sampling and random sampling can achieve the same optimal convergence rate. We give the proofs of Theorems 9 and 10 in Sect. A.2.

Distributed Filtered Hyperinterpolation on Manifolds
In this section, we describe the distributed learning by filtered hyperinterpolation for clean data with deterministic and random sampling's. Distributed data sets We say a large data set D is distributively stored in m local servers if for j = 1, . . . , m, m ≥ 2, the jth server contains a subset D j of D, and there is no common data between any pair of servers, that is, D j ∩ D j = ∅ for j = j , and D = ∪ m j=1 D j . The data sets D 1 , . . . , D m are called distributed data sets of D. In this case, the filtered hyperinterpolation V D,n which needs access to the entire data set D is infeasible. Instead, in this section, we construct a distributed filtered hyperinterpolation for the distributed data sets {D j } m j=1 of D by the divide and conquer strategy [22].
be a data set on manifold M. The distributed filtered hyperinterpolation (DFH) for distributed data sets {D j } m j=1 of D is a synthesized estimator of local estimators V D j ,n , j = 1, 2, . . . , m, each of which is the filtered hyperinterpolation for noisy data D j : where for j = 1, . . . , m, the local estimator is a filtered hyperinterpolation on D j : For noiseless data sets The synthesis here is a process when the local estimators communicate to a central processor to produce the global estimator V (m) D,n . The weight in the sum of (15) for each local server is proportional to the amount of data used in the server.

Distributed Filtered hyperinterpolation for Clean Data
The synthesis here is a process when the local estimators communicate to a central processor to produce the global estimator V (m) D,n . Like the non-distributed case, we start with the case of deterministic sampling. The following theorem shows that the distributed filtered hyperinterpolation V (m) D,n has a similar approximation performance as the non-distributed V D,n when the number of local servers is not too large as compared with the amount of data.
Theorem 11 (DFH for clean and deterministic data) Let d ≥ 2 and 1 ≤ p ≤ ∞, r > d/ p, n ∈ N, m ≥ 2 and D * a clean data set satisfying (12). Let {D * j } m j=1 be m distributed data sets of D * satisfying min j=1,...,m |D * j | ≥ |D * | d 2r +d . Let H be a filter given by Definition 1 with κ ≥ d + 1. For j = 1, . . . , m, the data set D * j on the jth server satisfies that Q D * j is a positive quadrature rule exact for polynomials of degree up to 3n − 1.
with constant c 5 only dependent on d, r and the filter H and its smoothness κ.
Theorem 11 illustrates that the distributed filtered hyperinterpolation for the deterministic samples has a slightly slower approximation rate as the non-distributed case of Theorem 6, where the latter processes all the distributed data sets in a single server.
2r +d imply the existence of the quadrature rule exact for polynomials of degree n for the points set of each local server D * j .
The distributed filtered hyperinterpolation with random sampling is a weighted average of individual non-distributed filtered hyperinterpolations on local servers, where each weight is in proportion to the amount of the data used by the corresponding local server. To be precise: on manifold M for m ≤ 2, and their points are random samples. Denote each subset by , and we assign the weights to D j to make a quadrature rule for random samples that satisfy Definition 5 and Lemma 6. Then, the distributed filtered hyperinterpolation V (m) D,n for D with random samples is defined by Definition 6, where each V D j ,n uses the For random samples, we first consider the clean data case. Let V D * j ,n be the nondistributed filtered hyperinterpolation for clean data D * j with random sampling points. We define the global estimator V (m) D * ,n as (15). In the following theorem, we show that the approximation error for V (m) D * ,n on d-manifold converges at rate |D * | − 2r 2r +d where r is the smoothness of the target function. It has the same convergence rate as the distributed case for deterministic samples in Theorem 11.  where c , c depend only on μ(M), d, r, and the filter H and its smoothness r , and constant C from Lemma 6.
The proofs of Theorems 11 and 12 are deferred to Sect. A.3.

Distributed Filtered Hyperinterpolation for Noisy Data
In this subsection, we describe the distributed learning by filtered hyperinterpolation for noisy data with deterministic and random sampling. As shown in the following theorem, we prove that the distributed filtered hyperinterpolation V where M is the upper bound for the noise as given by (12), and c 1 , c 2 , c 3 , c 5 are constants depending only on d, r , the filter H and its smoothness r , where c 1 is the constant c in Lemma 3, and c 2 , c 3 are from Assumption 8.
The distributed filtered hyperinterpolation for deterministic sampling has the same order |D| − 2r 2r +d of the approximation error as compared to the non-distributed case in Theorem 9. Thus, appropriately distributing data to local servers, the divide-and-conquer strategy does not reduce the approximation capability of filtered hyperinterpolation. We will see that it is also true when the sampling is random.

Remark 7
Suppose each server takes the same number of data. With less than |D| When the data D is noisy with random sampling points, the distributed V (m) D,n in (15) has the same approximation rate as the non-distributed case in Theorem 10.
Theorem 14 (DFH for noisy data and random samples) Let d ≥ 2, r > d/2, m ≥ 2 and D a noisy data set satisfying (12). The sampling points are i.i.d. random points on M with distribution μ in (11). If the target function f * ∈ W r 2 (M), n |D| 1/(2r +d) and min j=1,...,m |D j | ≥ |D| d+τ 2r +d for some τ in (0, 2r ), then the distributed filtered hyperinterpolation V (m) D,n for noisy D, as given by Definition 7, has the approximation error , where c , c and c depend only on μ(M), d, r , and the filter H and its smoothness r , and constant C from Lemma 6.
From Theorems 11 to 14, the distributed filtered hyperinterpolations for all cases (in terms of data type and sampling type) always have a lower rate than the non-distributed cases of Theorems 6 and 7. But the noise added to the data also deteriorates the rate of the non-distributed cases of Theorems 9 and 10, which leads to the same rate as the distributed cases.

Examples and Numerical Evaluation
We illustrate the notions and filtered hyperinterpolation for single and multiple servers on the 2-d mathematical torus T 2 . The torus T 2 can be parameterized by the product of unit circles S 1 × S 1 and is equivalent to [−π, π] 2 . Denote L 2 (T 2 ) the L 2 space on T 2 with the Lebesgue measure. On the manifold T 2 , the Laplacian is the Laplace-Beltrami operator with eigenfunctions { 1 2π exp(i k · x)} k∈Z 2 of x ∈ T 2 and eigenvalues {|k| 2 } k∈Z 2 , where i := √ −1 is the imaginary unit, k ·x = k 1 x 1 +k 2 x 2 and |k| := k 2 1 + k 2 2 . Here k = (k 1 , k 2 ) and x = (x 1 , x 2 ). The space of polynomials of degree n is Π n := span{ 1 2π e i k·x : |k| ≤ n}. For 1 ≤ p ≤ ∞, let L p (T d ) be the L p space with respect to the normalized Lebesgue measure dx on T 2 .
For our illustration, we define the filter H by the piece-wise polynomial function with H (t) = 1 for 0 ≤ t ≤ 1; for t ∈ (1, 2); and H (t) = 0 for t ≥ 2. Then, H is in C 5 (R + ) and satisfies Definition 1. Figure 2 shows the plot of this filter. This particular filter has been used in previous works for the sphere, see [31,[37][38][39]. We observe that the filter, which is constant 1 over [0, 1], enables the filtered approximation and filtered hyperinterpolation of degree n (as given below) to reproduce polynomials with degree up to n on T 2 . The finite support [0, 2] of H makes the filtered hyperinterpolation a polynomial of degree up to 2n − 1. The middle polynomial over [1,2] which is sufficiently smooth at the two ends modifies the Fourier coefficients from degree n + 1 to 2n − 1 and makes the resulting filtered hyperinterpolation a near best approximator, as shown by Theorem 6.
With the filter (18), the filtered kernel on T 2 with filter H is, for n ∈ N and x, y ∈ T 2 , As the support of filter H is [0, 2], the summation over k is constrained to |k| ≤ 2n −1.
The filtered approximation for f ∈ L 2 (T 2 ) is then As corresponds to Definition 4, this is the ideal approximation of degree n, which is hard to compute as it would require integrating the unknown target function.
To construct a non-distributed filtered hyperinterpolation on T 2 , we consider N = 9n 2 0 points x j,l = (2 jπ/(3n 0 ), 2lπ/(3n 0 )). For these we can use the quadrature rule Q D = {(w j,l , x j,l ) : j, l = 0, 1, . . . , 3n 0 − 1} with N = 9n 2 0 equal weights w j,l ≡ (2π) 2 /N . The quadrature rule Q D is exact for polynomials of degree n. To satisfy the conditions of Theorems 9 and 13, we can let hyperparameter n 0 = n. In general, the quadrature weights need to be constructed depending on the location of the input points x j,l . We have been able to show the existence of such weights for randomly sampled inputs (see Lemma 6); however, the explicit construction in such cases is yet to be explored.
Consider a noisy data set D = {(x j,l , y j,l ) : j, l = 0, 1, . . . , n 0 − 1} with y j,l = f * (x j,l ) + j,l , j, l = 0, 1, . . . , n 0 − 1, and f * ∈ C(M). Here f * is the ideal (noiseless) target function. The non-distributed filtered hyperinterpolation of degree This is our construction to obtain an approximation. It corresponds to Definition 4. The summation index k runs over a ball of radius 2n − 1 due to the compact support of the filter H . Thus, V D,n (x) is fully discrete and computable. By Theorem 9, the approximation error of V D,n for f * ∈ W r 2 (T 2 ), r > 1, has convergence rate at least of order |D| −r /(r +1) as n 0 (controlling the number of data points) and n (controlling the degree of the approximation) increase. In particular, if f * is a basis element in { 1 2π e ix·k }, then r = ∞, since a polynomial is infinitely smooth. In practice, we use the real part as the approximation and discard the complex part. Note that the amount of data, 9n 2 , determines the degree of the polynomials. In this example, the diffusion polynomials of degree n are given by sums of eigenfunctions with momentum vector k = (k 1 , k 2 ) in the same grid defining the data.
The corresponding distributed filtered hyperinterpolation with m servers is where D j is the data set on the jth local server, for j = 1, . . . , m. This corresponds to Definition 6. By Theorem 13, the approximation error of the distributed strategy V (m)

D,n
with m servers for f * ∈ W r 2 (T 2 ), r > 1, is at least of order |D| −r /(r +1) provided the number of servers used satisfies m ≤ |D| r r +d/2 . When the points of the data set D are randomly distributed and satisfy the condition of Lemma 6, the non-distributed filtered hyperinterpolation remains the same as (19) with the points x j,l replaced by the set of random points Λ D . But the local estimator V D j ,n in the distributed filtered hyperinterpolation uses the modified weights w * j,l := (2π) 2 /N , if 0≤k,l≤n−1 |w j,l | 2 ≤ 2/m, 0, otherwise, in the place of (19). For our illustration, we use the Wendland function on the torus as the target function: where φ(u) is the one-dimensional Wendland function φ(u) := (1 − u) 8 + (32u 3 + 25u 2 + 8u + 1), and x c = (0, 0) is the center, see [44,47]. We show in the right part of Figure 2 the Wendland function in (20) which is in C 6 (T 2 ). We generate noisy data set by adding Gaussian white noise at a particular noise level to the values from the Wendland function. Figure 3 shows the L 2 squared errors of both training and generalization for the approximation by non-distributed filtered hyperinterpolation on noisy data from the Wendland function (20) on T 2 , with six levels of noise from 0 to 0.1. The degree n for the approximation is up to 40, and the sample size is N = (3n) 2 . The right part shows a few examples of the noisy training data, all at a noise level of 0.01, and the corresponding learned functions. For noiseless case, the training and generalization errors both converge to zero rapidly at a rate of approximately V D,n − f * L 2 (M) ∼ N −4 . This is consistent with the theoretical upper bound N −3 given in Theorem 6 where s ≥ 6 and d = 2. The slightly higher rate N −4 is due to that the φ(x) may have a higher smoothness. For noisy data, the convergence of the error stops at a particular degree. The convergence rate is higher when the noise level is smaller. The mean squared error on the training data converges to a value close to the square of the noise level, which indicates that the trained function is filtering out the noise. For both noisy and noiseless cases, the generalization error is slightly lower than the training error. Also, the result has consistent stability over repetitions in all cases. Figure 4 shows the L 2 squared errors for distributed filtered hyperinterpolation. We also generate the data from Wendland function. For this experiment, we partition the data set equally into m = 4 servers. The ith server computes a filtered hyperinterpolation on the data D i which are defined on interleaved grids of the form mod(x j,l + s i , 2π), where s i is a shift number between (0, 2π) and s i are distinct for different subsets D i . The quadrature rule Q D i utilizes equal weights as the nondistributed case. The distributed filtered hyperinterpolation combines the results from all servers, which has similar approximation behavior as the non-distributed case. If using noisy data in training, the approximation error has saturation after a particular degree; while with noiseless data, the error decays to zero all through the degree. We observe here that the generalization error has a more significant gap with training error as compared to the non-distributed case, which may be partly due to the distributed strategy (on multiple servers) are adopted. These experiments show consistent results as the theory in previous sections.

Discussion
Rates of convergence In Table 1, we compare the theoretical convergence rates of the non-distributed and distributed filtered hyperinterpolation in noiseless and noisy cases, as obtained in the previous sections. It shows that the filtered hyperinterpolation for clean data can achieve an optimal convergence rate N −r /d in both non-distributed and distributed cases and both deterministic and random sampling cases. For noisy data, the non-distributed filtered hyperinterpolation has a slightly lower approximation rate at N −r /(r +d/2) , r > d/2, which in the limiting case r → d/2 becomes the optimal rate N −r /d . The distributed strategy preserves the convergence rate N −r /(r +d/2) of the non-distributed filtered hyperinterpolation for noisy data, provided that the number of data N increases sufficiently fast with the number of servers, but the condition of for the data set D j on M, satisfying Definition 7. 4 end 5 for each server j = 1, . . . , m in parallel do 6 Consider the filtered kernel functions K n,H (·, with filter H and degree specification n as shown in Def. 2.

7
Compute the approximation V D j ,n (·) = |D j | i=1 w i y i K n,H (·, x i ) using the output data y i ∈ D j and the quadrature rule w i ∈ Q j as shown in Def. 4. the number of servers in the deterministic sampling case is weaker than the random sampling case. Implementation and complexity We already illustrated the computation of the method in Sect. 5. A summary of the implementation is shown in Algorithm 1. 1 In deterministic sampling, we start with some given input data and a suitable quadrature rule for those input values. In the random sampling case, we begin with the data, which in theory is only assumed sampled from some distribution, and then construct a suitable quadrature rule. There are various details to consider for the implementation. First, we need to choose a filter H , which should be sufficiently smooth depending on the dimension of the manifold M. The support of the filter will constrain the degree of the polynomials in the approximation. Second, we need a quadrature rule. Once a quadrature rule has been determined for the input data on the manifold, it can be applied to any output data. For important families of manifolds and configurations of points, quadrature rules are available from the literature. For instance, on the torus, cubes [12,34]; sphere: Gauss-Legendre, spherical design [2,11,19,45]; Graph: its nodes. The practical computation of quadrature rules for general types of data (or random input data) is an interesting problem in its own right, which has yet to be developed in more detail. Once a quadrature rule is available, the time complexity for Algorithm 1 is O max j=1,...,m |D| d 2r +d |D j | . If |D j | are all equal, the time complexity becomes O |D| d 2r +d +1 /m . Final remarks We have provided the first complete theoretical foundation for distributed learning on manifolds by filtered hyperinterpolation. One appealing aspect of filtered hyperinterpolation is that it comes with strong theoretical guarantees on the error, which apply to the population error or generalization error. Obtaining accurate bounds of this kind with neural networks is an active topic of research (which needs to incorporate not only the theoretical capacity of the neural network but also implicit regularization effects from the parameter initialization and optimization procedures). In filtered hyperinterpolation, once the data and the corresponding approximation degree are given, the approximating function is computed in closed form, meaning that we do not require parameter optimization. Also, filtered hyperinterpolation is a method that allows us to tune the model complexity directly in terms of the amount of available data in a principled way. As we observe in numerical experiments, the population error often is better than the training error. An interpretation is that this method imposes priors in terms of the polynomial degree and thus it is able to filter out noise. The method incorporates the geometry of the input space through the basis functions which are utilized to construct the approximations. Here, the basis functions are eigenfunctions of the Laplace-Beltrami operator on the manifold. It also contributes to the interpretability of the approximations, which live in polynomial spaces for which we have a good intuition. On the downside, to obtain the approximating function, the method relies on numerical integration techniques, in particular, quadrature rules, which is non-trivial in general to obtain. For general Riemannian manifolds, we can use the eigenvalues and eigenvectors of the discrete version of the Laplacian to approximate the Laplace-Beltrami operator, where the sampling points can estimate the discrete Laplacian, see, e.g., [7,13,33].
Funding Open Access funding enabled and organized by Projekt DEAL.

Conflict of interest The authors declare that they have no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

A Proofs
The appendices contain the proofs of the theorems in Sects. 3.1, 3.2, 4.1 and 4.2 in turn.
Proof (Theorem 4) Let P ∈ Π n . By the linearity of V n,H and Lemma 5, which, as P is an arbitrary polynomial in Π n , together with Theorem 3 gives thus completing the proof.
We go to prove Theorem 6, for which we need some lemmas as given below. The following theorem shows a Marcinkiewicz-Zygmund inequality for a quadrature rule on M.
Then, for all 1 ≤ p 1 < ∞ and > n, where c 1 depends only on d, p 0 and c 0 .
Remark 9 Dai [9] proved Lemma 7 when M is the unit sphere S d , d ≥ 1.
The proof of Lemma 7 relies on the following lemma of Filbir and Mhaskar [14], which shows that the sum of the weights, the corresponding nodes of which lie in the region B (x 0 , β, β + α), is bounded by a constant multiple of the measure of this region.

Lemma 8 Let d ≥ 1 and let M be a d-dimensional compact Riemannian manifold.
be a positive quadrature rule on M satisfying (22) for some 1 ≤ p 0 < ∞, c 0 > 0 and n ∈ N 0 . Then, for β ≥ 0, α ≥ 1/n and x 0 ∈ M, where the constant c depends only on d. Let Lemma 8 implies the following estimate for a quadrature rule.

Lemma 9 Let d ≥ 1 and let M be a d-dimensional compact Riemannian manifold. Let
be a quadrature rule on M satisfying (22) for some 1 ≤ p 0 < ∞, c 0 > 0 and n ∈ N 0 . Let A n (θ ) be given by (24). Then, for ≥ n, where the constant c depends only on d.
Proof Let x ∈ M. Since M is compact, M is bounded, i.e., there exists 0 < r < ∞ such that M ⊆ B (x, r ). Using Lemma 8, where the last inequality uses Assumption 1 and μ B (x, k/n, Proof (Lemma 7) For 1 ≤ p 1 < ∞, using (21) and Hölder's inequality gives, for P ∈ Π n and x ∈ M, Lemma 3 shows that the second integral of the filtered kernel on the right-hand side is bounded. This with (25) gives where the constant c depends only on d, p 1 , H and κ. Summing over quadrature nodes, we then obtain by Lemmas 9 and 2 that where the constant c depends only on d, p 1 , H and κ.
The proof of optimal-order error for filtered hyperinterpolation utilizes its decomposition by framelets on manifolds [35,38,43]. Given H ∈ C κ (R + ), κ ≥ 1, we define recursively the contributions of levels j ∈ N 0 for f ∈ L p (M) by The following lemma shows that U j ( f ) forms a decomposition of f ∈ L p (M), and it gives an upper bound of the where the constant c depends only on d, p, s, H and κ.
This with lim J →∞ E 2 J −1 ( f ) p = 0 gives (27). For f ∈ W s p (M) and j ∈ N, Theorem 5 with (26) gives where the constant c depends only on d, p, s, H and κ.
Proof (Theorem 6) For p = ∞, Lemma 5 with the linearity of V D * ,n shows that for q ∈ Π n , where using standard arguments, Taking the minimum over q ∈ Π n of the right-hand side of (29) with Lemma 4 gives Since the quadrature rule Q D is exact for degree 3n − 1, the condition of Lemma 7 is satisfied for p 0 = 2, then (23) holds for p 1 = 1. This with (30) and Lemma 3 gives Thus, We next consider for p ∈ [1, ∞). Given n ≥ 0, let m be the integer satisfying To bound the right-hand side of the last inequality in (31), we need the following estimate.
where the constant c depends only on d, p, H and κ. For p = 1 and j ≥ m, where the penultimate inequality uses Lemma 3 and the last uses Lemma 7 with p 1 = 1. For 1 < p < ∞ and j ≥ m, by Hölder's inequality, Using Lemma 7 with p 1 = 1 and p 1 = p and Lemma 3 then gives which proves (32) for p ∈ (1, ∞). It follows from (31) and (32) that for f ∈ W s p (M), where the second inequality uses (28), and where since n 2 m and s > d/ p, thus proving (10).
Proof (Theorem 7) For the quadrature rule for random samples , we define four events, as follows. Let Ω D * be the event such that is a quadrature rule exact for polynomials in Π d n associated with the measure ν and Ξ c D * the complement event of Ξ D * . Then, by Lemma 6, We write Under the event Ω c D * , by Definition 7, we obtain that V D * ,n = 0. Then, by (33), This together with (33) gives To bound A D * ,n,1 , we observe that when the event Ω D * ∩ Ξ D * takes place, is a set of positive weights for quadrature rule Q |D * |,n . We then obtain from Theorem 6 and f * ∈ W r 2 (M) with r > d/2 that On the other hand, under the event Ω D * ∩ Ξ c D * , by Cauchy-Schwarz inequality, This with Lemma 3 and (33) gives Then, with (35), (34) and cn d+τ ≤ |D * | ≤ c n 2d , τ ∈ (0, d], with where c 5 and c depend only on μ(M), c 1 , d, r , and the filter H and its smoothness κ, and C from Lemma 6.

Lemma 11
For n ∈ N 0 and 0 < p < q ≤ ∞, where the constant c depends only on d, p, q.
We need the following concentration inequality, Lemma 12, established by [46]. Let F be a subset of a metric space. For ε > 0, the covering number N (F, ε) for F is the minimal natural integer such that F can be covered by balls of radius ε, see [8,49].

Lemma 12
Let G be a set of functions on a product space X ×Y with Borel probability measure ρ. For every g ∈ G, if |g − Eg| ≤ B almost everywhere and E(g 2 ) ≤c(Eg) α for some B ≥ 0, 0 ≤ α ≤ 1 andc ≥ 0. Then, for any ε > 0, where the expectation Eg is taken on the product space X × Y with respect to ρ.
The third one is a covering number estimate for Banach space, as given by [50].

Lemma 13 Let B be a finite-dimensional Banach space. Let B R be the closed ball of radius R centered at origin given by B
Let X be a finite dimensional vector space endowed with norm · X , and Z ⊂ X * be a finite set. We say that Z is a norm generating set for X if the mapping T Z : X → R |Z| defined by T Z (x) = (z(x)) z∈Z is injective. We call T Z sampling operator. Let W := T Z (X ) be the range of T Z , then the injectivity of T Z implies that T −1 Z : W → X exists. Denote by · R |Z| the norm of R |Z| , and · R |Z| * the dual norm on R |Z| * for · R |Z| . We equip the space W with the induced norm, and let T −1 Z := T −1 Z W →X be the operator norm. In addition, let K + be the positive cone of R |Z| which is the set of all (r z ) z∈Z ∈ R |Z| such that r z ≥ 0. Then, the following lemma [26] holds.

Lemma 14
Let Z be a norm generating set for X , with T Z the corresponding sampling operator. If L ∈ X * with L X * ≤ A, then there exist positive numbers {a z } z∈Z , depending only on L such that for every x ∈ X , Proof (Lemma 6) For p = 1, 2, without loss of generality, we prove Lemma 6 for P n ∈ Π d n satisfying P n L p,ν (M) = A for some constant A > 0. For arbitrary P n ∈ Π d n with P n L p,ν (M) = A, it follows from (11) and Lemma 11 that We need to estimate the above covering number N (G p , ε) for p = 1, 2. To this end, we let G 1 := {P n ∈ Π d n : P n L p,ν (M) = A} and G 2 := {P n ∈ Π d 2n : P n L p,ν (M) = A} . By definition, N (G 1 , ε) ≤ N (G 1 , ε) and N (G 2 , ε) = N (G 2 , ε), where for p = 2, we have used |P n | 2 ∈ Π d 2n . Then, by Lemma 13, for p = 1, 2, we obtain the upper bound where we have used the estimate dim G p ≤ ( pn) d . Let ε = A p /4. As N /n 2d > c for sufficiently large constant c, with confidence 1 − 2 exp −C N /n d , there holds From this, we then obtain (42) holds with probability at least 1 − 2 exp −C N /n d . We now apply (42) with p = 2 and Lemma 14 to prove Lemma 6. In Lemma 14, we take X = Π d n , P n X = P n L 2,ν (M) , and Z the set of point evaluation The operator T Z is then the restriction map P n → P n | X N and It follows from (42) that with confidence at least We let L be the functional By Hölder's inequality, y X * ≤ 1. Lemma 14 then shows that there exists a set of real numbers holds with confidence at least 1 − 2 exp −C 3 N /n d , subject to Finally, we use the second assertion of Lemma 14 and (42) with p = 1 to prove the positivity of w i,n . Since 1 ∈ Π d n , we have v 0 := 1| X N is an interior point of K + . For P n ∈ Π d n , T Z P n = P n | X N is in W ∩ K + if and only if P n (x i ) ≥ 0 for all x i ∈ X N . For arbitrary P n satisfying P n (x i ) ≥ 0, x i ∈ X N , define ξ i (P n ) = P n (x i ). From Lemma 11 and (11), we obtain for i = 1, . . . , N , Applying Lemma 12 with B = 2C 1 c 4 n d P n L 1,ν (M) ,c =C 1 c 4 n d P n 2 L 1,ν (M) and α = 0 to the set {P n : P n ∈ Π d n , P n L 1,ν (M) = A}, using Lemma 13, we obtain for any ε > 0, Let ε = A/4 = 1 4 P n L 1,ν (M) . We then obtain that with confidence 1 − 2 exp −C N /n d , This and (42) imply that for P n which satisfies that P n (x i ) ≥ 0 for all x i ∈ X N , the inequality holds with confidence 1 − 4 exp −C N /n d with the constant C depending only oñ C 3 and c 4 , then, for arbitrary P n ∈ Π d n satisfying P n (x i ) ≥ 0, x i ∈ X N . Lemma 14 then implies w i,n ≥ 0 with confidence 1 − 4 exp −C N /n d , thus completing the proof of Theorem 6.
Proof (Theorem 10) For the quadrature rule for random samples in Definition 7, we define four events, as follows. Let Ω D be the event such that is a quadrature rule exact for polynomials in Π d n associated with the measure ν and Ξ c D the complement event of Ξ D . Then, by Lemma 6, To estimate the approximation error, we write Under the event Ω c D , by Definition 7, V D,n = 0. Then, by (43), Next, we estimate the first term of RHS of (44) when the event Ω D takes place. Let be the set of points of data D. By the independence between and Λ D and E{ i } = 0 for i = 1, . . . , |D|, This and (48) and (47) give Putting the above estimate and (45) into (44), we obtain Taking account of n |D| 1 2r +d and r > d/2, we then have Thus, 2r +d with C 3 a constant independent of |D|, thus completing the proof.

A.3 Proofs for Section 4.1
Proof (Theorem 11) By Definition 6, for f * ∈ W r 2 (M), By min j=1,...,m |D * j | ≥ |D * | 2d 2r +d and Assumption 8, there exists a quadrature rule for each local server D * j that is exact for polynomials of degree 3n − 1 for n satisfying To prove Theorem 13, we need the following modified version of [16,Proposition 4].

Lemma 15 For V (m)
D,n in Definition 6 with quadrature rule given by (13), there holds where V D * j ,n is given by (9).
Proof Due to (15) and m j=1 .

Taking expectations gives
Here, Then, .
We will use the following lemma to prove Theorems 12 and 14, which can be obtained similarly as Lemma 15.

Lemma 16 For the distributed filtered hyperinterpolation V
Proof (Theorem 12) For each j = 1, . . . , m, Then, by Lemma 16, we need only to estimate E V D * j ,n − f * 2 L 2 (M) . To this end, we can use the similar proof of the non-distributed case of Theorem 10 to prove its distributed counterpart V (m) D * ,n . For each j = 1, . . . , m, we let Ω D * j be the event such that the sum of the quadrature weights i=1 w 2 i,n,D * j ≤ 2/|D * j |, and Ω c D * j the complement of Ω D * j . We write By (49) with D = D j , These give By min 1≤ j≤m |D * j | ≥ |D * | d+τ 2r +d , n |D * | 1 2r +d and 2r > d, 0 < τ < 2r , Then, by Lemma 16 and (53), for constants c , c dependent only on μ(M), d, r , and the filter H and its smoothness r , and constant C from Lemma 6.

A.4 Proofs for Section 4.2
Proof (Theorem 13) By Lemma 15, we only need to estimate the bounds of E V D j ,n − f * 2 L 2 (M) and V D * j ,n − f * 2 L 2 (M) . By min j=1,...,m |D j | ≥ |D| 2d/(2r +d) and D j and Assumption 8, there exists a quadrature rule for each local server which is exact for polynomials of degree 3n − 1 for n satisfying c 3 6 |D| 1/(2r +d) ≤ n ≤ Here, C 6 = 2 2r +1 · 3 2r c 2 5 c −2r To bound E{V D j ,n } − f * 2 L 2 (M) , we use (46) and Jensen's inequality to obtain We now use the similar proof as Theorem 10 to prove the error bound of distributed filtered hyperinterpolation V (m) D,n . For each j = 1, . . . , m, we let Ω D j be the event such that the sum of the quadrature weights i=1 w 2 i,n,D j ≤ 2/|D j |, and Ω c D j the complement of Ω D j . We write By ( Using (57) and (59), and Lemma 16, we then obtain for c , c and c dependent only on μ(M), d, r , and the filter H and its smoothness r , and constant C from Lemma 6. Ball with center x and radius α in manifold ρ(x, y)

B Table of notations
Distance between points x, y ∈ M induced by Riemannian metric Distributed filtered hyperinterpolation with degree n and m servers for noiseless data D * in (16)