[Numpy-discussion] numpy.mean still broken for largefloat32arrays

[Sturla Molden]

Sturla Molden <sturla.molden at gmail.com> wrote:
> Sebastian Berg <sebastian at sipsolutions.net> wrote:
> 
>> Yes, it is much more complicated and incompatible with naive ufuncs if
>> you want your memory access to be optimized. And optimizing that is very
>> much worth it speed wise...
> 
> Why? Couldn't we just copy the data chunk-wise to a temporary buffer of say
> 2**13 numbers and then reduce that? I don't see why we need another
> iterator for that.

I am sorry if this is a stupid suggestion. My knowledge of how NumPy ufuncs
works could have been better.

Sturla