wolffd@0: function D = distanceFull(W, X)
wolffd@0: 
wolffd@0:     D = PsdToEdm(X' * W * X);
wolffd@0: end