[Numpy-discussion] Optimizing mean(axis=0) on a 3D array

[Martin Spacek]

Tim Hochberg wrote:

> I'm actually surprised that the take version is faster than my original 
> version since it makes a big  ol' copy. I guess this is an indication 
> that indexing is more expensive than I realize. That's why nothing beats 
> measuring!

Actually, your original version is just as fast as the take() version. 
Both are about 9X faster than numpy.mean() on my system. I prefer the 
take() version because you only have to pass a single argument to 
mean_accum()

Martin