Implementation of a Gibbs-Metropolis sampling algorithm in CUDA