sampling runs without crashing.  possible problems in velocity update.
fixed indexing bug in sampler
HMC seems to work.  log likelihood improves during learning.