Daniel@0: function C = myintersect(A,B)
Daniel@0: % MYINTERSECT Intersection of two sets of positive integers (much faster than built-in intersect)
Daniel@0: % C = myintersect(A,B)
Daniel@0: 
Daniel@0: A = A(:)'; B = B(:)';
Daniel@0: 
Daniel@0: if isempty(A)
Daniel@0:   ma = 0;
Daniel@0: else
Daniel@0:   ma = max(A);
Daniel@0: end
Daniel@0: 
Daniel@0: if isempty(B)
Daniel@0:   mb = 0;
Daniel@0: else
Daniel@0:   mb = max(B);
Daniel@0: end
Daniel@0: 
Daniel@0: if ma==0 | mb==0
Daniel@0:   C = [];
Daniel@0: else
Daniel@0:   %bits = sparse(1, max(ma,mb));
Daniel@0:   bits = zeros(1, max(ma,mb));
Daniel@0:   bits(A) = 1;
Daniel@0:   C = B(logical(bits(B)));  
Daniel@0: end
Daniel@0: 
Daniel@0: %sum( bitget( bitand( cliquesb(i), cliquesb(j) ), 1:52 ) );