Incomplete factorization-based preconditionings for solving the Helmholtz equation

Preconditioning techniques based on incomplete factorization of matrices are investigated, to solve highly indefinite complex-symmetric linear systems. A novel preconditioning is introduced. The real part of the matrix is made positive definite, or less indefinite, by adding properly defined perturbations to the diagonal entries, while the imaginary part is unaltered. The resulting preconditioning matrix, which is obtained by applying standard methods to the perturbed complex matrix, turns out to perform significantly better than classical incomplete factorization schemes. For realistic values of the GMRES restart parameter, spectacular reduction of iteration counts is observed. A theoretical spectral analysis is provided, in which the spectrum of the preconditioner applied to indefinite matrix is related to the spectrum of the same preconditioner applied to a Stieltjes matrix extracted from the indefinite matrix. Results of numerical experiments are reported, which display the efficiency of the new preconditioning. Copyright © 2001 John Wiley & Sons, Ltd.

[1]  Fast iterative solvers for finite element analysis in general and shell analysis in particular , 1996 .

[2]  Chun-Hua Guo,et al.  Incomplete block factorization preconditioning for indefinite elliptic problems , 1999, Numerische Mathematik.

[3]  H. V. D. Vorst,et al.  Iterative solution methods for certain sparse linear systems with a non-symmetric matrix arising from PDE-problems☆ , 1981 .

[4]  R. Beauwens.  Modified incomplete factorization strategies , 1991 .

[5]  Richard Barrett,et al.  Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods , 1994, Other Titles in Applied Mathematics.

[6]  Dianne P. O'Leary,et al.  Efficient Iterative Solution of the Three-Dimensional Helmholtz Equation , 1998 .

[7]  Charles R. Johnson,et al.  Matrix analysis , 1985, Statistical Inference for Engineers and Data Scientists.

[8]  Charbel Farhat,et al.  Residual-Free Bubbles for the Helmholtz Equation , 1996 .

[9]  T. Manteuffel.  An incomplete factorization technique for positive definite linear systems , 1980 .

[10]  D. R. Fokkema,et al.  BICGSTAB( L ) FOR LINEAR EQUATIONS INVOLVING UNSYMMETRIC MATRICES WITH COMPLEX , 1993 .

[11]  Antonini Macedo,et al.  Two Level Finite Element Method and its Application to the Helmholtz Equation , 1997 .

[12]  Magolu Monga-Made.  Taking Advantage of the Potentialities of Dynamically Modified Block Incomplete Factorizations , 1998 .

[13]  Jack Dongarra,et al.  Numerical Linear Algebra for High-Performance Computers , 1998 .

[14]  Robert J. Plemmons,et al.  Nonnegative Matrices in the Mathematical Sciences , 1979, Classics in Applied Mathematics.

[15]  Robert Beauwens,et al.  HIGH-PERFORMANCE PCG SOLVERS FOR FEM STRUCTURAL ANALYSIS , 1996 .

[16]  J. H. Wilkinson.  The algebraic eigenvalue problem , 1966 .

[17]  D. Kershaw.  The incomplete Cholesky—conjugate gradient method for the iterative solution of systems of linear equations , 1978 .

[18]  Ivo Babuška,et al.  A Generalized Finite Element Method for solving the Helmholtz equation in two dimensions with minimal pollution , 1995 .

[19]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[20]  Chun-Hua Guo,et al.  Incomplete block factorization preconditioning for linear systems arising in the numerical solution of the Helmholtz equation , 1996 .

[21]  Elisabeth Larsson,et al.  Iterative Solution of the Helmholtz Equation by a Second-Order Method , 1999, SIAM J. Matrix Anal. Appl..

[22]  Thomas J. R. Hughes,et al.  Recent developments in finite element methods for structural acoustics , 1996 .

[23]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[24]  Magolu Monga-Made.  Ordering strategies for modified block incomplete factorizations , 1995 .

[25]  Oliver G. Ernst,et al.  A finite-element capacitance matrix method for exterior Helmholtz problems , 1996 .

[26]  Roland W. Freund,et al.  Conjugate Gradient-Type Methods for Linear Systems with Complex Symmetric Coefficient Matrices , 1992, SIAM J. Sci. Comput..

[27]  Jinchao Xu,et al.  The analysis of multigrid algorithms for nonsymmetric and indefinite elliptic problems , 1988 .

[28]  Jörg Liesen,et al.  Computable Convergence Bounds for GMRES , 2000, SIAM J. Matrix Anal. Appl..

[29]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[30]  O. Axelsson.  Iterative solution methods , 1995 .

[31]  Yvan Notay,et al.  DRIC: A dynamic version of the RIC method , 1994, Numer. Linear Algebra Appl..

[32]  J. Meijerink,et al.  An iterative solution method for linear systems of which the coefficient matrix is a symmetric -matrix , 1977 .

[33]  O. Axelsson,et al.  On the eigenvalue distribution of a class of preconditioning methods , 1986 .