Описание подпрограммы BCPOL (DBCPOL)

Подпро1рамма находит минимум функции «переменных ([Л6], ч.2, стр.150). Имеет вызов для расчётов с одинарной точностью ;

CALL BCPOL (FUNCT, npar, AR, ibtype, ZNIZWRK, ZVRXWRK, ftol, maxfcn, DTX, fvalue).

с двойной ;

CALL DBCPOL (FUNCT, npar, AR, ibtype, ZNIZWRK, ZVRXWRK, ftol, maxfcn, DTX, fvalue).

Если пользователь вызывает подпрограмму для вычислений с двойной точностью, то все вещественные аргументы должны быть описаны как REAL (8). Формальные аргументы означают:

FUNCTподпрограмма пользователя; вычисляет значение функции цели fcl по аргументам, передаваемым ей в массиве AR. В вызывающей программе должна быть описана оператором EXTERNAL. Оформляется пользователем как.

SUBROUTINE FUNCT (npar, AR, fcl), где ВХОД:

nparколичесгвово независимых переменных оптимизируемой функции;

AR — массив размерности npar, содержащий независимые переменные, для которых следует вычислить функцию цели; функция не должна изменять оптимизируемые переменные;

ВЫХОД:

fcl — значение функции цели, вычисленное по переменным в AR. Продолжение описания формальных аргументов пп (BCPOL) DBCPOL: ВХОД:

nparколичествово независимых переменных оптимизируемой функции;

ARмассив размерности npar, содержащий начальное значение оптимизируемых независимых переменных;

ibtype — скаляр, задающий вид ограничений на оптимизируемые переменные; оказывает следующие действия: ibtype = 0 — все ограничения задаются пользователем; ibtype = 1 — все переменные неотрицательные; ibtype = 2- все переменные неположительные;

ibtype = 3 — пользователь задаёт ограничения на первую переменную;

все другие переменные будут иметь те же ограничения; ZNIZWRK-массив размера праг, содержащий нижние границы переменных; ZVRXWRK-массив размера праг, содержащий верхние границы переменных; ftol — используется для выработки критерия сходимости. Первый критерий сходимости: вычисления останавливаются, когда относительная ошибка значений функции меньше, чем ftol, т. е. когда f (worst)-f (best) < ftol*(l+ABS (f (best))), где f (worst) и f (best) — соответственно значения функции в текущих худшей и лучшей точках многогранника. Второй критерий сходимости: вычисления останавливаются, когда стандартное отклонение значений функции в п+1 текущих точках меньше, чем ftol. Если подпрограмма завершает вычисления преждевременно, следует выполнить новый расчёт с меньшим значением ftol; maxfcn — на входе — максимально допустимое число обращений к функции FUNCT для вычисления целевой функции; на выходе — реальное число обращений к функции FUNCT;

ВЫХОД:

DTX — массив размерности праг, содержащий оптимальные переменные (найденное значение оптимизируемых независимых переменных); fvalue — скаляр, содержащий значение целевой функции в точке DTX (минимум функции цели).