Vector operations for accelerating expensive Bayesian computations - a tutorial guide

DOI10.1214/21-BA1265MaRDI QIDQ6202922zbMATH OpenOpenAlexFDO

Authors David J. Warne, S. A. Sisson, C. C. Drovandi

Publication date 27 February 2024

Published in Bayesian Analysis (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1902.09046, https://projecteuclid.org/journals/bayesian-analysis/volume-17/issue-2/Vector-Operations-for-Accelerating-Expensive-Bayesian-Computations-A-Tutorial/10.1214/21-BA1265.full

zbMATH Keywords

sequential Monte Carloapproximate Bayesian computationsingle instruction multiple datavectorisationadvanced vector extensionsweakly informative priors

Mathematics Subject Classification ID

Statistics (62-XX) Applied statistics (educational aspects) (97K80)

Abstract: Many applications in Bayesian statistics are extremely computationally intensive. However, they are often inherently parallel, making them prime targets for modern massively parallel processors. Multi-core and distributed computing is widely applied in the Bayesian community, however, very little attention has been given to fine-grain parallelisation using single instruction multiple data (SIMD) operations that are available on most modern commodity CPUs and is the basis of GPGPU computing. In this work, we practically demonstrate, using standard programming libraries, the utility of the SIMD approach for several topical Bayesian applications. We show that SIMD can improve the floating point arithmetic performance resulting in up to $6imes$ improvement in serial algorithm performance. Importantly, these improvements are multiplicative to any gains achieved through multi-core processing. We illustrate the potential of SIMD for accelerating Bayesian computations and provide the reader with techniques for exploiting modern massively parallel processing environments using standard tools.

Cites work

• scientific article; zbMATH DE number 236854 (Why is no real title available?)
• A comparative review of dimension reduction methods in approximate Bayesian computation
• A sequential particle filter method for static models
• Approximation of Bayesian predictive \(p\)-values with regression ABC
• Bad environments, good environments: a non-Gaussian asymmetric volatility model
• Bayesian computation: a summary of the current state, and samples backwards and forwards
• Bayesian parametric bootstrap for models with intractable likelihoods
• Continuous Markov processes and stochastic equations
• GPU-accelerated Gibbs sampling: a case study of the horseshoe probit model
• Generalized autoregressive conditional heteroscedasticity
• Julia: a fresh approach to numerical computing
• Massive Parallelization Boosts Big Bayesian Multidimensional Scaling
• Multilevel Monte Carlo in approximate Bayesian computation
• Multilevel rejection sampling for approximate Bayesian computation
• On the convergence of adaptive sequential Monte Carlo methods
• Particle Markov Chain Monte Carlo Methods
• Prior distributions for variance parameters in hierarchical models (Comment on article by Browne and Draper)
• Recalibration: a post-processing method for approximate Bayesian computation
• SIMD parallel MCMC sampling with applications for big-data Bayesian analytics
• Sequential Monte Carlo Samplers
• Sequential Monte Carlo samplers with independent Markov chain Monte Carlo proposals
• Sequential Monte Carlo without likelihoods
• The rate of convergence for approximate Bayesian computation
• Weak informativity and the information in one prior relative to another

Cited in

(1)

• Rapid Bayesian Inference for Expensive Stochastic Models

This page was built for publication: Vector operations for accelerating expensive Bayesian computations - a tutorial guide

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6202922)