On Tue, Jun 1, 2010 at 6:56 AM, Keith Goodman <kwgood...@gmail.com> wrote: > On Tue, Jun 1, 2010 at 6:47 AM, Neal Becker <ndbeck...@gmail.com> wrote: >> Not sure what to call this. >> >> Any suggestion on computing the vector: >> >> sum(u[i*M:i*M+N]) for i in range (len(u)/M) > > How about a cumsum and then a loop to take the differences of the > desired indices of the cumsum? Might be faster to convert the cumsum > to a list (looks like you are interested in 1d?) for use in the loop. > Converting the loop to cython would be the next step.
c = cumsum c.take(idx_left, axis=0) - c.take(idx_right, axis=0) _______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion