Changeset 1166 for trunk/packages/optimizer/src

Timestamp:

Sep 27, 2008, 1:58:13 AM (16 years ago)

Author:

liveletlive

Message:

Changes Made: Vanilla restart added, no modifications can be made to existing parameters in the middle of a run. Only mutation creates a new population on restart.
New configuration options added: MutationandCrossover? option. amount to be mutated made configurable. configration for (not)allowing duplicate strings added. app-qdot and simple modified to keep pace with these changes.

Location:

trunk/packages/optimizer/src

Files:

• plugin_pgapack.c (modified) (8 diffs)
• rp_optimizer.c (modified) (1 diff)
• rp_optimizer.h (modified) (2 diffs)
• rp_optimizer_tcl.c (modified) (4 diffs)

trunk/packages/optimizer/src/plugin_pgapack.c

@@ -22 +22 @@
     int popSize;         /* population size <=> PGASetPopSize() */
     int popRepl;         /* replacement <=> PGASetPopReplacementType() */
+    int numReplPerPop;    /* number of new strings created per population, the rest are the best strings from the previous population*/
     int stpcriteria;     /*stoppage criteria <=> PGASetStoppingRuleType()*/
     int randnumseed;  /*Random Number Seed <=> PGASetRandomSeed()*/
     double mutnrate;     /*Mutation Rate <=> PGASetMutationProb()*/
+    double mutnValue;     /*use this value while mutating*/
     double crossovrate;  /*Crossover Rate <=> PGASetCrossoverProb();*/
+    int allowdup;        /*Allow duplicate strings in the population or not*/
+    int mutnandcrossover;/*By default strings that do not undergo crossover undergo mutation, this option allows strings to crossover and be mutated*/
 } PgapackData;
 
@@ -58 +62 @@
   {"-operation", &RpOption_Oper, Rp_Offset(PgapackData,operation)},
   {"-poprepl", &RpOption_PopRepl, Rp_Offset(PgapackData,popRepl)},
+  {"-numReplPerPop",RP_OPTION_INT,Rp_Offset(PgapackData,numReplPerPop)},
   {"-popsize", RP_OPTION_INT, Rp_Offset(PgapackData,popSize)},
   {"-mutnrate",RP_OPTION_DOUBLE,Rp_Offset(PgapackData,mutnrate)},
+  {"-mutnValue",RP_OPTION_DOUBLE,Rp_Offset(PgapackData,mutnValue)},
   {"-crossovrate",RP_OPTION_DOUBLE,Rp_Offset(PgapackData,crossovrate)},
   {"-randnumseed",RP_OPTION_INT,Rp_Offset(PgapackData,randnumseed)},
   {"-stpcriteria",&RpOption_StpCriteria,Rp_Offset(PgapackData,stpcriteria)},
+  {"-allowdup",RP_OPTION_BOOLEAN,Rp_Offset(PgapackData,allowdup)},
+  {"-mutnandcrossover",RP_OPTION_BOOLEAN,Rp_Offset(PgapackData,mutnandcrossover)},
   {NULL, NULL, 0}
 };
@@ -104 +112 @@
     dataPtr->popRepl = PGA_POPREPL_BEST;
     dataPtr->popSize = 200;
+    dataPtr->numReplPerPop = (dataPtr->popSize)/10; /*10% replaced by default, change to whatever value you need*/
     dataPtr->crossovrate = 0.85;
     dataPtr->mutnrate = 0.05; /*by default in PGAPack 1/stringlength*/
+    dataPtr->mutnValue = 0.01;/*value of this number will be changed by plus/minus hundredth of its current value*/
     dataPtr->randnumseed = 1; /*should be a number greater than one, PGAPack requires it*/
     dataPtr->stpcriteria = PGA_STOP_NOCHANGE;
+    dataPtr->allowdup = PGA_FALSE; /*Do not allow duplicate strings by default*/
+    dataPtr->mutnandcrossover = PGA_FALSE;/*do not allow mutation and crossover to take place on the same string by default*/
     return (ClientData)dataPtr;
 }
 
 int pgapack_abort = 0;
+int pgapack_restart_user_action = 0;
 
 /*
@@ -136 +149 @@
     pgapack_abort = 0;          /* FALSE */
     PGASetAbortVar(&pgapack_abort);
+    PGASetRestartUserAction(&pgapack_restart_user_action);
 
     ctx = PGACreate(&argc, argv, PGA_DATATYPE_USER, envPtr->numParams,
@@ -145 +159 @@
     PGASetStoppingRuleType(ctx, dataPtr->stpcriteria);
     PGASetMutationProb(ctx,dataPtr->mutnrate);
+    PGASetMutationRealValue(ctx,dataPtr->mutnValue);
     PGASetCrossoverProb(ctx,dataPtr->crossovrate);
     PGASetRandomSeed(ctx,dataPtr->randnumseed);
     PGASetCrossoverType(ctx, PGA_CROSSOVER_UNIFORM);
+    PGASetNoDuplicatesFlag(ctx,dataPtr->allowdup);
+    PGASetMutationAndCrossoverFlag(ctx,dataPtr->mutnandcrossover);
+    PGASetNumReplaceValue(ctx,dataPtr->numReplPerPop);
 
 
@@ -342 +360 @@
 
     int n, ival,tempmr;
+    double mutnVal;
     RpOptimEnv *envPtr;
     RpOptimParam *paramPtr;
@@ -358 +377 @@
                 if(numPtr->mutnrate!=PARAM_NUM_UNSPEC_MUTN_RATE){
                         tempmr = numPtr->mutnrate;
+                        mutnVal = numPtr->mutnValue;
                 }else{
                         tempmr = mr;
+                        mutnVal = PGAGetMutationRealValue(ctx);
                 }
                 if (PGARandomFlip(ctx, tempmr)) {
@@ -367 +388 @@
                         /* bump the value up/down a little, randomly */
                         if (PGARandomFlip(ctx, 0.5)) {
-                            paramPtr[n].value.dval += 0.1*paramPtr[n].value.dval;
+                            paramPtr[n].value.dval += mutnVal*paramPtr[n].value.dval; /*Made the mutation amount configurable*/
                         } else {
-                            paramPtr[n].value.dval -= 0.1*paramPtr[n].value.dval;
+                            paramPtr[n].value.dval -= mutnVal*paramPtr[n].value.dval;
                         }
                         /* make sure the resulting value is still in bounds */

trunk/packages/optimizer/src/rp_optimizer.c

@@ -115 +115 @@
     numPtr->max = DBL_MAX;
     numPtr->mutnrate = PARAM_NUM_UNSPEC_MUTN_RATE; /*Unspecified by default, should result in global mutn rate being applied*/
+    numPtr->mutnValue = 0.01; /* value of this number will be changed by plus/minus hundredth of its current value*/
     numPtr->randdist = RAND_NUMBER_DIST_GAUSSIAN ;/*Gaussian by default*/
     numPtr->stddev = 1; 

trunk/packages/optimizer/src/rp_optimizer.h

@@ -48 +48 @@
  */
 typedef enum {
-    RP_OPTIM_SUCCESS=0, RP_OPTIM_UNKNOWN, RP_OPTIM_FAILURE, RP_OPTIM_ABORTED
+    RP_OPTIM_SUCCESS=0, RP_OPTIM_UNKNOWN, RP_OPTIM_FAILURE, RP_OPTIM_ABORTED, RP_OPTIM_RESTARTED
 } RpOptimStatus;
 
@@ -107 +107 @@
     double max;                     /* optimization constraint: max value */
     double mutnrate;                    /* independently sets mutation rate for each parameter*/
+    double mutnValue;               /* on mutation bumps up/down the value by val (+-) mutnValue*val */
     int randdist;                       /* gaussian or uniform distribution*/
     int strictmin;                                      /* whether a strict min is to be applied for gauss. rand. numbers*/

trunk/packages/optimizer/src/rp_optimizer_tcl.c

@@ -17 +17 @@
 
 extern int pgapack_abort;
+extern int pgapack_restart_user_action;
 
 /*
@@ -57 +58 @@
   {"-max", RP_OPTION_DOUBLE, Rp_Offset(RpOptimParamNumber,max)},
   {"-mutnrate",RP_OPTION_DOUBLE, Rp_Offset(RpOptimParamNumber,mutnrate)},
+  {"-mutnValue",RP_OPTION_DOUBLE, Rp_Offset(RpOptimParamNumber,mutnValue)},
   {"-randdist",&RpOption_RandDist,Rp_Offset(RpOptimParamNumber,randdist)},
   {"-strictmin",RP_OPTION_BOOLEAN,Rp_Offset(RpOptimParamNumber,strictmin)},
@@ -83 +85 @@
 #endif
 
-extern int pgapack_abort;
+
 extern void PGARuntimeDataTableInit();
 extern void PGARuntimeDataTableDeInit();
@@ -394 +396 @@
         pgapack_abort = value;
         return TCL_OK;
-    } else if (*option == 'g' && strcmp(option,"get") == 0) {
+    }else if (*option == 'r' && strcmp(option,"restart") == 0){
+                int value;
+        if (objc < 3) {
+            Tcl_WrongNumArgs(interp, 1, objv, "restart bool");
+            return TCL_ERROR;
+        }
+                if (Tcl_GetBooleanFromObj(interp, objv[2], &value) != TCL_OK) {
+                return TCL_ERROR;
+                }
+                pgapack_restart_user_action = value;
+                return TCL_OK;    
+        }else if (*option == 'g' && strcmp(option,"get") == 0) {
         /*
          * OPTION:  get ?globPattern? ?-option?