demo_prank_unique.f90 Source File

     program demo_prank_unique
     ! ranks array, removing duplicates
     use M_orderpack, only : prank_unique
     implicit none
     character(len=*),parameter :: g='(*(g0,1x))'
     integer,allocatable :: INVALS(:)
     integer,allocatable :: irngt(:)
     integer :: nord
     !
     write(*,g)'If enough values are unique, will return NORD indices'
     INVALS=[10,5,7,1,4,5,6,8,9,10,1]
     nord=5
     call printme()
     !
     write(*,g)'If not enough values are unique, will change NORD'
     INVALS=[-1,0,-1,0,-1,0,-1]
     nord=5
     if(allocated(irngt))deallocate(irngt)
     allocate(irngt(nord))
     call printme()
     !
     contains
     subroutine printme()
        write(*,g)'ORIGINAL:',INVALS
        write(*,g)'NUMBER OF INDICES TO SORT:',nord
        if(allocated(irngt))deallocate(irngt)
        allocate(irngt(nord))
        call prank_unique(INVALS,irngt,nord)
        write(*,g)'NUMBER OF INDICES RETURNED:',nord
        write(*,g)'RETURNED INDICES:',irngt(:nord)
        write(*,g)nord,'SMALLEST UNIQUE VALUES:',INVALS(irngt(:nord))
     end subroutine
     end program demo_prank_unique