This module contains procedures and generic interfaces for finding the specific array index whose element has the largest value smaller than the input value in arrays of various types.
More...
This module contains procedures and generic interfaces for finding the specific array index whose element has the largest value smaller than the input value in arrays of various types.
The input search array must be in ascending order.
- Benchmarks:
Benchmark :: The runtime performance of getBin with and without isLess
external optional input argument. ⛓
4 use iso_fortran_env,
only:
error_unit
15 integer(IK) :: fileUnit
16 integer(IK) ,
parameter :: NSIZE
= 20_IK
17 integer(IK) ,
parameter :: NBENCH
= 2_IK
18 integer(IK) :: arraySize(NSIZE)
19 logical(LK) :: dummy
= .true._LK
20 real(RKG) ,
allocatable :: Array(:)
23 type(bench_type) :: bench(NBENCH)
24 integer(IK) :: index, bin
25 real(RK) :: MeanTime(
2)
= 0._RK
27 bench(
1)
= bench_type(name
= SK_
"default" , exec
= default , overhead
= setOverhead)
28 bench(
2)
= bench_type(name
= SK_
"isLessArg", exec
= isLessArg , overhead
= setOverhead)
30 arraySize
= [(
2_IK**(isize
+1), isize
= 1_IK, NSIZE )]
32 write(
*,
"(*(g0,:,' '))")
33 write(
*,
"(*(g0,:,' '))")
"isLessArg vs. default"
34 write(
*,
"(*(g0,:,' '))")
36 open(newunit
= fileUnit, file
= "main.out", status
= "replace")
38 write(fileUnit,
"(*(g0,:,', '))")
"arraySize", (bench(i)
%name, i
= 1, NBENCH)
40 loopOverArraySize:
do isize
= 1, NSIZE
42 write(
*,
"(*(g0,:,' '))")
"Benchmarking with size", arraySize(isize)
44 call random_number(Limit)
46 Array
= getLinSpace(Limit(
1), Limit(
2), count
= arraySize(isize))
49 bench(i)
%timing
= bench(i)
%getTiming(minsec
= 0.03_RK)
50 MeanTime(i)
= MeanTime(i)
+ bench(i)
%timing
%mean
54 bench(i)
%timing
= bench(i)
%getTiming(minsec
= 0.03_RK)
55 MeanTime(i)
= MeanTime(i)
+ bench(i)
%timing
%mean
58 write(fileUnit,
"(*(g0,:,', '))") arraySize(isize), MeanTime
/ 2
60 end do loopOverArraySize
61 write(
*,
"(*(g0,:,' '))") dummy
62 write(
*,
"(*(g0,:,' '))")
72 subroutine setOverhead()
77 subroutine initialize()
83 dummy
= dummy
.and. bin
== index
86 subroutine isLessArg()
89 bin
= getBin(Array, value, isLess
= isLess)
100 pure function isLess(value, segment)
result(less)
102 integer(IK) ,
intent(in) ::
value, segment
104 less
= value
< segment
Generate and return bin, the index of the element of the input ascending-ordered array,...
Generate count evenly spaced points over the interval [x1, x2] if x1 < x2, or [x2,...
Generate and return an object of type timing_type containing the benchmark timing information and sta...
Return a uniform random scalar or contiguous array of arbitrary rank of randomly uniformly distribute...
Generate and return an array of size two containing the minimum and maximum of the two input values i...
This module contains procedures and generic interfaces for finding the specific array index whose ele...
This module contains procedures and generic interfaces for generating arrays with linear or logarithm...
This module contains abstract interfaces and types that facilitate benchmarking of different procedur...
This module contains classes and procedures for computing various statistical quantities related to t...
This module defines the relevant Fortran kind type-parameters frequently used in the ParaMonte librar...
integer, parameter RK
The default real kind in the ParaMonte library: real64 in Fortran, c_double in C-Fortran Interoperati...
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...
integer, parameter SK
The default character kind in the ParaMonte library: kind("a") in Fortran, c_char in C-Fortran Intero...
integer, parameter RKS
The single-precision real kind in Fortran mode. On most platforms, this is an 32-bit real kind.
This module contains procedures and generic interfaces for finding the minimum and maximum of two inp...
This is the class for creating benchmark and performance-profiling objects.
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
Postprocessing of the benchmark output ⛓
3import matplotlib.pyplot
as plt
8dirname = os.path.basename(os.getcwd())
12df = pd.read_csv(
"main.out", delimiter =
", ")
13colnames = list(df.columns.values)
19ax = plt.figure(figsize = 1.25 * np.array([6.4,4.6]), dpi = 200)
22for colname
in colnames[1:]:
23 plt.plot( df[colnames[0]].values
28plt.xticks(fontsize = fontsize)
29plt.yticks(fontsize = fontsize)
30ax.set_xlabel(colnames[0], fontsize = fontsize)
31ax.set_ylabel(
"Runtime [ seconds ]", fontsize = fontsize)
32ax.set_title(
" vs. ".join(colnames[1:])+
"\nLower is better.", fontsize = fontsize)
36plt.grid(visible =
True, which =
"both", axis =
"both", color =
"0.85", linestyle =
"-")
37ax.tick_params(axis =
"y", which =
"minor")
38ax.tick_params(axis =
"x", which =
"minor")
39ax.legend ( colnames[1:]
46plt.savefig(
"benchmark." + dirname +
".runtime.png")
52ax = plt.figure(figsize = 1.25 * np.array([6.4,4.6]), dpi = 200)
55plt.plot( df[colnames[0]].values
56 , np.ones(len(df[colnames[0]].values))
61for colname
in colnames[2:]:
62 plt.plot( df[colnames[0]].values
63 , df[colname].values / df[colnames[1]].values
67plt.xticks(fontsize = fontsize)
68plt.yticks(fontsize = fontsize)
69ax.set_xlabel(colnames[0], fontsize = fontsize)
70ax.set_ylabel(
"Runtime compared to {}".format(colnames[1]), fontsize = fontsize)
71ax.set_title(
"Runtime Ratio Comparison. Lower means faster.\nLower than 1 means faster than {}().".format(colnames[1]), fontsize = fontsize)
75plt.grid(visible =
True, which =
"both", axis =
"both", color =
"0.85", linestyle =
"-")
76ax.tick_params(axis =
"y", which =
"minor")
77ax.tick_params(axis =
"x", which =
"minor")
78ax.legend ( colnames[1:]
85plt.savefig(
"benchmark." + dirname +
".runtime.ratio.png")
Visualization of the benchmark output ⛓
Benchmark moral ⛓
- The procedures under the generic interface getBin take an optional
isLess()
external-function input argument which can perform custom user-defined comparisons between the input value
and array
elements.
By default, the comparison is the <
operator.
- The presence of an external
isLess
argument without inlining does appear to deteriorate the runtime performance of getBin.
However, the exact amount of this penalty appears to significantly depend on the computing platform and the ability of the compiler to inline the external function.
- Test:
- test_pm_arraySearch
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:
- Fatemeh Bagheri, Wednesday 12:20 AM, October 13, 2021, Dallas, TX