I cleaned up code...most debug comments can be disabled.
Energy conservation is good, state space resoration is good (upon rejection).
There are now two options (velocity_options).  The reversal option is still
under develoopment, as kinetic energy is not conserved upon reversal.
The leapfrog integrator is called, requiring context creation, but
we are running through many trials, so it might be okay.