simple code to make sure samplers converge to right distribution for Gaussian
Define momentum corruption rate in terms of alpha rather than beta