Bayesian optimization for computationally extensive probability distributions