1 files changed, 18 insertions, 0 deletions
diff --git a/sourcecodes/bnt-master/KPMtools/myismember.m b/sourcecodes/bnt-master/KPMtools/myismember.m
new file mode 100644
index 00000000..d339bb42
--- /dev/null
+++ b/sourcecodes/bnt-master/KPMtools/myismember.m
@@ -0,0 +1,18 @@
+function p = myismember(a,A)
+% MYISMEMBER Is 'a' an element of a set of positive integers? (much faster than built-in ismember)
+% p = myismember(a,A)
+
+%if isempty(A) | a < min(A) | a > max(A)  % slow
+
+if length(A)==0
+  p = 0;
+elseif a < min(A)
+  p = 0;
+elseif a > max(A)
+  p = 0;
+else
+  bits = zeros(1, max(A));
+  bits(A) = 1;
+  p = bits(a);
+end
+p = logical(p);