Function to calculate the kernel mean embedding for to distributional data
sets. It uses the empirical approximation for the integral
$$\int_{\mathcal X} \int_{\mathcal Y} K(x, y) d P_X d Q_Y $$ for a given
kernel \(K(\cdot, \cdot)\). Currently only supports radial basis function
kernel for fast computation.