A Fast Spectral Subtractional Solver for Elliptic Equations

  • Authors:
  • Elena Braverman;Boris Epstein;Moshe Israeli;Amir Averbuch

  • Affiliations:
  • Department of Mathematics and Statistics, University of Calgary, 2500 University Drive N.W., Calgary, Alberta T2N 1N4, Canada. maelena@math.ucalgary.ca;The Academic College of Tel Aviv-Yaffo, Antokolsky Str. 4, Tel Aviv 64044, Israel;Technion-Israel Institute of Technology, Computer Science Department, Haifa 32000, Israel;School of Computer Sciences, Tel Aviv University, Tel Aviv 69978, Israel

  • Venue:
  • Journal of Scientific Computing
  • Year:
  • 2004

Quantified Score

Hi-index 0.01

Visualization

Abstract

The paper presents a fast subtractional spectral algorithm for the solution of the Poisson equation and the Helmholtz equation which does not require an extension of the original domain. It takes O(N2 log N) operations, where N is the number of collocation points in each direction. The method is based on the eigenfunction expansion of the right hand side with integration and the successive solution of the corresponding homogeneous equation using Modified Fourier Method. Both the right hand side and the boundary conditions are not assumed to have any periodicity properties. This algorithm is used as a preconditioner for the iterative solution of elliptic equations with non-constant coefficients. The procedure enjoys the following properties: fast convergence and high accuracy even when the computation employs a small number of collocation points. We also apply the basic solver to the solution of the Poisson equation in complex geometries.