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