13 Sorting

• PDA_DSORT (SLATEC/CAMSUN)
  Sort an array and optionally make the same interchanges in an auxiliary array.
• PDA_IPERM
  Forms the inverse of a permutation.
• PDA_QSAD
  Sort a DOUBLE PRECISION array into ascending order.
• PDA_QSAI
  Sort an INTEGER array into ascending order.
• PDA_QSAR
  Sort a REAL array into ascending order.
• PDA_QSDD
  Sort a DOUBLE PRECISION array into descending order.
• PDA_QSDI
  Sort a INTEGER array into descending order.
• PDA_QSDR
  Sort a REAL array into descending order.
• PDA_QSIAD
  Sort an array of pointers to access a DOUBLE PRECISION array in ascending order.
• PDA_QSIAI
  Sort an array of pointers to access a INTEGER array in ascending order.
• PDA_QSIAR
  Sort an array of pointers to access a REAL array in ascending order.
• PDA_QSIDD
  Sort an array of pointers to access a DOUBLE PRECISION array in descending order.
• PDA_QSIDI
  Sort an array of pointers to access a INTEGER array in descending order.
• PDA_QSIDR
  Sort an array of pointers to access a REAL array in descending order.
• PDA_RINP
  Reorders an array in place using a permutation index.
• PDA_SAAC
  Sorts the columns of a two dimensional array into ascending order.
• PDA_SAAR
  Sorts the rows of a two dimensional array into ascending order.

13.1 Replacing calls to M01DJF

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.