Abstract:Phylogenies are central to many research areas in biology and commonly estimated using likelihood based methods. Unfortunately, any likelihood based method, including Bayesian inference, can be restrictively slow for large datasets —with many taxa and/or many sites in the sequence alignment— or complex substitutions models. The primary limiting factor when using large datasets and/or complex models in probabilistic phylogenetic analyses is the likelihood calculation, which dominates the total computation time. To address this bottleneck, we incorporated the high-performance phylogenetic library BEAGLE into RevBayes, which enables multi-threading on multi-core CPUs and GPUs, as well as hardware specific vectorized instructions for faster likelihood calculations. Our new implementation of RevBayes+BEAGLE retains the flexibility and dynamic nature that users expect from vanilla RevBayes. Additionally, we implemented native parallelization within RevBayes without an external library using the message passing interface (MPI); RevBayes+MPI. We evaluated our new implementation of RevBayes+BEAGLE using multi-threading on CPUs and two different powerful GPUs (NVidia Titan V and NVIDIA A100) against our native implementation of RevBayes+MPI. We found good improvements in speedup when multiple cores were used, with up to 20-fold speedup when using multiple CPU cores and over 90-fold speedup when using multiple GPU cores. The improvement depended on the data type used, DNA or amino acids, and the size of the alignment, but less on the size of the tree. We additionally investigated the cost of rescaling partial likelihoods to avoid numerical underflow and showed that unnecessarily frequent and inefficient rescaling can increase runtimes up to 4-fold. Finally, we presented and compared a new approach to store partial likelihoods on branches instead of nodes which can speed up computations up to 1.7-times but comes at twice the memory requirements.

Parallelizing Computation of Expected Values in Recombinant Binomial Trees

Single MCMC Chain Parallelisation on Decision Trees

GPGPUs in computational finance: massive parallel computing for American style options

Gaussian Recombining Split Tree

A new parallel solver suited for arbitrary semilinear parabolic partial differential equations based on generalized random trees

Parallelizing MCMC with Random Partition Trees

Recombined Multinomial Tree Based on Saddle-Point Approximation and Its Application to Levy Models Options Pricing

Parallel Computing Method of Valuing for Multi-asset European Option.

Parallel American Monte Carlo

The Boyle–Romberg Trinomial Tree, a Highly Efficient Method for Double Barrier Option Pricing

Bayesian Phylogenetic Analysis on multi-core Compute Architectures: Implementation and evaluation of BEAGLE in RevBayes with MPI

A Communication-Efficient Parallel Algorithm for Decision Tree

Parallel Accelerated Custom Correlation Coefficient Calculations for Genomics Applications

Belief Propagation by Message Passing in Junction Trees: Computing Each Message Faster Using GPU Parallelization

Randomized Binomial Tree and Pricing of American-Style Options

Fast generation of implied volatility surface: Optimize the traditional numerical analysis and machine learning

Improving Tree Probability Estimation with Stochastic Optimization and Variance Reduction

Towards Automatic Parallelization of Tree Reductions in Dynamic Programming

Parallel Sampling of Decomposable Graphs Using Markov Chains on Junction Trees

RJHMC-Tree for Exploration of the Bayesian Decision Tree Posterior