Matrix Preconditioning Techniques and Applications

There seem to be at least two important issues for which wavelet-like expansions have already proven to work with great success, namely preconditioning linear systems stemming from Galerkin approximations for elliptic problems and compressing full stiffness matrices arising in connection with integral or pseudodifferential operators, to facilitate nearly optimal complexity algorithms for the solution of the corresponding discrete problems.
Wolfgang Dahmen, et al. Multiscale methods for pseudodifferential equations . Recent Advances in Wavelet Analysis (1994)
In the usual FEM setting, Schur complement methods from Chapter 7 perform the best if there is some kind of diagonal dominance. This chapter proposes two related and efficient iterative algorithms based on the wavelet formulation for solving an operator equation with conventional arithmetic. In the new wavelet setting, the stiffness matrix possesses the desirable properties suitable for using the Schur complements. The proposed algorithms utilize the Schur complements recursively; they only differ in how to use coarse levels to solve Schur complements equations. In the first algorithm, we precondition a Schur complement by using coarse levels while in the second we use approximate Schur complements to construct a preconditioner. We believe that our algorithms can be adapted to higher dimensional problems more easily than previous work in the subject. The material is organized in the following
Section 9.1 Introduction
Section 9.2 Wavelets telescopic splitting of an operator
Section 9.3 An exact Schur preconditioner with level-by-level wavelets
Section 9.4 An approximate preconditioner with level-by-level wavelets
Section 9.5 Some analysis and...