sidx = np.argsort(searchvals)
left_idx = np.searchsorted(searchvals, values, sorter=sidx, side='left')
right_idx = np.searchsorted(searchvals, values, sorter=sidx, side='right')
out = np.where(left_idx != right_idx)[0]
Массив Numpy int: поиск индексов нескольких целевых целых чисел