This is beautiful! Thank you Pauli.<br><br><div class="gmail_quote">On Sat, Dec 4, 2010 at 7:16 PM, Pauli Virtanen <span dir="ltr"><<a href="mailto:pav@iki.fi">pav@iki.fi</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
On Sat, 04 Dec 2010 19:00:43 +0000, Eleftherios Garyfallidis wrote:<br>
[clip]<br>
<div class="im">> I am using ndindex and then a for loop. Is there a better/faster way?<br>
<br>
</div>Yes:<br>
<br>
import numpy as np<br>
from numpy import newaxis<br>
<br>
x = np.zeros((200, 200, 200, 3))<br>
x[...,0] = np.arange(200)[:,newaxis,newaxis]<br>
x[...,1] = np.arange(200)[newaxis,:,newaxis]<br>
x[...,2] = np.arange(200)[newaxis,newaxis,:]<br>
x[1,3,2]<br>
# -> array([ 1.,  3.,  2.])<br>
<br>
Depending on what you use this array for, it's possible that you can<br>
avoid constructing it (and use broadcasting etc. instead).<br>
<font color="#888888"><br>
--<br>
Pauli Virtanen<br>
<br>
_______________________________________________<br>
NumPy-Discussion mailing list<br>
<a href="mailto:NumPy-Discussion@scipy.org">NumPy-Discussion@scipy.org</a><br>
<a href="http://mail.scipy.org/mailman/listinfo/numpy-discussion" target="_blank">http://mail.scipy.org/mailman/listinfo/numpy-discussion</a><br>
</font></blockquote></div><br>