Optimizing list.sort() by checking type in advance

Thanks for looking at this! That's why I spent months of my life (overall) devising a sequence of sorting algorithms for Python that reduced the number of comparisons needed. Yes, that's why I think this is so cool: for a couple dozen lines of code, we can get (at least for some cases, according to my questionable benchmarks) the kinds of massive improvements you had to use actual computer science to achieve (as opposed to mere hackery). Note that when Python's current sort was adopted in Java, they still kept a quicksort variant for "unboxed" builtin types. The adaptive merge sort incurs many overheads that often cost more than they save unless comparisons are in fact very expensive compared to the cost of pointer copying (and in Java comparison of unboxed types is cheap). Indeed, for native numeric types, where comparison is dirt cheap, quicksort generally runs faster than mergesort despite that the former does _more_ comparisons (because mergesort does so much more pointer-copying). Ya, I think this may be a good approach for floats: if the list is all floats, just copy all the floats into a seperate array, use the standard library quicksort, and then construct a sorted PyObject* array. Like maybe set up a struct { PyObject* payload, float key } type of deal. This wouldn't work for strings (unicode is scary), and probably not for ints (one would have to check that all the ints are within C long bounds). Though on the other hand perhaps this would be too expensive? I had considered something "like this" for Python 2, but didn't pursue it because comparison was defined between virtually any two types (34 < [1], etc), and people were careless about that (both by design and by accident). In Python 3, comparison "blows up" for absurdly mixed types, so specializing for homogeneously-typed lists is a more promising idea on the face of it. The comparisons needed to determine _whether_ a list's objects have a common type is just len(list)-1 C-level pointer comparisons, and so goes fast. So I expect that, when it applies, this would speed even sorting an already-ordered list with at least 2 elements. That's what my crude benchmarks indicate... when I applied my sort to a list of 1e7 ints with a float tacked on the end, my sort actually ended up being a bit faster over several trials (which I attribute to PyObject_RichCompare == Py_True being faster than PyObject_RichCompareBool == 1, apologies for any typos in that code). For a mixed-type list with at least 2 elements, it will always be pure loss. But (a) I expect such lists are uncommon (and especially uncommon in Python 3); and (b) a one-time scan doing C-level pointer comparisons until finding a mismatched type is bound to be a relatively tiny cost compared to the expense of all the "rich comparisons" that follow. So +1 from me on pursuing this. Elliot, please: - Keep this on python-ideas. python-dev is for current issues in Python development, not for speculating about changes. - Open an issue on the tracker: https://bugs.python.org/ OK - At least browse the info for developers: https://docs.python.org/devguide/ Ya, I'm working on setting this up as a patch in the hg repo as opposed to an extension module to make benchmarking cleaner/more sane. - Don't overlook Lib/test/sortperf.py. As is, it should be a good test of what your approach so far _doesn't_ help, since it sorts only lists of floats (& I don't think you're special-casing them). If the timing results it reports aren't significantly hurt (and I expect they won't be), then add specialization for floats too and gloat about the speedup :-) Ya, I mean they aren't special-cased, but homogenous lists of floats still fit in the tp->rich_compare case, which still bypasses the expensive PyObject_RichCompare. I'll guess I'll see when I implement this as a patch and can run it on sortperf.py. - I expect tuples will also be worth specializing (complex sort keys are often implemented as tuples). I'm not sure what you mean here... I'm looking at the types of lo.keys, not of saved_ob_item (I think I said that earlier in this thread by mistake actually). So if someone is passing tuples and using itemgetter to extract ints or strings or whatever, the current code will work fine; lo.keys will be scalar types. Unless I misunderstand you here. I mean, when would lo.keys actually be tuples? Nice start! :-) Thanks!

On Tue, Oct 11, 2016 at 2:29 PM, Elliot Gorokhovsky <elliot.gorokhovsky@gmail.com> wrote:
Not quite sure what you mean here. What is payload, what is key? Are you implying that the original float objects could be destroyed and replaced with others of equal value? Python (unlike insurance claims) guarantees that you get back the exact same object as you started with. ChrisA

Oh no, the idea here is just you would copy over the floats associated with the PyObject* and keep them in an array of such structs, so that we know which PyObject* are associated with which floats. Then after the standard library quicksort sorts them you would copy the PyObject* into the list. So you sort the PyObject* keyed by the floats. Anyway, I think the copying back and forth would probably be too expensive, it's just an idea. Also, I apologize for the formatting of my last email, I didn't realize Inbox would mess up the quoting like that. I'll ensure I use plain-text quotes from now on. On Mon, Oct 10, 2016 at 9:38 PM Chris Angelico <rosuav@gmail.com> wrote:

On Tue, Oct 11, 2016 at 2:41 PM, Elliot Gorokhovsky <elliot.gorokhovsky@gmail.com> wrote:
It also wouldn't work if you have more than one object with the same value.
Python's sort is stable, so the three elements of the list (being all equal) must remain in the same order. ChrisA

It would still be stable. You would copy over {{x,1.0},{y,1.0},{x,1.0}}, and as long as a stable sort is used you would get out the same array, using the cmp function left->key < right->key. Then you would go in order, copying back [x,y,x]. On Mon, Oct 10, 2016 at 9:49 PM Chris Angelico <rosuav@gmail.com> wrote:

On Mon, Oct 10, 2016 at 11:30 PM Elliot Gorokhovsky < elliot.gorokhovsky@gmail.com> wrote:
If someone wanted to sort, e.g., a table (likely a list of tuples) by multiple columns at once, they might pass the key function as `itemgetter(3, 4, 5)`, meaning to sort by "column" (actually item) 3, then columns 4 and then 5 as tiebreakers. This itemgetter will return a new tuple of three items, that tuple being the key to sort by. Since tuples sort by the first different item, in this theoretical example the result of sort() will be exactly what the user wanted: a table sorted by three columns at once. A practical example of such a use case is sorting by last name first and then by first name where two people have the same last name. Assuming a list of dicts in this case, the key function passed to sort() would simply be `itemgetter('lastname", "firstname")`, which returns a tuple of two items to use as the key. So yes, there are perfectly valid use cases for tuples as keys.

On Mon, Oct 10, 2016 at 11:29 PM, Elliot Gorokhovsky <elliot.gorokhovsky@gmail.com> wrote:
I happened onto a page talking about float radix sort, and thought of this thread. Here it is: http://stereopsis.com/radix.html The author claimed an 8x speedup, though the test was done nearly fifteen years ago. I was unsure about posting publicly, because it's not as if an even faster float sort would help decide whether specialized sorts are worth adding to CPython. I'm posting for history.

On Tue, Oct 11, 2016 at 2:29 PM, Elliot Gorokhovsky <elliot.gorokhovsky@gmail.com> wrote:
Not quite sure what you mean here. What is payload, what is key? Are you implying that the original float objects could be destroyed and replaced with others of equal value? Python (unlike insurance claims) guarantees that you get back the exact same object as you started with. ChrisA

Oh no, the idea here is just you would copy over the floats associated with the PyObject* and keep them in an array of such structs, so that we know which PyObject* are associated with which floats. Then after the standard library quicksort sorts them you would copy the PyObject* into the list. So you sort the PyObject* keyed by the floats. Anyway, I think the copying back and forth would probably be too expensive, it's just an idea. Also, I apologize for the formatting of my last email, I didn't realize Inbox would mess up the quoting like that. I'll ensure I use plain-text quotes from now on. On Mon, Oct 10, 2016 at 9:38 PM Chris Angelico <rosuav@gmail.com> wrote:

On Tue, Oct 11, 2016 at 2:41 PM, Elliot Gorokhovsky <elliot.gorokhovsky@gmail.com> wrote:
It also wouldn't work if you have more than one object with the same value.
Python's sort is stable, so the three elements of the list (being all equal) must remain in the same order. ChrisA

It would still be stable. You would copy over {{x,1.0},{y,1.0},{x,1.0}}, and as long as a stable sort is used you would get out the same array, using the cmp function left->key < right->key. Then you would go in order, copying back [x,y,x]. On Mon, Oct 10, 2016 at 9:49 PM Chris Angelico <rosuav@gmail.com> wrote:

On Mon, Oct 10, 2016 at 11:30 PM Elliot Gorokhovsky < elliot.gorokhovsky@gmail.com> wrote:
If someone wanted to sort, e.g., a table (likely a list of tuples) by multiple columns at once, they might pass the key function as `itemgetter(3, 4, 5)`, meaning to sort by "column" (actually item) 3, then columns 4 and then 5 as tiebreakers. This itemgetter will return a new tuple of three items, that tuple being the key to sort by. Since tuples sort by the first different item, in this theoretical example the result of sort() will be exactly what the user wanted: a table sorted by three columns at once. A practical example of such a use case is sorting by last name first and then by first name where two people have the same last name. Assuming a list of dicts in this case, the key function passed to sort() would simply be `itemgetter('lastname", "firstname")`, which returns a tuple of two items to use as the key. So yes, there are perfectly valid use cases for tuples as keys.

On Mon, Oct 10, 2016 at 11:29 PM, Elliot Gorokhovsky <elliot.gorokhovsky@gmail.com> wrote:
I happened onto a page talking about float radix sort, and thought of this thread. Here it is: http://stereopsis.com/radix.html The author claimed an 8x speedup, though the test was done nearly fifteen years ago. I was unsure about posting publicly, because it's not as if an even faster float sort would help decide whether specialized sorts are worth adding to CPython. I'm posting for history.
participants (5)
-
Chris Angelico
-
Elliot Gorokhovsky
-
Franklin? Lee
-
Greg Ewing
-
Jonathan Goble