The routines PDA_SAACD solves the same problem as M01DJF, except that the columns of the array are sorted using an index vector rather than ranks. The index vector can be used to re-order other arrays (via the PDA_RINPD routine) or it can be permuted into a rank using PDA_IPERM. You should also take care to make sure that necessary extra elements are available in M, IP and the new workspace LINK.
PDA [1ex