[Numpy-discussion] another broadcast/fancy indexing question

Mark.Miller mpmusu@cc.usu....
Tue Jul 10 14:48:17 CDT 2007


Just ran across something that doesn't quite make sense to me at the moment.

Here's some code:

 >>> numpy.__version__
'1.0.2'
 >>>
 >>> def f1(b,c):
	b=b.astype(int)
	c=c.astype(int)
	return b,c

 >>> b,c = numpy.fromfunction(f1,(5,5))
 >>> a=numpy.zeros((2,12,5,5),int)
 >>> a1=a[0]
 >>> a1[:,b,c].shape
(12, 5, 5)
 >>> a[0,:,b,c].shape
(5, 5, 12)   ###why does this not return (12,5,5)?
 >>>

So in a nutshell, it's not completely clear to me why these are 
returning arrays of different shapes.  Can someone shed some light?

Thanks,

-Mark


More information about the Numpy-discussion mailing list