SUBROUTINE DLASRT( ID, N, D, INFO) CHARACTER * 1 ID INTEGER N, INFO DOUBLE PRECISION D(*) SUBROUTINE DLASRT_64( ID, N, D, INFO) CHARACTER * 1 ID INTEGER*8 N, INFO DOUBLE PRECISION D(*)
SUBROUTINE LASRT( ID, [N], D, [INFO]) CHARACTER(LEN=1) :: ID INTEGER :: N, INFO REAL(8), DIMENSION(:) :: D SUBROUTINE LASRT_64( ID, [N], D, [INFO]) CHARACTER(LEN=1) :: ID INTEGER(8) :: N, INFO REAL(8), DIMENSION(:) :: D
void dlasrt(char id, int n, double *d, int *info);
void dlasrt_64(char id, long n, double *d, long *info);
Use Quick Sort, reverting to Insertion sort on arrays of
size <= 20. Dimension of STACK limits N to about 2**32.