nag_sort_intvec_search (m01nb) is based on Professor Niklaus Wirth's implementation of the Binary Search algorithm (see
Wirth (2004)), but with two modifications. First, if the sought-after item is less than the value of the first element of the array to be searched,
is returned. Second, if a value equal to the sought-after item is not found, the index of the immediate lower value is returned.
Not applicable.
The argument
valid should be used with caution. Set it to
false only if you are confident that the other arguments are correct, in particular that array
iv is in fact arranged in ascending order. If you wish to search the same array
iv many times, you are recommended to set
valid to
true on first call of
nag_sort_intvec_search (m01nb) and to
false on subsequent calls, in order to minimize the amount of time spent checking
iv, which may be significant if
iv is large.
This example reads a list of integer numbers and sought-after items and performs the search for these items.
function m01nb_example
fprintf('m01nb example results\n\n');
iv = [int64( 5) 6 11 12 13 13 21 23 ...
23 41 58 59 65 65 86 99];
item = [int64(21) 4 71 100];
m1 = [int64( 1) 1 1 5];
m2 = [int64(10) 16 16 16];
validate = true;
for j = 1:4
[result, ifail] = m01nb( ...
validate, iv, item(j), 'm1', m1(j), 'm2', m2(j));
validate = false;
fprintf('Search for %5d in index range [%2d:%2d]: index = %2d\n', ...
item(j), m1(j), m2(j), result);
end