A fast poisson solver of arbitrary order accuracy in rectangular regions

A. Averbuch, M. Israeli, L. Vozovoi

Research output: Contribution to journalArticlepeer-review

45 Scopus citations


In this paper we propose a direct method for the solution of the Poisson equation in rectangular regions. It has an arbitrary order accuracy and low CPU requirements which makes it practical for treating large-scale problems. The method is based on a pseudospectral Fourier approximation and a polynomial subtraction technique. Fast convergence of the Fourier series is achieved by removing the discontinuities at the corner points using polynomial subtraction functions. These functions have the same discontinuities at the corner points as the sought solution. In addition to this, they satisfy the Laplace equation so that the subtraction procedure does not generate nonperiodic, nonhomogeneous terms. The solution of a boundary value problem is obtained in a series form in O(N log N) floating point operations, where N2 is the number of grid nodes. Evaluating the solution at all N2 interior points requires O(N2 log N) operations.

Original languageEnglish
Pages (from-to)933-952
Number of pages20
JournalSIAM Journal of Scientific Computing
Issue number3
StatePublished - May 1998
Externally publishedYes


  • Boundary value problem
  • Corner discontinuities
  • Poisson equation
  • Polynomial subtraction
  • Rectangular region
  • Spectral method


Dive into the research topics of 'A fast poisson solver of arbitrary order accuracy in rectangular regions'. Together they form a unique fingerprint.

Cite this