Share your thoughts, 1 month free Claude Pro on usSee more
WorkDL logo mark

Affine invariant interacting Langevin dynamics for Bayesian inference

About

We propose a computational method (with acronym ALDI) for sampling from a given target distribution based on first-order (overdamped) Langevin dynamics which satisfies the property of affine invariance. The central idea of ALDI is to run an ensemble of particles with their empirical covariance serving as a preconditioner for their underlying Langevin dynamics. ALDI does not require taking the inverse or square root of the empirical covariance matrix, which enables application to high-dimensional sampling problems. The theoretical properties of ALDI are studied in terms of non-degeneracy and ergodicity. Furthermore, we study its connections to diffusion on Riemannian manifolds and Wasserstein gradient flows. Bayesian inference serves as a main application area for ALDI. In case of a forward problem with additive Gaussian measurement errors, ALDI allows for a gradient-free approximation in the spirit of the ensemble Kalman filter. A computational comparison between gradient-free and gradient-based ALDI is provided for a PDE constrained Bayesian inverse problem.

Alfredo Garbuno-Inigo, Nikolas N\"usken, Sebastian Reich• 2019

Related benchmarks

TaskDatasetResultRank
Synthetic Distribution SamplingGMM2
KSD0.93
7
Synthetic Distribution SamplingFunnel20
KSD28.66
7
Synthetic Distribution SamplingFunnel10
KSD20.34
7
Synthetic Distribution SamplingGMM5
KSD1.516
7
Synthetic Distribution SamplingGMM10
KSD2.185
7
Synthetic Distribution SamplingJoker
KSD9.286
7
Synthetic Distribution SamplingFunnel2
KSD9.008
7
Synthetic Distribution SamplingFunnel5
KSD14.43
7
Synthetic Distribution SamplingGMM20
KSD3.355
7
Synthetic Distribution SamplingHimmelblau
KSD14.33
6
Showing 10 of 10 rows

Other info

Follow for update