Текст подпрограммы и версий ( Фортран ) pa07r.zip pa07d.zip |
Тексты тестовых примеров ( Фортран ) tpa07r.zip tpa07d.zip |
Текст подпрограммы и версий ( Си ) pa07r_c.zip pa07d_c.zip |
Тексты тестовых примеров ( Си ) tpa07r_c.zip tpa07d_c.zip |
Текст подпрограммы и версий ( Паскаль ) pa07r_p.zip pa07e_p.zip |
Тексты тестовых примеров ( Паскаль ) tpa07r_p.zip tpa07e_p.zip |
Определение узлов упорядоченной по возрастанию или по убыванию неравномерной сетки, между которыми лежит заданная точка.
Пусть заданы узлы сетки x1, x2,..., xN, упорядоченные либо по возрастанию, либо по убыванию. Подпрограмма PA07R по заданной точке X вычисляет методом бисекции значение индекса J, при котором точка X лежит между узлами xJ и xJ + 1.
Если J = 0, то это означает, что точка X лежит левее узла x1. Если J = N, то это означает, что точка X лежит правее узла xN.
SUBROUTINE PA07R (XX, N, X, J)
Параметры
XX - | вещественный вектор длины N, содержащий узлы заданной сетки x1, x2,..., xN, упорядоченной либо по возрастанию, либо по убыванию; |
N - | количество узлов сетки (тип: целый); |
X - | заданная точка, местоположение которой ищется по отношению к узлам сетки (тип: вещественный); |
J - | целая переменная, значение которой полагается равной индексу, при котором точка X лежит между узлами XX (J) и XX (J + 1) . |
Версии
PA07D - | определение узлов упорядоченной по возрастанию или убыванию неравномерной сетки, между которыми лежит заданная точка. При этом параметры XX и X должны иметь тип DOUBLE PRECISION. |
Вызываемые подпрограммы: нет
Замечания по использованию: нет
DIMENSION XX(20) N = 20 R = 0.0 DO 1 I = 1, N XX(I) = R 1 R = R + 0.2 X = 1.7 CALL PA07R (XX, N, X, J) Результат: J = 9