A fast parallel algorithm for the Poisson equation on a disk Academic Article uri icon

abstract

  • A parallel algorithm for solving the Poisson equation with either Dirichlet or Neumann conditions is presented. The solver follows some of the principles introduced in a previous fast algorithm for evaluating singular integral transforms by Daripa et al. Here we present recursive relations in Fourier space together with fast Fourier transforms which lead to a fast and accurate algorithm for solving Poisson problems within a unit disk. The algorithm is highly parallelizable and our implementation is virtually architecture-independent. Theoretical estimates show good parallel scalability of the algorithm, and numerical results show the accuracy of the method for problems with sharp variations on inhomogeneous term. Finally, performance results for sequential and parallel implementations are presented. 2001 Academic Press.

published proceedings

  • JOURNAL OF COMPUTATIONAL PHYSICS

author list (cited authors)

  • Borges, L., & Daripa, P.

citation count

  • 24

complete list of authors

  • Borges, L||Daripa, P

publication date

  • May 2001