function c=Cvec(C)
%function c=Cvec(C)
%INPUT:   C n by n complex matrix
%OUTPUT:  c = [real(C(:);imag(C(:)]
c=[real(C(:));imag(C(:))];
end  % end of function
