Log In Sign Up

A two-level shifted Laplace Preconditioner for Helmholtz Problems: Field-of-values analysis and wavenumber-independent convergence

by   Luis García Ramos, et al.

One of the main tools for solving linear systems arising from the discretization of the Helmholtz equation is the shifted Laplace preconditioner, which results from the discretization of a perturbed Helmholtz problem -Δ u - (k^2 + i ε )u = f where 0 ≠ε∈R is an absorption parameter. In this work we revisit the idea of combining the shifted Laplace preconditioner with two-level deflation and apply it to Helmholtz problems discretized with linear finite elements. We use the convergence theory of GMRES based on the field of values to prove that GMRES applied to the two-level preconditioned system with a shift parameter ε∼ k^2 converges in a number of iterations independent of the wavenumber k,provided that the coarse mesh size H satisfies a condition of the form Hk^2≤ C for some constant C depending on the domain but independent of the wavenumber k. This behaviour is sharply different to the standalone shifted Laplacian, for which wavenumber-independent GMRES convergence has been established only under the condition that ε∼ k by [M.J. Gander, I.G. Graham and E.A. Spence, Numer. Math., 131 (2015), 567-614]. Finally, we present numerical evidence that wavenumber-independent convergence of GMRES also holds for pollution-free meshes, where the coarse mesh size satisfies Hk^3/2≤ C, and inexact coarse grid solves.


page 1

page 2

page 3

page 4


An averaged space-time discretization of the stochastic p-Laplace system

We study the stochastic p-Laplace system in a bounded domain. We propose...

Relaxed Kacanov scheme for the p-Laplacian with large p

We introduce a globally convergent relaxed Kacanov scheme for the comput...

Semi matrix-free twogrid shifted Laplacian preconditioner for the Helmholtz equation with near optimal shifts

Due to its significance in terms of wave phenomena a considerable effort...

The Hierarchical Subspace Iteration Method for Laplace–Beltrami Eigenproblems

Sparse eigenproblems are important for various applications in computer ...

Estimation and numerical validation of inf-sup constant for bilinear form (p, div u)

We give a derivation for the value of inf-sup constant for the bilinear ...

Another source of mesh dependence in topology optimization

The topology optimization community has regularly employed nonlinear pro...