On Fri, Nov 19, 2010 at 3:18 PM, Christopher Barker <Chris.Barker@noaa.gov> wrote:
On 11/19/10 11:19 AM, Keith Goodman wrote:
> On Fri, Nov 19, 2010 at 10:55 AM, Nathaniel Smith<njs@pobox.com>  wrote:
>> Why not make this a patch to numpy/scipy instead?
>
> My guess is that having separate underlying functions for each dtype,
> ndim, and axis would be a nightmare for a large project like Numpy.

True, but:

1) Having special-cases for the most common cases is not such a bad idea.

2) could one use some sort of templating approach to get all the dtypes
and such that you want?

3) as for number of dimensions, I don't think to would be to hard to
generalize that -- at least for contiguous arrays.


Note that the fmax/fmin versions can be sped up in the same way as sum.reduce was. Also, you should pass the flattened array to the routine for the axis=None case.

Chuck