Abstract
Preconditioning techniques for dense linear systems arising from singular boundary integral equations are described and analyzed. A particular class of approximate inverse based preconditioners related to the mesh neighbor methods is known to be efficient. This paper shows that it is an operator splitting preconditioner and clusters eigenvalues for the normal equation matrix thus ensuring a fast convergence of the conjugate gradient normal method. Clustering of the eigenvalues of the preconditioned matrix and fast convergence of the generalized minimal residual method are also observed. For the type of problems considered, we demonstrate a crucial connection between two essential features of eigenvalue clustering for a sparse preconditioner - approximate inversion for a small cluster radius and operator splitting for a small cluster size. Experimental results from several boundary integral equations are presented.
| Original language | English |
|---|---|
| Pages (from-to) | 1058-1078 |
| Number of pages | 21 |
| Journal | SIAM Journal on Matrix Analysis and Applications |
| Volume | 22 |
| Issue number | 4 |
| DOIs | |
| Publication status | Published - 1 Jan 2001 |
Keywords
- approximate inversion
- conjugate gradient normal method
- conjugate gradients
- GMRES
- least squares solution
- operator splitting
- preconditioning
- singular boundary elements