select Interface

public interface select

Use an in-place quick select on an array of numbers

Example usage

program selectTest
use variableKind
use m_strings, only: str
use m_allocate, only: allocate
use m_random, only: rngInteger,rngNormal
use m_select, only: select
real(r64),allocatable :: d1D(:)
integer(i32),allocatable :: i1D(:)
integer(i32) :: k
real(r64) :: dv
integer(i32) :: iv

write(*,'(a)') 'Select the kth smallest element from a 10000 length array of random double precision numbers'
call allocate(d1D, 10000)
call rngNormal(d1D)
k = (1+size(d1D))/2
dv = select(d1D,k)
write(*,'(a)') 'kth element? '//str(dv)
write(*,'(a)') 'Select the kth smallest element from a 10000 length array of random integers''
call allocate(i1D, 10000)
call rngInteger(i1D)
iv = select(i1D, k)
write(*,'(a)') 'kth element? '//str(iv)
end program


Subroutines

public subroutine quickSelect_i1D(this, k, res)

Interfaced with select()

Arguments

Type IntentOptional AttributesName
integer(kind=i32), intent(inout) :: this(:)

Array to choose kth smallest from

integer(kind=i32), intent(in) :: k

kth smallest element

integer(kind=i32) :: res

public subroutine quickSelect_id1D(this, k, res)

Interfaced with select()

Arguments

Type IntentOptional AttributesName
integer(kind=i64), intent(inout) :: this(:)

Array to choose kth smallest from

integer(kind=i32), intent(in) :: k

kth smallest element

integer(kind=i64) :: res

public subroutine quickSelect_r1D(this, k, res)

Interfaced with select()

Arguments

Type IntentOptional AttributesName
real(kind=r32), intent(inout) :: this(:)

Array to choose kth smallest from

integer(kind=i32), intent(in) :: k

kth smallest element

real(kind=r32) :: res

public subroutine quickSelect_d1D(this, k, res)

Interfaced with select()

Arguments

Type IntentOptional AttributesName
real(kind=r64), intent(inout) :: this(:)

Array to choose kth smallest from

integer(kind=i32), intent(in) :: k

kth smallest element

real(kind=r64) :: res

allocate allocate_test appendString arange argInsertionsort argMedianOf3 argPartition argSelect argSort array1D_test binarySearch checkIsOpen closeFile compact compare compare_d1 copy copy_test countEntries crossproduct cumprod cumsum dArgDynamicArray dArgDynamicArray_test daysInMonth daysInYear dDynamicArray dDynamicArray_test deallocate debug debug_s deleteFile diff eMsg fastTwoDiff fastTwoSum fErr fileExists fileIO_test geometricMean getExtension getFileSize getNFileLines getRandomSeed hasExtension hasNentries iachar1D iArgDynamicArray iArgDynamicArray_test idArgDynamicArray idArgDynamicArray_test idDynamicArray idDynamicArray_test iDynamicArray iDynamicArray_test ind2sub indexing_test insertionsort integerBin intervalSearch isLeapYear isOpen isSorted isString KdTree KdTree_test lowerCase maths_test mean median medianOf3 mErr msg norm1 normI openBinaryFile openFile partition partition3way partition3way_d1D prependString Prng Prng_test ProgressBar project random_test rArgDynamicArray rArgDynamicArray_test rDynamicArray rDynamicArray_test read1Dble read1Integer readline readline reallocate reallocate_test removeComments repeat replacedelim rngExponential rngExponential_d1 rngExponential_d1D rngExponential_d2D rngExponential_d3D rngExponential_unscaled_d1 rngInteger rngInteger_i1 rngInteger_i1D rngInteger_i2D rngInteger_i3D rngNormal rngNormal_d1 rngNormal_d1D rngNormal_d2D rngNormal_d3D rngUniform rngUniform_d1 rngUniform_d1D rngUniform_d2D rngUniform_d3D rngUniform_xorshift rngWeibull rngWeibull_d1 rngWeibull_d1D rngWeibull_d2D rngWeibull_d3D secondsToHMS select select_test setPrng setPrng_withSeed setPrng_WOseed shuffle simpleSearch skipFileLines sort sorting_test std str strings_test sub2ind swap test tester time_test timeInSeconds timeToInteger trimExtension trimmedmean twoDiff twoSum upperCase variance wMsg writeline