Changeset 17440 for branches/3055_SMSEMOA/HeuristicLab.Algorithms.SMSEMOA

Timestamp:

02/17/20 16:30:40 (5 years ago)

Author:

kyang

Message:

#3055: Fixed the bug of "no reference point" on real-world test problems

Location:

branches/3055_SMSEMOA/HeuristicLab.Algorithms.SMSEMOA/3.4

Files:

• SMSEMOAAlgorithm.cs (modified) (1 diff)
• SMSEMOAlgorithmBase.cs (modified) (4 diffs)

branches/3055_SMSEMOA/HeuristicLab.Algorithms.SMSEMOA/3.4/SMSEMOAAlgorithm.cs

@@ -103 +103 @@
       var analyzer = Analyzer;
       var rand = RandomParameter.Value;
-      int lambda = 1;            // the size of offspring 
+      //int lambda = 1;            // the size of offspring 
+      var lambda = Lambda.Value;
       //int indexOffspring = 0;    // the index of offspring to be generated
 

branches/3055_SMSEMOA/HeuristicLab.Algorithms.SMSEMOA/3.4/SMSEMOAlgorithmBase.cs

@@ -101 +101 @@
     private DoubleArray ReferencePoint {
       get {
-        var problem = (MultiObjectiveTestFunctionProblem)Problem;
-        return problem.ReferencePoint;
+        if (Problem is MultiObjectiveTestFunctionProblem) {
+          var problem = (MultiObjectiveTestFunctionProblem)Problem;
+          return problem.ReferencePoint;
+        }
+        else {
+        return null;
+        }
       }
     }
@@ -250 +255 @@
 
     #region constructors
-    public SMSEMOAAlgorithmBase() { 
+    public SMSEMOAAlgorithmBase() {
       // Add or define or specify the parameters that may be use in SMS-EMOA.   
       // ***("Name", "Description", "Value")
@@ -480 +485 @@
         var lastLayer = pf2.Last();
 
-        try { // TODO: This can be use for dynamic reference point strategy later. Kaifeng , 02/2020
-          // smetric = Hypervolume.Calculate(lastLayer.Select(x => x.Item2), Enumerable.Repeat(11d, NadirPoint.Length).ToArray(), maximization);
-          smetric = Hypervolume.Calculate(lastLayer.Select(x => x.Item2), ReferencePoint.ToArray(), maximization);
-        }
-        catch {
-          smetric = int.MinValue;
-        }
-
-        var indices = Enumerable.Range(0, lastLayer.Count()).ToList();
-
-        for (int ii = 0; ii < lastLayer.Count(); ++ii) {
+        // TODO: This can be use for dynamic reference point strategy later. Kaifeng , 02/2020
+        // smetric = Hypervolume.Calculate(lastLayer.Select(x => x.Item2), Enumerable.Repeat(11d, NadirPoint.Length).ToArray(), maximization);
+
+
+        var reference = Enumerable.Repeat(double.MaxValue, maximization.Length).ToArray();
+        if (ReferencePoint is null) {     // KF, 20200217 -- fix no reference point on real-world applications. If No reference points in Algorithms, use 1.1 \times max objective values as the reference point
+          for (int i = 0; i < reference.Length; i++) {
+            reference[i] = 1.1 * reference[i];
+            if (reference[i] > 10000) {
+              reference[i] = 9999;          // set a upper bound for the reference point
+            }
+          }
+        }
+        else {
+          reference = ReferencePoint.ToArray();
+        }
+        //var reference = ReferencePoint.ToArray();
+        var nondominated = NonDominatedSelect.GetDominatingVectors(lastLayer.Select(x => x.Item2), reference, maximization, false);
+        smetric = nondominated.Any() ? Hypervolume.Calculate(nondominated, reference, maximization) : int.MinValue;
+
+        for (int ii = 0; ii < lastLayer.Count; ++ii) {
           try { // TODO: This can be use for dynamic reference point strategy later. Kaifeng , 02/2020
             // tempHV = Hypervolume.Calculate(indices.Where(idx => idx != ii).Select(idx => lastLayer[idx].Item2), Enumerable.Repeat(11d, NadirPoint.Length).ToArray(), maximization);
-            tempHV = Hypervolume.Calculate(indices.Where(idx => idx != ii).Select(idx => lastLayer[idx].Item2), ReferencePoint.ToArray(), maximization);
+            tempHV = Hypervolume.Calculate(Enumerable.Range(0, lastLayer.Count).Where(idx => idx != ii).Select(idx => lastLayer[idx].Item2), reference, maximization);
           }
           catch {
@@ -537 +552 @@
       //               True:  maximization problem 
       //               False: minimization problem
-      hypervolumes[0, 0] = Hypervolume.Calculate(pf.Select(x => x.Item2), Enumerable.Repeat(11d, NadirPoint.Length).ToArray(), maximization);
+      var reference = Enumerable.Repeat(double.MaxValue, maximization.Length).ToArray();
+      if (ReferencePoint is null) {     // KF, 20200217 -- fix no reference point on real-world applications. If No reference points in Algorithms, use 1.1 \times max objective values as the reference point
+        for (int i = 0; i < reference.Length; i++) {
+          reference[i] = 1.1 * reference[i];
+          if (reference[i] > 10000) {
+            reference[i] = 9999;          // set a upper bound for the reference point
+          }
+        }
+      }
+      else {
+        reference = ReferencePoint.ToArray();
+      }
+      //var reference = ReferencePoint.ToArray() ?? Enumerable.Repeat(double.MaxValue, maximization.Length).ToArray();
+      //var reference = ReferencePoint.ToArray();
+      var nondominated = NonDominatedSelect.GetDominatingVectors(pf.Select(x => x.Item2), reference, maximization, false);
+      hypervolumes[0, 0] = nondominated.Any() ? Hypervolume.Calculate(nondominated, reference, maximization) : int.MinValue;
+
+      //hypervolumes[0, 0] = Hypervolume.Calculate(pf.Select(x => x.Item2), reference, maximization);
       hypervolumes[0, 1] = pf.Count;
       Console.WriteLine("Current HV is", hypervolumes[0, 0]);