Generate and return the rank
th smallest value in the input array by first sorting its elements in ascending order (optionally only between the specified indices [lb, ub]
).
By default, the array will be sorted in ascending order, unless the user-specific isSorted()
external function is supplied to define a custom sorting criterion.
- Parameters
-
[in] | array | : The input contiguous array of shape (:) of either
-
type css_type (scalar string container of default kind SK) or
-
type css_pdt (scalar string container of kind any supported by the processor (e.g., SK, SKA, SKD , or SKU)) or
-
type
character of kind any supported by the processor (e.g., SK, SKA, SKD , or SKU), or
-
type
logical of kind any supported by the processor (e.g., LK), or
-
type
integer of kind any supported by the processor (e.g., IK, IK8, IK16, IK32, or IK64), or
-
type
complex of kind any supported by the processor (e.g., CK, CK32, CK64, or CK128), or
-
type
real of kind any supported by the processor (e.g., RK, RK32, RK64, or RK128),
or,
-
a scalar assumed-length
character of kind any supported by the processor (e.g., SK, SKA, SKD , or SKU),
containing the potentially unsorted sequence of values that is to be sorted and whose rank th smallest value is to be reported as the output selection . |
[in] | rank | : The input integer of default kind IK, representing the index of the rank th smallest value in the input array or the rank th ordered value according to the user-supplied external logical function isSorted() . |
| isSorted | : The external user-specified function that takes two input scalar arguments of the same type and kind as the input array .
It returns a scalar logical of default kind LK that is .true. if the first input scalar argument is sorted with respect to the second input argument according to the user-defined condition within isSorted , otherwise, it is .false. .
If array is a Fortran string (i.e., an assumed-length scalar character ), then both input arguments to isSorted() are character(1,SK) of default kind SK.
The following illustrates the generic interface of isSorted() , function isSorted(a,b) result (sorted)
TYPE(KIND) , intent(in) :: a, b
logical(LK) :: sorted
end function
This module defines the relevant Fortran kind type-parameters frequently used in the ParaMonte librar...
integer, parameter LK The default logical kind in the ParaMonte library: kind(.true.) in Fortran, kind(....
integer, parameter IK The default integer kind in the ParaMonte library: int32 in Fortran, c_int32_t in C-Fortran Interoper...
where TYPE(KIND) is the same as the type and kind of the input argument array , which can be one of the following.
character(*, SK), intent(in) :: a, b
character(1, SK), intent(in) :: a, b
type(css_type) , intent(in) :: a, b
type(css_pdt) , intent(in) :: a, b
integer(IK) , intent(in) :: a, b
real(RK) , intent(in) :: a, b
This module contains the derived types for generating allocatable containers of scalar,...
This is the css_pdt parameterized type for generating instances of container of scalar of string obje...
This is the css_type type for generating instances of container of scalar of string objects.
where the specified kind type parameters (SK , IK , LK , CK , RK ) can refer to any of the supported kinds by the processor.
This user-defined equivalence check is extremely useful where a user-defined sorting criterion other than simple ascending order is needed, for example, when the case-sensitivity of an input string or array of strings is irrelevant or when sorting of the absolute values matters excluding the signs of the numbers, or when descending order is desired.
In such cases, user can define a custom sorting condition within the user-defined external function isSorted to achieve the goal.
(optional, the default sorting condition is ascending order, that is a < b .) |
[in] | lb | : The optional input scalar integer of default kind IK representing the lower bound of the segment of array below which the array elements are assumed to be in ascending order, such that only the elements starting with and beyond lb require sorting.
If the initial segment of the input array is ordered, specifying this lower bound will lead to better performance of the algorithm.
(optional, default = 1_IK ) |
[in] | ub | : The optional input scalar integer of default kind IK representing the upper bound of the segment of array above which the array elements are assumed to be in ascending order, such that only the elements ending with and below ub require sorting.
If the final segment of the input array is ordered, specifying this upper bound will lead to better performance of the algorithm.
(optional, default = size(array, kind = IK) ) |
- Returns
selection
: The output scalar of the same type and kind as array
containing the rank
th smallest (or ordered) value in the input array
.
Possible calling interfaces ⛓
selection
= getSelected(array, rank, isSorted, lb
= lb, ub
= ub)
Generate and return the rankth smallest value in the input array by first sorting its elements in asc...
This module contains procedures and generic interfaces for selecting the th smallest element in unsor...
- Warning
- If
lb
is present as an input argument, its value must be larger than 0
and less than ub
or its default value.
Furthermore, the specified rank
must be larger than or equal to the specified value for lb
.
If ub
is present as an input argument, its value must be less than or equal to the length of array
and larger than ub
or its default value.
Furthermore, the specified rank
must be smaller than or equal to the specified value for ub
.
These conditions are verified only if the library is built with the preprocessor macro CHECK_ENABLED=1
.
-
The functions under this generic interface are
impure
when the external
input argument isSorted()
is present
.
-
The
pure
procedure(s) documented herein become impure
when the ParaMonte library is compiled with preprocessor macro CHECK_ENABLED=1
.
By default, these procedures are pure
in release
build and impure
in debug
and testing
builds.
- See also
- setSelected
setRankDense
setRankOrdinal
setRankFractional
setRankStandard
setRankModified
setSorted
setSorted
getRemoved
setReversed
getReplaced
setReplaced
setSplit
Example usage ⛓
10 type(display_type) :: disp
14 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
15 call disp%show(
"!Select the `rank`th smallest element in the input `integer` Array.")
16 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
21 integer(IK) :: rank, lenArray
22 integer,
parameter :: TKG
= merge(IKG, IK, IKG
> 0)
23 integer(TKG),
allocatable :: array(:)
24 integer(TKG) :: selection
26 call disp%show(
"lenArray = getUnifRand(5_IK, 9_IK)")
28 call disp%show(
"rank = getUnifRand(1_IK, lenArray)")
32 call disp%show(
"array = getUnifRand(-9, +9, lenArray)")
36 call disp%show(
"selection = getSelected(array, rank)")
47 integer(IK) :: rank, lenArray
48 integer,
parameter :: TKG
= merge(IKG, IK, IKG
> 0)
49 integer(TKG),
allocatable :: array(:)
50 integer(TKG) :: selection
52 call disp%show(
"lenArray = getUnifRand(5_IK, 9_IK)")
54 call disp%show(
"rank = getUnifRand(1_IK, lenArray)")
58 call disp%show(
"array = getUnifRand(-9, +9, lenArray)")
62 call disp%show(
"selection = getSelected(array, rank)")
73 integer(IK) :: rank, lenArray
74 integer,
parameter :: TKG
= merge(IKG, IK, IKG
> 0)
75 integer(TKG),
allocatable :: array(:)
76 integer(TKG) :: selection
78 call disp%show(
"lenArray = getUnifRand(5_IK, 9_IK)")
80 call disp%show(
"rank = getUnifRand(1_IK, lenArray)")
84 call disp%show(
"array = getUnifRand(-9, +9, lenArray)")
88 call disp%show(
"selection = getSelected(array, rank)")
99 integer(IK) :: rank, lenArray
100 integer,
parameter :: TKG
= merge(IKG, IK, IKG
> 0)
101 integer(TKG),
allocatable :: array(:)
102 integer(TKG) :: selection
104 call disp%show(
"lenArray = getUnifRand(5_IK, 9_IK)")
106 call disp%show(
"rank = getUnifRand(1_IK, lenArray)")
110 call disp%show(
"array = getUnifRand(-9, +9, lenArray)")
114 call disp%show(
"selection = getSelected(array, rank)")
124 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
125 call disp%show(
"!Select the `rank`th smallest element in the input `real` Array.")
126 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
131 integer(IK) :: rank, lenArray
132 integer,
parameter :: TKG
= merge(RKG, RK, RKG
> 0)
133 real(TKG),
allocatable :: array(:)
134 real(TKG) :: selection
136 call disp%show(
"lenArray = getUnifRand(5_IK, 9_IK)")
138 call disp%show(
"rank = getUnifRand(1_IK, lenArray)")
142 call disp%show(
"array = getUnifRand(-9, +9, lenArray)")
146 call disp%show(
"selection = getSelected(array, rank)")
157 integer(IK) :: rank, lenArray
158 integer,
parameter :: TKG
= merge(RKG, RK, RKG
> 0)
159 real(TKG),
allocatable :: array(:)
160 real(TKG) :: selection
162 call disp%show(
"lenArray = getUnifRand(5_IK, 9_IK)")
164 call disp%show(
"rank = getUnifRand(1_IK, lenArray)")
168 call disp%show(
"array = getUnifRand(-9, +9, lenArray)")
172 call disp%show(
"selection = getSelected(array, rank)")
183 integer(IK) :: rank, lenArray
184 integer,
parameter :: TKG
= merge(RKG, RK, RKG
> 0)
185 real(TKG),
allocatable :: array(:)
186 real(TKG) :: selection
188 call disp%show(
"lenArray = getUnifRand(5_IK, 9_IK)")
190 call disp%show(
"rank = getUnifRand(1_IK, lenArray)")
194 call disp%show(
"array = getUnifRand(-9, +9, lenArray)")
198 call disp%show(
"selection = getSelected(array, rank)")
208 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
209 call disp%show(
"!Select the `rank`th smallest element in the input Fortran `string` (character) array.")
210 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
215 character(
2, TKG) :: selection
216 character(
2, TKG),
allocatable :: array(:)
217 integer(IK) :: rank, lenArray
219 call disp%show(
"lenArray = getUnifRand(5_IK, 9_IK)")
221 call disp%show(
"rank = getUnifRand(1_IK, lenArray)")
225 call disp%show(
"array = getUnifRand('AA', 'ZZ', lenArray)")
228 call disp%show( array , deliml
= SK_
"""" )
229 call disp%show(
"selection = getSelected(array, rank)")
234 call disp%show( array , deliml
= SK_
"""" )
239 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
240 call disp%show(
"!Select the `rank`th smallest element in the input Fortran `string` (character) scalar.")
241 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
246 character(
1, TKG) :: selection
247 character(:, TKG),
allocatable :: array
248 integer(IK) :: rank, lenArray
250 call disp%show(
"lenArray = getUnifRand(5_IK, 9_IK)")
252 call disp%show(
"rank = getUnifRand(1_IK, lenArray)")
256 call disp%show(
"array = getUnifRand(repeat('A', lenArray), repeat('Z', lenArray))")
257 array
= getUnifRand(
repeat(
'A', lenArray),
repeat(
'Z', lenArray))
259 call disp%show( array , deliml
= SK_
"""" )
260 call disp%show(
"selection = getSelected(array, rank)")
265 call disp%show( array , deliml
= SK_
"""" )
270 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
271 call disp%show(
"!Select the `rank`th smallest element in the input array of containers of Fortran `string` (character) scalars.")
272 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
277 integer(IK) :: rank, lenArray, iell
278 type(css),
allocatable :: array(:)
279 type(css) :: selection
281 call disp%show(
"lenArray = getUnifRand(5_IK, 9_IK)")
283 call disp%show(
"rank = getUnifRand(1_IK, lenArray)")
287 call disp%show(
"array = [(css(getUnifRand('A', 'Z', getUnifRand(1_IK, 9_IK))), iell = 1, lenArray)]")
290 call disp%show( array , deliml
= SK_
"""" )
291 call disp%show(
"selection = getSelected(array, rank)")
296 call disp%show( array , deliml
= SK_
"""" )
301 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
302 call disp%show(
"!Select according to an input user-defined comparison function.")
303 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
308 integer(IK) :: rank, lenArray, i
309 integer(TKG),
allocatable :: array(:)
310 integer(TKG) :: selection
312 call disp%show(
"lenArray = getUnifRand(5_IK, 9_IK)")
314 call disp%show(
"rank = getUnifRand(1_IK, lenArray)")
318 call disp%show(
"array = int([((-2)**i, i = 1, lenArray)], kind = TKG)")
319 array
= int([((
-2)
**i, i
= 1, lenArray)],
kind = TKG)
322 call disp%show(
"!Select the `rank`th largest element via an input custom-designed `isSorted()` function.")
323 call disp%show(
"selection = getSelected(array, rank, isSorted_IK)")
334 integer(IK) :: rank, lenArray, i
335 real(TKG),
allocatable :: array(:)
336 real(TKG) :: selection
338 call disp%show(
"lenArray = getUnifRand(5_IK, 9_IK)")
340 call disp%show(
"rank = getUnifRand(1_IK, lenArray)")
344 call disp%show(
"array = int([((-2)**i, i = 1, lenArray)], kind = TKG)")
345 array
= int([((
-2)
**i, i
= 1, lenArray)],
kind = TKG)
348 call disp%show(
"!Select the `rank`th smallest element solely based on the magnitude of numbers using a custom comparison function.")
349 call disp%show(
"selection = getSelected(array, rank, isSorted_RK)")
360 integer(IK) :: rank, lenArray, i
361 character(:,TKG),
allocatable :: array
362 character(
1,TKG) :: selection
364 call disp%show(
"lenArray = getUnifRand(5_IK, 9_IK)")
366 call disp%show(
"rank = getUnifRand(1_IK, lenArray)")
370 call disp%show(
"array = 'ParaMonte'")
373 call disp%show( array , deliml
= SK_
"""" )
374 call disp%show(
"!Select the `rank`th smallest element with case-sensitivity (default behavior).")
375 call disp%show(
"selection = getSelected(array, rank)")
378 call disp%show( selection , deliml
= SK_
"""" )
380 call disp%show( array , deliml
= SK_
"""" )
386 integer(IK) :: rank, lenArray, i
387 character(:,TKG),
allocatable :: array
388 character(
1,TKG) :: selection
390 call disp%show(
"lenArray = getUnifRand(5_IK, 9_IK)")
392 call disp%show(
"rank = getUnifRand(1_IK, lenArray)")
396 call disp%show(
"array = 'ParaMonte'")
399 call disp%show( array , deliml
= SK_
"""" )
400 call disp%show(
"!Select the `rank`th smallest element WITHOUT case-sensitivity via a custom-designed input comparison function.")
401 call disp%show(
"selection = getSelected(array, rank, isSorted_SK)")
404 call disp%show( selection , deliml
= SK_
"""" )
406 call disp%show( array , deliml
= SK_
"""" )
411 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
412 call disp%show(
"!Expedite the selection process for a partially sorted array via optional arguments `lb` or `ub`.")
413 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
417 integer(IK),
allocatable :: array(:)
418 integer(IK) :: selection
423 call disp%show(
"array = [-3, -1, -2, 1, 2, 3] ! all elements after index `3` are sorted (`ub = 3_IK`).")
424 array
= [
-3,
-1,
-2,
1,
2,
3]
425 call disp%show(
"selection = getSelected(array, rank, ub = 3_IK) ! all elements after index `3` are sorted (`ub = 3_IK`).")
436 function isSorted_IK(a,b)
result(isSorted)
438 integer(IK),
intent(in) :: a, b
439 logical(LK) :: isSorted
443 function isSorted_RK(a,b)
result(isSorted)
445 integer(RK),
intent(in) :: a, b
446 logical(LK) :: isSorted
447 isSorted
= abs(a)
< abs(b)
450 function isSorted_SK(a,b)
result(isSorted)
453 character(
1, SK),
intent(in) :: a, b
454 logical(LK) :: isSorted
Generate and return a scalar or a contiguous array of rank 1 of length s1 of randomly uniformly distr...
This is a generic method of the derived type display_type with pass attribute.
This is a generic method of the derived type display_type with pass attribute.
Generate and return the input string where the uppercase English alphabets are all converted to lower...
This module contains classes and procedures for computing various statistical quantities related to t...
This module contains classes and procedures for input/output (IO) or generic display operations on st...
type(display_type) disp
This is a scalar module variable an object of type display_type for general display.
integer, parameter RK
The default real kind in the ParaMonte library: real64 in Fortran, c_double in C-Fortran Interoperati...
integer, parameter IKS
The single-precision integer kind in Fortran mode. On most platforms, this is a 32-bit integer kind.
integer, parameter IKL
The scalar integer constant of intrinsic default kind, representing the lowest range integer kind typ...
integer, parameter IKH
The scalar integer constant of intrinsic default kind, representing the highest range integer kind ty...
integer, parameter RKD
The double precision real kind in Fortran mode. On most platforms, this is an 64-bit real kind.
integer, parameter IK8
The integer kind for an 8-bits container.
integer, parameter IKD
The double precision integer kind in Fortran mode. On most platforms, this is a 64-bit integer kind.
integer, parameter SK
The default character kind in the ParaMonte library: kind("a") in Fortran, c_char in C-Fortran Intero...
integer, parameter RKH
The scalar integer constant of intrinsic default kind, representing the highest-precision real kind t...
integer, parameter RKS
The single-precision real kind in Fortran mode. On most platforms, this is an 32-bit real kind.
This module contains the uncommon and hardly representable ASCII characters as well as procedures for...
Generate and return an object of type display_type.
Example Unix compile command via Intel ifort
compiler ⛓
3ifort -fpp -standard-semantics -O3 -Wl,-rpath,../../../lib -I../../../inc main.F90 ../../../lib/libparamonte* -o main.exe
Example Windows Batch compile command via Intel ifort
compiler ⛓
2set PATH=..\..\..\lib;%PATH%
3ifort /fpp /standard-semantics /O3 /I:..\..\..\include main.F90 ..\..\..\lib\libparamonte*.lib /exe:main.exe
Example Unix / MinGW compile command via GNU gfortran
compiler ⛓
3gfortran -cpp -ffree-line-length-none -O3 -Wl,-rpath,../../../lib -I../../../inc main.F90 ../../../lib/libparamonte* -o main.exe
Example output ⛓
27+3,
-5,
+3,
-5,
-7,
-8,
-4,
+7,
+0
32+3,
-5,
+3,
-5,
-7,
-8,
-4,
+7,
+0
41-1,
+1,
+4,
+3,
+0,
+7,
+4,
+5,
+6
46-1,
+1,
+4,
+3,
+0,
+7,
+4,
+5,
+6
74-8.00000000,
+5.00000000,
-9.00000000,
+0.00000000,
+6.00000000,
+5.00000000
79-8.00000000,
+5.00000000,
-9.00000000,
+0.00000000,
+6.00000000,
+5.00000000
88-5.0000000000000000,
+2.0000000000000000,
-4.0000000000000000,
+4.0000000000000000,
+3.0000000000000000,
-8.0000000000000000,
+6.0000000000000000,
+1.0000000000000000,
-1.0000000000000000
93-5.0000000000000000,
+2.0000000000000000,
-4.0000000000000000,
+4.0000000000000000,
+3.0000000000000000,
-8.0000000000000000,
+6.0000000000000000,
+1.0000000000000000,
-1.0000000000000000
102-9.00000000000000000000000000000000000,
+8.00000000000000000000000000000000000,
+8.00000000000000000000000000000000000,
-8.00000000000000000000000000000000000,
-2.00000000000000000000000000000000000,
+6.00000000000000000000000000000000000,
+9.00000000000000000000000000000000000
105-9.00000000000000000000000000000000000
107-9.00000000000000000000000000000000000,
+8.00000000000000000000000000000000000,
+8.00000000000000000000000000000000000,
-8.00000000000000000000000000000000000,
-2.00000000000000000000000000000000000,
+6.00000000000000000000000000000000000,
+9.00000000000000000000000000000000000
121"CU",
"SU",
"ZL",
"NJ",
"BT",
"SF",
"NP",
"AK",
"GA"
126"CU",
"SU",
"ZL",
"NJ",
"BT",
"SF",
"NP",
"AK",
"GA"
138array
= getUnifRand(
repeat(
'A', lenArray),
repeat(
'Z', lenArray))
159"L",
"N",
"H",
"F",
"K",
"V",
"C",
"Y",
"V",
"Y",
"H",
"H",
"H",
"G",
"S",
"V",
"Y",
"P",
"I",
"B",
"M",
"D",
"S",
"Y",
"X",
"A",
"G",
"F",
"V"
164"L",
"N",
"H",
"F",
"K",
"V",
"C",
"Y",
"V",
"Y",
"H",
"H",
"H",
"G",
"S",
"V",
"Y",
"P",
"I",
"B",
"M",
"D",
"S",
"Y",
"X",
"A",
"G",
"F",
"V"
176array
= int([((
-2)
**i, i
= 1, lenArray)],
kind = TKG)
178-2,
+4,
-8,
+16,
-32,
+64
184-2,
+4,
-8,
+16,
-32,
+64
191array
= int([((
-2)
**i, i
= 1, lenArray)],
kind = TKG)
193-2.0000000000000000,
+4.0000000000000000,
-8.0000000000000000,
+16.000000000000000,
-32.000000000000000
199-2.0000000000000000,
+4.0000000000000000,
-8.0000000000000000,
+16.000000000000000,
-32.000000000000000
238array
= [
-3,
-1,
-2,
1,
2,
3]
243-3,
-1,
-2,
+1,
+2,
+3
- Test:
- test_pm_arraySelect
- Todo:
- Low Priority: This generic interface can be extended to higher-rank input objects.
Final Remarks ⛓
If you believe this algorithm or its documentation can be improved, we appreciate your contribution and help to edit this page's documentation and source file on GitHub.
For details on the naming abbreviations, see this page.
For details on the naming conventions, see this page.
This software is distributed under the MIT license with additional terms outlined below.
-
If you use any parts or concepts from this library to any extent, please acknowledge the usage by citing the relevant publications of the ParaMonte library.
-
If you regenerate any parts/ideas from this library in a programming environment other than those currently supported by this ParaMonte library (i.e., other than C, C++, Fortran, MATLAB, Python, R), please also ask the end users to cite this original ParaMonte library.
This software is available to the public under a highly permissive license.
Help us justify its continued development and maintenance by acknowledging its benefit to society, distributing it, and contributing to it.
- Copyright
- Computational Data Science Lab
- Author:
- Amir Shahmoradi, September 1, 2017, 12:00 AM, Institute for Computational Engineering and Sciences (ICES), The University of Texas Austin
Definition at line 546 of file pm_arraySelect.F90.