Changeset 1272 for trunk/packages/optimizer

Timestamp:

Feb 4, 2009, 10:27:21 AM (15 years ago)

Author:

liveletlive

Message:

Committing additions made for including time of execution as a stoppage criteria.
also committing some changes to pgapack.h required to make these changes in stoppage criteria.
This code is yet untested, and will be tested once some more test benches are added.

Location:

trunk/packages/optimizer/src/pgapack/pgapack

Files:

• include/pgapack.h (modified) (5 diffs)
• source/create.c (modified) (1 diff)
• source/parallel.c (modified) (1 diff)
• source/stop.c (modified) (10 diffs)

trunk/packages/optimizer/src/pgapack/pgapack/include/pgapack.h

@@ -137 +137 @@
 #define PGA_STOP_NOCHANGE       2    /* Stop: no change in best string    */
 #define PGA_STOP_TOOSIMILAR     4    /* Stop: homogeneous population      */
+#define PGA_STOP_AV_FITNESS     5    /* Stop: if within a certain % of target av.fitness per pop.*/
+#define PGA_STOP_BEST_FITNESS   6    /* Stop: if best fitness in a pop is within % of target best fitness*/
+#define PGA_STOP_VARIANCE       7        /* Stop: if variance of fitness in pop is within % of target variance */
+#define PGA_STOP_TIMEELAPSED    8        /* Stop: if elapsed time >= target*/
 
 /*****************************************
@@ -145 +149 @@
 #define PGA_CROSSOVER_UNIFORM   3    /* Uniform   crossover                */
 #define PGA_CROSSOVER_SBX               4    /* Simulated Binary crossover                 */
-
+#define PGA_CROSSOVER_TRIANGULAR 5   /* Triangular Crossover*/
 /*****************************************
 *            SELECTION                   *
@@ -252 +256 @@
     int PopSize;             /* Number of strings to use                  */
     int StringLen;           /* string lengths                            */
+    double TgtFitnessVal;                /* Target Fitness, best or average depending on stpg criteria*/
+    double FitnessTol;           /* Min Fitness Tolerance to be achieved for stoppage*/ 
+    double TgtFitnessVar;       /* Target Fitness Variance for stoppage */
+    double VarTol;                       /* Tolerance for variance to be achieved for stoppage*/        
+    double TgtElapsedTime;   /* Max Execution time of GA before stoppage */
     int StoppingRule;        /* Termination Criteria                      */
     int MaxIter;             /* Maximum number of iterations to run       */
@@ -338 +347 @@
      double Online;
      double Average;
+     double Variance;
      double Best;
      time_t starttime;
@@ -799 +809 @@
 void PGASetMaxNoChangeValue(PGAContext *ctx, int max_no_change);
 void PGASetMaxSimilarityValue(PGAContext *ctx, int max_similarity);
+void PGASetTgtFitnessVal(PGAContext *ctx, double tgt_fitness_val);
+void PGASetFitnessTol(PGAContext *ctx,double fitness_tol);
+void PGASetTgtFitnessVariance(PGAContext *ctx,double tgt_fitness_var);
+void PGASetVarTol(PGAContext *ctx,double var_tol);
+void PGASetTgtElapsedTime(PGAContext *ctx,double tgt_elapsed_time);
 
 /*****************************************

trunk/packages/optimizer/src/pgapack/pgapack/source/create.c

@@ -459 +459 @@
     if ( ctx->ga.MaxSimilarity      == PGA_UNINITIALIZED_INT)
          ctx->ga.MaxSimilarity       = 95;
+         
+    if (ctx->ga.TgtFitnessVal           == PGA_UNINITIALIZED_DOUBLE){
+         ctx->ga.TgtFitnessVal          = 0;
+    }
+    
+    if (ctx->ga.FitnessTol                      == PGA_UNINITIALIZED_DOUBLE){
+         ctx->ga.FitnessTol                     = 0;
+    }
+    
+    if (ctx->ga.VarTol                          == PGA_UNINITIALIZED_DOUBLE){
+         ctx->ga.VarTol                         = 0;
+    }
+    
+    if (ctx->ga.TgtFitnessVar           == PGA_UNINITIALIZED_DOUBLE){
+         ctx->ga.TgtFitnessVar          = 0;
+    }
+    
+    if (ctx->ga.TgtElapsedTime          == PGA_UNINITIALIZED_DOUBLE){
+         ctx->ga.TgtElapsedTime         = 0;
+    }     
 
     if ( ctx->ga.NumReplace         == PGA_UNINITIALIZED_INT)

trunk/packages/optimizer/src/pgapack/pgapack/source/parallel.c

@@ -70 +70 @@
 #include "pgapack.h"
 #include <setjmp.h>
-#include <time.h>
 
 extern jmp_buf pgapack_jmpbuf;

trunk/packages/optimizer/src/pgapack/pgapack/source/stop.c

@@ -71 +71 @@
 
 #include "pgapack.h"
+#include <time.h>
 
 int *pgapack_abortPtr = NULL;
@@ -144 +145 @@
 {
     int done = PGA_FALSE;
-    double diff;
+    double diff, process_execution_time;
+    time_t endtime;
 
     PGADebugEntered("PGACheckStoppingConditions");
@@ -162 +164 @@
         if(((ctx->ga.StoppingRule & PGA_STOP_AV_FITNESS) == PGA_STOP_AV_FITNESS)){
                 if(ctx->ga.TgtFitnessVal == PGA_UNINITIALIZED_DOUBLE){
-                        fprintf(stderr,"Uninitialized Value 'Target Average'");
+                        PGAError(ctx,"Uninitialized Value 'Target Average'",PGA_FATAL, PGA_DOUBLE, (void *)&(ctx->ga.TgtFitnessVal));
                         done |= PGA_TRUE;
                 }else{
                         diff = ((ctx->ga.TgtFitnessVal - ctx->rep.Average)/(ctx->ga.TgtFitnessVal))*100;
-                        if(fabs(diff) <= ctx->ga.FitnessTol){
+                        if(fabs(diff) <= fabs((double)ctx->ga.FitnessTol)){
                                 done |= PGA_TRUE;
                         }
@@ -174 +176 @@
         if(((ctx->ga.StoppingRule & PGA_STOP_BEST_FITNESS) == PGA_STOP_BEST_FITNESS)){
                 if(ctx->ga.TgtFitnessVal == PGA_UNINITIALIZED_DOUBLE){
-                        fprintf(stderr,"Uninitialized Value 'Target Best Fitness'");
+                        PGAError(ctx,"Undefined Value 'Target Best Fitness'", PGA_FATAL, PGA_DOUBLE, (void *) &(ctx->ga.TgtFitnessVal));
                         done |= PGA_TRUE;
                 }else{
                         diff =(((double)ctx->ga.TgtFitnessVal -(double)ctx->rep.Best)/((double)ctx->ga.TgtFitnessVal))*100;
-                        if(fabs(diff) <= ctx->ga.FitnessTol){
+                        if(fabs(diff) <= fabs((double)ctx->ga.FitnessTol)){
                                 done |= PGA_TRUE;
                         }
@@ -186 +188 @@
         if(((ctx->ga.StoppingRule & PGA_STOP_VARIANCE) == PGA_STOP_VARIANCE)){
                 if(ctx->ga.TgtFitnessVar == PGA_UNINITIALIZED_DOUBLE){
-                        fprintf(stderr,"Uninitialized Value 'Target Fitness Variance'");
+                        PGAError(ctx,"Undefined Value 'Target Fitness Variance'",PGA_FATAL,PGA_DOUBLE, (void *)&(ctx->ga.TgtFitnessVar));
                         done |= PGA_TRUE;
                 }else{
                         diff =(((double)ctx->ga.TgtFitnessVar -(double)ctx->rep.Variance)/((double)ctx->ga.TgtFitnessVar))*100;
-                        if(fabs(diff) <= ctx->ga.VarTol){
+                        if(fabs(diff) <= fabs((double)ctx->ga.VarTol)){
                                 done |= PGA_TRUE;
                         }
@@ -197 +199 @@
         
         if(((ctx->ga.StoppingRule & PGA_STOP_TIMEELAPSED) == PGA_STOP_TIMEELAPSED)){
-                //TODO :GTG Add code for setting initial time and current time and comparing
-                done |= PGA_TRUE;
+                endtime = time(NULL);
+                process_execution_time = difftime(endtime, (time_t)ctx->rep.starttime);
+                if(ctx->ga.TgtElapsedTime == PGA_UNINITIALIZED_DOUBLE){
+                        PGAError(ctx,"Undefined target time of execution",PGA_FATAL,PGA_DOUBLE,(void *)&(ctx->ga.TgtElapsedTime));
+                }
+                if(((double)process_execution_time/60) >= ((double)ctx->ga.TgtElapsedTime)){
+                        done |= PGA_TRUE;
+                }
         }
 
@@ -482 +490 @@
         PGADebugEntered("PGASetFitnessTol");
     PGAFailIfSetUp("PGASetFitnessTol");
+    if(fitness_tol < 0){
+        PGAError(ctx,"Fitness tolerance invalid. Should be a positive quantity.\n",PGA_FATAL,PGA_DOUBLE,(void *)&fitness_tol);
+    }
         ctx->ga.FitnessTol = fitness_tol;
         PGADebugExited("PGASetFitnessTol");
@@ -509 +520 @@
         PGADebugEntered("PGASetTgtFitnessVariance");
     PGAFailIfSetUp("PGASetTgtFitnessVariance");
+    if(tgt_fitness_var < 0){
+        PGAError(ctx,"Fitness variance invalid. Should be a positive quantity.\n",PGA_FATAL,PGA_DOUBLE,(void *)&tgt_fitness_var);
+    }
         ctx->ga.TgtFitnessVar = tgt_fitness_var;
         PGADebugExited("PGASetTgtFitnessVariance");
@@ -534 +548 @@
         PGADebugEntered("PGASetVarTol");
     PGAFailIfSetUp("PGASetVarTol");
+    if(var_tol < 0){
+        PGAError(ctx,"Variance Tolerance Invalid, should be positive quantity", PGA_FATAL, PGA_DOUBLE, (void *) &var_tol);
+    }
         ctx->ga.VarTol = var_tol;
         PGADebugExited("PGASetVarTol");
@@ -560 +577 @@
         PGADebugEntered("PGASetTgtElapsedTime");
     PGAFailIfSetUp("PGASetTgtElapsedTime");
+    if(tgt_elapsed_time < 0){
+        PGAError(ctx,"Invalid Time of execution, cannot be a negative quantity", PGA_FATAL, PGA_DOUBLE, (void*) &tgt_elapsed_time);
+    }
         ctx->ga.TgtElapsedTime = tgt_elapsed_time;
         PGADebugExited("PGASetTgtElapsedTime");