インデックス作成を行うための for ループがあります。
for (int i=0; i<N; i++){
a[i] = b[c[i]]
}
cは対象のインデックスであり、 are int *、 whilebおよびaarefloat *および操作された値です。
しかし、これには長い時間がかかります (それほど長くかかることはありません)。BLAS/LAPLACK/etc にある可能性が高いベクトル化バージョンが欲しいです。
を探していnested_indexing(float * output_vector, float * input_vector, int * input_indices)ます。
docsを調べてみましたが、何も見つかりませんでした。