Use ps_point for HMC temporaries.  Introduced NUTS (not sure NUTS algorithm is working yet)