[Numpy-discussion] single precision routines in NumPy?

I would like a single precision version of 'interp' in the Numeric Core. (I want such a routine since I'm operating on huge single precision arrays, that I don't want promoted to double precision.) I've written such a routine, but Paul Dubois and I are discussing the best way of integrating it into the core. One solution is to simply add a new function 'interpf' to arrayfnsmodule.c . Another solution is to add a typecode=Float option to interp. Any opinions on how this single precision version be handled? -- Joe VanAndel National Center for Atmospheric Research http://www.atd.ucar.edu/~vanandel/ Internet: vanandel@ucar.edu