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