Bayesian inference for the multivariate skew-normal model: A population Monte Carlo approach