Это может быть сделано в O (n) с использованием расширенного назначения:
my_vector = np.array([9., 6., 21., 17., 12.])
my_indices = np.array([1., 4., 2., 0., 3.])
my_result = np.empty_like(my_vector)
my_result[my_indices.astype(int)] = my_vector
my_result
# array([17., 9., 21., 12., 6.])
Сроки на более крупном примере:
from timeit import timeit
a = np.random.random(1000000)
b = np.random.permutation(1000000).astype(float)
def fargsort():
return(a[b.argsort()])
def fassign():
c = np.empty_like(a)
c[b.astype(int)] = a
return c
np.all(fargsort() == fassign())
# True
timeit(fargsort,number=10)
# 1.2764090860000579
timeit(fassign,number=10)
# 0.17262099700747058