(2018) relax the above constraint by penalizing the Wasserstein distance between q(z) and p(z).
it is computed by the Maximum Mean Discrepancy (MMD), defined as where P (z) and Q(z) are cumulative density functions.