move pivot storage to function
sort after wrapping kth around axis
fix NPY_SEPARATE_COMPILATION=0
support iterative partitioning for multidimensional arrays and add tests
documentation update