Biomolecular electrostatics using a fast multipole BEM on up to 512 GPUs and a billion unknowns
arXiv:1007.4591 · doi:10.1016/j.cpc.2011.02.013
Abstract
We present teraflop-scale calculations of biomolecular electrostatics enabled by the combination of algorithmic and hardware acceleration. The algorithmic acceleration is achieved with the fast multipole method (FMM) in conjunction with a boundary element method (BEM) formulation of the continuum electrostatic model, as well as the BIBEE approximation to BEM. The hardware acceleration is achieved through graphics processors, GPUs. We demonstrate the power of our algorithms and software for the calculation of the electrostatic interactions between biological molecules in solution. The applications demonstrated include the electrostatics of protein--drug binding and several multi-million atom systems consisting of hundreds to thousands of copies of lysozyme molecules. The parallel scalability of the software was studied in a cluster at the Nagasaki Advanced Computing Center, using 128 nodes, each with 4 GPUs. Delicate tuning has resulted in strong scaling with parallel efficiency of 0.8 for 256 and 0.5 for 512 GPUs. The largest application run, with over 20 million atoms and one billion unknowns, required only one minute on 512 GPUs. We are currently adapting our BEM software to solve the linearized Poisson-Boltzmann equation for dilute ionic solutions, and it is also designed to be flexible enough to be extended for a variety of integral equation problems, ranging from Poisson problems to Helmholtz problems in electromagnetics and acoustics to high Reynolds number flow.
References in corpus (1)
Cited by in corpus (19)
- Petascale turbulence simulation using a highly parallel fast multipole method on GPUs
- A Tuned and Scalable Fast Multipole Method as a Preeminent Algorithm for Exascale Systems
- A biomolecular electrostatics solver using Python, GPUs and boundary elements that can handle solvent-filled cavities and Stern layers
- Hierarchical N-body simulations with auto-tuning for heterogeneous systems
- Long range correction for multi-site Lennard-Jones models and planar interfaces
- A GPU-accelerated Direct-sum Boundary Integral Poisson-Boltzmann Solver
- A GPU-Accelerated Fast Summation Method Based on Barycentric Lagrange Interpolation and Dual Tree Traversal
- Pipelined Iterative Solvers with Kernel Fusion for Graphics Processing Units
- An O(N) and parallel approach to integral problems by a kernel-independent fast multipole method: Application to polarization and magnetization of interacting particles
- Parallel Higher-order Boundary Integral Electrostatics Computation on Molecular Surfaces with Curved Triangulation
- Parallelization of the inverse fast multipole method with an application to boundary element method
- Fast Multipole Method as a Matrix-Free Hierarchical Low-Rank Approximation
- A pilgrimage to gravity on GPUs
- Generalizing The Mean Spherical Approximation as a Multiscale, Nonlinear Boundary Condition at the Solute--Solvent Interface
- An Enhanced Spectral Boundary Integral Method for Modeling Highly Nonlinear Water Waves in Variable Depth
- Communication Reducing Algorithms for Distributed Hierarchical N-Body Problems with Boundary Distributions
- Inexact Krylov iterations and relaxation strategies with fast-multipole boundary element method
- Simple and efficient GPU parallelization of existing H-Matrix accelerated BEM code
- A GPU-Parallelized Interpolation-Based Fast Multipole Method for the Relativistic Space-Charge Field Calculation