Computational scales of Sobolev norms with application to preconditioning

This paper provides a framework for developing computationally efficient multilevel preconditioners and representations for Sobolev norms. Specifically, given a Hilbert space V and a nested sequence of subspaces V 1 C V 2 C... C V, we construct operators which are spectrally equivalent to those of the form A = Σ k μ k (Q k - Q k-1 ). Here μ k , k=1,2,..., are positive numbers and Q k is the orthogonal projector onto V k with Q 0 = 0. We first present abstract results which show when A is spectrally equivalent to a similarly constructed operator A defined in terms of an approximation Q k of Q k , for k = 1,2,... We show that these results lead to efficient preconditioners for discretizations of differential and pscudo-differential operators of positive and negative order. These results extend to sums of operators. For example, singularly perturbed problems such as I -eΔ can be preconditioned uniformly independently of the parameter e. We also show how to precondition an operator which results from Tikhonov regularization of a problem with noisy data. Finally, we describe how the technique provides computationally efficient bounded discrete extensions which have applications to domain decomposition.

[1]  Panayot S. Vassilevski,et al.  Stabilizing the Hierarchical Basis by Approximate Wavelets II: Implementation and Numerical Results , 1998, SIAM J. Sci. Comput..

[2]  Joseph E. Pasciak,et al.  Iterative schemes for nonsymmetric and indefinite elliptic boundary value problems , 1993 .

[3]  J. Pasciak,et al.  Parallel multilevel preconditioners , 1990 .

[4]  Philippe G. Ciarlet,et al.  The finite element method for elliptic problems , 2002, Classics in applied mathematics.

[5]  Wolfgang Dahmen,et al.  Local Decomposition of Refinable Spaces and Wavelets , 1996 .

[6]  Ulrich Tautenhahn,et al.  Error Estimates for Regularization Methods in Hilbert Scales , 1996 .

[7]  X. Saint-Raymond Non-unicité pour certains problèmes de Cauchy complexes non linéaires du premier ordre , 1984 .

[8]  Panayot S. Vassilevski,et al.  Stabilizing the Hierarchical Basis by Approximate Wavelets, I: Theory , 1997, Numer. Linear Algebra Appl..

[9]  Rob Stevenson,et al.  A robust hierarchical basis preconditioner on general meshes , 1997 .

[10]  P. Grisvard Elliptic Problems in Nonsmooth Domains , 1985 .

[11]  Arnd Meyer,et al.  Hierarchical Extension Operators and Local Multigrid Methods in Domain Decomposition Preconditioners , 1994 .

[12]  Peter Oswald,et al.  Multilevel Finite Element Approximation , 1994 .

[13]  Joseph E. Pasciak,et al.  New estimates for multilevel algorithms including the V-cycle , 1993 .

[14]  Panayot S. Vassilevski,et al.  Stabilizing the Hierarchical Basis by Approximate Wavelets, I: Theory , 1997 .

[15]  F. Natterer Error bounds for tikhonov regularization in hilbert scales , 1984 .

[16]  Rob Stevenson Piecewise linear (pre-)wavelets on non-uniform meshes , 1998 .

[17]  Sergey V. Nepomnyaschikh Optimal Multilevel Extension Operators , 1995 .