function H=sHMat(v)
%function H=sHMat(v)
%Input:  v = [v1;v2;d]  vi in R^(t(n)-n)  d in R^n
% Output:  H  Hermitian
n=round(sqrt(length(v)));
if n^2-length(v)>0
	disp('WARNING!!! catastrophic error; sHMat index is wrong')
end
tn = n*(n+1)/2; %number of variables in formulation
tnm = tn-n;
H=usHMat(v(1:tnm),n) +i*uskewHMat(v(tnm+1:2*tnm),n) +diag(v(2*tnm+1:end));
