Increment Numpy multi-d array with repeated indices

Question

I\'m interested in the multi-dimensional case of Increment Numpy array with repeated indices.

I have an N-dimensional array and a set N index arrays, who\'s values I

Answer 1

How about this:

import numpy as np
a = np.zeros((2,3,4))
i = np.array([0,0,1])
j = np.array([0,0,1])
k = np.array([0,0,3])

ijk = np.vstack((i,j,k)).T
H,edge = np.histogramdd(ijk,bins=a.shape)
a += H  

Answer 2

I don't know if there is an easier solution with direct array indexing, but this works:

for x,y,z in zip(i,j,k):
    a[x,y,z] +=1