On the computation of Wasserstein barycenters