On Fri, Jun 11, 2010 at 6:56 AM, Sturla Molden <stu...@molden.no> wrote:
> > Also about array iterators in NumPy's C base (i.e. for doing something > along an axis): we don't need those. There is a different way of coding > which leads to faster code. > > 1. Collect an array of pointers to each subarray (e.g. using > std::vector<dtype*> or dtype**) > 2. Dispatch on the pointer array... Do you have the code for this ? That's something I wanted to do, but never took the time to do. Faster generic iterator would be nice, but very hard to do in general. > Another thing I did when reimplementing lfilter was "copy-in copy-out" > for strided arrays. What is copy-in copy out ? I am not familiar with this term ? David _______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion