subroutine test_median
use M_orderpack__median, only : median
implicit none
real,allocatable :: xdont(:)
real(kind=dp),allocatable :: ddont(:)
integer,allocatable :: idont(:)
call unit_check_start('median', '-library orderpack') ! start tests
xdont=[80.0,70.0,20.0,10.0,1000.0]
call unit_check('median', median(xdont).eq.70.0, 'real median',median(xdont),70.0)
!
idont=[11, 22, 33, 44, 55, 66, 77, 88]
call unit_check('median', median(idont).eq.49, 'integer median',median(idont),49)
!
ddont=[11.0d0,77.0d0,22.0d0,66.0d0,33.0d0,88.0d0]
call unit_check('median', median(ddont).eq.49.5d0, 'doubleprecision median',median(ddont),49.5)
!
call unit_check_done('median',msg='test completed')
end subroutine test_median