Return NumPy values sorted along a selected axis together with the indices that produce the same ordering.
