On 17/04/2008, Zachary Pincus wrote: > But even if indices = array, one still needs to do something like: > for index in indices: histogram[index] += 1 > > Which is slow in python and fast in C. >
I haven't tried this, but if you want the sum in C you could do for x in unique(indices): histogram[x] = (indices==x).sum() Of course, this just replaces an O(N log N) algorithm by an O(N * M) (M is the number of bins), which is only going to be worth it for very small M. Peter _______________________________________________ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion