Introduction In modern Bayesian statistics, we often face posterior distributions that are difficult to compute. Let be prior density and be likelihood. The standard approach to compute posterior is to use MCMC (like Metropolis-Hastings, Gibbs sampling and HMC).