function d = mhdist0(X)
% Mahalanobis distance for N(0,I)
% returns the square distance

d = diag(X*X'); % for example

% equivalently, apply norm(...,2) for each row
return
