code: def search_indices(values, source):
      
      orig_indices = source.argsort()
      return orig_indices[np.searchsorted(source[orig_indices], values)]
docs: Given a set of values returns the indices of each of those values
      in the source array.