I would like to know if there is any way to make the following operation faster.
def test(): shape=(200,200,200,3) refinds = np.ndindex(shape[:3]) reftmp=np.zeros(shape) for ijk_t in refinds: i,j,k = ijk_t reftmp[i,j,k,0]=i reftmp[i,j,k,1]=j reftmp[i,j,k,2]=k
%timeit test() 1 loops, best of 3: 19.5 s per loop
I am using ndindex and then a for loop. Is there a better/faster way?
Thank you, Eleftherios