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

[Martin Spacek]

Martin Spacek wrote:

> 
> 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()

I forgot to mention that all my indices are, for now, sorted. I just 
tried shuffling them (as you did), but I still get the same 9x 
improvement in speed, so I don't know why you only get a 4x improvement 
on your system.

Martin