|
Title
|
|
|
|
Local Fourier analysis of the complex shifted Laplacian preconditioner for Helmholtz problems
|
|
|
Author
|
|
|
|
|
|
|
Abstract
|
|
|
|
| In this paper, we solve the Helmholtz equation with multigrid preconditioned Krylov subspace methods. The class of shifted Laplacian preconditioners is known to significantly speed up Krylov convergence. However, these preconditioners have a parameter inline image, a measure of the complex shift. Because of contradictory requirements for the multigrid and Krylov convergence, the choice of this shift parameter can be a bottleneck in applying the method. In this paper, we propose a wavenumber-dependent minimal complex shift parameter, which is predicted by a rigorous k-grid local Fourier analysis (LFA) of the multigrid scheme. We claim that, given any (regionally constant) wavenumber, this minimal complex shift parameter provides the reader with a parameter choice that leads to efficient Krylov convergence. Numerical experiments in one and two spatial dimensions validate the theoretical results. It appears that the proposed complex shift is both the minimal requirement for a multigrid V-cycle to converge and being near optimal in terms of Krylov iteration count. |
|
|
|
Language
|
|
|
|
English
|
|
|
Source (journal)
|
|
|
|
Numerical linear algebra with applications. - Chichester
|
|
|
Publication
|
|
|
|
Chichester
:
Wiley
,
2013
|
|
|
ISSN
|
|
|
|
1070-5325
|
|
|
DOI
|
|
|
|
10.1002/NLA.1881
|
|
|
Volume/pages
|
|
|
|
20
:4
(2013)
, p. 575-597
|
|
|
ISI
|
|
|
|
000322031700004
|
|
|
Full text (Publisher's DOI)
|
|
|
|
|
|
|
Full text (open access)
|
|
|
|
|
|