remove_implicit_linear.m

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

function [CEnew, dEnew, CInew, dInew] = remove_implicit_linear(CE, dE, CI, dI)
% remove implicit linearity in the linear constraints
    
    m = size(dI,1);
    d = size(dE,1);

    H.A = [CI; CE];
    H.B = [dI; dE];
    H.lin = m+1: m+d;
    
    r = cddmex('implicit_linear',H);
    
    if ~isempty(r)
        selector = true(size(CI,1),1);
        selector(r) = false;
        CInew = CI(selector, :);
        dInew = dI(selector, :);
        
        Rimp = [CI(r, :) dI(r, :);
                CE dE];
        [R, jb] = rref(Rimp);
        Radd = Rimp(jb, :);
        CEnew = Radd(:,1:end-1);
        dEnew = Radd(:,end);
        
    else
        CInew = CI;
        dInew = dI;
        CEnew = CE;
        dEnew = dE;
    end
    
end