Python で (for ループなしで) 高速な方法を見つけて、配列の繰り返しインデックスを割り当てたいと考えています。これは、for ループを使用した望ましい結果です。
import numpy as np
a=np.arange(9, dtype=np.float64).reshape((3,3))
# The array indices: [2,3,4] are identical.
Px = np.uint64(np.array([0,1,1,1,2]))
Py = np.uint64(np.array([0,0,0,0,0]))
# The array to be added at the array indices (may also contain random numbers).
x = np.array([.1,.1,.1,.1,.1])
for m in np.arange(len(x)):
a[Px[m]][Py[m]] += x
print a
%[[ 0.1 1. 2.]
%[ 3.3 4. 5.]
%[ 6.1 7. 8.]]
インデックスに追加しようとするx
と、明らかに同じ結果が得られません (3.3 対 3.1):a
Px,Py
a[Px,Py] += x
print a
%[[ 0.1 1. 2.]
%[ 3.1 4. 5.]
%[ 6.1 7. 8.]]
numpyでこれを行う方法はありますか? ありがとう。