Utility functions in Bayesian experimental design are usually based on the posterior distribution. When the posterior is found by simulation. it must be sampled from for each future dataset drawn from the prior predictive distribution. Many thousands of posterior distributions are often required. https://www.remtavares.com/