Есть ли библиотека Nim, содержащая реализацию `argsort` или оболочку?
Я ищу версию argsort
Например, существует в NumPy или в Fortran. Есть ли реализация argsort
в Ним... или доступны Ним в какой-то библиотеке? Кажется немного удивительным, что это отсутствует.
ОБНОВИТЬ
Следующее, кажется, работает для argsort
:
proc argsort[T](a : T) : seq[int] =
result = toSeq(0..a.len - 1)
sort(result, proc (i, j: int): int = cmp(a[i], a[j]))
Предположительно, однако, это может быть более эффективно написано изначально и избегая указателя функции....
0 ответов
Он еще не стал шустрым, но у Arraymancer есть argsort
Вот
import arraymancer,algorithm,sequtils
proc argsort*[T](t: Tensor[T], order = SortOrder.Ascending): Tensor[int] =
## Returns the indices which would sort `t`. Useful to apply the same sorting to
## multiple tensors based on the order of the tensor `t`.
##
## Sorts the raw underlying data!
# TODO: should we clone `t` so that if `t` is a view we don't access the whole
# data?
assert t.rank == 1, "Only 1D tensors can be sorted at the moment!"
proc cmpIdxTup(x, y: (T, int)): int = system.cmp(x[0], y[0])
# make a tuple of input & indices
var tups = zip(toOpenArray(t.storage.Fdata, 0, t.size - 1),
toSeq(0 ..< t.size))
# sort by custom sort proc
tups.sort(cmp = cmpIdxTup, order = order)
result = newTensorUninit[int](t.size)
for i in 0 ..< t.size:
result[i] = tups[i][1]
let x = @[3,9,4,1,5].toTensor()
echo x.argsort()
производит:
Tensor[system.int] of shape [5]" on backend "Cpu"
3 0 2 4 1