A fast Fourier transform based direct solver for the Helmholtz problem

DOI10.1002/NLA.2283MaRDI QIDQ6306552arXivWikidataFDO

Publication date 11 September 2018

Direct numerical methods for linear systems and matrix inversion (65F05) Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation (35J05) Boundary value problems for second-order elliptic equations (35J25) Fourier and Fourier-Stieltjes transforms and other transforms of Fourier type (42A38) Numerical solution of discretized equations for boundary value problems involving PDEs (65N22) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30) Theoretical approximation in context of PDEs (35A35) Numerical methods for discrete and fast Fourier transforms (65T50)

Abstract: This paper is devoted to the efficient numerical solution of the Helmholtz equation in a two- or three-dimensional rectangular domain with an absorbing boundary condition (ABC). The Helmholtz problem is discretized by standard bilinear and trilinear finite elements on an orthogonal mesh yielding a separable system of linear equations. The main key to high performance is to employ the Fast Fourier transform (FFT) within a fast direct solver to solve the large separable systems. The computational complexity of the proposed FFT based direct solver is O(N log N) operations. Numerical results for both two- and three-dimensional problems are presented confirming the efficiency of the method discussed.

This page was built for publication: A fast Fourier transform based direct solver for the Helmholtz problem

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