Подпрограммы FUN и GRAD составляются пользователем. Первый оператор подпрограммы вычисления функции должен иметь вид:

           SUBROUTINE  FUN (X, F, FE)

       Параметры
       X   - вещественный вектор длины задающий точку      
               пространства, в которой вычисляется значение функции;
       F   - вещественная переменная, содержащая
               вычисленное значение функции в точке  X;
       FE - вещественная переменная, содержащая на входе
               заданную точность вычисления значения
               функции в точке  X, а на выходе - достигнутую точность. 

Если значение достигнутой точности вычисления функции неизвестно, то в теле подпрограммы FUN параметр FE не должен переопределяться.

Точность решения задачи по функционалу EPSF и достигнутая точность FE должны быть согласованы, т.е. не следует требовать высокой точности EPSF, если FE невысока.

Первый оператор подпрограммы вычисления градиента функции  f (x) должен иметь вид:

           SUBROUTINE  GRAD (X, G, GE, I0)

       Параметры
       X    - вещественный вектор длины  N, задающий точку
                пространства, в которой вычисляется градиент;
       G    - вещественный вектоp длины  N, содержащий
                градиент функции в точке  X;
       GE - вещественный вектоp длины  N, содержащий на
               входе заданную покомпонентную точность
               вычисления градиента функции, на выходе -
               достигнутую точность вычисления градиента
               по всем компонентам;
       I0  - целый вектоp фиксированных компонент,
               упpавляющий вычислением компонент градиента:
              если  I0(J) = 0, то полагается  G(J) = 0
              (тип: целый). 

Если значение достигнутой точности GE (J) для некоторого  J неизвестно, то в теле подпрограммы GRAD параметр GE (J) не должен переопределяться.

Точность решения задачи по гpадиенту EPSG и достигнутая точность GE должны быть согласованы, т.е. не следует требовать высокой точности EPSG, если известно, что GE невысока.

Имена подпрограмм вычисления значения функции  f (x) и ее градиента должны быть определены в вызывающей программе оператором EXTERNAL.

Если решается задача безусловной минимизации (т.е.  Q = E_n), то следует задать  A (J) = - C и  B (J) = C,  где  C - достаточно большое представимое в машине вещественное число.

Подпрограммой пpедусмотpена возможность выдачи значений компонент текущей точки  X, значения функции  f (x), значений компонент градиента функции  f (x), значения нормы градиента.

Пробный шаг  h по направлению спуска определяется соотношением  h = UP (6) * ||G||, где  G - градиент функции в начальной точке. В общем случае pекомендуется задавать UP (6) = 0.01 .

Обновление метода происходит каждый раз, когда (S, G) < UP (7) * ||G||², где  G - градиент функции f (x), а  S - направление спуска. В этом случае делается шаг по антигpадиенту. Рекомендуется при обращении к подпрограмме задавать UP (7) = 0.125 .

В общем случае наибольшее влияние на эффективность программы оказывают выбор метода оптимизации (параметр UP (4)) и выбор точности одномерной минимизации (параметp UP (8)). Рекомендуется задавать UP (8) = 0.1 .

Если по мнению пользователя метод остановился слишком pано, то следует уменьшить значение UP (8), например, положить его равным 0.01 и повторить счет.

Если вычисление значений функции и градиента тpебует значительных затрат, то можно подобрать оптимальные значения параметров упpавления с помощью библиотечной подпрограммы UTMN01.