[Numpy-discussion] Openmp support (was numpy's future (1.1 and beyond): which direction(s) ?)

[Charles R Harris]

On Sat, Mar 22, 2008 at 12:01 PM, James Philbin <philbinj at gmail.com> wrote:

> >  gcc keeps advancing autovectorization.  Is manual vectorization worth
> the
> >  trouble?
>
> Well, the way that the ufuncs are written at the moment,
> -ftree-vectorize will never kick in due to the non-constant strides.
> To get this to work, one has to special case out unary strides. Even
> with constant strides -ftree-vectorize often produces sub-optimal code
> as it has to make very conservative assumptions about the content of
> variables (it can do better if -fstrict-aliasing is used, but I think


Numpy would die with strict-aliasing because it relies on casting char
pointers to various types. It might be possible to use void pointers, but
that would be a major do over and it isn't clear what the strides and
offsets, currently in char units, would become.

Chuck
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/numpy-discussion/attachments/20080322/f0b60473/attachment.html>