The Galerkin method for solving problem (1.1) numerically is to select a finitedimensional
|
|
- Alberta Booth
- 6 years ago
- Views:
Transcription
1 COMPRESSON OF ANSOTROPC TENSOR-PRODUCT DSCRETZATONS S KNAPEK D MÜNCHEN STKNAPEK}@YAHOODE Abstract This paper is concerned with the discretization and compression of integral operators t continues the work of [9] Based on the framework of tensor-product biorthogonal wavelet bases we construct compression schemes for operator adapted sparse grid type discretizations These discretization and compression schemes preserve the approximation order of the standard full-grid spaces We give detailed information on the cases in which the curse of dimension can be broken We show that our compression schemes are optimal up to logarithmic factors AMS subject classification 4A7, 4A25, 4A30, 4A65, 45L05, 45L0, 65D99, 65N2, 65N30, 65N38, 65N55 Key words biorthogonal wavelets, compression, optimized approximation spaces, sparse grids, hyperbolic cross approximation, integral equations ntroduction Consider a variational problem: Given f H s, s R, find u H s such that () a(u, v) = (f, v) v H s The Galerkin method for solving problem () numerically is to select a finitedimensional subspace V J from H s L 2 and to solve the variational problem therein Fixing a basis of V J finally leads to a linear system of equations (2) A J x J = b J of dimension dim(v J ) The work count to solve (2) to a prescribed accuracy can often be made proportional to the number of non-zero entries in A J, see eg [5, 2] For conventional spaces that belong to uniform full grids with subdivision rate two, this number is at least of the order dim(v J ) = O(2 nj ), ie it grows exponentially with the dimension n This is particularly problematic in the case of integral operators with global kernel K, (3) Au(x) := K(x, y)u(y)dy, where a naive Galerkin discretization on a uniform full grid with O(2 nj ) unknowns leads to a discrete operator with O(2 2nJ ) entries This makes matrix vector multiplications, as they are needed in iterative methods, prohibitively expensive for large n or large J Approaches to reduce the work count for approximating the solution of () are either based on the smoothness properties of the solution u (and eventually the kernel K), ie on decay properties in Fourier-space, or on the decay properties of the kernel K in physical space The original idea of the first approach applied to () was to exploit smoothness properties of the solution u This leads to the replacement of the full grid approximation space by a so-called sparse grid space (that also appeared under the names hyperbolic cross approximation or boolean blending schemes) [, 8, 0, 3, 4, 5, 24, 25, 26, 27] t has been shown that under some assumptions on the approximation and smoothness properties of the underlying anisotropic tensor-product basis functions and provided
2 that certain additional regularity assumptions on the solution u of () are fulfilled, the resulting approximation spaces exhibit the same order of approximation as the full grid space, while having significantly less dimension This approach has been applied successfully to differential [4, 5, 27] as well as integral operators [5, 6] The full potential of the sparse grid approach for integral operators is revealed when it is applied to the kernel K, exploiting K s smoothness properties This leads to a significant further reduction of the work count, see [2, 8, 9, 22], and hence, reduces the need for additional compression schemes on top of these discretizations, compare section 3 Another well known strategy to reduce the cost of matrix vector multiplications is to exploit eventual decay properties of the kernel in physical space Such decay properties do hold for pseudo-differential operators and hence for example for the single and double layer potential operator as well as for classical differential operators [20] For isotropic biorthogonal wavelets it has been shown that, with a sufficient number of vanishing moments, most entries in the corresponding stiffness matrices can be replaced by zero without destroying the order of approximation [3, 6, 7, 9, 23] Then for isotropic tensor-product wavelets on a full grid the number of entries in the stiffness matrix after compression is only O(J k 2 nj ) with some k R + or even O(2 nj ) [23] Corresponding investigations concerning anisotropic tensor-product discretizations seem still to be missing, especially together with sparse grid type discretizations n this paper we continue the work of [9] and develop compression schemes for anisotropic tensor-product type discretisations of (3) We show that in many cases a suitable discretization scheme renders additional compression unnecessary We identify the cases, where additional compression is worthwhile For the relevant cases we we develop a suitable compression scheme along the lines of [23] and show how optimal compression (up to logarithmic factors) can be obtained The remainder of the paper is as follows Section 2 introduces the necessary notation and summarizes the basic facts about biorthogonal wavelet bases, tensor-product spaces, norm equivalences and the smoothness classes we consider (they are certain intersections of classes of functions with dominating ed derivatives, see (2) below) n Section 3 we discuss anisotropic tensor-product discretizations for integral operators including full and sparse grid type discretizations n Section 4 we introduce and discuss our compression scheme 2 Preliminaries Multi-indices (vectors) are written boldface, for example j = (j,, j n ) nequalities like l t or l 0 are to be understood componentwise We write x y if there exist C, C 2 independent of any parameters x or y may depend on such that C y x C 2 y n the rest of the paper C denotes a generic constant which may depend on the smoothness assumptions and on the dimension n of the problem Moreover, for j, k N n we write j := i n ( j i ), j := n i= j i, (j, k) := ( j, k ) and (j, k) := j + k With dist(x, y) we denote the Euclidian distance between x and y 2 Sobolev spaces Let us denote by H t ( n ), t R, a scale of standard Sobolev spaces on n := [0, ] n and by L 2 ( n ) the space of L 2 -integrable functions on n When t < 0, H t ( n ) is defined as the dual of H t ( n ), ie, H t := (H t ) (incorporating eventual boundary conditions) Now we fix the smoothness assumptions we consider Let t R, l R + 0, := (,, ) and e i := (0,, 0,, 0,, 0) the i-th unit-vector in R n Then we define (2) H t,l (n ) := H t+le ( n ) H t+le n ( n ), 2
3 where H k (n ) := H k () H k n () for k R n This definition includes the class of functions with dominating ed derivative H t,0 (n ) as well as the standard isotropic Sobolev spaces H t ( n ) = H 0,t (n ) Note that H t,l is a ture of classes of functions with dominating ed derivative H t (T n ) = H t,0 (T n ) and standard Sobolev spaces H t (T n ) = H 0,t (T n ) For additional information on these spaces, specifically on the definition via Fourier transform, the treatment of boundary conditions and the connections between spaces of bounded ed derivative and these anisotropic tensor-product spaces see [5, 7, 8] n the rest of the paper we will drop n in H t,l (n ) and write H t,l 22 Biorthogonal wavelet bases The approximation spaces considered here are anisotropic tensor-products of univariate function spaces We start from a onedimensional multi-resolution analysis L 2 = j 0 V j and we assume that the complement spaces W j = V j V j of V j in V j are spanned by some multiscale basis functions such that we have W j = spanψ, k τ j }, where τ j is some index set defined from the subdivision rate of successive refinement levels We will consider dyadic refinement throughout the paper Moreover, we assume that ψ, k τ j } forms a Riesz-basis of W j and that there exists a dual system ψ, k τ j, j N 0 }, such that ψ, ψ j k L 2 = δ jj δ kk, j, j N 0, k τ j, k τ j holds We assume ψ L 2 = and diam(supp(ψ )) 2 j n the following let ψ and ψ have N and Ñ vanishing moments, respectively Moreover we write W j = span ψ, k τ j } for the complement spaces spanned by the dual wavelets For the higher-dimensional case n >, let j N n 0 be given Consider the tensorproduct partition with uniform step size 2 j i into the i-th coordinate direction By W j and W j we denote the corresponding tensor-product function spaces W j := W j W jn and W j := W j W jn A basis of W j is then given by k τj ψ (x) = ψ jk (x ) ψ jnk n (x n )}, with τ j = τ j τ jn To simplify things we assume in the following the validity of the norm equivalences (22) u 2 H t,l j t, t + l ( r, r) and 2 2t j +2l j w j 2 L 2, for u = j w j H t,l, w j W j, (23) u 2 H t,l j 2 2t j +2l j w j 2 L 2, for u = j w j H t,l, w j W j, t, t + l ( r, r), see [5] Such two-sided estimates can be inferred from the validity of norm equivalences in one dimension which in turn can be inferred from direct estimates (estimates of Jackson type) and inverse estimates (Bernstein inequalities) for the primal and the dual wavelets as a consequence of approximation theory in Sobolev spaces together with interpolation and duality arguments, see [5, 5, 2] Note that for t = 0 or l = 0 we regain the (standard) norm equivalences for the isotropic Sobolev space H l and the Sobolev space with dominating ed derivative H t,0 The different factors 22t j and 2 2l j in these equivalences reflect the different smoothness requirements 3 Anisotropic tensor-product discretization Given an index set N 2n 0 we consider the approximation space V := W j W k (j,k) 3
4 and the corresponding Galerkin stiffness matrix alm for (j, k) and (l, m) τ (A ) lm :=, 0 else, with a lm := Aψ jl, ψ km L 2 and τ := τ j τ k Here we assume that eventual boundary conditions are implemented into the definition of the primal and dual wavelets The corresponding integral operator is given by A u(x) := K (x, y)u(y)dy with K (x, y) := a lm ψjl (x) ψ km (y) (j,k) (l,m) τ Typical approximation spaces are defined via the index sets (3) J = (j, k) N 2n 0 : j, k J}, J = (j, k) N 2n 0 : j J and k J} or J = (j, k) N 2n 0 : j, k J}, J N, leading to the full grid space with dimension O(2 2nJ ), the tensor product of two sparse grid spaces in n dimensions with dimension O(2 2J J 2n 2 ) and the sparse grid space in 2n dimensions with dimension O(2 J J 2n ), respectively Note that the dimension of the sparse grid spaces compares favorably with the dimension of the full grid space n the rest of the paper we assume l, t, p, q 0 and l + t s (p + q) Then the approximation error when using the full grid approximation space induced by the index set (3) is (32) (A A J )u H s C 2 (s l t)j K H t,l u H p,q for u Hp,q, K Ht,l, compare [9] From the results in [9] we also obtain the following Theorem Theorem We assume that the norm equivalences (22) and (23) hold Let the parameters r, r from (22) and (23) fulfill r < s < r and p [0, r), p + q [0, r), t [0, r), t + l [0, r) Let the index set J (s, t, l, p, q) be defined by (33) J(s,t,l,p,q) := (j, k) N 2n 0 : s j + l j, k + t j, k + q k + p k (l + t s)j} We denote the corresponding integral operator by A J (s,t,l,p,q) Moreover, let (A A J (s,t,l,p,q))u H s, u H p,q, K Ht,l Then (34) (A A J (s,t,l,p,q))u H s C 2 (s l t)j K H t,l u H p,q That is, the order of approximation of the full grid approximation space is preserved when using the approximation space V J (s, t, l, p, q) induced by the index set (33) 4
5 The number of elements in the approximation V J (s,t,l,p,q) space is given by (35) O(2 J ) for l s < 0, s l t nj dim(v J (s,t,l,p,q)) = O(2 s l nt ) for 0 < l s < q + np, O(2 2nJ s l t s l 2nt q np ) for 0 q + np < l s n the extremal cases l s = 0 and 0 < l s = q + np there appear additional logarithmic terms These estimates follow directly from general estimates in [] Combining (34) with (35) we see that the number of non-zero elements N in the stiffness matrix A J (s,t,l,p,q) needed to obtain (A A J (s,t,l,p,q))u H s ɛ, ɛ R +, is bounded by (36) O(ɛ s l t ) for l s < 0, s l N O(ɛ n t ) for 0 < l s < q + np, 2 s l q p 2t O(ɛ n ) for 0 q + np < l s Again, in the extremal cases l s = 0 and 0 < l s = q + np there appear additional logarithmic terms These estimates need to be compared with (37) N = O(ɛ 2n s l t ) for the full grid approximation space Hence, for a kernel with dominating ed smoothness, ie t > 0, we have the following results: n the case l s < 0 the number of non-zero elements in the resulting stiffness matrix does not depend on the dimension n n the case 0 l s < q + np there appears a slight n-dependence, which is further increasing in the case 0 q + np < l s n the case t = 0, ie for a kernel with isotropic smoothness, the above approximation scheme in most cases still significantly reduces the number of non-zero entries in the stiffness matrix compared to the full grid scheme Consider for example the case t = p = 0, 0 < l s < q, ie, the case of fully isotropic smoothness Then the dependence on the dimension n of the number of non-zero elements is reduced from 2n in the full grid case to n Therefore, in many cases the discretization scheme induced by the index set (33) renders additional compression schemes unnecessary Only in the case 0 q + np < l s there can be a relatively strong dependence on the dimension n reflected by the exponent 2n in the third estimates of (35) and (36) Here additional compression may be worthwhile The rest of the paper is concerned with the question whether the dependence on n in this case can be reduced significantly via additional compression 4 Compression An additional reduction of the cost of evaluating (3) approximately may be obtained by exploiting the decay properties of the kernel K in physical space The idea is again to drop small entries in the stiffness matrix without destroying the order of approximation Here one usually assumes specific decay properties of the Schwarz-kernel of the pseudodifferential operator under consideration For singular kernels these are typically of the form (4) C = C αβ < : x α y β K(x, y) C x y f(α,β), f(α, β) > 0, x y 5
6 Here f is typically of the form f(α, β) = α + β + m with some m R + Then estimates of the size of the entries a lm of the Galerkin stiffness matrix are obtained by expansions of the Schwarz-kernel in a polynomial basis together with the cancellation properties of the primal wavelets A Taylor expansion of the kernel together with the cancellation properties of the primal wavelets together with (4) shows (42) Aψ jl, ψ km L 2 C 2 l,k (N+ 3 2 ) dist(supp(ψ jl ), supp(ψ km )) f(n+,n+), for dist(supp(ψ jl ), supp(ψ km )) > 0 Given an additional index set τ comp with A comp τ we can replace the stiffness matrix A by (A comp (A ) lm for (j, k) and (l, m) τ comp, ) lm = 0 for (j, k) and (l, m) τ comp The corresponding integral operator is given by A comp u(x) := K comp (x, y)u(y)dy with K comp (x, y) := (j,k) (l,m) τ comp a lm ψjl (x) ψ km (y) The remaining task is to define the τ comp in such a way that the order of approximation does not deteriorate because of this compression The following result tells us how to define the index sets τ comp in order to balance the overall complexity and the error after compression Theorem 2 Let the assumptions of Theorem hold Moreover, let the decay property (42) hold Let (A J (s,t,l,p,q) A comp J (s,t,l,p,q) )u Hs Here J (s, t, l, p, q) is again given by (33) For (j, k) J (s, t, l, p, q) let We define n (t+l s)j+s j q k p k (N+) j,k B := J f(n+,n+) 2 f(n+,n+) (43) τ comp := (l, m) τ : dist(supp(ψ jl ), supp(ψ km )) B } Then (A A comp J (s,t,l,p,q) )u H s C 2(s l t)j u H p,q Proof: For shortness we use the abbreviations for J (s, t, l, p, q), A for A J (s,t,l,p,q) and analogously A comp for A comp J (s,t,l,p,q) Let u := u, ψ L 2 be the unique wavelet coefficients of u L 2, ie, u = j=0 k τ j u ψ We use the biorthogonality 6
7 between the dual and the primal wavelets and the stability of the dual wavelets in H s and the primal wavelets in H p,q, ie the validity of (22) and (23) t holds (A A comp )u 2 H = a s lm u km ψjl (x) 2 H s = j (23) j = j j l τ j k:(j,k) l τ j 2 2s j l τ j 2 2s j l τ j 2 2s j C (22) (j,k) (j,k) m:(l,m) τ \τ comp k:(j,k) k:(j,k) k:(j,k) (l,m) τ \τ comp k:(l,m) τ \τ comp m:(l,m) τ \τ comp m:(l,m) τ \τ comp ( k N n (l,m) τ \τ comp a lm u km ψjl (x) 2 H s a lm u km 2 2 q k p k a lm 2 q k +p k u km 2 2q k 2p k a 2 lm m τ k 2 2q k +2p k u 2 km 2 2s j 2q k 2p k a 2 lm u 2 H p,q With the definition of τ comp together with (42) we then obtain (for shortness we use the abbreviation f = f(n +, N + )) (A A comp )u 2 H s C ( (j,k) ) 2 2s j 2q k 2p k 2 2 j,k (N+ 3 2 ) ) 2f dist(supp(ψ jl), supp(ψ km)) u 2 H p,q (l,m) τ, dist(supp(ψ jl ),supp(ψ km ))>B ( C 2 2s j 2q k 2p k 2 2 j,k (N+ 3 2 ) B 2f (j,k) (l,m) τ, dist(supp(ψ jl ),supp(ψ km ))>B ) u 2 H p,q Now we apply (l,m) τ = τ 2 j,k and plug in the definition of B We obtain (A A comp )u 2 H C 2 2s j 2q k 2p k 2 2 j,k (N+) B 2f s u 2 H p,q (j,k) C 2 2(s l t)j J 2n 7 (j,k) u 2 H p,q 2
8 (44) C 2 2(s l t)j u 2 H p,q Now we combine (34) from Theorem with (44) and obtain (A A comp )u H s (A A )u H s + (A A comp )u H s C 2 (s l t)j K H t,l u H p,q + C 2(s l t)j u H p,q C 2 (s l t)j u H p,q Theorem 2 shows that the order of approximation of the full grid approximation space is preserved Note that the resulting stiffness matrix A J (s,t,l,p,q) will in general be non-symmetric A symmetric operator may be derived by taking the imum of B and B kj in the definition of τ comp The following Lemma gives an upper estimate for the number of remaining non-zero elements in the compressed stiffness matrix Lemma Let 0 q + np < l s Then the number of non-zero elements N in the stiffness matrix A comp J (s,t,l,p,q) is bounded by N CJ n(2+ n f ) dim ( ) n f V (N+ t l) J (s,t,l,p,q) CJ n(2+ n f ) 2 nj s l t s l 2nt q np 2( n (N+ t l)) f Specifically for f(n +, N + ) = 2n(N + + m), m R +, it holds N O(J n ) 2 nj s l t s l 2nt q np(+ N++m) t+l+m Proof: For shortness we use again the abbreviations f = f(n +, N + ) as well as J := J (s, t, l, p, q) and V J := V J (s,t,l,p,q) for the index set (33) and the corresponding approximation space The number of non-zero elements in the stiffness matrix A comp J is given by N = (43) = (45) (j,k) J (l,m) τ comp (j,k) J (l,m) τ :dist(supp(ψ jl ),supp(ψ km )) B } Now, for every fixed index pair (j, l) the number of pairs (k, m) with is bounded by C n i= dist(supp(ψ jl ), supp(ψ km )) B (, 2 j i + B 2 ki ) = C C µ, ν: µ ν=,,n}, µ ν= n i= (, 2 j i 2 ki + B 2 ki ) (, 2 j i 2 ) (, B ki 2 ) ki i µ 8 i ν
9 (46) C µ, ν: µ ν=,,n}, µ ν= (, 2 j i 2 ) (, B ki 2 ) ki i µ i ν Because of (47) and similarly (48) i µ l ν (, 2 ji 2 k i ) C µ:µ µ} 2 i µ (ki ji) (, B 2 k ) C l ν:ν ν} 2 l ν k l B ν we obtain N = (j,k) J (l,m) τ :dist(supp(ψ jl ),supp(ψ km )) B } (46) C (j,k) J 2 j (47),(48) C C (j,k) J 2 j (j,k) J CJ 2n (j,k) J µ, ν: µ ν=,,n}, µ ν= µ,ν:µ ν,,n}, µ ν= µ,ν:µ ν,,n}, µ ν= µ,ν:µ ν,,n}, µ ν= Now we plug in the definition of B and obtain with K := (j,k) J µ,ν:µ ν,,n}, µ ν= (, 2 ji 2 k ) i i µ l ν 2 i µ (k i j i ) )2 } i ν k i B ν } 2 ( i µ ν ki)+( i µ ji) B ν N CJ n(2+ n f ) 2 K k i + i µ ν i µ } 2 ( i µ ν ki)+( i µ ji) B ν j i (, B 2 k l + ν } f [(t + l s)j + s j q k p k (N + ) j, k ] The imum in µ, ν is obtained for µ =, ν =,, n} Hence K (j,k) J j, k + n f [(t + l s)j + s j q k p k (N + ) j, k ] 9 ) }
10 ( nf ) (N + t l) j, k (j,k) J + n } f [( t l) j, k + s j q k p k + (t + l s)j] (j,k) J ( nf (N + t l) ) j, k + n } f [s j q k p k t j, k l j, k + (t + l s)j] Using the definition of J in (33) we obtain K (j,k) J ( nf (N + t l) ) j, k and hence together with the third inequality in (35) N CJ n(2+ n f ) dim(v J ) n f (N+ t l) CJ n(2+ n f ) 2 nj s l t s l 2nt q np 2( n (N+ t l)) f Lemma shows that for N + t + l there will be no compression effect However, for N + > t+l there will be a significant reduction of the number of non-zero entries in the stiffness matrix Combining Theorem 2 with the estimate for the number of non-zero elements in the stiffness matrix in Lemma we obtain the following result Theorem 3 number of non-zero elements in the stiffness matrix A comp J (s,t,l,p,q) Under the assumptions of Theorem 2 and 0 q + np < l s the needed to obtain (A A comp J (s,t,l,p,q) )u H s ɛ, ɛ R+, is bounded by (49) ( N O ɛ 2 s l q n ) ( n p 2t f (N+ t l)) ln(ɛ ) n2 Specifically for f(n +, N + ) = 2n(N + + m) with m R + it holds (40) ( N O ɛ s l q p 2t n (+ t+l+m N++m ) ln(ɛ ) n ) A comparison of (49) with the last estimate in (36) shows that compression reduces the quadratic exponential factor in the right hand side of (36) to 2 2 n f (N + t l)) Additionally, (40) shows that for t + l + m/n + + m 0 we obtain optimal compression up to a logarithmic factor Note however, that generalizations to more general geometries are not as easy as in the isotropic case REFERENCES [] K Babenko, Approximation by trigonometric polynomials in a certain class of periodic functions of several variables, Dokl Akad Nauk SSSR 32 (960),
11 [2] G Baszenski, N-th order polynomial spline blending, in Multivariate approximation theory, W Schemp, K Zeller (eds), Birkhäuser, Basel, 985, [3] G Beylkin, R Coifman, V Rokhlin, Fast wavelet transforms and numerical algorithms, Comm Pure and Appl Math 44 (99), 4-83 [4] H-J Bungartz, M Griebel, D Röschke, C Zenger, Pointwise convergence of the combination technique for Laplace s equation, East-West J of Numerical Mathematics (2) (994), 2-45 [5] W Dahmen, Stability of Multiscale Transformations, Journal of Fourier Analysis and Applications 2 (996), [6] W Dahmen, Wavelet and Multiscale Methods for Operator equations, Acta Numerica (997), [7] W Dahmen, S Prößdorf, R Schneider, Wavelet approximation methods for pseudodifferential equations : Matrix compression and fast solution, Advances in computational Mathematics (993), [8] F-J Delvos, d-variate Boolean interpolation, J Approx Theory 34 (982), 99-4 [9] R A DeVore, B Jawerth, V Popov, Compression of Wavelet decompositions, Amer J Math 4 (992), [0] Din Zung, The approximation of classes of periodic functions of many variables, Russion Math Surveys 38 (983), 7-8 [] Din Zung, Number of integral points in a certain set and the approximation of functions of several variables, Math Notes of the Academy of Sciences of the USSR 36(4) (984), [2] K Frank, S Heinrich, S Pereverzev, nformation complexity of multivariate Fredholm equations in sobolev classes, J Complexity 2 (996), 7-34 [3] E M Galeev, Approximation by fourier sums of functions with several bounded derivatives, Math Notes 23 (978), 09-7 [4] W J Gordon, Blending-function methods for bivariate and multivariate interpolation and approximation, SAM J Numerical Analysis 8 (97), [5] M Griebel, S Knapek, Optimized tensor-product approximation spaces, Constructive Approximation 6(4) (2000), [6] M Griebel, P Oswald, T Schiekofer, Sparse grids for boundary integral equations, Numer Math 83(2) (999), [7] R Hochmuth, S Knapek, G Zumbusch, Tensor-products of Sobolev spaces and Applications, SFB-256 Bericht Nr 685, Universität Bonn, 2000 [8] S Knapek, Hyperbolic cross approximation of integral operators with smooth kernel, SFB-256 Bericht Nr 665, Universität Bonn, 2000 [9] S Knapek, F Koster ntegral operators on sparse grids, SAM J Numer Anal 39(5) (2002), [20] Y Meyer, R Coifman Wavelets: Calderon-Zymund and multilinear operators, Cambridge studies in advanced mathematics 48, Cambridge University Press (997) [2] P Oswald, On discrete norm estimates related to multilevel preconditioners in the finite element method, in: Constructive Theory of Functions, KG vanov, P Petrushev, B Sendov, (eds), Proc nt Conf Varna, 99, Bulg Acad Sci, Sofia, 992, [22] S V Pereverzev, On the complexity of finding solutions of Fredholm equations of the second kind with differentiable kernels,, Ukrain Mat Zh 40 (988), 84-9 [23] R Schneider, Multiskalen- und Wavelet-Kompression: Analysisbasierte Methoden zur effizienten Lösung großer vollbesetzter Gleichungssysteme, Advances in Numerical Analysis, Teubner, Stuttgart, 998 [24] S A Smolyak, Quadrature and interpolation formulas for tensor products of certain classes of functions, Dokl Akad Nauk SSSR 4 (963), [25] V N Temlyakov, Approximation of periodic functions, Nova Science, New York, 993 [26] J F Traub, G W Wasilkowski, H Wozniakowski, nformation-based Complexity, Academic Press, New York, 988 [27] C Zenger, Sparse grids, Parallel Algorithms for Partial Differential Equations, Notes on Num Fluid Mech 3, W Hackbusch (ed), Vieweg, Braunschweig, 99
On the convergence of the combination technique
Wegelerstraße 6 53115 Bonn Germany phone +49 228 73-3427 fax +49 228 73-7527 www.ins.uni-bonn.de M. Griebel, H. Harbrecht On the convergence of the combination technique INS Preprint No. 1304 2013 On
More informationComputation of operators in wavelet coordinates
Computation of operators in wavelet coordinates Tsogtgerel Gantumur and Rob Stevenson Department of Mathematics Utrecht University Tsogtgerel Gantumur - Computation of operators in wavelet coordinates
More informationDISCRETE CDF 9/7 WAVELET TRANSFORM FOR FINITE-LENGTH SIGNALS
DISCRETE CDF 9/7 WAVELET TRANSFORM FOR FINITE-LENGTH SIGNALS D. Černá, V. Finěk Department of Mathematics and Didactics of Mathematics, Technical University in Liberec Abstract Wavelets and a discrete
More informationSchemes. Philipp Keding AT15, San Antonio, TX. Quarklet Frames in Adaptive Numerical. Schemes. Philipp Keding. Philipps-University Marburg
AT15, San Antonio, TX 5-23-2016 Joint work with Stephan Dahlke and Thorsten Raasch Let H be a Hilbert space with its dual H. We are interested in the solution u H of the operator equation Lu = f, with
More informationN-Widths and ε-dimensions for high-dimensional approximations
N-Widths and ε-dimensions for high-dimensional approximations Dinh Dũng a, Tino Ullrich b a Vietnam National University, Hanoi, Information Technology Institute 144 Xuan Thuy, Hanoi, Vietnam dinhzung@gmail.com
More informationConstruction of wavelets. Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam
Construction of wavelets Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam Contents Stability of biorthogonal wavelets. Examples on IR, (0, 1), and (0, 1) n. General domains
More informationL. Levaggi A. Tabacco WAVELETS ON THE INTERVAL AND RELATED TOPICS
Rend. Sem. Mat. Univ. Pol. Torino Vol. 57, 1999) L. Levaggi A. Tabacco WAVELETS ON THE INTERVAL AND RELATED TOPICS Abstract. We use an abstract framework to obtain a multilevel decomposition of a variety
More informationMax-Planck-Institut fur Mathematik in den Naturwissenschaften Leipzig H 2 -matrix approximation of integral operators by interpolation by Wolfgang Hackbusch and Steen Borm Preprint no.: 04 200 H 2 -Matrix
More informationCubature formulas for function spaces with moderate smoothness
Cubature formulas for function spaces with moderate smoothness Michael Gnewuch, René Lindloh, Reinhold Schneider, and Anand Srivastav Institut für Informatik Christian-Albrechts-Universität zu Kiel Christian-Albrechts-Platz
More informationNonlinear Approximation Schemes Associated With Nonseparable Wavelet Bi-frames
Nonlinear Approximation Schemes Associated With Nonseparable Wavelet Bi-frames Martin Ehler September 26, 2008 Abstract In the present paper, we study nonlinear approximation properties of multivariate
More informationApproximation of High-Dimensional Rank One Tensors
Approximation of High-Dimensional Rank One Tensors Markus Bachmayr, Wolfgang Dahmen, Ronald DeVore, and Lars Grasedyck March 14, 2013 Abstract Many real world problems are high-dimensional in that their
More informationMatrix assembly by low rank tensor approximation
Matrix assembly by low rank tensor approximation Felix Scholz 13.02.2017 References Angelos Mantzaflaris, Bert Juettler, Boris Khoromskij, and Ulrich Langer. Matrix generation in isogeometric analysis
More informationSparse Tensor Galerkin Discretizations for First Order Transport Problems
Sparse Tensor Galerkin Discretizations for First Order Transport Problems Ch. Schwab R. Hiptmair, E. Fonn, K. Grella, G. Widmer ETH Zürich, Seminar for Applied Mathematics IMA WS Novel Discretization Methods
More informationWavelet bases on a triangle
Wavelet bases on a triangle N. Ajmi, A. Jouini and P.G. Lemarié-Rieusset Abdellatif.jouini@fst.rnu.tn neyla.ajmi@laposte.net Département de Mathématiques Faculté des Sciences de Tunis 1060 Campus-Tunis,
More informationSpace-Time Adaptive Wavelet Methods for Parabolic Evolution Problems
Space-Time Adaptive Wavelet Methods for Parabolic Evolution Problems Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam joint work with Christoph Schwab (ETH, Zürich), Nabi
More informationError estimates for two-dimensional cross approximation
Journal of Approximation Theory 162 (21) 1685 17 www.elsevier.com/locate/jat Error estimates for two-dimensional cross approximation Jan Schneider Max Planck Institute for Mathematics in the Sciences,
More informationINDUSTRIAL MATHEMATICS INSTITUTE. B.S. Kashin and V.N. Temlyakov. IMI Preprint Series. Department of Mathematics University of South Carolina
INDUSTRIAL MATHEMATICS INSTITUTE 2007:08 A remark on compressed sensing B.S. Kashin and V.N. Temlyakov IMI Preprint Series Department of Mathematics University of South Carolina A remark on compressed
More informationMRA Frame Wavelets with Certain Regularities Associated with the Refinable Generators of Shift Invariant Spaces
Chapter 6 MRA Frame Wavelets with Certain Regularities Associated with the Refinable Generators of Shift Invariant Spaces Tian-Xiao He Department of Mathematics and Computer Science Illinois Wesleyan University,
More informationBesov regularity for operator equations on patchwise smooth manifolds
on patchwise smooth manifolds Markus Weimar Philipps-University Marburg Joint work with Stephan Dahlke (PU Marburg) Mecklenburger Workshop Approximationsmethoden und schnelle Algorithmen Hasenwinkel, March
More informationOn the Numerical Evaluation of Fractional Sobolev Norms. Carsten Burstedde. no. 268
On the Numerical Evaluation of Fractional Sobolev Norms Carsten Burstedde no. 268 Diese Arbeit ist mit Unterstützung des von der Deutschen Forschungsgemeinschaft getragenen Sonderforschungsbereiches 611
More informationON THE INFLUENCE OF THE WAVENUMBER ON COMPRESSION IN A WAVELET BOUNDARY ELEMENT METHOD FOR THE HELMHOLTZ EQUATION
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 4, Number 1, Pages 48 62 c 2007 Institute for Scientific Computing and Information ON THE INFLUENCE OF THE WAVENUMBER ON COMPRESSION IN A
More informationfor compression of Boundary Integral Operators. Steven Paul Nixon B.Sc.
Theory and Applications of the Multiwavelets for compression of Boundary Integral Operators. Steven Paul Nixon B.Sc. Institute for Materials Research School of Computing, Science & Engineering, University
More informationWavelets, Multiresolution Analysis and Fast Numerical Algorithms. G. Beylkin 1
A draft of INRIA lectures, May 1991 Wavelets, Multiresolution Analysis Fast Numerical Algorithms G. Beylkin 1 I Introduction These lectures are devoted to fast numerical algorithms their relation to a
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 informationOn the fast algorithm for multiplication of functions in the wavelet bases
Published in Proceedings of the International Conference Wavelets and Applications, Toulouse, 1992; Y. Meyer and S. Roques, edt., Editions Frontieres, 1993 On the fast algorithm for multiplication of functions
More informationThe tensor algebra of power series spaces
The tensor algebra of power series spaces Dietmar Vogt Abstract The linear isomorphism type of the tensor algebra T (E) of Fréchet spaces and, in particular, of power series spaces is studied. While for
More informationAALBORG UNIVERSITY. Compactly supported curvelet type systems. Kenneth N. Rasmussen and Morten Nielsen. R November 2010
AALBORG UNIVERSITY Compactly supported curvelet type systems by Kenneth N Rasmussen and Morten Nielsen R-2010-16 November 2010 Department of Mathematical Sciences Aalborg University Fredrik Bajers Vej
More informationSparse Grids in a Nutshell
Sparse Grids in a Nutshell Jochen Garcke Abstract The technique of sparse grids allows to overcome the curse of dimensionality, which prevents the use of classical numerical discretization schemes in more
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationAn Introduction to Hierachical (H ) Rank and TT Rank of Tensors with Examples
An Introduction to Hierachical (H ) Rank and TT Rank of Tensors with Examples Lars Grasedyck and Wolfgang Hackbusch Bericht Nr. 329 August 2011 Key words: MSC: hierarchical Tucker tensor rank tensor approximation
More informationUniversal Scalable Robust Solvers from Computational Information Games and fast eigenspace adapted Multiresolution Analysis.
Universal Scalable Robust Solvers from Computational Information Games and fast eigenspace adapted Multiresolution Analysis Houman Owhadi Joint work with Clint Scovel IPAM Apr 3, 2017 DARPA EQUiPS / AFOSR
More informationWavelet Bases of the Interval: A New Approach
Int. Journal of Math. Analysis, Vol. 1, 2007, no. 21, 1019-1030 Wavelet Bases of the Interval: A New Approach Khaled Melkemi Department of Mathematics University of Biskra, Algeria kmelkemi@yahoo.fr Zouhir
More informationThe adaptive tensor product wavelet scheme: Sparse matrices and the application to singularly perturbed problems
IMA Journal of Numerical Analysis () Page of 8 doi:.93/imanum/drn The adaptive tensor product wavelet scheme: Sparse matrices and the application to singularly perturbed problems NABI CHEGINI AND ROB STEVENSON
More informationA Tutorial on Wavelets and their Applications. Martin J. Mohlenkamp
A Tutorial on Wavelets and their Applications Martin J. Mohlenkamp University of Colorado at Boulder Department of Applied Mathematics mjm@colorado.edu This tutorial is designed for people with little
More informationSMALL SUPPORT SPLINE RIESZ WAVELETS IN LOW DIMENSIONS BIN HAN, QUN MO, AND ZUOWEI SHEN
SMALL SUPPORT SPLINE RIESZ WAVELETS IN LOW DIMENSIONS BIN HAN, QUN MO, AND ZUOWEI SHEN Abstract. In [B. Han and Z. Shen, SIAM J. Math. Anal., 38 (006), 530 556], a family of univariate short support Riesz
More informationA New Scheme for the Tensor Representation
J Fourier Anal Appl (2009) 15: 706 722 DOI 10.1007/s00041-009-9094-9 A New Scheme for the Tensor Representation W. Hackbusch S. Kühn Received: 18 December 2008 / Revised: 29 June 2009 / Published online:
More informationQuarkonial frames of wavelet type - Stability, approximation and compression properties
Quarkonial frames of wavelet type - Stability, approximation and compression properties Stephan Dahlke 1 Peter Oswald 2 Thorsten Raasch 3 ESI Workshop Wavelet methods in scientific computing Vienna, November
More informationPapers Published. Edward Neuman. Department of Mathematics, Southern Illinois University, Carbondale. and
Papers Published Edward Neuman Department of Mathematics, Southern Illinois University, Carbondale and Mathematical Research Institute, 144 Hawthorn Hollow, Carbondale, IL, 62903, USA 134. On Yang means
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 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 informationList of the scientific publications of prof. D.Sc. Kamen Ganchev Ivanov, Institute for Mathematics and Informatics BAS (August 2016)
List of the scientific publications of prof. D.Sc. Kamen Ganchev Ivanov, Institute for Mathematics and Informatics BAS (August 2016) 1. K. G. Ivanov. New estimates of errors of quadrature formulae, formulae
More informationB. Vedel Joint work with P.Abry (ENS Lyon), S.Roux (ENS Lyon), M. Clausel LJK, Grenoble),
Hyperbolic wavelet analysis of textures : global regularity and multifractal formalism B. Vedel Joint work with P.Abry (ENS Lyon), S.Roux (ENS Lyon), M. Clausel LJK, Grenoble), S.Jaffard(Créteil) Harmonic
More informationSharp estimates for a class of hyperbolic pseudo-differential equations
Results in Math., 41 (2002), 361-368. Sharp estimates for a class of hyperbolic pseudo-differential equations Michael Ruzhansky Abstract In this paper we consider the Cauchy problem for a class of hyperbolic
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 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 informationCOMPACTLY SUPPORTED ORTHONORMAL COMPLEX WAVELETS WITH DILATION 4 AND SYMMETRY
COMPACTLY SUPPORTED ORTHONORMAL COMPLEX WAVELETS WITH DILATION 4 AND SYMMETRY BIN HAN AND HUI JI Abstract. In this paper, we provide a family of compactly supported orthonormal complex wavelets with dilation
More informationDimension Adaptive Tensor Product Quadrature
Dimension Adaptive Tensor Product Quadrature Thomas Gerstner and Michael Griebel Department for Applied Mathematics University of Bonn {gerstner,griebel}@iam.uni-bonn.de Abstract We consider the numerical
More informationAn Adaptive Wavelet Viscosity Method for Hyperbolic Conservation Laws
An Adaptive Wavelet Viscosity Method for Hyperbolic Conservation Laws Daniel Castaño Díez, 1 Max Gunzburger, 2 Angela Kunoth 3 1 EMBL Heidelberg, Meyerhofstraße 1, 69117 Heidelberg, Germany 2 School of
More informationFitting multidimensional data using gradient penalties and combination techniques
1 Fitting multidimensional data using gradient penalties and combination techniques Jochen Garcke and Markus Hegland Mathematical Sciences Institute Australian National University Canberra ACT 0200 jochen.garcke,markus.hegland@anu.edu.au
More informationADAPTIVE MULTISCALE MOMENT METHOD FOR SOLVING TWO-DIMENSIONAL FREDHOLM INTEGRAL EQUATION OF THE FIRST KIND
Progress In Electromagnetics Research, PIER 21, 173 201, 1999 ADAPTIVE MULTISCALE MOMENT METHOD FOR SOLVING TWO-DIMENSIONAL FREDHOLM INTEGRAL EQUATION OF THE FIRST KIND C. Su Dept. of Applied Math. Northwestern
More information524 Jan-Olov Stromberg practice, this imposes strong restrictions both on N and on d. The current state of approximation theory is essentially useless
Doc. Math. J. DMV 523 Computation with Wavelets in Higher Dimensions Jan-Olov Stromberg 1 Abstract. In dimension d, a lattice grid of size N has N d points. The representation of a function by, for instance,
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 informationApplications of Polyspline Wavelets to Astronomical Image Analysis
VIRTUAL OBSERVATORY: Plate Content Digitization, Archive Mining & Image Sequence Processing edited by M. Tsvetkov, V. Golev, F. Murtagh, and R. Molina, Heron Press, Sofia, 25 Applications of Polyspline
More informationConstruction of Multivariate Compactly Supported Orthonormal Wavelets
Construction of Multivariate Compactly Supported Orthonormal Wavelets Ming-Jun Lai Department of Mathematics The University of Georgia Athens, GA 30602 April 30, 2004 Dedicated to Professor Charles A.
More informationA Recursive Trust-Region Method for Non-Convex Constrained Minimization
A Recursive Trust-Region Method for Non-Convex Constrained Minimization Christian Groß 1 and Rolf Krause 1 Institute for Numerical Simulation, University of Bonn. {gross,krause}@ins.uni-bonn.de 1 Introduction
More informationFast Sparse Spectral Methods for Higher Dimensional PDEs
Fast Sparse Spectral Methods for Higher Dimensional PDEs Jie Shen Purdue University Collaborators: Li-Lian Wang, Haijun Yu and Alexander Alekseenko Research supported by AFOSR and NSF ICERM workshop, June
More informationVandermonde Matrices for Intersection Points of Curves
This preprint is in final form. It appeared in: Jaén J. Approx. (009) 1, 67 81 Vandermonde Matrices for Intersection Points of Curves Hakop Hakopian, Kurt Jetter and Georg Zimmermann Abstract We reconsider
More informationGreedy bases are best for m-term approximation
Greedy bases are best for m-term approximation Witold Bednorz Institute of Mathematics, Warsaw University Abstract We study the approximation of a subset K in a Banach space X by choosing first basis B
More informationWe simply compute: for v = x i e i, bilinearity of B implies that Q B (v) = B(v, v) is given by xi x j B(e i, e j ) =
Math 395. Quadratic spaces over R 1. Algebraic preliminaries Let V be a vector space over a field F. Recall that a quadratic form on V is a map Q : V F such that Q(cv) = c 2 Q(v) for all v V and c F, and
More informationSpline Interpolation on Sparse Grids
Spline Interpolation on Sparse Grids Winfried Sickel a and Tino Ullrich b a Mathematical Institute, Friedrich-Schiller-University, D-07737 Jena, Germany; b Hausdorff-Center for Mathematics, 53115 Bonn,
More information2014:05 Incremental Greedy Algorithm and its Applications in Numerical Integration. V. Temlyakov
INTERDISCIPLINARY MATHEMATICS INSTITUTE 2014:05 Incremental Greedy Algorithm and its Applications in Numerical Integration V. Temlyakov IMI PREPRINT SERIES COLLEGE OF ARTS AND SCIENCES UNIVERSITY OF SOUTH
More informationHow much should we rely on Besov spaces as a framework for the mathematical study of images?
How much should we rely on Besov spaces as a framework for the mathematical study of images? C. Sinan Güntürk Princeton University, Program in Applied and Computational Mathematics Abstract Relations between
More informationRiesz Bases of Wavelets and Applications to Numerical Solutions of Elliptic Equations
Riesz Bases of Wavelets and Applications to Numerical Solutions of Elliptic Equations Rong-Qing Jia and Wei Zhao Department of Mathematical and Statistical Sciences University of Alberta Edmonton, Canada
More information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationConstruction of Biorthogonal B-spline Type Wavelet Sequences with Certain Regularities
Illinois Wesleyan University From the SelectedWorks of Tian-Xiao He 007 Construction of Biorthogonal B-spline Type Wavelet Sequences with Certain Regularities Tian-Xiao He, Illinois Wesleyan University
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationNonstationary Subdivision Schemes and Totally Positive Refinable Functions
Nonstationary Subdivision Schemes and Totally Positive Refinable Functions Laura Gori and Francesca Pitolli December, 2007 Abstract In this paper we construct a class of totally positive refinable functions,
More informationarxiv:math/ v1 [math.fa] 5 Aug 2005
arxiv:math/0508104v1 [math.fa] 5 Aug 2005 G-frames and G-Riesz Bases Wenchang Sun Department of Mathematics and LPMC, Nankai University, Tianjin 300071, China Email: sunwch@nankai.edu.cn June 28, 2005
More informationNEW CONSTRUCTIONS OF PIECEWISE-CONSTANT WAVELETS
NEW CONSTRUCTIONS OF PIECEWISE-CONSTANT WAVELETS YOUNGMI HUR AND AMOS RON Abstract. The classical Haar wavelet system of L 2 (R n ) is commonly considered to be very local in space. We introduce and study
More informationAdaptive Wavelet Algorithms
Adaptive Wavelet Algorithms for solving operator equations Tsogtgerel Gantumur August 2006 Contents Notations and acronyms vii 1 Introduction 1 1.1 Background.............................. 1 1.2 Thesis
More informationMultiscale Frame-based Kernels for Image Registration
Multiscale Frame-based Kernels for Image Registration Ming Zhen, Tan National University of Singapore 22 July, 16 Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image
More informationEXAMPLES OF REFINABLE COMPONENTWISE POLYNOMIALS
EXAMPLES OF REFINABLE COMPONENTWISE POLYNOMIALS NING BI, BIN HAN, AND ZUOWEI SHEN Abstract. This short note presents four examples of compactly supported symmetric refinable componentwise polynomial functions:
More informationQUADRATURES INVOLVING POLYNOMIALS AND DAUBECHIES' WAVELETS *
QUADRATURES INVOLVING POLYNOMIALS AND DAUBECHIES' WAVELETS * WEI-CHANG SHANN AND JANN-CHANG YAN Abstract. Scaling equations are used to derive formulae of quadratures involving polynomials and scaling/wavelet
More informationAdaptive Control with Multiresolution Bases
Adaptive Control with Multiresolution Bases Christophe P. Bernard Centre de Mathématiques Appliquées Ecole Polytechnique 91128 Palaiseau cedex, FRANCE bernard@cmapx.polytechnique.fr Jean-Jacques E. Slotine
More informationA Nodal Spline Collocation Method for the Solution of Cauchy Singular Integral Equations 1
European Society of Computational Methods in Sciences and Engineering (ESCMSE) Journal of Numerical Analysis, Industrial and Applied Mathematics (JNAIAM) vol. 3, no. 3-4, 2008, pp. 211-220 ISSN 1790 8140
More informationAdaptive Wavelet BEM for Boundary Integral Equations: Theory and Numerical Experiments
Adaptive Wavelet BEM for Boundary Integral Equations: Theory and Numerical Experiments S. Dahlke, H. Harbrecht, M. Utzinger, M. Weimar October 7, 06 Abstract In this paper, we are concerned with the numerical
More informationAdaptive sparse grids
ANZIAM J. 44 (E) ppc335 C353, 2003 C335 Adaptive sparse grids M. Hegland (Received 1 June 2001) Abstract Sparse grids, as studied by Zenger and Griebel in the last 10 years have been very successful in
More informationDecoupling of modes for the elastic wave equation in media of limited smoothness
Decoupling of modes for the elastic wave equation in media of limited smoothness Valeriy Brytik 1 Maarten V. de Hoop 1 Hart F. Smith Gunther Uhlmann September, 1 1 Center for Computational and Applied
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 informationQuadrature Prefilters for the Discrete Wavelet Transform. Bruce R. Johnson. James L. Kinsey. Abstract
Quadrature Prefilters for the Discrete Wavelet Transform Bruce R. Johnson James L. Kinsey Abstract Discrepancies between the Discrete Wavelet Transform and the coefficients of the Wavelet Series are known
More informationDomain decomposition methods via boundary integral equations
Domain decomposition methods via boundary integral equations G. C. Hsiao a O. Steinbach b W. L. Wendland b a Department of Mathematical Sciences, University of Delaware, Newark, Delaware 19716, USA. E
More informationOrthogonal Polynomials, Quadratures & Sparse-Grid Methods for Probability Integrals
1/33 Orthogonal Polynomials, Quadratures & Sparse-Grid Methods for Probability Integrals Dr. Abebe Geletu May, 2010 Technische Universität Ilmenau, Institut für Automatisierungs- und Systemtechnik Fachgebiet
More informationThe Mortar Wavelet Method Silvia Bertoluzza Valerie Perrier y October 29, 1999 Abstract This paper deals with the construction of wavelet approximatio
The Mortar Wavelet Method Silvia Bertoluzza Valerie Perrier y October 9, 1999 Abstract This paper deals with the construction of wavelet approximation spaces, in the framework of the Mortar method. We
More informationTensor Sparsity and Near-Minimal Rank Approximation for High-Dimensional PDEs
Tensor Sparsity and Near-Minimal Rank Approximation for High-Dimensional PDEs Wolfgang Dahmen, RWTH Aachen Collaborators: Markus Bachmayr, Ron DeVore, Lars Grasedyck, Endre Süli Paris, Oct.11, 2013 W.
More informationTHE DUAL FORM OF THE APPROXIMATION PROPERTY FOR A BANACH SPACE AND A SUBSPACE. In memory of A. Pe lczyński
THE DUAL FORM OF THE APPROXIMATION PROPERTY FOR A BANACH SPACE AND A SUBSPACE T. FIGIEL AND W. B. JOHNSON Abstract. Given a Banach space X and a subspace Y, the pair (X, Y ) is said to have the approximation
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 informationClassical Fourier Analysis
Loukas Grafakos Classical Fourier Analysis Second Edition 4y Springer 1 IP Spaces and Interpolation 1 1.1 V and Weak IP 1 1.1.1 The Distribution Function 2 1.1.2 Convergence in Measure 5 1.1.3 A First
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 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 informationAlgorithms for Scientific Computing
Algorithms for Scientific Computing Hierarchical Methods and Sparse Grids d-dimensional Hierarchical Basis Michael Bader Technical University of Munich Summer 208 Intermezzo/ Big Picture : Archimedes Quadrature
More informationOn pointwise estimates for maximal and singular integral operators by A.K. LERNER (Odessa)
On pointwise estimates for maximal and singular integral operators by A.K. LERNER (Odessa) Abstract. We prove two pointwise estimates relating some classical maximal and singular integral operators. In
More informationBAND-LIMITED REFINABLE FUNCTIONS FOR WAVELETS AND FRAMELETS
BAND-LIMITED REFINABLE FUNCTIONS FOR WAVELETS AND FRAMELETS WEIQIANG CHEN AND SAY SONG GOH DEPARTMENT OF MATHEMATICS NATIONAL UNIVERSITY OF SINGAPORE 10 KENT RIDGE CRESCENT, SINGAPORE 119260 REPUBLIC OF
More informationWavelets and modular inequalities in variable L p spaces
Wavelets and modular inequalities in variable L p spaces Mitsuo Izuki July 14, 2007 Abstract The aim of this paper is to characterize variable L p spaces L p( ) (R n ) using wavelets with proper smoothness
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 informationTensor networks, TT (Matrix Product States) and Hierarchical Tucker decomposition
Tensor networks, TT (Matrix Product States) and Hierarchical Tucker decomposition R. Schneider (TUB Matheon) John von Neumann Lecture TU Munich, 2012 Setting - Tensors V ν := R n, H d = H := d ν=1 V ν
More informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
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 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 informationA Multiresolution Strategy for Reduction of Elliptic PDEs and Eigenvalue Problems
APPLIED AND COMPUTATIONAL HARMONIC ANALYSIS 5, 129 155 (1998) ARTICLE NO. HA970226 A Multiresolution Strategy for Reduction of Elliptic PDEs and Eigenvalue Problems Gregory Beylkin* and Nicholas Coult
More informationNumerical differentiation by means of Legendre polynomials in the presence of square summable noise
www.oeaw.ac.at Numerical differentiation by means of Legendre polynomials in the presence of square summable noise S. Lu, V. Naumova, S. Pereverzyev RICAM-Report 2012-15 www.ricam.oeaw.ac.at Numerical
More information