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