Algorithms for mean-field variational inference via polyhedral optimization in the Wasserstein space