  penm.ubYx{3} = 100;
% 
% *Starting point* 
%
% is an optional input. When set, it should be ideally feasible. For a matrix
% variable |Y{i}|, the starting point |Yinit{i}| must be a symmetric matrix
% of the same non-zero pattern as |Y{i}|.
%
%   penm.xinit=rand(penm.Nx,1);
%   penm.Yinit{1} = 3.*eye(n);
% 
% The starting point might be automatically modified within solve()
% * if barrier (for example, penm.lbxbar) is used and the point is not
%   feasible w.r.t. the box constraints (or very 