April 30, 2008
4:29 a.m.
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