function v=Hvec(H)
%function H=HMat(v,n)
%  same as sHvec but with multiplication by sqrt 2
%Input:  v = [v1;v2;d]  vi in R^(t(n)-n)  d in R^n
% Output:  H  Hermitian
%    
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=[RH(svind);IH(svind);diag(H)];
