Daniel@0: function p = myismember(a,A)
Daniel@0: % MYISMEMBER Is 'a' an element of a set of positive integers? (much faster than built-in ismember)
Daniel@0: % p = myismember(a,A)
Daniel@0: 
Daniel@0: %if isempty(A) | a < min(A) | a > max(A)  % slow
Daniel@0: 
Daniel@0: if length(A)==0
Daniel@0:   p = 0;
Daniel@0: elseif a < min(A)
Daniel@0:   p = 0;
Daniel@0: elseif a > max(A)
Daniel@0:   p = 0;
Daniel@0: else
Daniel@0:   bits = zeros(1, max(A));
Daniel@0:   bits(A) = 1;
Daniel@0:   p = bits(a);
Daniel@0: end
Daniel@0: p = logical(p);