TestConPolyLogicalExactLogic.m

clear all
 
% dim=7;
% c1= [logical(randi([0 1],dim,1)) ];
% % for i=1:4
% %     gI1{i} =logical(randi([0 1],dim,1));
% % end
% %{randi([0 1],dim,1)};
% numOfInd1 = 3;
% for i=1:numOfInd1
%    g1{i} = randi([0 1],dim,1);
% end
% E1 = [1 0 1;1 1 0];
% id1 = [1 ;2];
% % g1 ={};
% % E1=[];
% 
% c2= [logical(randi([0 1],dim,1)) ];
% % for i=1:4
% %     gI2{i} =logical(randi([0 1],dim,1));
% % end
% for i=1:2
%     g2{i} = logical(randi([0 1],dim,1));
% end
% E2 =  [1 0;0 1];
% id2 = [3; 4];
% % g2 = {};
% % E2=[];
% 
% 
% 
% 
% 
% pZ1 = logicalPolyZonotope(c1,g1,E1,id1);
% pZ2 = logicalPolyZonotope(c2,g2,E2,id2);


% pp1= logical([0 1 0 round(rand);0 1 1 round(rand); 1 0 1 round(rand);0 1 1 round(rand); 0 0 0 round(rand); 1 1 1 round(rand)]);
% pp2= logical([0 0 1 0 round(rand);0 1 1 1 round(rand);0 1 0 0 round(rand); 0 0 1 1 round(rand); 0 0 0 0 round(rand); 1 1 1 1 round(rand)]);
% pp3= logical([0 1 0 round(rand);1 0 1 round(rand); 0 0 1 round(rand);1 0 1 round(rand); 0 0 0 round(rand); 1 1 1 round(rand)]);
% 
% plz1 = logicalPolyZonotope.enclosePoints(pp1);
% plz2 = logicalPolyZonotope.enclosePoints(pp2);
% plz3 = logicalPolyZonotope.enclosePoints(pp3);
% 
% evaluate(plz1);
% evaluate(plz2);
% 
% pc1=logicalConPolyZonotope(plz1.c,plz1.G,plz1.E,[],[],[],plz1.id);
% pc2=logicalConPolyZonotope(plz2.c,plz2.G,plz2.E,[],[],[],plz2.id);
% pc3=logicalConPolyZonotope(plz3.c,plz3.G,plz3.E,[],[],[],plz3.id);
% 
% pZ1=intersect(pc1,pc2);
% pZ2=intersect(pc1,pc3);


pp1= logical([0 1 0 round(rand);0 1 1 round(rand); 1 0 1 round(rand);0 1 1 round(rand); 0 0 0 round(rand); 1 1 1 round(rand)]);
pp2= logical([0 0 1 0 round(rand);0 1 1 1 round(rand);0 1 0 0 round(rand); 0 0 1 1 round(rand); 0 0 0 0 round(rand); 1 1 1 1 round(rand)]);
pp3= logical([0 1 0 round(rand);1 0 1 round(rand); 0 0 1 round(rand);1 0 1 round(rand); 0 0 0 round(rand); 1 1 1 round(rand)]);

plz1 = logicalPolyZonotope.enclosePoints(pp1);
plz2 = logicalPolyZonotope.enclosePoints(pp2);
plz3 = logicalPolyZonotope.enclosePoints(pp3);

evaluate(plz1);
evaluate(plz2);

pc1=logicalConPolyZonotope(plz1.c,plz1.G,plz1.E,[],[],[],plz1.id);
pc2=logicalConPolyZonotope(plz2.c,plz2.G,plz2.E,[],[],[],plz2.id);
pc3=logicalConPolyZonotope(plz3.c,plz3.G,plz3.E,[],[],[],plz3.id);

pZ1=intersect(pc1,pc2);
pZ2=intersect(pc1,pc3);
pZ3=intersect(pc2,pc3);



%plot(Z1,[1 2 3],'r*');
% manually
p1=evaluate(pZ1);
p2=evaluate(pZ2);
[~,p1Col]=size(p1);
[~,p2Col]=size(p2);
ExactXorPoints=[];
ExactAndPoints=[];
ExactNandPoints=[];
ExactOrPoints=[];
ExactAndORFormulaPoints =[] ;
tic
for i =1:p1Col
        ExactXorPoints = [ExactXorPoints, xor(p1(:,i),p1(:,i))];
        ExactAndPoints = [ExactAndPoints, and(p1(:,i),p1(:,i))];
        ExactNandPoints =[ExactNandPoints, nand(p1(:,i),p1(:,i))];
        ExactOrPoints =[ExactOrPoints, or(p1(:,i),p1(:,i))];
end

for i =1:p1Col
    for j=1:p2Col
        ExactAndORFormulaPoints = [ExactAndORFormulaPoints, or(p1(:,i),and(not(p1(:,i)),p2(:,j)))];
    end
end
exeTime=toc;


ExactXorPoints = unique(ExactXorPoints','rows')';
ExactAndPoints = unique(ExactAndPoints','rows')';
ExactNandPoints = unique(ExactNandPoints','rows')';
ExactOrPoints = unique(ExactOrPoints','rows')';
ExactAndORFormulaPoints = unique(ExactAndORFormulaPoints','rows')';

pZ1ExactXorpZ1 = exactXor(pZ1,pZ1);
pZ1ExactAndpZ1 = exactAnd(pZ1,pZ1);
pZ1ExactOrpZ1 = exactOr(pZ1,pZ1);
pZ1ExactNandpZ1 = exactNand(pZ1,pZ1);

tempAnd= and(not(pZ1),pZ2);
%tempAnd.id = [1;2;1;2];
pZ1ExactAndORpZ2 = exactOr(pZ1,tempAnd);


pZ1XorpZ1Points= evaluate(pZ1ExactXorpZ1);
pZ1AndpZ1Points= evaluate(pZ1ExactAndpZ1);
pZ1OrpZ1Points= evaluate(pZ1ExactOrpZ1);
pZ1NandpZ1Points= evaluate(pZ1ExactNandpZ1);
pZ1ExactAndORpZ2Points= evaluate(pZ1ExactAndORpZ2);

%% checks
if isequal(ExactXorPoints,pZ1XorpZ1Points) && isequal(ExactNandPoints,pZ1NandpZ1Points) && isequal(ExactAndPoints,pZ1AndpZ1Points) %&& isequal(ExactAndORFormulaPoints,pZ1ExactAndORpZ2Points)
    disp('All tests passed')
else
    disp('Test problem(s)')
end


if ~isequal(ExactXorPoints,pZ1XorpZ1Points)
    disp('Exact Xor problem')
end

if ~isequal(ExactNandPoints,pZ1NandpZ1Points)
    disp('Exact Nand problem')
end

if ~isequal(ExactAndPoints,pZ1AndpZ1Points)
    disp('Exact And problem')
end

if ~isequal(ExactOrPoints,pZ1OrpZ1Points)
    disp('Exact Or problem')
end

% if ~isequal(ExactAndORFormulaPoints,pZ1ExactAndORpZ2Points)
%     disp('exactOr(pZ1,and(pZ1,pZ2)) problem')
% end