An efficient high-order compact finite difference method for the Helmholtz equation

Abstract

This paper is devoted to applying the sixth-order compact finite difference approach to the Helmholtz equation. Instead of using matrix inversion, a discrete sinusoidal transform is used as a quick solver to solve the discretized system resulted from the compact finite difference method. Through this way, the computational costs of the method with large numbers of nodes are greatly reduced. The efficiency and accuracy of the scheme are investigated by solving some illustrative examples, having the exact solutions.