Fourier analysis of the CGMN method for solving the Helmholtz equation

arXiv1210.2644MaRDI QIDQ6236243FDOQ6236243

Publication date: 9 October 2012

Abstract: The Helmholtz equation arises in many applications, such as seismic and medical imaging. These application are characterized by the need to propagate many wavelengths through an inhomogeneous medium. The typical size of the problems in 3D applications precludes the use of direct factorization to solve the equation and hence iterative methods are used in practice. For higher wavenumbers, the system becomes increasingly indefinite and thus good preconditioners need to be constructed. In this note we consider an accelerated Kazcmarz method (CGMN) and present an expression for the resulting iteration matrix. This iteration matrix can be used to analyze the convergence of the CGMN method. In particular, we present a Fourier analysis for the method applied to the 1D Helmholtz equation. This analysis suggests an optimal choice of the relaxation parameter. Finally, we present some numerical experiments.

Mathematics Subject Classification ID

Preconditioners for iterative methods (65F08) Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation (35J05)

This page was built for publication: Fourier analysis of the CGMN method for solving the Helmholtz equation

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