You [mentioned](https://discourse.julialang.org/t/mcmc-landscape/25654/40) having a toy implementation for RHMC. Is the code available somewhere? I would like to experiment with it using [higher order AD for matrix expressions](http://www.matrixcalculus.org/matrixcalculus.pdf).