Changeset 10279

Timestamp:

01/05/14 01:30:28 (10 years ago)

Author:

ascheibe

Message:

#1886 added calculation of the ratio between the convex hull volume and the max. hypercube volume

Location:

branches/HeuristicLab.Analysis.AlgorithmBehavior/HeuristicLab.Analysis.AlgorithmBehavior.Analyzers/3.3

Files:

• ConvexHullMeasures.cs (modified) (1 diff)
• RunCollectionModifiers/RealVectorConvexHullModifier.cs (modified) (3 diffs)

branches/HeuristicLab.Analysis.AlgorithmBehavior/HeuristicLab.Analysis.AlgorithmBehavior.Analyzers/3.3/ConvexHullMeasures.cs

@@ -161 +161 @@
       return maxDist;
     }
+
+    private static List<double[]> GenerateHyperCube(double[] bounds, int dim) {
+      if (bounds.Length != 2)
+        throw new ArgumentException("Only bounds of length 2 are supported!", "bounds");
+
+      int numPoints = (int)Math.Pow(2, dim);
+
+      var hyperCube = new List<double[]>(numPoints);
+      var binaryVectors = new List<string>(numPoints);
+      for (int i = 0; i < numPoints; i++) {
+        string binary = Convert.ToString(i, 2);
+        binaryVectors.Add(binary);
+      }
+
+      for (int i = 0; i < numPoints; i++) {
+        for (int j = 0; j < dim - binaryVectors[i].Length; j++) {
+          //TODO: do this properly
+          binaryVectors[i] = "0" + binaryVectors[i];
+        }
+      }
+
+      foreach (var binaryVector in binaryVectors) {
+        double[] point = new double[dim];
+        for (int i = 0; i < binaryVector.Length; i++) {
+          if (binaryVector[i] == '0') {
+            point[i] = bounds[0];
+          } else {
+            point[i] = bounds[1];
+          }
+        }
+        hyperCube.Add(point);
+      }
+      return hyperCube;
+    }
+
+    public static double CalculateHypercubeVolume(double[] bounds, int dim) {
+      return CalculateVolume(GenerateHyperCube(bounds, dim));
+    }
   }
 }

branches/HeuristicLab.Analysis.AlgorithmBehavior/HeuristicLab.Analysis.AlgorithmBehavior.Analyzers/3.3/RunCollectionModifiers/RealVectorConvexHullModifier.cs

@@ -73 +73 @@
         DataTable centroidMotionTable = new DataTable(CentroidMotionName);
         DoubleValue overallVolume = new DoubleValue();
+        DoubleValue overallVolumeRatio = new DoubleValue();
 
         DataRow dtVolumeRow = new DataRow("Volume");
@@ -112 +113 @@
         if (completeHull != null && completeHull.Any() && completeHull.First().Length < completeHull.Count) {
           overallVolume.Value = ConvexHullMeasures.CalculateVolume(completeHull);
+          var boundsMatrix = run.Parameters["Bounds"] as DoubleMatrix;
+          var bounds = new double[] { boundsMatrix[0, 0], boundsMatrix[0, 1] };
+          double hyperHullVolume = ConvexHullMeasures.CalculateHypercubeVolume(bounds, BestSolutionParameter.Value.Length);
+          overallVolumeRatio.Value = overallVolume.Value / hyperHullVolume;
         } else {
           overallVolume.Value = double.NaN;
+          overallVolumeRatio.Value = double.NaN;
         }
 
@@ -121 +127 @@
 
         run.Results["Overall volume"] = overallVolume;
+        run.Results["Overall volume ratio"] = overallVolumeRatio;
         run.Results["Overall diameter"] = new DoubleValue(ConvexHullMeasures.CalculateMaxDiameter(completeHull));
         run.Results["Convex hull volume"] = volDataTable;