Sampling with Riemannian Hamiltonian Monte Carlo in a Constrained Space

Neural Information Processing Systems 

Our algorithm incorporates constraints into the Riemannian version of Hamiltonian Monte Carlo and maintains sparsity.