Overlapping additive Schwarz preconditioners for interpolation on the unit sphere by spherical radial basis functions
|
|
- Rudolph Cross
- 5 years ago
- Views:
Transcription
1 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 decomposition technique for the interpolation problem on the unit sphere using spherical basis functions. 1 Introduction When solving the interpolation problem on the unit sphere based on scattered data, we observe a very ill-conditioned linear system. In this paper, we introduce overlapping additive Schwarz preconditioners for the interpolation problem using spherical basis functions in which the native space can be identified with a Sobolev space. Let S n be the unit sphere in R n+1, n = 1,,..., and suppose that X = {x 1,...,x N } is a set of scattered points lying on S n. Given real numbers f i, i = 1,...,N, we want to find a smooth function u defined on S n which interpolates the data, namely, u(x i ) = f i, i = 1,...,N. (1.1) The problem (1.1) where the underlying domain is the sphere arises in many areas including, e.g., geodesy and earth science in which the sphere is taken as a model for the earth. Even though in practice measured data usually contain noise which means that approximation of data is necessary, for simplicity, in this paper we consider only the interpolation problem, with scattered data collected from satellites. A review paper by Fasshauer and Schumaker [5] discusses available methods for interpolation of scattered data on the sphere. Two main methods are the one that uses spherical splines [1], and another that uses spherical radial basis functions [13, 14]. We follow the second approach in this paper. School of Mathematics and Statistics, The University of New South Wales, Sydney 05, Australia. qlegia@unsw.edu.au, Thanh.Tran@unsw.edu.au, I.Sloan@unsw.edu.au 1
2 It is perceived both theoretically [10] and numerically (Section 7) that the matrices arising from this interpolation problem are very ill-conditioned. We shall present an algorithm which uses spherical radial basis functions to construct the interpolant and an additive Schwarz preconditioner to accelerate the solution process. We shall give an estimate for the condition number of the preconditioned system. Even though domain decomposition methods have been extensively studied for finite-element and boundary-element methods, not much has been done for meshless methods using radial basis functions. For the interpolation problem in R n using radial basis functions, the idea of dividing the scattered data set into smaller subsets for the purpose of defining the Schwarz alternating algorithm is proposed in []. Work on applying the multiplicative Schwarz alternating algorithm using spherical splines is also carried out in [7], but in that work the data points are not scattered, and again the Schwarz method is not used as a preconditioner. None of these papers uses domain decomposition methods as preconditioners to be solved with the conjugate-gradient method, and none studies the condition numbers of the preconditioned systems. The purpose of this paper is to fill this gap. We shall study the condition number of the system preconditioned by an overlapping additive Schwarz method. In [9] we studied the use of additive Schwarz preconditioners for elliptic partial differential equations on the sphere, and proved a bound for the condition number. Similar approach will be used in the present paper for the interpolation problem. Preliminaries In this section, we will review spherical harmonics, function spaces on the unit Euclidean sphere S n R n+1, and spherical basis functions..1 Spherical harmonics Spherical harmonics are the restriction of homogeneous harmonic polynomials in R n+1 to the unit sphere S n. We denote an orthonormal (with respect to the L (S n ) inner product) basis for the spherical harmonics of degree l by {Y l,k : k = 1,...,N(n,l)}, l = 0, 1,..., where N(n,l) is the dimension of the space of all spherical harmonics of degree l; the values of N(n,l) are (see [1]): N(n, 0) = 1 and N(n,l) = (l + n 1)Γ(l + n 1) Γ(l + 1)Γ(n) for l 1. The asymptotic behaviour of N(n,l) for fixed n and increasing l is O(l n 1 ). The spherical harmonics {Y l,k : l = 0, 1,...;k = 1,...,N(n,l)} form a complete orthonormal basis for L (S n ). Correspondingly, for a given function f L (S n ), we
3 define its Fourier coefficients by f l,k = f(x)y l,k (x)ds(x), S n where ds is the surface measure of the sphere S n, and represent f as a Fourier series, f = l=0 N(n,l) k=1 f l,k Y l,k, in which the equal sign is understood in the L (S n ) sense. The addition formula for spherical harmonics of the same degree l (see [1]) is N(n,l) k=0 Y l,k (x)y l,k (y) = 1 ω n N(n,l)P l (n + 1;x y), (.1) where P l (n + 1;t) is the normalized Legendre polynomial of degree l in R n+1 and ω n is the surface area of the unit sphere S n. Recall from [1] that P l (n + 1; 1) = 1 and +1 1 P l (n + 1;t)P k (n + 1;t)(1 t ) (n )/ dt = ω n ω n 1 N(n,l) δ l,k, (.) where ω n 1 is the surface area of the sphere S n 1, and δ l,k is the Kronecker delta. Spherical harmonics of degree l are eigenfunctions of the Laplace Beltrami operator on S n, with eigenvalues λ l, with λ l = l(l + n 1), (.3) i.e., Y l,k = λ l Y l,k.. Sobolev spaces on S n For a given s 0, the Sobolev space H s (S n ) on the unit sphere is defined in terms of the eigenvalues of the Laplace-Beltrami operator (see [11]), N(n,l) H s (S n ) := f L (S n ) : (1 + λ l ) s f l,k <. (.4) l=0 k=1 The norm of a function f in this space is defined to be N(n,l) f H s (S n ) = (1 + λ l ) s f l,k l=0 k=1 3 1/.
4 Sobolev spaces on S n can also be defined using local charts (see [11]). Here we use a specific atlas of charts, as in [8]. Let a spherical cap of radius α centered at p S n be defined as C(p,α) := {x S n : θ(p,x) < α}, (.5) where θ(p,x) = cos 1 (p x) is the geodesic distance between two points x,p S n. Let ˆn and ŝ denote the north and south poles of S n, respectively. Then a simple cover for the sphere is provided by U 1 = C(ˆn,θ 0 ) and U = C(ŝ,θ 0 ), where θ 0 (π/, π/3). (.6) The stereographic projection σˆn of the punctured sphere S n \ {ˆn} onto R n is defined as a mapping that maps x S n \ {ˆn} to the intersection of the equatorial hyperplane {z = 0} and the extended line that passes through x and ˆn. The stereographic projection σŝ based on ŝ can be defined analogously. We set ψ 1 = 1 tan(θ 0 /) σ ŝ U1 and ψ = 1 tan(θ 0 /) σˆn U, (.7) so that ψ k, k = 1,, maps U k onto B(0, 1), the unit ball in R n. We conclude that A = {U k,ψ k } k=1 is a C atlas of covering coordinate charts for the sphere. It is known (see [15]) that the stereographic coordinate charts {ψ k } k=1 as defined in (.7) map spherical caps to Euclidean balls, but in general concentric spherical caps are not mapped to concentric Euclidean balls. The projection ψ k, for k = 1,, does not distort too much the geodesic distance between two points x,y S n, as shown in [6]. With the atlas so defined, we define the map π k which takes a real-valued function g with compact support in U k into a real-valued function on R n by { g ψ 1 k (x), if x B(0, 1), π k (g)(x) = 0, otherwise. Let {χ k : S n R} k=1 be a partition of unity subordinated to the atlas, i.e., a pair of non-negative infinitely differentiable functions χ k on S n with compact support in U k, such that k χ k = 1. For any function f : S n R, we can use the partition of unity to write f = (χ k f), where (χ k f)(p) = χ k (p)f(p), p S n. k=1 The Sobolev space H s (S n ) is defined to be the set {f L (S n ) : π k (χ k f) H s (R n ) for k = 1, }, which is equipped with the norm ( 1/ f H s (S n ) = π k (χ k f) H s (R )). (.8) n k=1 This H s (S n ) norm is equivalent to the H s (S n ) norm given in Section.1 (see [11]). 4
5 .3 Interpolation between Sobolev spaces As outlined in [17], two Banach spaces A 0 and A 1 are said to be an interpolation pair (A 0,A 1 ) if there is a continuous inclusion A 1 A 0. For any f A 0 and t > 0, we define the K-functional as K(t,f) = inf g A 1 ( f g A 0 + t g A 1 ) 1/. (.9) For τ (0, 1) the interpolation space A τ = (A 0,A 1 ) τ is defined to be the Banach subspace of A 0 for which the following norm is finite: ( ( ) ) 1/ K(t,f) f τ = K(t,f)t τ dt L ((0, ),dt/t) =. t τ t We now apply these results to the Sobolev spaces H k (Ω), where k is a non-negative integer and Ω is a bounded open set in R n. f H k (Ω) = 0 α k 0 D α f L (Ω) The fractional Sobolev space H k+δ (Ω) (for 0 < δ < 1) is defined as follows: H k+δ (Ω) = (H k (Ω),H k+1 (Ω)) δ = {f H k (Ω) : f H k+δ (Ω) < }, 1/. where by (.9) and K(t,f) = ( ( ) ) 1/ K(t,f) dt f H k+δ (Ω) = t δ t 0 inf ( f g g H k+1 H (Ω) k (Ω) + t g H k+1 (Ω) )1/. 3 Interpolation using spherical radial basis functions 3.1 Positive definite kernels In this section, we will review necessary background on positive definite kernels on the unit sphere and spherical basis functions. A real-valued kernel Φ in C(S n S n ) is termed positive definite on S n if Φ(x,y) = Φ(y,x) and if for every finite set of distinct points X = {x 1,...,x N } on S n, the symmetric N N matrix A with entries A i,j = Φ(x i,x j ) is positive semi-definite. If 5
6 the matrix A is positive definite then Φ is called a strictly positive definite kernel (see [16, ]). Let φ be a univariate function defined on [ 1, 1] which can be expanded in terms of Legendre polynomials as φ(t) = 1 ω n l=0 N(n,l) φ(l)p l (n + 1;t), (3.1) where +1 φ(l) = ω n 1 φ(t)p l (n + 1;t)(1 t ) (n )/ dt. (3.) 1 Due to the addition formula (.1), a kernel Φ defined by Φ(x,y) = φ(x y) (3.3) can be represented as Φ(x,y) = l=0 N(n,l) k=1 φ(l)y l,k (x)y l,k (y). (3.4) In [4], a complete characterization of strictly positive definite kernels is established: the kernel Φ is strictly positive definite if and only if φ(l) 0 for all l 0 and φ(l) > 0 for infinitely many even values of l and infinitely many odd values of l; see also [16] and []. The native space N φ associated with the kernel Φ is defined as If N φ := f L (S n ) : l=0 N(n,l) k=1 f l,k φ(l) <. (3.5) c 1 (1 + λ l ) τ φ(l) c (1 + λ l ) τ, (3.6) where c 1,c > 0 are some constants and τ > n/, the native space N φ can be identified with the Sobolev space H τ (S n ) defined in (.4). Henceforth, the condition (3.6) is shortened to φ(l) (1 + λ l ) τ. 3. Spherical radial basis functions With this kernel Φ, we can now establish a set of spherical basis functions (SBFs) {Φ 1,...,Φ N } associated with a set X = {x 1,...,x N } of scattered and distinct points on S n, where Φ i (x) := Φ(x i,x). 6
7 The finite dimensional space spanned by these SBFs is denoted by V X : V X := span {Φ i : i = 1,...,N}. (3.7) We note that the SBFs Φ i, i = 1,...,N, depend only on the geodesic distance between the points x and x i. The set X is characterized by its mesh norm h X and separation radius q X defined as h X := sup y S n min x i X θ(x i,y), and q X := 1 min i j θ(x i,x j ). 3.3 Spherical basis functions derived from compactly supported radial basis functions In the following we define a positive definite kernel Φ from a univariate function φ satisfying (3.6) by using Wendland s compactly supported radial basis functions [18]. For any non-negative integer j, let { (1 r) j, 0 < r 1, ρ j (r) = 0, r > 1, and let Iρ j (r) = We define, for any non-negative integer m, r sρ j (s)ds, r 0. ρ n+1,m (r) = I m ρ m+ n+1 +1(r). The radial basis function Ψ n+1,m is defined in R n+1 as Ψ n+1,m (x) = ρ n+1,m ( x ), where denotes the Euclidean norm in R n+1. It is shown in [19, Corollary.3] that Ψ n+1,m C m (R n+1 ). For any given N and any set of N pairwise distinct points {x 1,...,x N } in R n+1, the matrix [Ψ n+1,m (x i x j )] N i, (3.8) is positive definite; see [0, Theorem 9.13]. Since x y = x y for any x,y S n, the kernel Φ defined by (3.3) with is related to the above radial basis function Ψ n+1,m by φ(t) = ρ n+1,m ( t) (3.9) Φ(x,y) = Ψ n+1,m (x y), x,y S n. 7
8 Since the matrix (3.8) is positive definite, Φ is a strictly positive definite kernel on the sphere S n. Moreover, the asymptotic behaviour of the Fourier coefficients φ(l) is, see [14, Proposition 4.6], φ(l) = O(l m n 1 ). Using (.3), we deduce that φ satisfies (3.6) with τ = m + (n + 1)/. (3.10) 3.4 Interpolation problem as a variational problem The kernel Φ being strictly positive definite, the interpolation problem in V X of scattered points using spherical basis functionsis always solvable. Given a function f whose values f(x j ) for j = 1,...,N are known, the interpolant I X f of f is defined as a linear combination of the SBFs which satisfies (I X f)(x j ) = f(x j ) for all j = 1,...,N. It is observed in [10] that the matrix A with entries A i,j = Φ i (x j ) arising from this interpolation problem is ill-conditioned. More fully, it is shown there that the least eigenvalue of the matrix A depends on the separation radius q X of the set X, which can be very small for a large set of scattered data, and also on the smoothness of the kernel Φ. (The smoother the kernel the smaller the least eigenvalue of A.) In this section, we will show how the interpolation problem can be written in the variational problem in terms of inner product in the native space. First we show that Φ(x,y) is the reproducing kernel in the native space, i.e., for any function f in the native space, we have f(x) = f, Φ(,x) φ. To see that, notice that and hence f, Φ(,x) φ = l=0 N(n,l) k=1 Φ(,x) l,k = φ(l)y l,k (x) f l,k φ(l)yl,k (x) φ(l) = l=0 N(n,l) k=1 f l,k Y l,k (x) = f(x). Suppose f is a function is the native space, and we seek the spherical basis interpolant I X f, which is a linear combination of Φ(x j, ), j = 1,...,N, so that I X f(x j ) = f(x j ), then by the property of the reproducing kernel Φ, let Φ j ( ) = Φ(x j, ), it is equivalent to I X f, Φ j φ = f, Φ j φ. Suppose V X is the finite dimensional space spanned by {Φ j : j = 1,...,N}, the interpolation problem using spherical basis functions can be written as the following problem Find I X f V X such that I X f,v φ = f,v φ v V X. (3.11) 8
9 Let A be the matrix with entries A i,j = Φ(x i,x j ) = Φ(x i, ), Φ(x j, ) φ. The problem will reduce to the problem of solving the following linear system Ac = f, where the entries of the vector f is given as f = [f j ] N in which f j = f, Φ j φ = f(x j ) for j = 1,...,N. We have the following Pythagoras theorem for the interpolation in the native space norm: I X f φ + I X f f φ = f φ. (3.1) From (3.1) we can deduce that the interpolation operator is bounded in the φ norm. 4 Additive Schwarz preconditioners 4.1 Additive Schwarz operator Additive Schwarz methods provide fast solutions to equation (3.11) by solving, at the same time, problems of smaller size. Let the space V be decomposed as V X = V V J, (4.1) where V j, j = 0,...,J are subspaces of V X, and let P j : V X V j, j = 0,...,J, be projections defined by If we define P j v,w φ = v,w φ v V X, w V j. (4.) P := P P J, (4.3) then the additive Schwarz method for equation (3.11) consists in solving, by an iterative method, the equation Pu X = g, (4.4) where the righthand side is given by g = J j=0 g j, with g j V j being solutions of g j,w φ = f,w φ, for any w V j. (4.5) A practical method to solve (4.4) is the conjugate gradient method; the additive Schwarz method (see Section 6) can be viewed as a preconditioned conjugate gradient method. 9
10 4. Bound on the condition number Bounds for λ min (P) and λ max (P), the minimum and maximum eigenvalues of the additive Schwarz operator P can be obtained by using the following lemma, see [1]. Lemma 4.1 (i) Assume that there exists a constant c 1 > 0 such that, for any u V X satisfying u = J j=0 u j with u j V j for j = 0,...,J the following inequality holds. Then u,u φ c 1 u j,u j φ j=0 λ max (P) c 1. (ii) Assume that there exists a constant c > 0 such that any u V X has a decomposition u = J j=0 u j satisfying u j,u j φ c u,u φ. j=0 Then λ min (P) c 1. 5 Additive Schwarz method for interpolation on the unit sphere In this section, we will present our main theoretical result of the paper, namely the estimate for condition number of the additive Schwarz method for the interpolation problem using spherical basis functions. Let α and β be fixed numbers satisfying 0 < α < β < π/, and let X 0 := {p j : j = 1,...,J} be a subset of X such that the cap C(p j,β) is wholly in U 1 or U (or both) of the atlas A and J S n = C(p j,α). (5.1) Assume that the support of Φ(p, ), which is a spherical cap centered at p, has radius γ. (In the case of SBF constructed from Wendland s radial basis functions, γ = π/3). For j = 1,...,J, the subset X j is defined as X j := {x k X : θ(x k,p j ) α}. (5.) 10
11 The sets X j may have different numbers of elements and may overlap each other. Because of (5.1) it is clear that X is decomposed into J overlapping subsets {X j : j = 1,...,J} of discrete points such that X = J X j. We define V j = V Xj, j = 0,...,J, i.e. V j = span {Φ k : x k X j }, so that V X = V V J. The Schwarz operator P is then defined by (4.) and (4.3). Functions in V j have supports in Γ j, where Γ j := C(p j,α + γ), j = 1,...,J. We assume that: Assumption 5.1 We can partition the index set {1,...,J} into M (for 1 M J) sets J m, for 1 m M such that if i,j J m and i j then Γ i Γ j =. The partitioning problem mentioned in Assumption 5.1 is related to the graph colouring problem [3]. We can define an undirected graph G = (V,E) in which the set of vertices V = {ν 1,...,ν J } is identified with the set of caps Γ j, and E is the set of edges, where if Γ i Γ j then there is an edge between ν i and ν j. A partition satisfying Assumption 5.1 is equivalent to a colouring of the vertices of G so that adjacent vertices have different colours. The minimal number of colours needed is called the chromatic number of G, and is denoted by δ(g). In general, it is difficult to determine the chromatic number of a graph. However, it is easy to see that δ(g) ω(g), where ω(g) is the maximal order of a complete subgraph of G, that is, it is the maximal number of vertices all of which are mutually connected. In terms of the caps, every point on the sphere S n lies in at most M 1 = ω(g) spherical caps Γ j. An upper bound of δ(g) is given in [3, Theorem 3, Chapter 5]: when G is neither a complete graph nor an odd cycle, then δ(g) (G), with (G) being the maximal degree of G. In terms of our spherical caps, each cap Γ j intersects at most M = (G) other caps. Therefore, for a given set X 0 and parameters α, γ, we can compute the lower bound M 1 and the upper bound M so that 1 M 1 = ω(g) M M = (G) J. (5.3) 11
12 5.1 A bound for λ max Lemma 5. There exists a positive constant c independent of the set X such that for any u V X satisfying u = J j=0 u j with u j V j for j = 0,...,J, u,u φ cm u j,u j φ. j=0 Proof. Using the inequality a + b ( a + b ), we have u H τ (S n ) u 0 H τ (S n ) + u j. H τ (S n ) Let k = τ, from the definition of the Sobolev norm (.8), u j H k (S n ) = π 1 (χ 1 u j ) H k (R n ) + π (χ u j ) H k (R n ) (5.4) Now, from the fact that u j V j together with Assumption 5.1 we can partition the index set {1,...,J} to M sets of indices J m so that if i,j J m then supp u i supp u j =. Then, in this proof only, let g j = π 1 (χ 1 u j ), by using the Cauchy-Schwarz inequality, we have g j H k (R n ) M = g j m=1 j J m H k (R n ) cm M m=1 j J m g j Since the supports of g i and g j are disjoint for i,j J m, i j, Thus, g j Similarly, we obtain g j H k (R n ) H k+1 (R n ) j J m g j M M H k (R n ) M m=1 M m=1 = j J m g j H k (R n ). j J m g j H k (R n ) = M j J m g j H k+1 (R n ) = M H k (R n ). (5.5) g j H k (R n ). (5.6) g j H k+1 (R n ). (5.7) 1
13 Using the definition of the K-functional (.9) and (5.6) (5.7), we conclude that g j M g j H τ (R ). n H τ (R n ) Hence, by using similar arguments for π (χ u j ), we conclude ( ) u j cm π 1 (χ 1 u j ) H τ (R n ) + π (χ u j ) H τ (R n ) H τ (S n ) = cm u j H τ (S ). n Therefore, u H τ (S n ) = u j j=0 H τ (S n ) cm u j H τ (S ). n Using the fact that u,u φ u H τ (S n ) we obtain the result. j=0 5. A bound for λ min We proceed as in [9, Lemma 5.3], in which the bilinear form a(, ) is re-defined by a(, ) :=, φ and the norm of a linear operator T is defined by T φ = sup{ Tv φ : v V X and v φ 1}. Lemma 5.3 For any u V X there exist u j V j, j = 0,...,J, satisfying u = J j=0 u j and ( ) J u j,u j φ 1 + (1 Q u,u φ ) φ, j=0 where Q = Q J Q 1, in which Q i is the orthogonal projection from V X to V with respect to, φ. The above lemma and Lemma 4.1 yeild the following estimate for the minimum eigenvalue of P: ( ) 1 J λ min (P) 1 + (1 Q. (5.8) φ ) This estimate is not an optimal bound for λ min (P), as can be seen from Table 1. In that table C 1 = P 1 (I Q) 1 φ and C j = P j Q j 1...Q 1 (I Q) 1 φ,j =,...,J. The norms of the operators were computed by using their matrix representations. 13
14 N q X cos(α) cos(β) J λ 1 min 1 + J C j 1 + J/(1 Q φ ) 1344 π/ π/ π/ Main result Table 1: Upper bound for λ 1 min (P) Theorem 5.4 The condition number of the additive Schwarz operator P is bounded by ( ) J κ(p) cm 1 + (1 Q, φ ) where c is a constant independent of M,J and the set X, and Q is defined as in An overlapping additive Schwarz algorithm Since the Earth is rotating around its own axis, and the satellite traverses from near the North pole to near the South pole then back to the North pole in an elliptical path, data along the track of satellite forms of a sequence of discrete points, see Figure 1. The projection of the satellite position onto the surface of the earth is the ground track, or subsatellite track, and this path over a period of time is a measure of the coverage of the globe from which data can be acquired. Suppose we number the scattered data following the satellite track as {x 1,x,...,x N } =: X. The algorithm to partition X in the form X = X 0 X 1 X J can be described as the following. (1) Select α (0,π/3) and β [α,π]. () The first center is p 1 = x 1 X. 14
15 (3) Define X 1 := {x k X : θ(x k,p 1 ) α}. (4) Suppose X j 1, for j > 1, has been selected around center p j 1. The next center p j is chosen from X \{p 1,...,p j 1 } such that θ(p j 1,p j ) β and θ(p k,p j ) α for k = 1,...,j. (5) The subset X j is defined as X j := {x k X : θ(x k,p j ) α}. (6) Repeat (4) and (5) until every point in X is in at least one X k. (7) Define X 0 = {p 1,...,p J }. In the following, we describe the preconditioned conjugate gradient method based on the decomposition of the scattered data X = X 0... X J. For k = 0,...,J, let A k be the restriction of the matrix A onto each subspace V k, and let R k be matrix that transform the basis of V to the basis for V k. A pseudocode. INPUT Input the scattered set X on the sphere, the right-hand side f, and the desired accuracy ǫ. SETUP (1) Partition the scattered set X into X 0... X J. () The residual vector r = [f(x j )] N. (3) The pseudo-residual vector p = 0. (4) The initial solution vector s = 0. (5) Set the iteration counter iter = 0. ITERATIVE SOLUTION (1) while r > ǫ () for to J (3) p = p + R T j A 1 j R j r. 15
16 Figure 1: Global scattered MAGSAT satellite data Figure : Domain decomposition based on satellite data 16
17 (4) end for (5) p = p + R T 0 A 1 0 R 0 r (6) If iter > 0 then set ζ 0 = ζ 1. (7) Set ζ 1 = p r. (8) iter = iter + 1. (9) If iter = 1 then define p 1 = p else p 1 = p + (ζ 1 /ζ 0 )p 1. (10) Update the residual vector (11) Update the solution vector (1) end while r = r s = s + r p p 1 Ap 1 Ap 1. r p p 1 Ap 1 p 1. 7 Numerical results In this section, we present numerical experiments on S based on globally scattered data extracted from a very large data set collected by NASA s satellite MAGSAT. Given a positive real number q, the set of scattered point X is extracted along the satellite track from the original data set so that the separation radius q X q. The set X is constructed after a two-stage thinning process: (1) Points are taken along the satellite track so that the geodesic distance between two sucessive points is q. () Points from stage (1) are re-selected so that the separation radius q X q. The number of points and their q X of each data set are listed in Table. The local supported spherical basis functions induced by compactly supported radial basis functions introduced by Wendland [18] for R 3 are used. The spherical basis interpolation I X f = j α jφ(x j, ), in which the kernel Φ is defined by where ρ 3,m (r) are given in Table 3. Φ(x,y) = ρ 3,m ( x y), 17
18 N q X π/ π/ π/ π/ π/ π/600 Table : Separation radius q X of different data sets m ρ 3,m (r) τ 1 (1 r) 4 +(4r + 1).5 (1 r) 6 +(35r + 18r + 3) (1 r) 8 +(3r 3 + 5r + 8r + 1) 4.5 Table 3: Wendland s RBFs In the numerical experiment, the basis functions used is Φ(x,y) = ρ 3,3 ( x y). The right hand side are values f(x j ) with f(x) = exp(x 1 + x + x 3 ). The stopping criteria for the conjugate gradient method is Aα f f References [1] P. Alfeld, M. Neamtu, and L. L. Schumaker. Fitting scattered data on sphere-like surfaces using spherical splines. J. Comput. Appl. Math., 73:5 43, [] R. K. Beatson, W. A. Light, and S. Billings. Fast solution of the radial basis function interpolation equations domain decomposition method. SIAM J. Sci. Comput., : , 000. [3] B. Bollobás. Modern Graph Theory. Springer Verlag, New York, [4] D. Chen, V. A. Menegatto, and X. Sun. A necessary and sufficient condition for strictly positive definite functions on spheres. Proc. Amer. Math. Soc., 131: , 003. [5] G. E. Fasshauer and L. L. Schumaker. Scattered data fitting on the sphere. In M. Dahlen, T. Lyche, and L. L. Schumaker, editors, Mathematical Methods for Curves and Surfaces II, pages , Nashville, Vanderbilt University Press. 18
19 N cos α cos β M 1 M J λ min λ max κ(p) CPU iter Table 4: Preconditioned systems [6] Q. T. Le Gia, F. J. Narcowich, J. D. Ward, and H. Wendland. Continuous and discrete least square error approximation by radial basis functions on spheres. J. Approx. Theory, to appear. [7] K. Hesse. Domain Decomposition Methods in Multiscale Geopotential Determination from SST and SGG. PhD thesis, University of Kaiserslautern, Germany, 00. [8] S. Hubbert and T. M. Morton. A Duchon framework for the sphere. J. Approx. Theory, 19:8 57, 004. [9] Q. T. Le Gia, I. H. Sloan, and T. Tran. Overlapping additive Schwarz preconditioners for elliptic PDEs on the unit sphere. Math. Comp., 008. (to appear). 19
20 N λ min λ max κ(a) CPU iter E E E E E E E E E E E E E E+04 0.E Table 5: Unpreconditioned systems [10] J. Levesley, Z. Luo, and X. Sun. Norm estimates of interpolation matrices and their inverses associated with strictly positive definite functions. Proc. Amer. Math. Soc., 17:17 134, [11] J. L. Lions and E. Magenes. Non-Homogeneous Boundary Value Problems and Applications I. Springer-Verlag, New York, 197. [1] C. Müller. Spherical Harmonics, volume 17 of Lecture Notes in Mathematics. Springer-Verlag, Berlin, [13] F. J. Narcowich, X. Sun, J. D. Ward, and H. Wendland. Direct and inverse Sobolev error estimates for scattered data interpolation via spherical basis functions. Found. Comput. Math., 7: , 007. [14] F. J. Narcowich and J. D. Ward. Scattered data interpolation on spheres: error estimates and locally supported basis functions. SIAM J. Math. Anal., 33: , 00. [15] J. G. Ratcliffe. Foundations of Hyperbolic Manifolds. Springer, New York, [16] I. J. Schoenberg. Positive definite function on spheres. Duke Math. J., 9:96 108, 194. [17] H. Triebel. Interpolation Theory, Function Spaces and Differential Operators. North-Holland, Amsterdam, [18] H. Wendland. Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree. Advances in Comp. Math., 4: , [19] H. Wendland. Error estimates for interpolation by compactly supported radial basis functions of minimal degree. J. Approx. Theory, 93:58 7, [0] H. Wendland. Scattered Data Approximation. Cambridge University Press, Cambridge,
21 [1] O. B. Widlund. Optimal iterative refinement methods. In T. F. Chan, Roland Glowinski, Jacques Périaux, and O. B. Widlund, editors, Proc. of the Second International Symposium on Domain Decomposition Methods for Partial Differential Equations, pages , Philadelphia, PA, SIAM. [] Y. Xu and E. W. Cheney. Strictly positive definite functions on spheres. Proc. Amer. Math. Soc., 116: ,
Boundary 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 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 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 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 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 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 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 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 informationAPPROXIMATION OF LINEAR PARTIAL DIFFERENTIAL EQUATIONS ON SPHERES. A Dissertation QUOC THONG LE GIA
APPROIMATION OF LINEAR PARTIAL DIFFERENTIAL EQUATIONS ON SPHERES A Dissertation by QUOC THONG LE GIA Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the requirements
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationDomain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions
Domain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions Ernst P. Stephan 1, Matthias Maischak 2, and Thanh Tran 3 1 Institut für Angewandte Mathematik, Leibniz
More informationTHE L 2 -HODGE THEORY AND REPRESENTATION ON R n
THE L 2 -HODGE THEORY AND REPRESENTATION ON R n BAISHENG YAN Abstract. We present an elementary L 2 -Hodge theory on whole R n based on the minimization principle of the calculus of variations and some
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 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 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 informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
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 informationMATH 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 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 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 informationAnalysis in weighted spaces : preliminary version
Analysis in weighted spaces : preliminary version Frank Pacard To cite this version: Frank Pacard. Analysis in weighted spaces : preliminary version. 3rd cycle. Téhéran (Iran, 2006, pp.75.
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 informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationThe Closed Form Reproducing Polynomial Particle Shape Functions for Meshfree Particle Methods
The Closed Form Reproducing Polynomial Particle Shape Functions for Meshfree Particle Methods by Hae-Soo Oh Department of Mathematics, University of North Carolina at Charlotte, Charlotte, NC 28223 June
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 informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
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 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 informationApproximation Theory on Manifolds
ATHEATICAL and COPUTATIONAL ETHODS Approximation Theory on anifolds JOSE ARTINEZ-ORALES Universidad Nacional Autónoma de éxico Instituto de atemáticas A.P. 273, Admon. de correos #3C.P. 62251 Cuernavaca,
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 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 informationMath 350 Fall 2011 Notes about inner product spaces. In this notes we state and prove some important properties of inner product spaces.
Math 350 Fall 2011 Notes about inner product spaces In this notes we state and prove some important properties of inner product spaces. First, recall the dot product on R n : if x, y R n, say x = (x 1,...,
More informationSplines which are piecewise solutions of polyharmonic equation
Splines which are piecewise solutions of polyharmonic equation Ognyan Kounchev March 25, 2006 Abstract This paper appeared in Proceedings of the Conference Curves and Surfaces, Chamonix, 1993 1 Introduction
More informationOptimization and Optimal Control in Banach Spaces
Optimization and Optimal Control in Banach Spaces Bernhard Schmitzer October 19, 2017 1 Convex non-smooth optimization with proximal operators Remark 1.1 (Motivation). Convex optimization: easier to solve,
More informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
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 informationSolutions to Problem Set 5 for , Fall 2007
Solutions to Problem Set 5 for 18.101, Fall 2007 1 Exercise 1 Solution For the counterexample, let us consider M = (0, + ) and let us take V = on M. x Let W be the vector field on M that is identically
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 informationOn Riesz-Fischer sequences and lower frame bounds
On Riesz-Fischer sequences and lower frame bounds P. Casazza, O. Christensen, S. Li, A. Lindner Abstract We investigate the consequences of the lower frame condition and the lower Riesz basis condition
More informationEXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018
EXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018 While these notes are under construction, I expect there will be many typos. The main reference for this is volume 1 of Hörmander, The analysis of liner
More informationOn uniqueness in the inverse conductivity problem with local data
On uniqueness in the inverse conductivity problem with local data Victor Isakov June 21, 2006 1 Introduction The inverse condictivity problem with many boundary measurements consists of recovery of conductivity
More informationWave equation on manifolds and finite speed of propagation
Wave equation on manifolds and finite speed of propagation Ethan Y. Jaffe Let M be a Riemannian manifold (without boundary), and let be the (negative of) the Laplace-Beltrami operator. In this note, we
More informationNew constructions of domain decomposition methods for systems of PDEs
New constructions of domain decomposition methods for systems of PDEs Nouvelles constructions de méthodes de décomposition de domaine pour des systèmes d équations aux dérivées partielles V. Dolean?? F.
More informationVectors in Function Spaces
Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also
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 informationMathematics 530. Practice Problems. n + 1 }
Department of Mathematical Sciences University of Delaware Prof. T. Angell October 19, 2015 Mathematics 530 Practice Problems 1. Recall that an indifference relation on a partially ordered set is defined
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
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 informationProjected Surface Finite Elements for Elliptic Equations
Available at http://pvamu.edu/aam Appl. Appl. Math. IN: 1932-9466 Vol. 8, Issue 1 (June 2013), pp. 16 33 Applications and Applied Mathematics: An International Journal (AAM) Projected urface Finite Elements
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 informationADJOINTS, ABSOLUTE VALUES AND POLAR DECOMPOSITIONS
J. OPERATOR THEORY 44(2000), 243 254 c Copyright by Theta, 2000 ADJOINTS, ABSOLUTE VALUES AND POLAR DECOMPOSITIONS DOUGLAS BRIDGES, FRED RICHMAN and PETER SCHUSTER Communicated by William B. Arveson Abstract.
More information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
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 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 informationPICARD S THEOREM STEFAN FRIEDL
PICARD S THEOREM STEFAN FRIEDL Abstract. We give a summary for the proof of Picard s Theorem. The proof is for the most part an excerpt of [F]. 1. Introduction Definition. Let U C be an open subset. A
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationDerivatives of Harmonic Bergman and Bloch Functions on the Ball
Journal of Mathematical Analysis and Applications 26, 1 123 (21) doi:1.16/jmaa.2.7438, available online at http://www.idealibrary.com on Derivatives of Harmonic ergman and loch Functions on the all oo
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 informationRecall that any inner product space V has an associated norm defined by
Hilbert Spaces Recall that any inner product space V has an associated norm defined by v = v v. Thus an inner product space can be viewed as a special kind of normed vector space. In particular every inner
More informationThe oblique derivative problem for general elliptic systems in Lipschitz domains
M. MITREA The oblique derivative problem for general elliptic systems in Lipschitz domains Let M be a smooth, oriented, connected, compact, boundaryless manifold of real dimension m, and let T M and T
More informationEigenvalue comparisons in graph theory
Eigenvalue comparisons in graph theory Gregory T. Quenell July 1994 1 Introduction A standard technique for estimating the eigenvalues of the Laplacian on a compact Riemannian manifold M with bounded curvature
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationarxiv: v1 [math.na] 11 Jul 2011
Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients arxiv:07.260v [math.na] Jul 20 Blanca Ayuso De Dios, Michael Holst 2, Yunrong Zhu 2, and Ludmil Zikatanov
More informationDeviation Measures and Normals of Convex Bodies
Beiträge zur Algebra und Geometrie Contributions to Algebra Geometry Volume 45 (2004), No. 1, 155-167. Deviation Measures Normals of Convex Bodies Dedicated to Professor August Florian on the occasion
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationAn homotopy method for exact tracking of nonlinear nonminimum phase systems: the example of the spherical inverted pendulum
9 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June -, 9 FrA.5 An homotopy method for exact tracking of nonlinear nonminimum phase systems: the example of the spherical inverted
More informationBIHARMONIC WAVE MAPS INTO SPHERES
BIHARMONIC WAVE MAPS INTO SPHERES SEBASTIAN HERR, TOBIAS LAMM, AND ROLAND SCHNAUBELT Abstract. A global weak solution of the biharmonic wave map equation in the energy space for spherical targets is constructed.
More informationPost-processing of solutions of incompressible Navier Stokes equations on rotating spheres
ANZIAM J. 50 (CTAC2008) pp.c90 C106, 2008 C90 Post-processing of solutions of incompressible Navier Stokes equations on rotating spheres M. Ganesh 1 Q. T. Le Gia 2 (Received 14 August 2008; revised 03
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationThe Nearest Doubly Stochastic Matrix to a Real Matrix with the same First Moment
he Nearest Doubly Stochastic Matrix to a Real Matrix with the same First Moment William Glunt 1, homas L. Hayden 2 and Robert Reams 2 1 Department of Mathematics and Computer Science, Austin Peay State
More informationYour first day at work MATH 806 (Fall 2015)
Your first day at work MATH 806 (Fall 2015) 1. Let X be a set (with no particular algebraic structure). A function d : X X R is called a metric on X (and then X is called a metric space) when d satisfies
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationMULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday.
MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* DOUGLAS N ARNOLD, RICHARD S FALK, and RAGNAR WINTHER Dedicated to Professor Jim Douglas, Jr on the occasion of his seventieth birthday Abstract
More informationMatrix square root and interpolation spaces
Matrix square root and interpolation spaces Mario Arioli and Daniel Loghin m.arioli@rl.ac.uk STFC-Rutherford Appleton Laboratory, and University of Birmingham Sparse Days,CERFACS, Toulouse, 2008 p.1/30
More informationRemarks on the Rademacher-Menshov Theorem
Remarks on the Rademacher-Menshov Theorem Christopher Meaney Abstract We describe Salem s proof of the Rademacher-Menshov Theorem, which shows that one constant works for all orthogonal expansions in all
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 informationThe following definition is fundamental.
1. Some Basics from Linear Algebra With these notes, I will try and clarify certain topics that I only quickly mention in class. First and foremost, I will assume that you are familiar with many basic
More informationNUMERICAL RADIUS OF A HOLOMORPHIC MAPPING
Geometric Complex Analysis edited by Junjiro Noguchi et al. World Scientific, Singapore, 1995 pp.1 7 NUMERICAL RADIUS OF A HOLOMORPHIC MAPPING YUN SUNG CHOI Department of Mathematics Pohang University
More informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
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 informationFriedrich symmetric systems
viii CHAPTER 8 Friedrich symmetric systems In this chapter, we describe a theory due to Friedrich [13] for positive symmetric systems, which gives the existence and uniqueness of weak solutions of boundary
More informationDomain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions
Domain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions Bernhard Hientzsch Courant Institute of Mathematical Sciences, New York University, 51 Mercer Street, New
More informationMINIMAL DEGREE UNIVARIATE PIECEWISE POLYNOMIALS WITH PRESCRIBED SOBOLEV REGULARITY
MINIMAL DEGREE UNIVARIATE PIECEWISE POLYNOMIALS WITH PRESCRIBED SOBOLEV REGULARITY Amal Al-Rashdan & Michael J. Johnson* Department of Mathematics Kuwait University P.O. Box: 5969 Safat 136 Kuwait yohnson1963@hotmail.com*
More informationConnected sum constructions in geometry and nonlinear analysis. Frank Pacard
Connected sum constructions in geometry and nonlinear analysis Frank Pacard January 19, 2008 2 Contents 1 Laplace-Beltrami operator 5 1.1 Definition................................ 5 1.2 Spectrum of the
More informationCharacterization of half-radial matrices
Characterization of half-radial matrices Iveta Hnětynková, Petr Tichý Faculty of Mathematics and Physics, Charles University, Sokolovská 83, Prague 8, Czech Republic Abstract Numerical radius r(a) is the
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 informationLEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.
LEGENDRE POLYNOMIALS AND APPLICATIONS We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.. Legendre equation: series solutions The Legendre equation is
More informationA Computational Approach to Study a Logistic Equation
Communications in MathematicalAnalysis Volume 1, Number 2, pp. 75 84, 2006 ISSN 0973-3841 2006 Research India Publications A Computational Approach to Study a Logistic Equation G. A. Afrouzi and S. Khademloo
More informationThere are two things that are particularly nice about the first basis
Orthogonality and the Gram-Schmidt Process In Chapter 4, we spent a great deal of time studying the problem of finding a basis for a vector space We know that a basis for a vector space can potentially
More informationMEAN CURVATURE FLOW OF ENTIRE GRAPHS EVOLVING AWAY FROM THE HEAT FLOW
MEAN CURVATURE FLOW OF ENTIRE GRAPHS EVOLVING AWAY FROM THE HEAT FLOW GREGORY DRUGAN AND XUAN HIEN NGUYEN Abstract. We present two initial graphs over the entire R n, n 2 for which the mean curvature flow
More informationStatistics 612: L p spaces, metrics on spaces of probabilites, and connections to estimation
Statistics 62: L p spaces, metrics on spaces of probabilites, and connections to estimation Moulinath Banerjee December 6, 2006 L p spaces and Hilbert spaces We first formally define L p spaces. Consider
More informationNEW BOUNDS FOR TRUNCATION-TYPE ERRORS ON REGULAR SAMPLING EXPANSIONS
NEW BOUNDS FOR TRUNCATION-TYPE ERRORS ON REGULAR SAMPLING EXPANSIONS Nikolaos D. Atreas Department of Mathematics, Aristotle University of Thessaloniki, 54006, Greece, e-mail:natreas@auth.gr Abstract We
More informationOn the bang-bang property of time optimal controls for infinite dimensional linear systems
On the bang-bang property of time optimal controls for infinite dimensional linear systems Marius Tucsnak Université de Lorraine Paris, 6 janvier 2012 Notation and problem statement (I) Notation: X (the
More informationComplex geometrical optics solutions for Lipschitz conductivities
Rev. Mat. Iberoamericana 19 (2003), 57 72 Complex geometrical optics solutions for Lipschitz conductivities Lassi Päivärinta, Alexander Panchenko and Gunther Uhlmann Abstract We prove the existence of
More informationUSING FUNCTIONAL ANALYSIS AND SOBOLEV SPACES TO SOLVE POISSON S EQUATION
USING FUNCTIONAL ANALYSIS AND SOBOLEV SPACES TO SOLVE POISSON S EQUATION YI WANG Abstract. We study Banach and Hilbert spaces with an eye towards defining weak solutions to elliptic PDE. Using Lax-Milgram
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 informationTHE STOKES SYSTEM R.E. SHOWALTER
THE STOKES SYSTEM R.E. SHOWALTER Contents 1. Stokes System 1 Stokes System 2 2. The Weak Solution of the Stokes System 3 3. The Strong Solution 4 4. The Normal Trace 6 5. The Mixed Problem 7 6. The Navier-Stokes
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 information290 J.M. Carnicer, J.M. Pe~na basis (u 1 ; : : : ; u n ) consisting of minimally supported elements, yet also has a basis (v 1 ; : : : ; v n ) which f
Numer. Math. 67: 289{301 (1994) Numerische Mathematik c Springer-Verlag 1994 Electronic Edition Least supported bases and local linear independence J.M. Carnicer, J.M. Pe~na? Departamento de Matematica
More information