intersect2.m

来自「CheckMate is a MATLAB-based tool for mod」· M 代码 · 共 26 行

function CON = intersect2(CONfeas,CONadd, varargin)
    
    h1 = CONfeas;
    h2 = CONadd;
    
    H.A = [h1.CI; h2.CI];
    H.B = [h1.DI; h2.DI];
    
    linstart = size(H.B,1) + 1;
    
    H.A = [H.A; h1.CE; h2.CE];
    H.B = [H.B; h1.DE; h2.DE];
    linend = size(H.B, 1);
    
    if (linend >= linstart)
        H.lin = linstart: linend;
    end
    
    h = cddmex('reduce_h', H);
    
    CE = h.A(h.lin,:);
    DE = h.B(h.lin,:);
    CI = h.A;
    DI = h.B;
    CON = linearicon(CE,DE,CI,DI);