[Numpy-discussion] dot/tensordot limitations

[Travis E. Oliphant]

Nadav Horesh wrote:
> You open here a Pandora box:
> What should I do if I want to run an operation on elements of an array which are not in the library. The usual answers are either use more memory ( (A*B).sum(axis=1) ), or more time ([dot(A[i],B[i]) for i in len(A)]).
>
> Would your issue can be rephrase like this:
> A way to iterate over array in the C level, where the function/operation is defined in the C level, without the overhead of python?
>   
My plan for this is the general function listed on the NumPy Trac area 
along with a weave-like kernel creation from Python code.

http://www.scipy.org/scipy/numpy/wiki/GeneralLoopingFunctions

I'd love to get time to do this or mentor someone else in doing it.

-Travis