[Numpy-discussion] Problems with Numexpr and discontiguous arrays

Ivan Vilata i Balaguer ivilata at carabos.com
Thu Oct 5 03:42:05 CDT 2006

En/na Tim Hochberg ha escrit::

> 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). [...]

I noticed something strange in those statements when implementing
support for strings, and I must confess that I didn't grasp their
meaning, so I implemented it a little differently for strings::

    #define s1    ((char   *)x1 + j*params.memsteps[arg1])

That seemed to work, but it might not be right (though I tested a bit),
and certainly it may not be efficient enough.  Here you have my previous
patches if you want to have a look at how I (try to) do it:



	Ivan Vilata i Balaguer   >qo<   http://www.carabos.com/
	       Cárabos Coop. V.  V  V   Enjoy Data

-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 309 bytes
Desc: OpenPGP digital signature
Url : http://projects.scipy.org/pipermail/numpy-discussion/attachments/20061005/47422051/attachment.bin 

More information about the Numpy-discussion mailing list