A multigrid-based shifted Laplacian preconditioner for a fourth-order Helmholtz discretization

N. Umetani, S. P. MacLachlan, C. W. Oosterlee*

*Corresponding author for this work

Research output: Contribution to journalArticleAcademicpeer-review

Abstract

In this paper, an iterative solution method for a fourth-order accurate discretization of the Helmholtz equation is presented. The method is a generalization of that presented in (SIAM J. Sci. Comput. 2006; 27:1471-1492), where multigrid was employed as a preconditioner for a Krylov subspace iterative method. The multigrid preconditioner is based on the solution of a second Helmholtz operator with a complex-valued shift. In particular, we compare preconditioners based on a point-wise Jacobi smoother with those using an ILU(0) smoother, we compare using the prolongation operator developed by de Zeeuw in (J. Comput. Appl. Math. 1990; 33:1-27) with interpolation operators based on algebraic multigrid principles, and we compare the performance of the Krylov subspace method Bi-conjugate gradient stabilized with the recently introduced induced dimension reduction method, IDR(s). These three improvements are combined to yield an efficient solver for heterogeneous problems.

Original languageEnglish
Pages (from-to)603-626
Number of pages24
JournalNumerical Linear Algebra with Applications
Volume16
Issue number8
DOIs
Publication statusPublished - Aug 2009
Externally publishedYes

Keywords

  • Algebraic multigrid
  • Complex-valued multigrid preconditioner
  • Helmholtz equation
  • ILU smoother
  • Non-constant wavenumber

Fingerprint

Dive into the research topics of 'A multigrid-based shifted Laplacian preconditioner for a fourth-order Helmholtz discretization'. Together they form a unique fingerprint.

Cite this