[Numpy-discussion] Possible to use numexpr with user made ufuncs/scipy ufuncs?

[Francesc Alted]

A Saturday 26 June 2010 19:17:43 Pauli Virtanen escrigué:
> But what if such an expression does not exist? For example, there is no
> finite closed form expression for the Bessel functions in terms of
> elementary functions. An accurate implementation is rather complicated,
> and usually piecewise defined.
> 
> For convenience reasons, it could be useful if one could do something like
> 
>     numexpr.evaluate("cos(iv(0, x))", functions=dict(iv=scipy.special.iv))
> 
> and this would be translated to numexpr bytecode that would make a Python
> function call to obtain "iv(0, x)" for each block of data required,
> assuming "iv" is a vectorized function. It's of course possible to
> precompute the value of "iv(0, x)", but this is extra hassle and requires
> additional memory.

Yeah, I can see the merit of implementing such a thing, mainly for avoiding 
additional memory consumption.

But again, the nice thing would be to implement such a special functions in 
terms of numexpr expressions so that the evaluation itself can be faster.  
Admittedly, that would take a bit more time.

Anyway, if someone comes with patches for implementing this, I'd glad to 
commit them.
 
-- 
Francesc Alted