A scalable Helmholtz solver combining the shifted Laplace preconditioner with multigrid deflation
A.H. Sheikh
D. Lahaye
C. Vuik
More Info
expand_more
Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.
Abstract
A Helmholtz solver whose convergence is parameter independent can be obtained by combining the shifted Laplace preconditioner with multigrid deflation. To proof this claim, we develop a Fourier analysis of a two-level variant of the algorithm proposed in [1]. In this algorithm those eigenvalues that prevent the shifted Laplace preconditioner from being scalable are removed by deflation using multigrid vectors. Our analysis shows that the spectrum of the two-grid operator consists of a cluster surrounded by a few outliers, yielding a number of outer Krylov subspace iterations that remains constant as the wave number increases. Our analysis furthermore shows that the imaginary part of the shift in the two-grid operator can be made arbitrarily large without affecting the convergence. This opens promising perspectives on obtaining a very good preconditioner at very low cost. Numerical tests for problems with constant and non-constant wave number illustrate our convergence theory.