starting to add CMC algorithm in
running constrained (non-hamiltonian) monte carlo successfully, but it's not producing the right answer