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