wolffd@0: function B = mk_unit_norm(A)
wolffd@0: % MK_UNIT_NORM Make each column be a unit norm vector
wolffd@0: % function B = mk_unit_norm(A)
wolffd@0: % 
wolffd@0: % We divide each column by its magnitude
wolffd@0: 
wolffd@0: 
wolffd@0: [nrows ncols] = size(A);
wolffd@0: s = sum(A.^2);
wolffd@0: ndx = find(s==0);
wolffd@0: s(ndx)=1; 
wolffd@0: B = A ./ repmat(sqrt(s), [nrows 1]);
wolffd@0: