Set operations for arrays based on sorting.
For floating point arrays, inaccurate results may appear due to usual round-off
and floating point comparison issues.
Speed could be gained in some operations by an implementation of
numpy.sort, that can provide directly the permutation vectors, thus avoiding
calls to numpy.argsort.
Original author: Robert Cimrman