wolffd@0: function W = initializeDiag(X)
wolffd@0: 
wolffd@0:     [d,n,m] = size(X);
wolffd@0: 
wolffd@0: %     W = ones(d, 1);
wolffd@0:     W = 1./std(X,1,2);
wolffd@0:     W(isinf(W)) = 1;
wolffd@0: end