Particle filter efficiency under limited communication

04/21/2019
by   Deborshee Sen, et al.
0

Sequential Monte Carlo (SMC) methods are typically not straightforward to implement on parallel architectures. This is because standard resampling schemes involve communication between all particles in the system. In this article, we consider the α-SMC algorithm, a generalised particle filter algorithm with limited communications between particles. The communication between different particles is controlled through the so-called α-matrices. We study the influence of the communication structure on the convergence and stability properties of the resulting algorithms. We prove that under standard assumptions, it is possible to use randomised communication structures where each particle only communicates with a small number of neighbouring particles while still having good mixing properties, and this ensures that the resulting algorithms are stable in time and converge at the usual Monte Carlo rate. A particularly simple approach to implement these ideas consists of choosing the α-matrices as the Markov transition matrices of random walks on Ramanujan graphs.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
02/04/2022

De-Sequentialized Monte Carlo: a parallel-in-time particle smoother

Particle smoothers are SMC (Sequential Monte Carlo) algorithms designed ...
research
06/18/2018

Stability of Conditional Sequential Monte Carlo

The particle Gibbs (PG) sampler is a Markov Chain Monte Carlo (MCMC) alg...
research
09/18/2019

Collective sampling through a Metropolis-Hastings like method: kinetic theory and numerical experiments

The classical Metropolis-Hastings algorithm provides a simple method to ...
research
05/05/2023

Counterexample to a transition probability formula for the ancestral process

We consider weighted particle systems in which new generations are re-sa...
research
02/08/2018

A local parallel communication algorithm for polydisperse rigid body dynamics

The simulation of large ensembles of particles is usually parallelized b...
research
11/13/2015

k-means: Fighting against Degeneracy in Sequential Monte Carlo with an Application to Tracking

For regular particle filter algorithm or Sequential Monte Carlo (SMC) me...
research
10/11/2021

Weak Convergence of Non-neutral Genealogies to Kingman's Coalescent

Interacting particle populations undergoing repeated mutation and fitnes...

Please sign up or login with your details

Forgot password? Click here to reset