> indices = argsort(a1) > ranks = zeros_like(indices) > ranks[indices] = arange(len(indices))
Doesn't answer your original question directly, but I only recently learned from this list that the following does the same as the above: ranks = a1.argsort().argsort() Will wonders never cease... So does ranks=a2.argsort(axis=0).argsort(axis=0) then do the trick? Zach _______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion
