Constrained Gaussian Wasserstein Optimal Transport with Commutative Covariance Matrices