TY - JOUR
T1 - Parallel hierarchical solvers and preconditioners for boundary element methods
AU - Grama, Ananth
AU - Kumar, Vipin
AU - Sameh, Ahmed
PY - 1998
Y1 - 1998
N2 - 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.
AB - 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.
KW - Barnes-Hut method
KW - Boundary element method
KW - Dense iterative solver
KW - Fast multipole method
KW - Hierarchical dense matrix-vector product
KW - Parallel treecode
KW - Preconditioning boundary element methods
UR - http://www.scopus.com/inward/record.url?scp=0032131422&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=0032131422&partnerID=8YFLogxK
U2 - 10.1137/S1064827596313322
DO - 10.1137/S1064827596313322
M3 - Article
AN - SCOPUS:0032131422
SN - 1064-8275
VL - 20
SP - 337
EP - 358
JO - SIAM Journal on Scientific Computing
JF - SIAM Journal on Scientific Computing
IS - 1
ER -