Abstract:SUMMARY Performance improvements in biomolecular simulations based on molecular dynamics (MD) codes are widely desired. Unfortunately, the factors, which allowed past performance improvements, particularly the microprocessor clock frequencies, are no longer increasing. Hence, novel software and hardware solutions are being explored for accelerating performance of widely used MD codes. In this paper, we describe our efforts on porting, optimizing and tuning of Large‐scale Atomic/Molecular Massively Parallel Simulator, a popular MD framework, on heterogeneous architectures: multi‐core processors with graphical processing unit (GPU) accelerators. Our implementation is based on accelerating the most computationally expensive non‐bonded interaction terms on the GPUs and overlapping the computation on the CPU and GPUs. This functionality is built on top of message passing interface that allows multi‐level parallelism to be extracted even at the workstation level with the multi‐core CPUs and allows extension of the implementation on GPU‐enabled clusters. We hypothesize that the optimal benefit of heterogeneous architectures for applications will come by utilizing all possible resources (for example, CPU‐cores and GPU devices on GPU‐enabled clusters). Benchmarks for a range of biomolecular system sizes are provided, and an analysis is performed on four generations of NVIDIA's GPU devices. On GPU‐enabled Linux clusters, by overlapping and pipelining computation and communication, we observe up to 10‐folds application acceleration in multi‐core and multi‐GPU environments illustrating significant performance improvements. Detailed analysis of the implementation is presented that allows identification of bottlenecks in algorithm, indicating that code optimization and improvements on GPUs could allow microsecond scale simulation throughput on workstations and inexpensive GPU clusters, putting widely desired biologically relevant simulation time‐scales within reach of a large user community. In order to systematically optimize simulation throughput and to enable performance prediction, we have developed a parameterized performance model that will allow developers and users to explore the performance potential of future heterogeneous systems for biological simulations. Copyright © 2012 John Wiley & Sons, Ltd.

Population genetic simulation: Benchmarking frameworks for non‐standard models of natural selection

Towards Simulation Optimization: An Examination of the Impact of Scaling on Coalescent and Forward Simulations

SimBit: A high performance, flexible and easy-to-use population genetic simulator

SLiM 2: Flexible, Interactive Forward Genetic Simulations

SLiM 3: Forward Genetic Simulations Beyond the Wright–Fisher Model

Clotho: addressing the scalability of forward time population genetic simulation

Exact and efficient phylodynamic simulation from arbitrarily large populations

Population size rescaling significantly biases outcomes of forward-in-time population genetic simulations

SLiM: Simulating Evolution with Selection and Linkage

Population genetics meets ecology: a guide to individual-based simulations in continuous landscapes

Simulation-Based Evaluation of Methods, Data Types, and Temporal Sampling Schemes for Detecting Recent Population Declines

Accurate inference of population history in the presence of background selection

Robust Forward Simulations of Recurrent Hitchhiking

Recent advances in stochastic simulation algorithms create opportunities to study new ecological and eco‐evolutionary problems

Genealogical processes of non-neutral population models under rapid mutation

Simultaneous Inference of Past Demography and Selection from the Ancestral Recombination Graph under the Beta Coalescent

AdmixSim 2: a forward-time simulator for modeling complex population admixture

Sampling through time and phylodynamic inference with coalescent and birth-death models

Understanding Past Population Dynamics: Bayesian Coalescent-Based Modeling with Covariates

Exact Simulation of Conditioned Wright-Fisher Models.

Performance modeling of microsecond scale biological molecular dynamics simulations on heterogeneous architectures