I have the following function:
l2norm_squared = function(X) Re(sum(Conj(t(X)) %*% X)) fast_mmd = function(X, Y, nBasis = 2^6, seed = NULL) { if (!is.nul