INTEGRATION BY RBF OVER THE SPHERE
|
|
- Thomasine Wilson
- 5 years ago
- Views:
Transcription
1 INTEGRATION BY RBF OER THE SPHERE ALISE SOMMARIA AND ROBERT S. WOMERSLEY Abstract. In this paper we consider numerical integration over the sphere by radial basis functions (RBF). After a brief introduction on RBF and spherical radial basis functions (SRBF), we show how to compute integrals of functions whose values are known at scattered data points. Numerical examples are given. Key words. Scattered data, numerical integration on spheres, radial basis functions. AMS subject classifications. 41A3, 65D3 1. Introduction. Radial Basis Functions (RBF) and Spherical Radial Basis Functions (SRBF) are well-established tools for approximating functions whose values are known on scattered data, respectively on R n and the unit sphere S n 1 {x R n : x 1} R n (the so-called unit n 1-sphere). Given the values of a function f at N centers X : {x i },...,N R n, one has to find scalars {λ i },...,N such that f(x) s f,x (x) : λ i φ(dist(x, x i )) (1.1) where φ : R + R is the given RBF/SRBF and dist is the Euclidean or the geodesic distance [4], [3] respectively. The function s f,x can be obtained, for instance, by interpolation of f at the centers X, independently of their geometry. Such a property is particularly appealing since it is not necessarily possible in the case of polynomial interpolation, where the distribution of the centers may imply the singularity of the corresponding andermonde matrix. By the Curtis-Mairhuber theorem [3, Theorem.3] there is no N dimensional space of continuous functions on the domain R n (having an interior point) that contains a unique interpolant for every set of centers X consisting of N data points. Consequently it is natural to choose a family of interpolants that depends on X, as in the case of RBF [7], [3]. A wide variety of RBF/SRBF have been proposed. Hardy introduced in 1971 the Multiquadrics (denoted by MQ) and Inverse Multiquadrics (IMQ) in conjunction with studies on topography, while Duchon in considered the Thin Plate Splines (TPS) and the Polyharmonic Plate Splines (see [3, 8.6]). The interest in RBF has increased during the last two decades, especially after the work in 198 by Franke on scattered data interpolation methods [1]. He showed that interpolation by multiquadrics and thin plate splines led to very good results, conjecturing also that the interpolation matrix that arises in the use of multiquadrics is invertible. This result was proved in 1986 by Micchelli [] who linked this problem to the work of Schoenberg on positive definite functions [3, Theorem 7.13]. In [] the analysis covers a wider class of RBF, including the Thin Plate Splines and the Gaussians (denoted by G). After these pioneering results, work blossomed in several directions: the discovery of RBF having compact support, fast methods for the School of Mathematics, University of New South Wales, Sydney NSW 5 (Australia) (alvise@maths.unsw.edu.au). School of Mathematics, University of New South Wales, Sydney NSW 5 (Australia) (R.Womersley@unsw.edu.au). 1
2 computation of RBF interpolants with thousand points (in spite of the high condition number of the interpolation matrix), error estimates and the properties of the interpolation matrix and choice of a good scaling parameter. Of course it is impossible to mention in detail all the research that has been done on this subject, and from this point of view we suggest to the interested reader the monographs [4], [5], [7], [3]. The problem of interpolation and approximation on scattered data of S n by SRBF has also been considered. Starting from the concept of positive definiteness on the sphere (also introduced by Schoenberg [3, Corollary 17.9]), several authors have studied the requirements under which the interpolation matrix is nonsingular, independently of the distribution of the centers X S n (see [6] and references therein). For examples of effective SRBF, see [1], [8], [11], [13]. If the values of the function f are known on scattered data over a square, some results concerning numerical cubature by RBF have been obtained in [3]. On the other hand, there has been an increasing interest in integrating functions on the unit sphere by cubature rules on special sets of points ([3], [8], [9]) as well as on scattered data on S n, caps or more generally compact subsets of caps (see [18], [19] also for various applications over S n ). The latter techniques allow for the construction of rules with positive weights and a certain polynomial degree of precision depending on the distribution of the points. In the framework of Galerkin methods, integrals of convolution type on S n have been studied in [11], [16] and []. The purpose of this paper is to show how the more often used RBF and SRBF are able to approximate integrals of the form If f(x)dµ(x) (1.) S where S R 3 is the unit -sphere, dµ(x) denotes the surface measure on S and f : S R a function whose properties will be specified later. We will show that once an RBF/SRBF approximation is at hand, the numerical cubature can be achieved easily and that the rule is optimal in the sense of Golomb-Weinberger [3, Section 13]. In section 6 numerical experiments illustrate the performance of the method and the size of the weights of the cubature rule.. Numerical integration on the sphere by RBF and SRBF. Let X {x 1,..., x N } S be a set of N centers and let f : S R be a continuous function defined on S. Suppose that you know an approximation s f,x which is a linear combination of radial basis functions and low order polynomials s f,x (x) : M 1 λ i φ( x x i /σ) + k k l γ k,l P k,l (x), x S. (.1) Here σ > is the scaling parameter [5], φ is a conditionally positive RBF of order M in R 3 [3], and P k,l (x) ξ k1 η k ζ k3, (.) where x (ξ, η, ζ) T and k k 1 + k + k 3, k 1, k, k 3, so P k,l is a homogenous polynomial of degree k in three variables. For some RBF, M so there is no low order polynomial component. Typically the function s f,x is obtained by interpolating f on the centers but other strategies may be used (for instance least squares methods [5], [8]).
3 We will restrict our attention to S to simplify the analysis but it is not difficult to extend the results to the n dimensional case. If for every i 1,..., N, the integral S φ( x x i /σ)dµ(x) (.3) exists, then integrating (.1) over S gives S f(x)dµ(x) + λ i φ( x x i /σ)dµ(x) S M 1 k k l Every rotation R : S S is an isometry, that is γ k,l S P k,l (x)dµ(x). (.4) φ( Rx Rx i /σ) φ( x x i /σ). (.5) Consider the rotation R that maps the center x i to the North Pole P (,, 1). As det(jac R 1 ) 1, the substitution x Rx gives φ( x x i /σ)dµ(x) S φ( Rx Rx i /σ)dµ(x) (.6) S φ( x P /σ) det(jac R 1 ) dµ(x ) S φ( x P /σ) dµ(x ). (.7) S As this holds for every center {x i },...,N, (.4) gives ( ) f(x)dµ(x) φ( x P /σ)dµ(x) + S S M 1 k λ i + γ k,l P k,l (x)dµ(x). (.8) S k The first integral on the right side can be easily evaluated for most functions φ. A point on the unit sphere S having cartesian coordinates (ξ, η, ζ) can also be described by spherical coordinates (ψ, θ) with l ξ cos ψ sin θ, η sin ψ sin θ, ζ cos θ. (.9) Here ψ is the azimuthal angle in the xy plane measured from the positive part of the x axis with ψ < π (the longitude) and θ the polar angle from the z axis with θ π (the colatitude). 3
4 RBF φ(r) φ( s) ds MQ (1 + r ) 1 3 IMQ (1 + r ) 1 (1 + s) 1 G exp( r ) exp( s) TPS r 1 log r 4 s log s 1 8 s Buhm. C r 4 log r 7 r r3 r s3 log s 3 18 s s + 1 s 6 Wend. C (1 r) s3/ + 1 s + s Wend. C (1 r) 4 8 +(4r + 1) 7 5 s s 5 5 s + s Wend. C 4 (1 r) 6 +(35r + 18r + 3) s s 4 18 s s 3 14 s + 3 s Table.1 Indefinite integrals of some RBF. Using this change of coordinates, remembering that the determinant of the Jacobian matrix of the transformation is simply sin θ and setting s : (1 cos θ)/σ, using x P (1 cos θ), (.1) we get I φ,σ φ ( x P /σ ) dµ(x) S π ( ) π (1 cos θ) φ sin θ dθ dψ 4/σ πσ σ φ ( s ) ds. (.11) For some RBF, an explicit expression for the indefinite integral φ ( s) ds are provided in Table I. For RBF with (compact) support [, 1] we used the notation r min{r, 1}, s min{s, 1} and s + : max(s, ). In particular for the Wendland [31] and Buhmann [3] RBF the upper limit of integration in the last integral in (.11) becomes min{4/σ, 1}. The approximation (.8) of the integral also requires I M M 1 k k l γ k,l S P k,l (x) dµ(x). (.1) 4
5 To this end, we observe that P k,l (x)dµ(x) S ξ k1 η k ζ k3 dµ(x) S π π π π ( π (cos ψ sin θ) k 1 (sin ψ sin θ) k (cos θ) k 3 dθ dψ cos k1 ψ sin k 1+k θ sin k ψ cos k3 θ dθ dψ ) ( π cos k 1 ψ sin k ψ dψ ) cos k 3 θ sin k1+k θ dθ. (.13) The Beta function β satisfies π/ sin m 1 θ cos m θ dθ 1 ( β m1 + 1, m ) + 1, (.14) and β(t 1, t ) β(t, t 1 ), so π sin m 1 θ cos m θ dθ π/ π/ π sin m 1 θ cos m θ dθ + sin m 1 θ cos m θ dθ π/ π/ sin m 1 θ cos m θ dθ + ( 1) m1 cos m1 θ sin m θ dθ 1 ( β m1 + 1, m ) ( ( 1)m1 m + 1 β, m ) ( 1 + ( 1) m 1 ) ( m1 + 1 β, m + 1 ). (.15) Similarly if m 1 and m are both even, then π ( sin m m θ cos m θ dθ β, m ) + 1, (.16) while π sin m 1 θ cos m θ dθ in all the other cases. Thus, from (.13), (.15) and (.16), one can easily compute S P k,l (x)dµ(x). The only RBF in Table.1 for which M > are the Multiquadrics (MQ) and the Thin Plate Splines (TPS), which are conditionally positive of order 1 and respectively. In these cases, when γ k,l are known, the computation of I M is particularly easy. For the Multiquadrics it is simply required to integrate the constant 1 on S, hence S P, (x)dµ(x) 4π. (.17) In the case of the Thin Plate Splines ξ dµ(x) η dµ(x) ζ dµ(x), (.18) S S S 5
6 SRBF Analytic Expression φ(α) Sph. Multiq. 1 + h hc Sph. Multiq. II (1 h ) (1 + h hc) 3/ Sph. Rec. Multiq. 1/ 1 + h hc Abel-Poisson Spline (1 h )/(1 + h hc) 3/ 1 Log. Spline log(1 + h/(1 h h h hc)) Sph. Spline I 1 a + b log((a + b)/b) 1 Sph. Spline II (4 3ab 3c + ((4 3c)c 1) log b b(3c 1) log(a + b)) Table. SRBF, with a c, b 1 c, c cos(α), α [, π] and h (, 1). which implies that S P, (x)dµ(x) 4π, For the integration by SRBF, suppose that f(x) S P 1,l (x)dµ(x), l 1,, 3. (.19) M 1 λ i φ(d(x, x i )) + k k l γ k,l Y k,l (x). (.) where Y k,l, l,..., k are spherical harmonics of degree k and d(x, y) arccos(x T y) is the geodesic distance (see [3], [9]). As d(rx, Ry) d(x, y) for any rotation R : S S then for any center x i X I φ,h φ(d(x, x i ))dµ(x) S φ( d(x, P) ) dµ(x). S (.1) Using spherical coordinates θ, ψ from (.9) with d(x, y) arccos(x T y) (so d(x, P) θ) we have I φ,h π π φ (θ) sin θ dθ dψ π π φ (θ) sin θ dθ. (.) Some of the more common SRBF are listed in Table., in which the argument α [, π] is the geodesic distance and h (, 1) is a spherical localization parameter. For more information about the Spherical Multiquadrics and Spherical Inverse Multiquadric see [8], for Spherical Multiquadrics II [1], while for Spherical Spline I, Spherical Spline II, Poisson Spline and Logarithmic Spline consider [1], [11] and [13]. Note that in [8] the argument of the SRBF is the geodesic distance d(x, y) while in [13] it is a function of c x T y cos d(x, y). As in the case of the RBF, we show in Table III that for most of the SRBF the indefinite integral φ (θ) sin θ dθ (.3) 6
7 SRBF Indefinite integrals φ(θ) sin θ dθ Sph. Multiq. (1 + h hc) 3/ /(3h) Sph. Multiq. II (1 h )(1 + h hc) 5/ /(5h) Sph. Rec. Multiq. 1 + h hc/h Abel-Poisson Spline (1 h )/(h 1 + h hc) Log. Spline Long Expression Sph. Spline I see (.4) Sph. Spline II see (.5) Table.3 Indefinite integrals of SRBF, with a c, b 1 c, c cos(θ), θ [, π] and h (, 1). is known explicitly. For the Log Spline and Spherical Spline I and II a short espression is not available. According to Maple and putting c cos(θ), a c, b 1 c, we get in the case of Spherical Spline I ln ( 1 + ) I φ,1 (c c 3) ln a while for Spherical Spline II ( a + a ) (ac c + a 1) + ln( c) + ln (.4) I φ,1 c + a 3 + 5ac 6 + ( ac + c + c ln(b) + c 3 ln(b)) c ln(b) 3 4 (c c + c 3 ) ln (a + 1 c) ln (a + ). (.5) To finish our analysis we simply observe that the integral of the spherical harmonics Y k,l over S, is always equal to, except for the case k l when, depending on the normalization of Y k,l, it is 4π. Summarizing, numerical integration based on approximation/interpolation by a linear combination of RBF/SRBF is particularly inexpensive when the coefficients {λ i }, {γ k,l } in (.1) or (.) are known. In practice, one has to compute only one definite integral where the indefinite integral is known explicitly. 3. Cubature and optimal recovery. The results in this section are an application to integration of a more general theory on linear functionals developed by Golomb-Weinberger [3, Section 13], in the framework of optimal recovery in Hilbert spaces. It is required that the space is endowed with a reproducing kernel, a tool that is well-known to researchers in approximation theory from 195 [3, Section 1]. Suppose that H is a real Hilbert space of functions f : R, where is a region containing at least one point. Furthermore let (, ) H and f H (f, f) 1/ H be respectively the inner product and the corresponding norm in H. A function K H : R is called a reproducing kernel for H if 1. K H (, x) H for all x ;. f(x) (f, K H (, x)) H for all f H and all x. Such a Hilbert space H is often referred to by the acronym RKHS. Among the properties of reproducing kernels we mention the following 1. K H (x, y) (K H (, y), K H (, y)) H for all x ;. K H (x, y) K H (y, x) for x, y ; 7
8 3. If f, f n H, n N are given such that f n converges to f in the Hilbert space norm then f n also converges pointwise to f. It is not difficult to show that every RKHS has exactly one reproducing kernel. Let us consider now the linear functional I : f f(x) dµ(x), f H (3.1) where µ(x) is a measure on. It is easy to prove that if KH (x, x) dµ(x) < +, (3.) then the functional I is also bounded because f(x) dµ(x) (f, K H (, x)) H dµ(x) f H K H (, x) H dµ(x) f H (KH (, x), K H (, x)) H dµ(x) f H KH (x, x) dµ(x). (3.3) By the Riesz theorem [3], there is a unique v I H such that As K H (, x) H, IK H (, x) gives If (f, v I ) H, for all f H. (3.4) K H (x, x)dµ(x ) (K H (, x), v I ) H v I (x), (3.5) v I K H (, x)dµ(x). (3.6) Given a fixed set of nodes X {x i },...,N, if the N by N matrix A defined by A i,j K H (x i, x j ), i, j 1,..., N is non singular, then there are N weights {w i,n },...,N (uniquely determined) such that K H (x, x j ) dµ(x) Note that by (3.7), the function w i,n K H (x i, x j ), j 1,..., N. (3.7) N s (int) I,X (x) : w i,n K H (x i, x) (3.8) is the interpolant of the representer v I on the set of nodes X. 8
9 At this point, if : span{k H (x i, )} then one can define the cubature rule Q (int) f : w i,n f(x i ) (f, s (int) I,X ) H (3.9) which can be interpreted as approximating (f, v I ) H by (f, s (int) I,X ) H. Observe that Q (int) f : N w i,n f(x i ) is meaningless when f L (), but L () is not a RKHS so what has been described above cannot be applied. In general, H is a RKHS if and only if the point evaluation functionals are continuous, i.e. δ x H for all x, where H is the dual space of H. This property suggests why (3.9) makes sense when H is a RKHS. The rule Q (int) is exact on, as for f N γ ik H (x i, ) f(x) dµ(x) γ i K H (x i, x) dµ(x) γ i j1 γ i N j1 w j,n K H (x i, x) dµ(x) w j,n K H (x i, x j ) γ i K H (x i, x j ) j1 w i,n f(x i ). (3.1) Notice that in (3.1) the weights {w i,n } are independent of the integrand f as they were obtained by interpolating the representer v I on the set of nodes X. Another approach consists in approximating If via the orthogonal projection P onto span{k H (x i, )} with respect to the inner product (, ) H of the RKHS H. If v I is the representer previously defined, then from (P f, g) H (f, P g) H for all f, g H (3.11) we can define another quadrature rule Q (proj) by Q (proj) f : (P f)(x) dµ(x) H (P f, v I ) H (P v I, f) H. (3.1) Maybe a little surprisingly, the quadrature rule Q (int) coincides with Q (proj) when applied to functions f H. To this purpose, we first notice that by the properties of P Q (proj) f (f, P v I ) H (P f, v I ) H, f H, (3.13) Q (int) f (f, s (int) I ) H (P f, s (int) I ) H, f H. (3.14) 9
10 Thus Q (proj) f Q (int) f on H, if and only if they coincide on, i.e. (K H (x i, ), v I ) (K H (x i, ), s (int) I ), i 1,..., N. (3.15) But (3.15) holds since, v I is the representer of I on H and Q (int) is exact on giving (K H (x i, ), v I ) K H (x i, ) dµ(x) (K H (x i, ), s (int) I ), i 1,..., N. (3.16) 4. An abstract setting. Denoting by H the dual space of H, a cubature rule Q w f : w i,n f(x i,n ) (4.1) with weights w {w i,n } is said optimal in the RKHS H if it minimizes the worst case error N E H (Q w ) : I w i,n δ xj,n (4.) H over all the rules with nodes {x i,n }. From (3.6) v I K H(, x)dµ(x) is the Riesz representer, so we get as in [1] E H (Q w ) sup f(x)dµ(x) w f H, f H 1 i,n f(x i,n ) ( ) sup f, K f H, f H 1 H (, x)dµ(x) w i,n K H (, x i,n ) H K H (, x)dµ(x) w i,n K H (, x i,n ) H v I w i,n K H (, x i,n ). (4.3) H This implies that the minimum of Q (opt) Q (opt) Q (proj) is attained by Q (int). (4.4) It is interesting at this point to see some examples of RKHS in which optimal recovery can be applied. It is well known that if ϕ : R is a symmetric (completely) positive definite kernel, then its associated native space K H N ϕ () is a Hilbert space with reproducing kernel ϕ (see, e.g. [3, p. 138]). By (3.7), we have that for a fixed set of nodes X (that are also centers of the radial basis function), the weights {w i,n } such that ϕ(x, x j ) dµ(x) w i,n ϕ(x i, x j ), j 1,..., N. (4.5) 1
11 This provides the cubature rule Q (int) f : w i,n f(x i ) (4.6) that is optimal with respect to the inner product of the RKHS N ϕ (). Notice that by (4.5) Q (int) is the rule of our numerical experiments since w i,n f(x i ) N λ j j1 w i,n w i,n ϕ(x i, x j ) N j1 λ j ϕ(x i, x j ) λ j j1 ϕ(x, x j ) dµ(x). The case of conditionally positive definite kernels is more complicated, but again it is possible to show that their native spaces N ϕ () are RKHS with respect to a known inner product dependent on ϕ (see, e.g. [3], p.146). For a detailed error analysis see [3] where upper bounds of the cubature error have been achieved in terms of mesh-norm and separation distance. 5. Numerical implementation. The results this section are probably known to experts in RBF, but in our opinion they are still useful for a general audience. In our numerical implementation, the coefficients {λ i }, {γ k,l } have been computed by interpolating the values {f i } at the centers {x i },...,N : R n. In the case of conditionally positive RBF of degree M one has to solve the linear system M 1 λ i φ( x j x i /σ) + k where x j S, j 1,..., N, plus the additional conditions k l γ k,l P k,l (x j ) f(x j ), j 1,..., N (5.1) λ i P k,l (x i ), k,..., M 1, l,..., k. (5.) The linear system defined by (5.1), (5.), can be rewritten as [ ] [ ] [ ] A Q λ f Q T γ (5.3) where λ (λ 1,..., λ N ) T, γ (γ,, γ 1,, γ 1,1, γ 1,,..., γ M 1,(M 1) ) T A i,j φ( x i x j /σ), i, j 1,..., N (5.4) Q i,j P k,l (x i ), i 1,..., N, k,..., M, l,..., k, j k + l + 1 (5.5) and solved by an appropriate numerical method [5]. The case of the SRBF is similar, with A i,j φ(d(x i, x j )), i, j 1,..., N (5.6) Q i,j Y k,l (x i ), i 1,..., N, j k + l + 1. (5.7) 11
12 where d(x, y) is the geodesic distance and Y k,l, k,..., M 1, l,..., k the spherical harmonics. By the uncertainty relation [6] the attainable error and the condition number κ( ) of the matrix [ A Q Q T ] (5.8) cannot both be kept small. As result, a large number of centres N in (.1) can lead to a better approximation but also to a very ill-conditioned interpolation matrix. To mitigate this problem several techniques have been proposed [], [5], [3, p. 53]. One of the more common approaches is due to Floater and Iske [9], which is used here for S although more general domains R n can be treated. The condition number of depends on the separation distance (including the factor means this is the packing radius on the sphere) ρ X 1 min x i X min x j X x j x i x i x j, (5.9) while the available error estimates (see [5], [6], [3]) depend on the mesh norm (in some papers also called fill-in distance or the covering radius on the sphere) h X, sup x min x x i. (5.1),...,N The ideal situation is to have a large separation distance ρ X and a small mesh norm h X,. A common criterion is the mesh ratio q X, : h X, ρ X 1, (5.11) which one wants to keep as small as possible. The given set of centers X {x i } is often decomposed into a hierarchy of nested subsets X 1 X... X L 1 X L X (5.1) in which, roughly speaking, the ratio q Xl, is small and at each level the interpolation problem has a small condition number. To reach this target, in [9], [14] different thinning algorithms have been proposed and analysed in connection to a variant of the so called k-center problem. Once the subsets {X l } l1,...,l are at hand, we can interpolate the data f Xl {f i,l } relative to the centers X l {x i,l } by the Multilevel Scheme described in [9], [14] and analysed in [1]. To be more precise, starting with l 1, at the l-th level one matches the error function f (s s l 1 ) (5.13) by a RBF (or SRBF) s l, as described respectively in (.1) and (.). 1
13 Denoting by g Xl the restriction of the function g : R to the set X l we can rewrite the sequence of interpolation problems (5.13) as s 1 X1 f X1 s X (f s 1 ) X... (5.14) L s L 1 XL 1 (f s k ) XL 1 s L XL As X X L it follows that that k1 L 1 (f s k ) XL. k1 (s s L ) X f X (5.15) i.e. s f,x s s L. Note that in general at level l of the scheme (5.14) we can choose the RBF scaling parameter σ l (or the spherical localization paramater h) as well as which RBF/SRBF to use. For a wide class of domains, Iske [15] suggests using at the first level globally supported Thin Plates Splines, and in the next levels compactly supported RBF with scaling parameters σ l monotonically decreasing. A procedure for computing a good σ l is also given. The multilevel interpolation method (5.14) can be easily used for our numerical integration purposes. This is a consequence of the fact that from (5.15) we can assume that hence in the domain R n f(x) dx f s X s XL (5.16) s 1 (x) dx s L (x) dx. (5.17) When S, one can of course use the technique described in the previous section for computing the L terms on the right hand side of (5.17). In our Matlab program, we have implemented these thinning algorithms and multilevel scheme to provide a good approximation of the function f in the centers X. The set X is a portion of the data from the orbit of a MAGSAT/NASA satellite and decomposed into subsets X l as established respectively by the separation distance based thinning algorithm [9]. As alternative one can use the progressive algorithm proposed and analysed in [14]. A sketch of our code is hence the following: 1. Set the value of the integral I equal to ;. Compute by thinning a hierarchical decomposition {X l } l1,...,l of the initial set of centers X; 3. At level l, once it has been decided the RBF/SRBF to use (and its parameters), compute the function s l by solving the linear system required by (5.14); 4. Compute I l S s l and put I I + I l ; 5. Iterate steps 3. and 4. of this algorithm until a stopping criterion is satisfied or l L. 13
14 First example: the Franke function f 1 RBF σ e e e I κ( ) 1 MQ. 6-E4 1-E 1-E6 +E7 +E7 G.5 7-E4 9-E3 7-E5 8+E7 +E7 IMQ.5 8-E4 1-E -E5 5+E5 +E7 Wend. C E4 1-E 8-E6 6+E4 +E7 TPS 1. 7-E4 1-E 9-E4 3+E5 +E7 SRBF h e e e I κ( ) 1 Abel-Poisson spline E4 9-E3 5-E5 1+E6 1+E3 Log. spline.85 7-E4 1-E 5-E6 +E6 +E3 Sph. Rec. Multiq E4 9-E3 -E5 1+E6 5+E Table 6.1 Numerical results for the Franke function f Numerical examples. We now show the behaviour of these cubature rules when applied to two well-known test functions restricted to the unit sphere S. Before passing to the examples, we briefly describe some technical details. The initial set of centers X T consists of 1 points obtained by thinning a selected portion of MAGSAT data X M [17]. In practice, we use a subset X containing only the first 1 points of X T that already provides a good compromise between separation distance ρ X and mesh norm h X,S. In our specific case h X,S.11, ρ X.37 and q X,S h X,S /ρ X 3.. We perform our cubature tests on two Franke functions [1] adapted by Renka to the three dimensional case [4], namely f 1 (x, y, z).75 exp ( (9x ) /4 (9y ) /4 (9z ) /4 ) and +.75 exp ( (9x + 1) /49 (9y + 1)/1 (9z + 1)/1 ) +.5 exp ( (9x 7) /4 (9y 3) /4 (9z 5) /4 ). exp ( (9x 4) (9y 7) (9z 5) ) (6.1) f (x, y, z) (1 + tanh ( 9x 9y + 9z)) /9. (6.) The value of If 1 computed by Maple to significant digits is If 1 f 1 (x) dµ(x) , (6.3) S while If f (x) dµ(x) 4π S 9. (6.4) We estimate the interpolation errors as follows. First we generate 3 random points X R with a uniform probability distribution and evaluate the RBF/SRBF interpolant s X as well as the test functions f 1, f, on X R and then we compute the relative errors e (f) : s X XR f XR f XR, e (f) : s X XR f XR f XR. (6.5) 14
15 Second example: the Franke function f RBF σ e e e I κ( ) 1 MQ E1 -E 1-E4 +E15 1+E1 G.5 -E1 3-E 8-E6 3+E1 +E1 IMQ.95 -E1 -E 9-E5 5+E15 4+E1 Wend. C E 7-E3 5-E4 1+E5 +E3 TPS 1. 9-E 9-E3 4-E4 3+E5 4+E SRBF h e e e I κ( ) 1 Abel-Poisson spline.35 -E -E1 5-E5 5+E15 5+E1 Log. spline.45 -E 1-E -E4 3+E15 3+E1 Sph. Rec. Multiq..4 -E -E 9-E5 5+E15 5+E1 Table 6. Numerical results for the Franke function f. The relative quadrature error is e I (f) f(x) dµ(x) Q (int) f S. (6.6) f(x) dµ(x) S where Q (int) f is one of the rules introduced in the previous sections. The norm condition number of the interpolation matrix is denoted by κ( ). alues of κ( ) larger than one over the machine precision, that is larger than 1 15 must be treated with care. In Tables 6.1 and 6., we consider only nearly optimal values of the scaling/localization parameters without showing what happens for general σ (or h). As is well-known for interpolation by RBF [5], a bad choice of σ may produce poor numerical results. For example, in the case of the function f 1 and σ, radial basis functions MQ, G and IMQ provide results whose errors e, e, e I are bigger than 1 1. Numerical experiments show that in general the TPS and the RBF with compact support (e.g. Wendland and the Buhmann) are more robust, since they provide good results even for non optimal σ (typically in parallel with not too high condition numbers). Table 6. also reports the values of 1, a fundamental parameter for analysing the sensitivity to perturbations of the RBF/SRBF process (see [3]). In Tables 6.3, we report the performance of the Multilevel scheme (5.14) when it is used to approximate If by the strategy described in (5.17). To this purpose, we consider interpolants s l (with l 1,... L) that are linear combinations of the (possibly scaled) Wendland function φ(r) (1 r) 4 +(4r +1). Furthermore, we take as centers the first, 4, 6, 8 and 1 points of the thinned set X, while for the sparsity of the interpolation matrix, we compute the percentage of nonzero entries. In particular we illustrate the differences between the stationary case, in which the ratio between the mesh-norm h X,S and scaling the parameter σ is constant (in our tests we set σ 4 h X,S ), and the non-stationary case in which we have choosen σ 1. The numerical results show the main advantages of the multilevel scheme in the stationary case not only for interpolation but also for cubature. With the exception of the first level, the linear systems that arise in (5.3) are sparser and also have better condition numbers, without affecting after these improvements the cubature error e I. Finally, in Table 6.4 we give information about the cubature weights w. At this 15
16 Multilevel method and the Franke function f Points σ κ( ) Sparsity e I 1 3+E 5% 8-E E3 5% 3-E E3 5% -E E3 5% 1-E E4 5% 6-E4 Points σ κ( ) Sparsity e I E 8% 8-E E 13% 3-E E 9% -E E 7% 1-E E 6% 5-E4 Table 6.3 Stationary and non-stationary multilevel method applied to f. stage we need to be careful, especially in the case of RBF/SRBF conditionally positive of order M >. We will use the bold face to represent vectors. With an obvious notation, the (possibly augmented) linear system (5.3) can be rewritten as A (aug) c (aug) f (aug). (6.7) On the other hand, knowing the weights w and denoting by (, ) R N in R N, one can write the scalar product If (w, f) R N. (6.8) Now we have to compute the weights w. If I (aug) is the vector ( ) I (aug) Iφ,σ e N b (6.9) where e N (1,..., 1) R N and b is the vector whose entries are the integrals of the polynomials required by the conditionally positive RBF/SRBF, then for a certain w (aug) R N+M we have If (c (aug), I (aug) ) R N+M (w(aug), f (aug) ) R N+M. (6.1) Consequently, by (6.7) and the fact that (A (aug) ) 1 is symmetric, we get (w (aug), f (aug) ) R N+M (c(aug), I (aug) ) R N+M ((A (aug) ) 1 f (aug), I (aug) ) R N+M (f (aug), (A (aug) ) 1 I (aug) ) R N+M. (6.11) Since (6.11) holds for any f (aug), we have w (aug) (A (aug) ) 1 I (aug). (6.1) 16
17 Weights RBF Min Max % (wi ) + (wi ) MQ -E 5-E 87% G -E 5-E 88% IMQ 3-E 6-E 86% Wend. C +1-E3 3-E 1% 1.6 TPS +4-E3 3-E 1% 1.6 SRBF Min Max % (wi ) + (wi ) Abel-Poisson spline 4-E 7-E 84% Log. spline 3-E 5-E 87% Sph. Rec. Multiq. 3-E 6-E 85% Table 6.4 Weights of the RBF/SRBF used in the numerical tests on the Franke function f. As the last M components of f (aug) are zero, we deduce that only the first N terms of w (aug) are relevant for the cubature. Thus the vector w of cubature weights for use in (6.8) is just the first N elements of the solution w (aug) to A (aug) w (aug) I (aug). The same results can alternatively be established working directly with the Riesz representer v I of the linear functional I. In Table 6.4, Min, Max are respectively the minimum and maximum weights, % is the percentage of positive terms in the vector w, (w i ) + and (w i ) are respectively the sum of positive and negative weights. One common stability criterion, (wi ) + (w i ), is satisfied by all the RBF/SRBF. In particular the percentage of positive weights is particularly high for the Wendland C and the Thin Plate Splines. 7. Acknowledgements. The first author is grateful to Dr. K. Hesse, Dr. T. Legia, Prof. I. Sloan, Prof. M. ianello and Prof. H. Wendland for useful discussions. Financial support of the Australian Research Council, Centre of Excellence for Mathematics and Statistics of Complex Systems is gratefully acknowledged. REFERENCES [1] B.J.C. Baxter and S. Hubbert, Radial basis functions for the sphere, Recent progress in multivariate approximation (Witten-Bommerholz, ), pp [] R. K. Beatson, W. A. Light, and S. Billings, Fast solution of the radial basis function interpolation equations: Domain decomposition methods, SIAM J. Sci. Comput., vol. (), no. 5, pp [3] M.D. Buhmann, A new class of radial basis functions with compact support, Math. Comp. 7, (1), no. 33, [4] M.D. Buhmann, Radial basis functions, Acta Numer., 9 (), pp [5] M.D. Buhmann, Radial basis functions: theory and implementations, Cambridge Monographs on Applied and Computational Mathematics, 1 (3). [6] D. Chen,.A. Menegatto and X. Sun, A necessary and sufficient condition for strictly positive definite functions on spheres, Proc. Amer. Math. Society, 131, (3), no. 9, pp [7] G.E. Fasshauer, Meshfree Methods, Handbook of Theoretical and Computational Nanotechnology, M. Rieth and W. Schommers (eds.), American Scientific Publishers, 5. [8] G.E. Fasshauer and L.L. Schumaker, Scattered Data Fitting on the Sphere, Mathematical Methods for Curves and Surfaces II, M. Daehlen, T. Lyche, and L. L. Schumaker (eds.), anderbilt University Press, 1998, pp
18 [9] M.S. Floater and A. Iske, Multistep scattered data interpolation using compactly supported radial basis functions, J. Comput. Appl. Math., 73 (1996), pp [1] R. Franke, Scattered data interpolation: tests of some methods, Math. Comp., 48 (198), pp [11] W. Freeden, T. Gervens, and M. Schreiner, Constructive Approximation on the Sphere (With Applications to Geomathematics), Clarendon Press, [1] K. Hesse and I.H. Sloan, Worst-Case errors in a Sobolev space setting for cubature over the sphere S, Bull. Aust. Math. Soc., 71 (5), pp [13] S. Hubbert, Radial Basis Function Interpolation on the Sphere, Ph.D. dissertation, (). [14] A. Iske, Multiresolution Methods in Scattered Data Modelling, Lecture Notes in Computational Science and Engineering, ol. 37, Springer-erlag, Heidelberg, 4. [15] A. Iske, Hierarchical Scattered Data Filtering for Multilevel Interpolation Schemes, Mathematical Methods for Curves and Surfaces: Oslo, T. Lyche and L. L. Schumaker (eds.), anderbilt University Press, Nashville, 1, pp [16] Q.T. Le Gia, Galerkin approximation for elliptic PDEs on spheres, J. Approx. Theory, 13 (4), pp [17] MAGSAT, MAGSAT-NASA/United States Geological Survey (USGS) effort to measure nearearth magnetic fields on a global basis. [18] H.N. Mhaskar, F.J. Narcowich and J.D. Ward, Spherical Marcinkiewicz-Zygmund inequalities and positive quadrature. Math. Comp., 7 (1), pp [19] H.N. Mhaskar, Local quadrature formulas on the sphere. J. Complexity, (4), no. 5, pp [] C.A. Micchelli, Interpolation of scattered data: distance matrices and conditionally positive definite functions, Constr. Approx., (1986), 11-. [1] F.J. Narcowich, R. Schaback and J.D. Ward, Multilevel interpolation and approximation. Appl. Comput. Harmon. Anal., 7 (1999), no. 3, pp [] F.J. Narcowich, J.D. Ward, Nonstationary wavelets on the m-sphere for scattered data. Appl. Comput. Harmon. Anal., 3 (1996), no. 4, pp [3] M. Reimer, Multivariate polynomial approximation. International Series of Numerical Mathematics, 144. Birkhuser erlag, Basel, 3. [4] R. J. Renka, Multivariate interpolation of large sets of scattered data., ACM Trans. Math. Software, 14 (1988), no., pp [5] S. Rippa, An algorithm for selecting a good value for the parameter c in a radial function interpolation, Adv. in Comput. Math., 11 (1999), pp [6] R. Schaback, Error estimates and condition numbers for radial basis function interpolation, Adv. in Comput. Math., 3 (1995), pp [7] R. Schaback, Lower bounds for norms of inverses of interpolation matrices for radial basis functions, J. Approx. Theory, 79 (1994), pp [8] N.J.A. Sloane, Tables of Sphere Packings a nd Spherical Codes, IEEE Trans. Information Theory, IT-7 (1981), pp , njas/. [9] I.H. Sloan and R. Womersley, How good can polynomial interpolation on the sphere be?, Adv. Comput. Math., 14 (1), no. 3, pp [3] A. Sommariva and M. ianello, Numerical cubature on scattered data by Radial Basis Functions, accepted for the publication by Computing. [31] H. Wendland, Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree, Adv. in Comput. Math., 4 (1995), pp [3] H. Wendland, Scattered Data Approximation, Cambridge Monographs on Applied and Computational Mathematics, 5. 18
MATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 9: Conditionally Positive Definite Radial Functions Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH
More informationInterpolation by Basis Functions of Different Scales and Shapes
Interpolation by Basis Functions of Different Scales and Shapes M. Bozzini, L. Lenarduzzi, M. Rossini and R. Schaback Abstract Under very mild additional assumptions, translates of conditionally positive
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 33: Adaptive Iteration Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 33 1 Outline 1 A
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 33: Adaptive Iteration Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 33 1 Outline 1 A
More informationNumerical cubature on scattered data by radial basis functions
Numerical cubature on scattered data by radial basis functions A. Sommariva, Sydney, and M. Vianello, Padua September 5, 25 Abstract We study cubature formulas on relatively small scattered samples in
More informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
More informationApproximation by Conditionally Positive Definite Functions with Finitely Many Centers
Approximation by Conditionally Positive Definite Functions with Finitely Many Centers Jungho Yoon Abstract. The theory of interpolation by using conditionally positive definite function provides optimal
More informationApplied and Computational Harmonic Analysis
Appl Comput Harmon Anal 32 (2012) 401 412 Contents lists available at ScienceDirect Applied and Computational Harmonic Analysis wwwelseviercom/locate/acha Multiscale approximation for functions in arbitrary
More informationDIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION
Meshless Methods in Science and Engineering - An International Conference Porto, 22 DIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION Robert Schaback Institut für Numerische und Angewandte Mathematik (NAM)
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 5: Completely Monotone and Multiply Monotone Functions Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH
More informationKernel B Splines and Interpolation
Kernel B Splines and Interpolation M. Bozzini, L. Lenarduzzi and R. Schaback February 6, 5 Abstract This paper applies divided differences to conditionally positive definite kernels in order to generate
More informationOverlapping additive Schwarz preconditioners for interpolation on the unit sphere by spherical radial basis functions
Overlapping additive Schwarz preconditioners for interpolation on the unit sphere by spherical radial basis functions Q. T. Le Gia T. Tran I. H. Sloan October, 008 Abstract We present an overlapping domain
More informationStudia Geophysica et Geodaetica. Meshless BEM and overlapping Schwarz preconditioners for exterior problems on spheroids
Meshless BEM and overlapping Schwarz preconditioners for exterior problems on spheroids Journal: Studia Geophysica et Geodaetica Manuscript ID: SGEG-0-00 Manuscript Type: Original Article Date Submitted
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 34: Improving the Condition Number of the Interpolation Matrix Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationAN ELEMENTARY PROOF OF THE OPTIMAL RECOVERY OF THE THIN PLATE SPLINE RADIAL BASIS FUNCTION
J. KSIAM Vol.19, No.4, 409 416, 2015 http://dx.doi.org/10.12941/jksiam.2015.19.409 AN ELEMENTARY PROOF OF THE OPTIMAL RECOVERY OF THE THIN PLATE SPLINE RADIAL BASIS FUNCTION MORAN KIM 1 AND CHOHONG MIN
More informationSpherical designs with close to the minimal number of points 1
Spherical designs with close to the minimal number of points Robert S. Womersley School of Mathematics and Statistics, University of New South Wales, Sydney, NSW, 05, Australia Abstract In this paper we
More informationSolving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels
Solving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels Y.C. Hon and R. Schaback April 9, Abstract This paper solves the Laplace equation u = on domains Ω R 3 by meshless collocation
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 2: Radial Basis Function Interpolation in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590
More informationD. Shepard, Shepard functions, late 1960s (application, surface modelling)
Chapter 1 Introduction 1.1 History and Outline Originally, the motivation for the basic meshfree approximation methods (radial basis functions and moving least squares methods) came from applications in
More informationMultiscale RBF collocation for solving PDEs on spheres
Multiscale RBF collocation for solving PDEs on spheres Q. T. Le Gia I. H. Sloan H. Wendland Received: date / Accepted: date Abstract In this paper, we discuss multiscale radial basis function collocation
More informationMultivariate Interpolation with Increasingly Flat Radial Basis Functions of Finite Smoothness
Multivariate Interpolation with Increasingly Flat Radial Basis Functions of Finite Smoothness Guohui Song John Riddle Gregory E. Fasshauer Fred J. Hickernell Abstract In this paper, we consider multivariate
More informationRKHS, Mercer s theorem, Unbounded domains, Frames and Wavelets Class 22, 2004 Tomaso Poggio and Sayan Mukherjee
RKHS, Mercer s theorem, Unbounded domains, Frames and Wavelets 9.520 Class 22, 2004 Tomaso Poggio and Sayan Mukherjee About this class Goal To introduce an alternate perspective of RKHS via integral operators
More informationNumerical solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method
Numerical solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method Y. Azari Keywords: Local RBF-based finite difference (LRBF-FD), Global RBF collocation, sine-gordon
More informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 7: Conditionally Positive Definite Functions Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 14: The Power Function and Native Space Error Estimates Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationA new stable basis for radial basis function interpolation
A new stable basis for radial basis function interpolation Stefano De Marchi and Gabriele Santin Department of Mathematics University of Padua (Italy) Abstract It is well-known that radial basis function
More informationScattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions
Chapter 3 Scattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions 3.1 Scattered Data Interpolation with Polynomial Precision Sometimes the assumption on the
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 40: Symmetric RBF Collocation in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 2 Part 3: Native Space for Positive Definite Kernels Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2014 fasshauer@iit.edu MATH
More informationData fitting by vector (V,f)-reproducing kernels
Data fitting by vector (V,f-reproducing kernels M-N. Benbourhim to appear in ESAIM.Proc 2007 Abstract In this paper we propose a constructive method to build vector reproducing kernels. We define the notion
More informationMeshfree Approximation Methods with MATLAB
Interdisciplinary Mathematical Sc Meshfree Approximation Methods with MATLAB Gregory E. Fasshauer Illinois Institute of Technology, USA Y f? World Scientific NEW JERSEY LONDON SINGAPORE BEIJING SHANGHAI
More informationKernel-based Approximation. Methods using MATLAB. Gregory Fasshauer. Interdisciplinary Mathematical Sciences. Michael McCourt.
SINGAPORE SHANGHAI Vol TAIPEI - Interdisciplinary Mathematical Sciences 19 Kernel-based Approximation Methods using MATLAB Gregory Fasshauer Illinois Institute of Technology, USA Michael McCourt University
More informationComputational Aspects of Radial Basis Function Approximation
Working title: Topics in Multivariate Approximation and Interpolation 1 K. Jetter et al., Editors c 2005 Elsevier B.V. All rights reserved Computational Aspects of Radial Basis Function Approximation Holger
More informationRadial Basis Functions I
Radial Basis Functions I Tom Lyche Centre of Mathematics for Applications, Department of Informatics, University of Oslo November 14, 2008 Today Reformulation of natural cubic spline interpolation Scattered
More informationSOLVING HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS IN SPHERICAL GEOMETRY WITH RADIAL BASIS FUNCTIONS
SOLVING HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS IN SPHERICAL GEOMETRY WITH RADIAL BASIS FUNCTIONS NATASHA FLYER 1. Introduction Mathematical modeling of space and climate phenomena generally requires
More informationStability constants for kernel-based interpolation processes
Dipartimento di Informatica Università degli Studi di Verona Rapporto di ricerca Research report 59 Stability constants for kernel-based interpolation processes Stefano De Marchi Robert Schaback Dipartimento
More informationQMC designs and covering of spheres by spherical caps
QMC designs and covering of spheres by spherical caps Ian H. Sloan University of New South Wales, Sydney, Australia Joint work with Johann Brauchart, Josef Dick, Ed Saff (Vanderbilt), Rob Womersley and
More informationLocal radial basis function approximation on the sphere
Local radial basis function approximation on the sphere Kerstin Hesse and Q. T. Le Gia School of Mathematics and Statistics, The University of New South Wales, Sydney NSW 05, Australia April 30, 007 Abstract
More informationStability of Kernel Based Interpolation
Stability of Kernel Based Interpolation Stefano De Marchi Department of Computer Science, University of Verona (Italy) Robert Schaback Institut für Numerische und Angewandte Mathematik, University of Göttingen
More informationNumerical Integration in Meshfree Methods
Numerical Integration in Meshfree Methods Pravin Madhavan New College University of Oxford A thesis submitted for the degree of Master of Science in Mathematical Modelling and Scientific Computing Trinity
More informationNumerical Integration over unit sphere by using spherical t-designs
Numerical Integration over unit sphere by using spherical t-designs Congpei An 1 1,Institute of Computational Sciences, Department of Mathematics, Jinan University Spherical Design and Numerical Analysis
More informationA Comparison between Solving Two Dimensional Integral Equations by the Traditional Collocation Method and Radial Basis Functions
Applied Mathematical Sciences, Vol. 5, 2011, no. 23, 1145-1152 A Comparison between Solving Two Dimensional Integral Equations by the Traditional Collocation Method and Radial Basis Functions Z. Avazzadeh
More informationNEEDLET APPROXIMATION FOR ISOTROPIC RANDOM FIELDS ON THE SPHERE
NEEDLET APPROXIMATION FOR ISOTROPIC RANDOM FIELDS ON THE SPHERE Quoc T. Le Gia Ian H. Sloan Yu Guang Wang Robert S. Womersley School of Mathematics and Statistics University of New South Wales, Australia
More informationStrictly Positive Definite Functions on a Real Inner Product Space
Strictly Positive Definite Functions on a Real Inner Product Space Allan Pinkus Abstract. If ft) = a kt k converges for all t IR with all coefficients a k 0, then the function f< x, y >) is positive definite
More informationRecent Results for Moving Least Squares Approximation
Recent Results for Moving Least Squares Approximation Gregory E. Fasshauer and Jack G. Zhang Abstract. We describe two experiments recently conducted with the approximate moving least squares (MLS) approximation
More informationDiscrete Projection Methods for Integral Equations
SUB Gttttingen 7 208 427 244 98 A 5141 Discrete Projection Methods for Integral Equations M.A. Golberg & C.S. Chen TM Computational Mechanics Publications Southampton UK and Boston USA Contents Sources
More informationSolving the Dirichlet-to-Neumann map on an oblate spheroid by a mesh-free method
Solving the Dirichlet-to-Neumann map on an oblate spheroid by a mesh-free method A. Costea a, Q.T. Le Gia b,, D. Pham c, E. P. Stephan a a Institut für Angewandte Mathematik and QUEST (Centre for Quantum
More informationA orthonormal basis for Radial Basis Function approximation
A orthonormal basis for Radial Basis Function approximation 9th ISAAC Congress Krakow, August 5-9, 2013 Gabriele Santin, joint work with S. De Marchi Department of Mathematics. Doctoral School in Mathematical
More informationCompression on the digital unit sphere
16th Conference on Applied Mathematics, Univ. of Central Oklahoma, Electronic Journal of Differential Equations, Conf. 07, 001, pp. 1 4. ISSN: 107-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu
More informationINFINITUDE OF MINIMALLY SUPPORTED TOTALLY INTERPOLATING BIORTHOGONAL MULTIWAVELET SYSTEMS WITH LOW APPROXIMATION ORDERS. Youngwoo Choi and Jaewon Jung
Korean J. Math. (0) No. pp. 7 6 http://dx.doi.org/0.68/kjm.0...7 INFINITUDE OF MINIMALLY SUPPORTED TOTALLY INTERPOLATING BIORTHOGONAL MULTIWAVELET SYSTEMS WITH LOW APPROXIMATION ORDERS Youngwoo Choi and
More informationConsistency Estimates for gfd Methods and Selection of Sets of Influence
Consistency Estimates for gfd Methods and Selection of Sets of Influence Oleg Davydov University of Giessen, Germany Localized Kernel-Based Meshless Methods for PDEs ICERM / Brown University 7 11 August
More informationUniform approximation on the sphere by least squares polynomials
Uniform approximation on the sphere by least squares polynomials Woula Themistoclakis Marc Van Barel August 7, 8 Abstract The paper concerns the uniform polynomial approximation of a function f, continuous
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 39: Non-Symmetric RBF Collocation in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter
More informationLeast Squares Approximation
Chapter 6 Least Squares Approximation As we saw in Chapter 5 we can interpret radial basis function interpolation as a constrained optimization problem. We now take this point of view again, but start
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 1 Part 3: Radial Basis Function Interpolation in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2014 fasshauer@iit.edu
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 6: Scattered Data Interpolation with Polynomial Precision Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationDual Bases and Discrete Reproducing Kernels: A Unified Framework for RBF and MLS Approximation
Dual Bases and Discrete Reproducing Kernels: A Unified Framework for RBF and MLS Approimation G. E. Fasshauer Abstract Moving least squares (MLS) and radial basis function (RBF) methods play a central
More informationScattered Data Approximation of Noisy Data via Iterated Moving Least Squares
Scattered Data Approximation o Noisy Data via Iterated Moving Least Squares Gregory E. Fasshauer and Jack G. Zhang Abstract. In this paper we ocus on two methods or multivariate approximation problems
More informationApproximately dual frames in Hilbert spaces and applications to Gabor frames
Approximately dual frames in Hilbert spaces and applications to Gabor frames Ole Christensen and Richard S. Laugesen October 22, 200 Abstract Approximately dual frames are studied in the Hilbert space
More informationKey words. Radial basis function, scattered data interpolation, hierarchical matrices, datasparse approximation, adaptive cross approximation
HIERARCHICAL MATRIX APPROXIMATION FOR KERNEL-BASED SCATTERED DATA INTERPOLATION ARMIN ISKE, SABINE LE BORNE, AND MICHAEL WENDE Abstract. Scattered data interpolation by radial kernel functions leads to
More informationarxiv: v1 [math.ca] 23 Oct 2018
A REMARK ON THE ARCSINE DISTRIBUTION AND THE HILBERT TRANSFORM arxiv:80.08v [math.ca] 3 Oct 08 RONALD R. COIFMAN AND STEFAN STEINERBERGER Abstract. We prove that if fx) ) /4 L,) and its Hilbert transform
More informationA Posteriori Error Bounds for Meshless Methods
A Posteriori Error Bounds for Meshless Methods Abstract R. Schaback, Göttingen 1 We show how to provide safe a posteriori error bounds for numerical solutions of well-posed operator equations using kernel
More informationKernel Method: Data Analysis with Positive Definite Kernels
Kernel Method: Data Analysis with Positive Definite Kernels 2. Positive Definite Kernel and Reproducing Kernel Hilbert Space Kenji Fukumizu The Institute of Statistical Mathematics. Graduate University
More informationA Tilt at TILFs. Rod Nillsen University of Wollongong. This talk is dedicated to Gary H. Meisters
A Tilt at TILFs Rod Nillsen University of Wollongong This talk is dedicated to Gary H. Meisters Abstract In this talk I will endeavour to give an overview of some aspects of the theory of Translation Invariant
More informationFast and accurate methods for the discretization of singular integral operators given on surfaces
Fast and accurate methods for the discretization of singular integral operators given on surfaces James Bremer University of California, Davis March 15, 2018 This is joint work with Zydrunas Gimbutas (NIST
More informationWavelets and Image Compression. Bradley J. Lucier
Wavelets and Image Compression Bradley J. Lucier Abstract. In this paper we present certain results about the compression of images using wavelets. We concentrate on the simplest case of the Haar decomposition
More information1. Introduction. A radial basis function (RBF) interpolant of multivariate data (x k, y k ), k = 1, 2,..., n takes the form
A NEW CLASS OF OSCILLATORY RADIAL BASIS FUNCTIONS BENGT FORNBERG, ELISABETH LARSSON, AND GRADY WRIGHT Abstract Radial basis functions RBFs form a primary tool for multivariate interpolation, and they are
More informationINVERSE AND SATURATION THEOREMS FOR RADIAL BASIS FUNCTION INTERPOLATION
MATHEMATICS OF COMPUTATION Volume 71, Number 238, Pages 669 681 S 0025-5718(01)01383-7 Article electronically published on November 28, 2001 INVERSE AND SATURATION THEOREMS FOR RADIAL BASIS FUNCTION INTERPOLATION
More informationOn the quasi- interpolation
On the quasi- interpolation WU Zongmin Shanghai Key Laboratory for Contemporary Applied Mathematics School of Mathematical Sciences, Fudan University zmwu@fudan.edu.cn http://homepage.fudan.edu.cn/ ZongminWU/
More informationIntegral Interpolation
Integral Interpolation R. K. Beatson and M. K. Langton Department of Mathematics and Statistics University of Canterbury Private Bag 4800, Christchurch 8140 New Zealand R.Beatson@math.canterbury.ac.nz
More informationReproducing Kernel Hilbert Spaces Class 03, 15 February 2006 Andrea Caponnetto
Reproducing Kernel Hilbert Spaces 9.520 Class 03, 15 February 2006 Andrea Caponnetto About this class Goal To introduce a particularly useful family of hypothesis spaces called Reproducing Kernel Hilbert
More informationD. Shepard, Shepard functions, late 1960s (application, surface modelling)
Chapter 1 Introduction 1.1 History and Outline Originally, the motivation for the basic meshfree approximation methods (radial basis functions and moving least squares methods) came from applications in
More information3. Some tools for the analysis of sequential strategies based on a Gaussian process prior
3. Some tools for the analysis of sequential strategies based on a Gaussian process prior E. Vazquez Computer experiments June 21-22, 2010, Paris 21 / 34 Function approximation with a Gaussian prior Aim:
More informationOn some positive definite functions
isid/ms/205/0 September 2, 205 http://www.isid.ac.in/ statmath/index.php?module=preprint On some positive definite functions Rajendra Bhatia and Tanvi Jain Indian Statistical Institute, Delhi Centre 7,
More informationReproducing Kernels of Generalized Sobolev Spaces via a Green Function Approach with Distributional Operators
Noname manuscript No. (will be inserted by the editor) Reproducing Kernels of Generalized Sobolev Spaces via a Green Function Approach with Distributional Operators Gregory E. Fasshauer Qi Ye Abstract
More informationOn Multivariate Newton Interpolation at Discrete Leja Points
On Multivariate Newton Interpolation at Discrete Leja Points L. Bos 1, S. De Marchi 2, A. Sommariva 2, M. Vianello 2 September 25, 2011 Abstract The basic LU factorization with row pivoting, applied to
More informationSpatio-Spectral Analysis on the Sphere Using Spatially Localized Spherical Harmonics Transform
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 60, NO. 3, MARCH 2012 1487 Spatio-Spectral Analysis on the Sphere Using Spatially Localized Spherical Harmonics Transform Zubair Khalid, Salman Durrani, Parastoo
More informationDecomposition of Riesz frames and wavelets into a finite union of linearly independent sets
Decomposition of Riesz frames and wavelets into a finite union of linearly independent sets Ole Christensen, Alexander M. Lindner Abstract We characterize Riesz frames and prove that every Riesz frame
More informationJUHA KINNUNEN. Harmonic Analysis
JUHA KINNUNEN Harmonic Analysis Department of Mathematics and Systems Analysis, Aalto University 27 Contents Calderón-Zygmund decomposition. Dyadic subcubes of a cube.........................2 Dyadic cubes
More informationMaximum norm estimates for energy-corrected finite element method
Maximum norm estimates for energy-corrected finite element method Piotr Swierczynski 1 and Barbara Wohlmuth 1 Technical University of Munich, Institute for Numerical Mathematics, piotr.swierczynski@ma.tum.de,
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationKernel approximation, elliptic PDE and boundary effects
Kernel approximation, elliptic PDE and boundary effects Thomas Hangelbroek University of Hawaii at Manoa Schloss Rauischholzhausen, 2016 work supported by: NSF DMS-1413726 Surface Spline Approximation
More informationA. Iske RADIAL BASIS FUNCTIONS: BASICS, ADVANCED TOPICS AND MESHFREE METHODS FOR TRANSPORT PROBLEMS
Rend. Sem. Mat. Univ. Pol. Torino Vol. 61, 3 (23) Splines and Radial Functions A. Iske RADIAL BASIS FUNCTIONS: BASICS, ADVANCED TOPICS AND MESHFREE METHODS FOR TRANSPORT PROBLEMS Abstract. This invited
More informationBiorthogonal Spline Type Wavelets
PERGAMON Computers and Mathematics with Applications 0 (00 1 0 www.elsevier.com/locate/camwa Biorthogonal Spline Type Wavelets Tian-Xiao He Department of Mathematics and Computer Science Illinois Wesleyan
More informationQuadrature using sparse grids on products of spheres
Quadrature using sparse grids on products of spheres Paul Leopardi Mathematical Sciences Institute, Australian National University. For presentation at 3rd Workshop on High-Dimensional Approximation UNSW,
More informationReproducing Kernels of Generalized Sobolev Spaces via a Green Function Approach with Differential Operators
Reproducing Kernels of Generalized Sobolev Spaces via a Green Function Approach with Differential Operators Qi Ye Abstract In this paper we introduce a generalization of the classical L 2 ( )-based Sobolev
More informationc 2005 Society for Industrial and Applied Mathematics
SIAM J. UMER. AAL. Vol. 43, o. 2, pp. 75 766 c 25 Society for Industrial and Applied Mathematics POLYOMIALS AD POTETIAL THEORY FOR GAUSSIA RADIAL BASIS FUCTIO ITERPOLATIO RODRIGO B. PLATTE AD TOBI A. DRISCOLL
More informationKatholieke Universiteit Leuven Department of Computer Science
Extensions of Fibonacci lattice rules Ronald Cools Dirk Nuyens Report TW 545, August 2009 Katholieke Universiteit Leuven Department of Computer Science Celestijnenlaan 200A B-3001 Heverlee (Belgium Extensions
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Accuracy and Optimality of RKHS Methods Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2014 fasshauer@iit.edu MATH 590 1 Outline 1 Introduction
More informationPositive Definite Kernels: Opportunities and Challenges
Positive Definite Kernels: Opportunities and Challenges Michael McCourt Department of Mathematical and Statistical Sciences University of Colorado, Denver CUNY Mathematics Seminar CUNY Graduate College
More informationRecent Developments in Numerical Methods for 4d-Var
Recent Developments in Numerical Methods for 4d-Var Mike Fisher Slide 1 Recent Developments Numerical Methods 4d-Var Slide 2 Outline Non-orthogonal wavelets on the sphere: - Motivation: Covariance Modelling
More informationReproducing Kernel Hilbert Spaces
Reproducing Kernel Hilbert Spaces Lorenzo Rosasco 9.520 Class 03 February 11, 2009 About this class Goal To introduce a particularly useful family of hypothesis spaces called Reproducing Kernel Hilbert
More informationLecture 2: Reconstruction and decomposition of vector fields on the sphere with applications
2013 Dolomites Research Week on Approximation : Reconstruction and decomposition of vector fields on the sphere with applications Grady B. Wright Boise State University What's the problem with vector fields
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationRadial basis functions topics in slides
Radial basis functions topics in 40 +1 slides Stefano De Marchi Department of Mathematics Tullio Levi-Civita University of Padova Napoli, 22nd March 2018 Outline 1 From splines to RBF 2 Error estimates,
More informationWendland Functions A C++ code to compute them
Wendland Functions A C++ code to compute them Carlos Argáez 1, Sigurdur Hafstein 1 and Peter Giesl 2 1 Faculty of Physical Sciences, University of Iceland, 107 Reykjavík, Iceland 2 Department of Mathematics,
More informationNumerische Mathematik
Numer. Math. (997) 76: 479 488 Numerische Mathematik c Springer-Verlag 997 Electronic Edition Exponential decay of C cubic splines vanishing at two symmetric points in each knot interval Sang Dong Kim,,
More informationREAL RENORMINGS ON COMPLEX BANACH SPACES
REAL RENORMINGS ON COMPLEX BANACH SPACES F. J. GARCÍA PACHECO AND A. MIRALLES Abstract. In this paper we provide two ways of obtaining real Banach spaces that cannot come from complex spaces. In concrete
More informationRecent developments in the dual receiprocity method using compactly supported radial basis functions
Recent developments in the dual receiprocity method using compactly supported radial basis functions C.S. Chen, M.A. Golberg and R.A. Schaback 3 Department of Mathematical Sciences University of Nevada,
More information