While writing a functional testing notebook for this sampler, I noticed very poor performance on the correlated Gaussian even with a large number of iterations, indicating a possible bug.
[x] Check equations for relativistic mass and kinetic energy.
[x] Check the distribution for momentum. In the paper they show that it has a hyperbolic distribution, while we are using the multivariate normal from standard Hamiltonian.
While writing a functional testing notebook for this sampler, I noticed very poor performance on the correlated Gaussian even with a large number of iterations, indicating a possible bug.