[Numpy-discussion] Problems with Numexpr and discontiguous arrays
Tim Hochberg
tim.hochberg at ieee.org
Wed Oct 4 13:19:08 EDT 2006
Ivan Vilata i Balaguer wrote:
> It seemed that discontiguous arrays worked OK in Numexpr since r1977 or
> so, but I have come across some alignment or striding problems which can
> be seen with the following code::
I looked at this just a little bit and clearly this bit from interp_body
cannot work in the presence of recor arrays:
//....
intp sf1 = sb1 / sizeof(double); \
//...
#define f1 ((double *)x1)[j*sf1]
There are clearly some assumptions that sb1 is evenly divisible by
sizeof(double). Blech!. This is likely my fault, and I expect it won't
be too horrible to fix, but I don't know that I'll have time immediately.
-tim
More information about the NumPy-Discussion
mailing list