Is there a way to handle numeric (or numarray or numpy) arrays as sets and compute efficiently their intersection, union, etc. ? I hope there is a faster way than s = array(set(A) & set(B)). Can this be done with masked arrays maybe ? I've never used them though and browsing through the docs didn't make clear if it is possible. Thanks, George