function   x=svec(X)
%  function   x=svec(X);
%  Use the isometry version
n=size(X,1);
E=ones(n);
indsu=find(triu(E,0));
indsuu=find(triu(E,1));
X(indsuu)=sqrt(2)*X(indsuu);
x=X(indsu);