Publication
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)
UAntwerpen
Faculty/Department
Research group
Project info
Simulation of image formation in X-ray phase contrast tomography
Publication type
Subject
Affiliation
Publications with a UAntwerp address
External links
Web of Science
Record
Identifier
Creation 05.09.2013
Last edited 05.01.2025
To cite this reference