On the Optimality of Shifted Laplacian in a Class of Polynomial Preconditioners for the Helmholtz Equation

DOI10.1007/978-3-319-28832-1_3zbMATH Open1366.65093arXiv1501.04445OpenAlexW2594348025MaRDI QIDQ5266543FDOQ5266543

Siegfried Cools, Wim Vanroose

Publication date: 16 June 2017

Published in: Modern Solvers for Helmholtz Problems (Search for Journal in Brave)

Abstract: This paper introduces and explores the class of expansion preconditioners EX(m) that forms a direct generalization to the classic complex shifted Laplace (CSL) preconditioner for Helmholtz problems. The construction of the EX(m) preconditioner is based upon a truncated Taylor series expansion of the original Helmholtz operator inverse. The expansion preconditioner is shown to significantly improve Krylov solver convergence rates for the Helmholtz problem for growing values of the number of series terms m. However, the addition of multiple terms in the expansion also increases the computational cost of applying the preconditioner. A thorough cost-benefit analysis of the addition of extra terms in the EX(m) preconditioner proves that the CSL or EX(1) preconditioner is the practically most efficient member of the expansion preconditioner class. Additionally, possible extensions to the expansion preconditioner class that further increase preconditioner efficiency are suggested.

Full work available at URL: https://arxiv.org/abs/1501.04445

zbMATH Keywords

convergenceKrylov subspace methodspreconditioningfinite difference methodHelmholtz equationnumerical experimentmultigrid methodsshifted Laplacian

Mathematics Subject Classification ID

Multigrid methods; domain decomposition for boundary value problems involving PDEs (65N55) Preconditioners for iterative methods (65F08) Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation (35J05) Finite difference methods for boundary value problems involving PDEs (65N06) Stability and convergence of numerical methods for boundary value problems involving PDEs (65N12)

Cites Work

• Title not available (Why is that?)
• GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems
• A perfectly matched layer for the absorption of electromagnetic waves
• A class of analytic perturbations for one-body Schrödinger Hamiltonians
• Title not available (Why is that?)
• Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems
• Multi-Level Adaptive Solutions to Boundary-Value Problems
• Title not available (Why is that?)
• Title not available (Why is that?)
• Why it is Difficult to Solve Helmholtz Problems with Classical Iterative Methods
• Krylov Subspace Methods
• Recent computational developments in Krylov subspace methods for linear systems
• A Multigrid Tutorial, Second Edition
• Spectral Analysis of the Discrete Helmholtz Operator Preconditioned with a Shifted Laplacian
• A Flexible Inner-Outer Preconditioned GMRES Algorithm
• The Numerical Solution of Laplace's Equation
• An iterative method for the Helmholtz equation
• Applying GMRES to the Helmholtz equation with shifted Laplacian preconditioning: What is the largest shift for which wavenumber-independent convergence is guaranteed?
• On accuracy conditions for the numerical computation of waves
• Comparison of multigrid and incomplete LU shifted-Laplace preconditioners for the inhomogeneous Helmholtz equation
• On a class of preconditioners for solving the Helmholtz equation
• A multigrid method enhanced by Krylov subspace iteration for discrete Helmholtz equations
• A Novel Multigrid Based Preconditioner For Heterogeneous Helmholtz Problems
• Preconditioning Helmholtz linear systems
• Incomplete factorization-based preconditionings for solving the Helmholtz equation
• Local Fourier analysis of the complex shifted Laplacian preconditioner for Helmholtz problems
• On the convergence of shifted Laplace preconditioner combined with multilevel deflation
• Polynomial Preconditioned GMRES and GMRES-DR
• Flexible Inner-Outer Krylov Subspace Methods
• On the indefinite Helmholtz equation: Complex stretched absorbing boundary layers, iterative analysis, and preconditioning
• On a multilevel Krylov method for the Helmholtz equation preconditioned by shifted Laplacian
• GMRES with multiple preconditioners
• Analyzing the wave number dependency of the convergence rate of a multigrid preconditioned Krylov method for the Helmholtz equation with an absorbing layer
• Two-Level preconditioned Krylov subspace methods for the solution of three-dimensional heterogeneous Helmholtz problems in seismics

Cited In (6)

• Contraction operator transformation for the complex heterogeneous Helmholtz equation
• Preconditioning the Helmholtz equation with the shifted Laplacian and Faber polynomials
• Symmetry and invariant bases in finite element exterior calculus
• Robust Superlinear Krylov Convergence for Complex Noncoercive Compact-Equivalent Operator Preconditioners
• Analysis of the Shifted Helmholtz Expansion Preconditioner for the Helmholtz Equation
• How Large a Shift is Needed in the Shifted Helmholtz Preconditioner for its Effective Inversion by Multigrid?

Uses Software

• Eigtool

This page was built for publication: On the Optimality of Shifted Laplacian in a Class of Polynomial Preconditioners for the Helmholtz Equation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5266543)