Parallel hierarchical solvers and preconditioners for boundary element methods

Ananth Grama, Vipin Kumar, Ahmed Sameh

Research output: Contribution to journalArticlepeer-review

34 Scopus citations

Abstract

The method of moments is an important tool for solving boundary integral equations arising in a variety of applications. It transforms the physical problem into a dense linear system. Due to the large number of variables and the associated computational requirements, these systems are solved iteratively using methods such as GMRES and CG and its variants. The core operation of these iterative solvers is the application of the system matrix to a vector. This requires 0(n2) operations and memory using accurate dense methods. The computational complexity can be reduced to O(n log n) and the memory requirement to ⊖(n) using hierarchical approximation techniques. The algorithmic speedup from approximation can be combined with parallelism to yield very fast dense solvers. In this paper, we present efficient parallel formulations of dense iterative solvers based on hierarchical approximations for solving potential integral equations of the first kind. We study the impact of various parameters on the accuracy and performance of the parallel solver. We demonstrate that our parallel formulation incurs minimal parallel processing overhead and scales up to a large number of processors. We present two preconditioning techniques for accelerating the convergence of the iterative solver. These techniques are based on an inner-outer scheme and a block-diagonal scheme based on a truncated Green's function. We present detailed experimental results on up to 256 processors of a Cray T3D. Our code achieves raw computational speeds of over 5 GFLOPS. When compared to the accurate solver, this corresponds to a speed of approximately 776 GFLOPS.

Original languageEnglish (US)
Pages (from-to)337-358
Number of pages22
JournalSIAM Journal of Scientific Computing
Volume20
Issue number1
DOIs
StatePublished - Jan 1 1998

Keywords

  • Barnes-Hut method
  • Boundary element method
  • Dense iterative solver
  • Fast multipole method
  • Hierarchical dense matrix-vector product
  • Parallel treecode
  • Preconditioning boundary element methods

Fingerprint Dive into the research topics of 'Parallel hierarchical solvers and preconditioners for boundary element methods'. Together they form a unique fingerprint.

Cite this