  A method for the introduction of second-order derivatives of the log
likelihood into HMC algorithms is introduced, which does not require the
Hessian to be evaluated at each leapfrog step but only at the start and end of
trajectories.
