wolffd@0: function W = initializeFull(X)
wolffd@0: 
wolffd@0:     [d,n,m] = size(X);
wolffd@0: 
wolffd@0:     W = zeros(d,d,m);
wolffd@0:     for i = 1:m
wolffd@0:         W(:,:,i) = diag(1./std(X(:,:,i),1,2));
wolffd@0:     end
wolffd@0:     W(isinf(W)) = 1;
wolffd@0: end