期刊
APPLIED NUMERICAL MATHEMATICS
卷 50, 期 3-4, 页码 409-425出版社
ELSEVIER SCIENCE BV
DOI: 10.1016/j.apnum.2004.01.009
关键词
Helmholtz equation; Krylov subspace; preconditioner
In 1983, a preconditioner was proposed [J. Comput. Phys. 49 (1983) 443] based on the Laplace operator for solving the discrete Helmholtz equation efficiently with CGNR. The preconditioner is especially effective for low wavenumber cases where the linear system is slightly indefinite. Laird [Preconditioned iterative solution of the 2D Helmholtz equation, First Year's Report, St. Hugh's College, Oxford, 2001] proposed a preconditioner where an extra term is added to the Laplace operator. This term is similar to the zeroth order term in the Helmholtz equation but with reversed sign. In this paper, both approaches are further generalized to a new class of preconditioners, the so-called shifted Laplace preconditioners of the form Deltaphi - ak(2)phi with alpha is an element of C. Numerical experiments for various wavenumbers indicate the effectiveness of the preconditioner. The preconditioner is evaluated in combination with GMRES, Bi-CGSTAB, and CGNR. (C) 2004 IMACS. Published by Elsevier B.V. All rights reserved.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据