function v=sHvectwo(H)
%function H=sHMattwo(v,n)
%% sHmat function without sqrt 2
% Input:  H  Hermitian
%Output:  v = [v1;v2;d]  vi in R^(t(n)-n)  d in R^n
%sqrt2=sqrt(2);    % not used!!
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([RH(svind);IH(svind);diag(H)]); % removes 0i
%v=real([sqrt2*RH(svind);sqrt2*IH(svind);diag(H)]); % removes 0i
