matrix_hmc
