Local Fourier Analysis of the Complex Shifted Laplacian preconditioner for Helmholtz problems

Abstract

In this paper we solve the Helmholtz equation with multigrid preconditioned Krylov subspace methods. The class of Shifted Laplacian preconditioners are known to significantly speed-up Krylov convergence. However, these preconditioners have a parameter beta, a measure of the complex shift. Due to 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, as well as being near-optimal in terms of Krylov iteration count.

arXiv comment20 pages

References in corpus (2)

• On the indefinite Helmholtz equation: complex stretched absorbing boundary layers, iterative analysis, and preconditioning
• Analyzing the wave number dependency of the convergence rate of a multigrid preconditioned Krylov method for the Helmholtz equation with an absorbing layer

Cited by in corpus (8)

• A dispersion minimizing scheme for the 3-D Helmholtz equation based on ray theory
• A new level-dependent coarsegrid correction scheme for indefinite Helmholtz problems
• Complex additive geometric multilevel solvers for Helmholtz equations on spacetrees
• An efficient multigrid calculation of the far field map for Helmholtz and Schrödinger equations
• The Helmholtz equation with uncertainties in the wavenumber
• A hybrid shifted Laplacian multigrid and domain decomposition preconditioner for the elastic Helmholtz equations
• Evolving Generalizable Multigrid-Based Helmholtz Preconditioners with Grammar-Guided Genetic Programming
• Regula falsi based automatic regularization method for PDE constrained optimization