function v=sHvec(H)
%function H=sHMat(v,n)
%Input:  v = [v1;v2;d]  vi in R^(t(n)-n)  d in R^n
% Output:  H  Hermitian
sqrt2=sqrt(2);
n=size(H,1);
%svind: linear indices of the **strictly** upper triangle by columns
svind = logical(triu(ones(n),1));
RH=real(H);
IH=imag(H);
v=real([sqrt2*RH(svind);sqrt2*IH(svind);diag(H)]); % removes 0i
