Размещение площадок по отношению к площадкам сетки
pointer = sorted(meshsites,sites)
Различные команды в этой панели инструментов должны определять индекс j, для которого данный x лежит в интервале [tj.. tj + 1], с (ti) данной последовательностью без возрастания, например последовательностью узлов. Эта работа выполняется sorted следующим образом.
pointer = sorted(meshsites,sites) - целочисленный вектор строки, j-ая запись равняется количеству записей в meshsites это ≤ ssites(j), с ssites вектор sort(sites). Таким образом, если оба meshsites и sites не повторяются, то
meshsites(pointer(j)) ≤ sites(j) < meshsites(pointer(j)+1)
с очевидными интерпретациями, когда
pointer(j) < 1 or length(meshsites) < pointer(j) + 1
В частности, имея pointer(j) < 1 тогда соответствует sites(j) строго слева от meshsites(1), при этом имея length(meshsites) < pointer(j)+1 тогда соответствует sites(j) на, или справа от, meshsites(end).
Заявление
sorted([1 1 1 2 2 3 3 3],[0:4])
будет генерировать выходные данные 0 3 5 8 8, как и заявление
sorted([3 2 1 1 3 2 3 1],[2 3 0 4 1])
Выходные данные индексирования из sort([meshsites(:).',sites(:).']) используется.