demo_sort_quick_compact.f90 Source File

       program demo_sort_quick_compact
       use M_sort, only : sort_quick_compact
       implicit none
       integer,parameter            :: isz=10000
       real                         :: rrin(isz)
       real                         :: rrout(isz)
       integer                      :: i
       write(*,*)'initializing array with ',isz,' random numbers'
       CALL RANDOM_NUMBER(rrin)
       rrin=rrin*450000.0
       write(*,*)'sort real array with sort_quick_compact(3f)'
       rrout=sort_quick_compact(rrin)
       write(*,*)'checking '
       do i=1,isz-1
          if(rrout(i).lt.rrout(i+1))then
             write(*,*)'Error in sorting reals', &
             & i,rrout(i),rrout(i+1)
          endif
       enddo
       write(*,*)'test of sort_quick_compact(3f) complete'
       end program demo_sort_quick_compact