Организация компьютерных экспериментов

****************************************************
*              Dis_Server                          *
*      Факторный отсеивающий эксперимент           *
****************************************************
Dis_Server_Results       MATRIX                ,2,2,2
      INITIAL Dis_Server_Results,UNSPECIFIED
Dis_Server_NextRunNumber EQU                    0
EXPERIMENT Dis_Server() BEGIN

     /* Наблюдение 1 */
     T1_ = 60;
     Koef = 0.5;
     Q_ = 300000;
     IF (StringCompare(DataType(Dis_Server_Results[1,1,1]),
               "UNSPECIFIED")'E'0)
     THEN BEGIN
/* Установить начальное значение переменной количества наблюдений */
     Dis_Server_NextRunNumber = 1;
/* Записать данные наблюдения и запустить процесс моделирования*/
     Dis_Server_GetResult();
     Dis_Server_Results[1,1,1] = N$ObrZap/N$KolZap;
     END;

     /* Наблюдение 2 */
     T1_ = 60;
     Koef = 0.5;
     Q_ = 700000;
     IF (StringCompare(DataType(Dis_Server_Results[1,1,2]),
          "UNSPECIFIED")'E'0)
     THEN BEGIN
/* Записать данные наблюдения и запустить процесс моделирования */
     Dis_Server_GetResult();
     Dis_Server_Results[1,1,2] = N$ObrZap/N$KolZap;
	 END;

     /* Наблюдения 3 - 7 для краткости пропущены */
     /* Наблюдение 8 */
     T1_ = 180;
     Koef = 1.5;
     Q_ = 700000;
     IF (StringCompare(DataType(Dis_Server_Results[2,2,2]),
               "UNSPECIFIED")'E'0)
     THEN BEGIN
/* Записать данные наблюдения и запустить процесс моделирования */
     Dis_Server_GetResult();
     Dis_Server_Results[2,2,2] = N$ObrZap/N$KolZap;
     END;
/* Эффекты смешивания в дробном факторном эксперименте */
     SE_Effects(Dis_Server_Results,"I");
END;

*******************************************************
*        Процедура запуска наблюдения                 *
*******************************************************
PROCEDURE Dis_Server_GetResult() BEGIN
/* Выполнить указанное число прогонов и записать результаты. */
/* Факторы для этого наблюдения уже были определены. */
TEMPORARY CurrentYield,ShowString,CommandString;
/* Вызов процедуры запуска */
     Dis_Server_Run(Dis_Server_NextRunNumber);
     CurrentYield = N$ObrZap/N$KolZap;
     ShowString = PolyCatenate("Run ",String(Dis_Server_NextRunNumber),
	      ". ", "" );
     ShowString = PolyCatenate(ShowString," Yield=",String(CurrentYield),
	      ". ");
     ShowString = PolyCatenate(ShowString," T1_=",String(T1_), ";" );
     ShowString = PolyCatenate(ShowString," Koef=",String(Koef), ";" );
     ShowString = PolyCatenate(ShowString," Q_=",String(Q_), ";" );
     CommandString = PolyCatenate("SHOW """,ShowString,
          """", "" );
     DoCommand(CommandString);
     Dis_Server_NextRunNumber = Dis_Server_NextRunNumber + 1;
     RETURN CurrentYield;
END;
*******************************************************
*             Процедура запуска                       *
*******************************************************
PROCEDURE Dis_Server_Run(Run_Number) BEGIN
     DoCommand("CLEAR OFF"); /* Использовать OFF для сохранения результата. */
/* Увеличьте число команд RMULT, если у вас большее число ГСЧ. */
/* Задать новые случайные числа всем потокам случайных чисел. */
     TEMPORARY CommandString;
/* Вычислить, прежде чем перейти к DoCommand. */
     CommandString = Catenate("RMULT ",Run_Number#111);
/* DoCommand контролирует строку в глобальном контексте. */
     DoCommand(CommandString);
/* Установить собственные условия наблюдения. */
     DoCommand("START 100,NP");  /* Пройти неустоявшийся режим. */
     DoCommand("RESET");         /* Начать период измерений. */
     DoCommand("START 9604,NP"); /* Провести моделирование. */
END;