Classifique invariável para numpy.argsort com várias dimensões

numpy.argsort docs state

Devoluções:
index_array: ndarray, int Matriz de índices que ordenam a ao longo do eixo especificado. Se a é unidimensional,a[index_array] produz um a classificado.

Como posso aplicar o resultado denumpy.argsort para uma matriz multidimensional para recuperar uma matriz classificada? (NÃO é apenas uma matriz 1-D ou 2-D; pode ser uma matriz N-dimensional em que N é conhecido apenas em tempo de execução)

>>> import numpy as np
>>> np.random.seed(123)
>>> A = np.random.randn(3,2)
>>> A
array([[-1.0856306 ,  0.99734545],
       [ 0.2829785 , -1.50629471],
       [-0.57860025,  1.65143654]])
>>> i=np.argsort(A,axis=-1)
>>> A[i]
array([[[-1.0856306 ,  0.99734545],
        [ 0.2829785 , -1.50629471]],

       [[ 0.2829785 , -1.50629471],
        [-1.0856306 ,  0.99734545]],

       [[-1.0856306 ,  0.99734545],
        [ 0.2829785 , -1.50629471]]])

Para mim, não é apenas uma questão de usarsort() em vez de; Eu tenho outra matrizB e eu quero pedirB usando os resultados denp.argsort(A) ao longo do eixo apropriado. Considere o seguinte exemplo:

>>> A = np.array([[3,2,1],[4,0,6]])
>>> B = np.array([[3,1,4],[1,5,9]])
>>> i = np.argsort(A,axis=-1)
>>> BsortA = ???             
# should result in [[4,1,3],[5,1,9]]
# so that corresponding elements of B and sort(A) stay together

Parece que essa funcionalidade éjá é uma solicitação de aprimoramento numpy.