Hamiltonian Monte Carlo on Symmetric and Homogeneous Spaces via Symplectic Reduction

by   Alessandro Barp, et al.

The Hamiltonian Monte Carlo method generates samples by introducing a mechanical system that explores the target density. For distributions on manifolds it is not always simple to perform the mechanics as a result of the lack of global coordinates, the constraints of the manifold, and the requirement to compute the geodesic flow. In this paper we explain how to construct the Hamiltonian system on naturally reductive homogeneous spaces using symplectic reduction, which lifts the HMC scheme to a matrix Lie group with global coordinates and constant metric. This provides a general framework that is applicable to many manifolds that arise in applications, such as hyperspheres, hyperbolic spaces, symmetric positive-definite matrices, Grassmannian, and Stiefel manifolds.


page 1

page 2

page 3

page 4


Note on the geodesic Monte Carlo

Geodesic Monte Carlo (gMC) comprises a powerful class of algorithms for ...

Hamiltonian Monte-Carlo for Orthogonal Matrices

We consider the problem of sampling from posterior distributions for Bay...

Simulating sticky particles: A Monte Carlo method to sample a Stratification

Many problems in materials science and biology involve particles interac...

Dynamically rescaled Hamiltonian Monte Carlo for Bayesian Hierarchical Models

Dynamically rescaled Hamiltonian Monte Carlo (DRHMC) is introduced as a ...

Hamiltonian Flow Simulation of Rare Events

Hamiltonian Flow Monte Carlo(HFMC) methods have been implemented in engi...

Optimization on manifolds: A symplectic approach

There has been great interest in using tools from dynamical systems and ...

Magnetic Manifold Hamiltonian Monte Carlo

Markov chain Monte Carlo (MCMC) algorithms offer various strategies for ...