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