Hi Wes,
I do not profess to be an expert, but I have been off loading a fair number
of loops to C from Python code and achieved significant improvements most
have been of the following form (which I have found to be the fastest):
size = *incomingArrayObj->dimensions;
r_dptr = PyArray_DATA(resultArray);
while(size--) {
r_dptr = result;
r_dptr++;
}
Where for multidimensional arrays r_dptr could be incremented by the number
of dims rather than just ++:
dims = PyArray_DIM(incomingArrayObj,1);
i have not however actually used PyArray_SETITEM so cannot comment on the
issue you are having.
Hanni
2009/2/11 Wes McKinney <[email protected]>
> Hello,
>
> I am writing some Cython code and have noted that the buffer interface
> offers very little speedup for PyObject arrays. In trying to rewrite the
> same code using the C API in Cython, I find I can't get PyArray_SETITEM to
> work, in a call like:
>
> PyArray_SETITEM(result, <void *> iterresult.dataptr, obj)
>
> where result is an ndarray of dtype object, and obj is a PyObject*.
>
> Anyone have some experience with this can offer pointers (no pun
> intended!)?
>
> Thanks,
> Wes
>
> _______________________________________________
> Numpy-discussion mailing list
> [email protected]
> http://projects.scipy.org/mailman/listinfo/numpy-discussion
>
>
_______________________________________________
Numpy-discussion mailing list
[email protected]
http://projects.scipy.org/mailman/listinfo/numpy-discussion
