ables |Y|. By default a        
% reciprocal penalty/barrier function is used which might cause that    
% the matrix variables are feasible only up to the given tolerance      
% (see option settings). It is possible to switch some or all of them   
% to logaritmic barrier which secures strict feasibility of |Y| during  
% the whole computation. The way to do it is similar to |lbxbar| and    
% |ubxbar|. All such matrix variable bounds need to be listed in arrays 
% penm.lbYbar and penm.ubYbar. Any 