wolffd@0: function D = distanceFullMKL(W, X)
wolffd@0: 
wolffd@0:     [d, n, m] = size(X);
wolffd@0: 
wolffd@0:     D = zeros(n);
wolffd@0:     parfor i = 1:m
wolffd@0:         D = D + PsdToEdm(X(:,:,i)' * W(:,:,i) * X(:,:,i));
wolffd@0:     end
wolffd@0: end