Related Experiment Video
Updated: Apr 6, 2026

Image-based Lagrangian Particle Tracking in Bed-load Experiments
Published on: July 20, 2017
Markov Chain Monte Carlo from Lagrangian Dynamics
Shiwei Lan1, Vasileios Stathopoulos2, Babak Shahbaba1
1Department of Statistics, University of California, Irvine, Irvine, CA 92697, USA.
Hamiltonian Monte Carlo (HMC) improves the computational e ciency of the Metropolis-Hastings algorithm by reducing its random walk behavior. Riemannian HMC (RHMC) further improves the performance of HMC by exploiting the geometric properties of the parameter space. However, the geometric integrator used for RHMC involves implicit equations that require fixed-point iterations. In some cases, the computational overhead for solving implicit equations undermines RHMC's benefits. In an attempt to circumvent this problem, we propose an explicit integrator that replaces the momentum variable in RHMC by velocity. We show that the resulting transformation is equivalent to transforming Riemannian Hamiltonian dynamics to Lagrangian dynamics. Experimental results suggests that our method improves RHMC's overall computational e ciency in the cases considered. All computer programs and data sets are available online (http://www.ics.uci.edu/~babaks/Site/Codes.html) in order to allow replication of the results reported in this paper.
Hamiltonian Monte Carlo (HMC) improves the computational e ciency of the Metropolis-Hastings algorithm by reducing its random walk behavior. Riemannian HMC (RHMC) further improves the performance of HMC by exploiting the geometric properties of the parameter space. However, the geometric integrator used for RHMC involves implicit equations that require fixed-point iterations. In some cases, the computational overhead for solving implicit equations undermines RHMC's benefits. In an attempt to circumvent this problem, we propose an explicit integrator that replaces the momentum variable in RHMC by velocity. We show that the resulting transformation is equivalent to transforming Riemannian Hamiltonian dynamics to Lagrangian dynamics. Experimental results suggests that our method improves RHMC's overall computational e ciency in the cases considered. All computer programs and data sets are available online (http://www.ics.uci.edu/~babaks/Site/Codes.html) in order to allow replication of the results reported in this paper.
Related Concept Videos
Principle of Linear Impulse and Momentum for a System of Particles
Notably, internal forces between particles, occurring in equal and opposite collinear pairs, cancel out and are not part of the equation of motion. This exclusion simplifies the...
Mechanistic Models: Compartment Models in Algorithms for Numerical Problem Solving
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...
Linear Momentum in Control Volume
Poisson's And Laplace's Equation
Eulerian and Lagrangian Flow Descriptions
The Eulerian method focuses on fixed points in space where fluid properties, such as velocity, pressure, and temperature, are observed as the fluid moves between these...
Principle of Linear Impulse and Momentum for a Single Particle
Delving...

