Changeset 17657 for branches/2825-NSGA3

Timestamp:

07/07/20 21:51:44 (4 years ago)

Author:

dleko

Message:

#2825 Implement recombination.

Location:

branches/2825-NSGA3/HeuristicLab.Algorithms.NSGA3/3.3

Files:

• HeuristicLab.Algorithms.NSGA3-3.3.csproj (modified) (2 diffs)
• NSGA3.cs (modified) (10 diffs)
• Plugin.cs (modified) (1 diff)
• SimulatedBinaryCrossover.cs (added)
• Solution.cs (modified) (2 diffs)

branches/2825-NSGA3/HeuristicLab.Algorithms.NSGA3/3.3/HeuristicLab.Algorithms.NSGA3-3.3.csproj

@@ -69 +69 @@
       <Private>False</Private>
     </Reference>
+    <Reference Include="HeuristicLab.Operators-3.3, Version=3.3.0.0, Culture=neutral, PublicKeyToken=ba48961d6f65dcec" />
     <Reference Include="HeuristicLab.Optimization-3.3, Version=3.3.0.0, Culture=neutral, PublicKeyToken=ba48961d6f65dcec, processorArchitecture=MSIL">
       <SpecificVersion>False</SpecificVersion>
@@ -112 +113 @@
     <Compile Include="Properties\AssemblyInfo.cs" />
     <Compile Include="ReferencePoint.cs" />
+    <Compile Include="SimulatedBinaryCrossover.cs" />
     <Compile Include="Solution.cs" />
     <Compile Include="Utility.cs" />

branches/2825-NSGA3/HeuristicLab.Algorithms.NSGA3/3.3/NSGA3.cs

@@ -1 +1 @@
 using System;
 using System.Collections.Generic;
+using System.Diagnostics;
 using System.Linq;
 using System.Threading;
@@ -64 +65 @@
         private const string SetSeedRandomlyName = "SetSeedRandomly";
         private const string PopulationSizeName = "PopulationSize";
-        private const string CrossoverProbabilityName = "CrossOverProbability";
+        private const string CrossoverProbabilityName = "CrossoverProbability";
+        private const string CrossoverContiguityName = "CrossoverContiguity";
         private const string MutationProbabilityName = "MutationProbability";
         private const string MaximumGenerationsName = "MaximumGenerations";
@@ -98 +100 @@
         }
 
+        private IFixedValueParameter<DoubleValue> CrossoverContiguityParameter
+        {
+            get { return (IFixedValueParameter<DoubleValue>)Parameters[CrossoverContiguityName]; }
+        }
+
         private IFixedValueParameter<PercentValue> MutationProbabilityParameter
         {
@@ -124 +131 @@
 
         public PercentValue CrossoverProbability => CrossoverProbabilityParameter.Value;
+
+        public DoubleValue CrossoverContiguity => CrossoverContiguityParameter.Value;
 
         public PercentValue MutationProbability => MutationProbabilityParameter.Value;
@@ -159 +168 @@
             Parameters.Add(new FixedValueParameter<IntValue>(SeedName, "The random seed used to initialize the new pseudo random number generator.", new IntValue(0)));
             Parameters.Add(new FixedValueParameter<BoolValue>(SetSeedRandomlyName, "True if the random seed should be set to a random value, otherwise false.", new BoolValue(true)));
-            Parameters.Add(new FixedValueParameter<IntValue>(PopulationSizeName, "The size of the population of Individuals.", new IntValue(100)));
+            Parameters.Add(new FixedValueParameter<IntValue>(PopulationSizeName, "The size of the population of Individuals.", new IntValue(200)));
             Parameters.Add(new FixedValueParameter<PercentValue>(CrossoverProbabilityName, "The probability that the crossover operator is applied on two parents.", new PercentValue(0.9)));
+            Parameters.Add(new FixedValueParameter<DoubleValue>(CrossoverContiguityName, "The contiguity value for the Simulated Binary Crossover that specifies how close a child should be to its parents (larger value means closer). The value must be greater than or equal than 0. Typical values are in the range [2;5]."));
             Parameters.Add(new FixedValueParameter<PercentValue>(MutationProbabilityName, "The probability that the mutation operator is applied on a Individual.", new PercentValue(0.05)));
             Parameters.Add(new FixedValueParameter<IntValue>(MaximumGenerationsName, "The maximum number of generations which should be processed.", new IntValue(1000)));
@@ -195 +205 @@
 
             InitFields();
+            InitResults();
             InitReferencePoints();
-            InitResults();
             Analyze();
         }
@@ -218 +228 @@
                 RealVector randomRealVector = new RealVector(Problem.Encoding.Length, random, minBound, maxBound);
 
-                solutions.Add(new Solution(StorableConstructorFlag.Default)
-                {
-                    Chromosome = randomRealVector
-                });
+                solutions.Add(new Solution(randomRealVector));
                 solutions[i].Fitness = Evaluate(solutions[i].Chromosome);
             }
@@ -236 +243 @@
         private void InitResults()
         {
-            Results.Add(new Result(GeneratedReferencePointsResultName, "The initially generated reference points", Utility.ConvertToDoubleMatrix(ReferencePoints)));
+            Results.Add(new Result(GeneratedReferencePointsResultName, "The initially generated reference points", new DoubleMatrix()));
             Results.Add(new Result(CurrentFrontResultName, "The Pareto Front", new DoubleMatrix()));
         }
@@ -288 +295 @@
 
             // Do non-dominated sort
-            var qualities = Utility.ToFitnessMatrix(solutions);
+            var qualities = Utility.ToFitnessMatrix(rt);
             // compute the pareto fronts using the DominationCalculator and discard the qualities
             // part in the inner tuples
@@ -578 +585 @@
         private List<Solution> Recombine(List<Solution> solutions)
         {
-            throw new NotImplementedException();
+            List<Solution> childSolutions = new List<Solution>();
+
+            for (int i = 0; i < solutions.Count; i += 2)
+            {
+                int parentIndex1 = random.Next(solutions.Count);
+                int parentIndex2 = random.Next(solutions.Count);
+                // ensure that the parents are not the same object
+                if (parentIndex1 == parentIndex2) parentIndex2 = (parentIndex2 + 1) % solutions.Count;
+                var parent1 = solutions[parentIndex1];
+                var parent2 = solutions[parentIndex2];
+
+                // Do crossover with crossoverProbabilty == 1 in order to guarantee that a crossover happens
+                var children = SimulatedBinaryCrossover.Apply(random, Problem.Encoding.Bounds, parent1.Chromosome, parent2.Chromosome, 1);
+                Debug.Assert(children != null);
+
+                var child1 = new Solution(children.Item1);
+                var child2 = new Solution(children.Item2);
+                child1.Fitness = Evaluate(child1.Chromosome);
+                child2.Fitness = Evaluate(child1.Chromosome);
+
+                childSolutions.Add(child1);
+                childSolutions.Add(child2);
+            }
+
+            return childSolutions;
         }
 
         private List<Solution> Mutate(List<Solution> solutions)
         {
-            throw new NotImplementedException();
+            foreach (var solution in solutions)
+            {
+                UniformOnePositionManipulator.Apply(random, solution.Chromosome, Problem.Encoding.Bounds);
+            }
+            return solutions;
         }
 

branches/2825-NSGA3/HeuristicLab.Algorithms.NSGA3/3.3/Plugin.cs

@@ -17 +17 @@
     [PluginDependency("HeuristicLab.Encodings.BinaryVectorEncoding", "3.3")]
     [PluginDependency("HeuristicLab.Encodings.RealVectorEncoding", "3.3")]
+    [PluginDependency("HeuristicLab.Operators", "3.3")]
     [PluginDependency("HeuristicLab.Optimization", "3.3")]
     [PluginDependency("HeuristicLab.Parameters", "3.3")]

branches/2825-NSGA3/HeuristicLab.Algorithms.NSGA3/3.3/Solution.cs

@@ -3 +3 @@
 using HeuristicLab.Common;
 using HeuristicLab.Encodings.RealVectorEncoding;
-using HeuristicLab.Optimization;
 
 namespace HeuristicLab.Algorithms.NSGA3
 {
     [StorableType("D35C9D79-77DF-4CF3-AB44-FCCC00E44C5F")]
-    public class Solution : Individual, IDeepCloneable
+    public class Solution : IDeepCloneable
     {
         // Chromosome
@@ -16 +15 @@
         public double[] Fitness { get; set; }
 
-        public Solution(StorableConstructorFlag _) : base(_)
+        public Solution(RealVector chromosome)
         {
+            Chromosome = chromosome;
+        }
+
+        public Solution(Solution solution, Cloner cloner)
+        {
+            Chromosome = cloner.Clone(solution.Chromosome);
+            Fitness = new double[solution.Fitness.Length];
+            Array.Copy(solution.Fitness, Fitness, solution.Fitness.Length);
         }
 
         public IDeepCloneable Clone(Cloner cloner)
         {
-            throw new NotImplementedException();
+            return new Solution(this, cloner);
         }
 
         public object Clone()
         {
-            throw new NotImplementedException();
-        }
-
-        public override TEncoding GetEncoding<TEncoding>()
-        {
-            throw new NotImplementedException();
-        }
-
-        public override Individual Copy()
-        {
-            throw new NotImplementedException();
+            return new Cloner().Clone(this);
         }
     }