OPTION, ECHO=FALSE;
OPTION, INFO=TRUE;

TITLE, STRING="OPAL FD MOGA Optimization, 2020-06-10";

REAL betax0 = 3.73099;
REAL betay0 = betax0 / 2.0;
REAL eps = 1E-6;
REAL sigmax = sqrt(betay0 * eps); // at the end of the FD
REAL sigmay = sqrt(betax0 * eps); // x and y trade sizes
REAL loc = 1.0;

dv0: DVAR, VARIABLE="QFX1_K1", LOWERBOUND=  0, UPPERBOUND= 15;
dv1: DVAR, VARIABLE="QDX1_K1", LOWERBOUND=-15, UPPERBOUND=  0;

drmsx:   OBJECTIVE, EXPR="fabs(statVariableAt('rms_x', ${loc}) - ${sigmax})";
drmsy:   OBJECTIVE, EXPR="fabs(statVariableAt('rms_y', ${loc}) - ${sigmay})";

OPTIMIZE, INPUT="tmpl/fdopt.tmpl", OBJECTIVES = {drmsx, drmsy},
          DVARS = {dv0, dv1}, INITIALPOPULATION=10, MAXGENERATIONS=1,
          NUM_IND_GEN=10, MUTATION_PROBABILITY=0.43,
          NUM_MASTERS=1, NUM_COWORKERS=1, SIMTMPDIR="simtmpdir",
          TEMPLATEDIR="tmpl", FIELDMAPDIR="fieldmaps", OUTPUT="optLinac",
          OUTDIR="results";

QUIT;