Changeset 15563 for branches/GeneralizedQAP/HeuristicLab.Problems.GeneralizedQuadraticAssignment.Algorithms/3.3/Evolutionary

Timestamp:

12/30/17 23:10:29 (7 years ago)

Author:

abeham

Message:

#1614:

• Added LAHC and pLAHC-s
• Changed all algorithms to update high frequency results only every second

Location:

branches/GeneralizedQAP/HeuristicLab.Problems.GeneralizedQuadraticAssignment.Algorithms/3.3/Evolutionary

Files:

• ESContext.cs (modified) (4 diffs)
• ESGQAPSolution.cs (modified) (3 diffs)
• EvolutionStrategy.cs (modified) (7 diffs)

branches/GeneralizedQAP/HeuristicLab.Problems.GeneralizedQuadraticAssignment.Algorithms/3.3/Evolutionary/ESContext.cs

@@ -45 +45 @@
     }
 
+    [Storable]
+    private IValueParameter<IRandom> normalRand;
+    public IRandom NormalRand {
+      get { return normalRand.Value; }
+      set { normalRand.Value = value; }
+    }
+
     public void ReplacePopulation(IEnumerable<ISingleObjectiveSolutionScope<ESGQAPSolution>> replacement) {
       Scope.SubScopes.Replace(replacement);
@@ -55 +62 @@
       problem = cloner.Clone(original.problem);
       bestSolution = cloner.Clone(original.bestSolution);
+      normalRand = cloner.Clone(original.normalRand);
     }
     public ESContext() : this("Evolution Strategy (GQAP) Context") { }
@@ -60 +68 @@
       Parameters.Add(problem = new ValueParameter<GQAP>("Problem"));
       Parameters.Add(bestSolution = new ValueParameter<ESGQAPSolution>("BestFoundSolution"));
+      Parameters.Add(normalRand = new ValueParameter<IRandom>("NormalRand"));
     }
 
@@ -74 +83 @@
       scope.Variables.Add(new Variable(name, code.Assignment));
       scope.Variables.Add(new Variable("Evaluation", code.Evaluation));
+      scope.Variables.Add(new Variable("StrategyParameter", code.sParam));
       return scope;
     }

branches/GeneralizedQAP/HeuristicLab.Problems.GeneralizedQuadraticAssignment.Algorithms/3.3/Evolutionary/ESGQAPSolution.cs

@@ -21 +21 @@
 
 using HeuristicLab.Common;
+using HeuristicLab.Data;
 using HeuristicLab.Encodings.IntegerVectorEncoding;
 using HeuristicLab.Persistence.Default.CompositeSerializers.Storable;
@@ -28 +29 @@
 
     [Storable]
-    private double sParam;
+    internal DoubleValue sParam;
     public double SParam {
-      get { return sParam; }
+      get { return sParam.Value; }
       set {
-        if (sParam == value) return;
-        sParam = value;
+        if (sParam.Value == value) return;
+        sParam.Value = value;
         OnPropertyChanged(nameof(SParam));
       }
@@ -42 +43 @@
     protected ESGQAPSolution(ESGQAPSolution original, Cloner cloner)
     : base(original, cloner) {
-      sParam = original.sParam;
+      sParam = cloner.Clone(original.sParam);
     }
     public ESGQAPSolution(IntegerVector assignment, Evaluation eval, double sParam)
       : base(assignment, eval) {
-      this.sParam = sParam;
+      this.sParam = new DoubleValue(sParam);
     }
     

branches/GeneralizedQAP/HeuristicLab.Problems.GeneralizedQuadraticAssignment.Algorithms/3.3/Evolutionary/EvolutionStrategy.cs

@@ -106 +106 @@
       base.Initialize(cancellationToken);
 
+      Context.NormalRand = new NormalDistributedRandom(Context.Random, 0, 1);
       Context.Problem = Problem;      
       Context.BestQuality = double.NaN;
       Context.BestSolution = null;
-
+      
       for (var m = 0; m < Mu; m++) {
         var assign = new IntegerVector(Problem.ProblemInstance.Demands.Length, Context.Random, 0, Problem.ProblemInstance.Capacities.Length);
@@ -115 +116 @@
         Context.EvaluatedSolutions++;
 
-        var ind = new ESGQAPSolution(assign, eval, 1.0 / assign.Length);
+        var ind = new ESGQAPSolution(assign, eval, Context.Random.NextDouble() * 2 - 1);
         var fit = Problem.ProblemInstance.ToSingleObjective(eval);
         Context.AddToPopulation(Context.ToScope(ind, fit));
@@ -133 +134 @@
 
     protected override void Run(CancellationToken cancellationToken) {
+      var lastUpdate = ExecutionTime;
+
       while (!StoppingCriterion()) {
         var nextGen = new List<ISingleObjectiveSolutionScope<ESGQAPSolution>>(Lambda);
@@ -141 +144 @@
           var offspring = (ESGQAPSolution)m.Solution.Clone();
           var count = Mutate(m, offspring);
-          offspring.SParam += ((1.0 / count) - offspring.SParam) / 10.0;
+          offspring.SParam += 0.7071 * Context.NormalRand.NextDouble(); //((1.0 / count) - offspring.SParam) / 10.0;
 
           offspring.Evaluation = Problem.ProblemInstance.Evaluate(offspring.Assignment);
@@ -162 +165 @@
 
         IResult result;
-        if (Results.TryGetValue("Iterations", out result))
-          ((IntValue)result.Value).Value = Context.Iterations;
-        else Results.Add(new Result("Iterations", new IntValue(Context.Iterations)));
-        if (Results.TryGetValue("EvaluatedSolutions", out result))
-          ((IntValue)result.Value).Value = Context.EvaluatedSolutions;
-        else Results.Add(new Result("EvaluatedSolutions", new IntValue(Context.EvaluatedSolutions)));
+        if (ExecutionTime - lastUpdate > TimeSpan.FromSeconds(1)) {
+          if (Results.TryGetValue("Iterations", out result))
+            ((IntValue)result.Value).Value = Context.Iterations;
+          else Results.Add(new Result("Iterations", new IntValue(Context.Iterations)));
+          if (Results.TryGetValue("EvaluatedSolutions", out result))
+            ((IntValue)result.Value).Value = Context.EvaluatedSolutions;
+          else Results.Add(new Result("EvaluatedSolutions", new IntValue(Context.EvaluatedSolutions)));
+          lastUpdate = ExecutionTime;
+        }
         if (Results.TryGetValue("BestQuality", out result))
           ((DoubleValue)result.Value).Value = Context.BestQuality;
@@ -180 +186 @@
         if (cancellationToken.IsCancellationRequested) break;
       }
+      IResult result2;
+      if (Results.TryGetValue("Iterations", out result2))
+        ((IntValue)result2.Value).Value = Context.Iterations;
+      else Results.Add(new Result("Iterations", new IntValue(Context.Iterations)));
+      if (Results.TryGetValue("EvaluatedSolutions", out result2))
+        ((IntValue)result2.Value).Value = Context.EvaluatedSolutions;
+      else Results.Add(new Result("EvaluatedSolutions", new IntValue(Context.EvaluatedSolutions)));
     }
 
     private int Mutate(ISingleObjectiveSolutionScope<ESGQAPSolution> m, ESGQAPSolution offspring) {
+      var stopProb = (Math.Tanh(m.Solution.SParam) + 1) / 2.0; // squash strategy parameter to ]0;1[
       var offspringFeasible = offspring.Evaluation.IsFeasible;
       double[] slack = null;
@@ -210 +224 @@
           offspring.Assignment[equip] = newLoc;
         }
-        if (Context.Random.NextDouble() < m.Solution.SParam) break;
+        if (Context.Random.NextDouble() < stopProb) break;
         count++;
       }