Changeset 10287 for branches/HeuristicLab.Analysis.AlgorithmBehavior/HeuristicLab.Analysis.AlgorithmBehavior.Analyzers/3.3/RunCollectionModifiers

Timestamp:

01/06/14 14:40:45 (11 years ago)

Author:

ascheibe

Message:

#1886 adapted convex hull modifier to use fractions for big volumes

File:

• branches/HeuristicLab.Analysis.AlgorithmBehavior/HeuristicLab.Analysis.AlgorithmBehavior.Analyzers/3.3/RunCollectionModifiers/RealVectorConvexHullModifier.cs (modified) (4 diffs)

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

@@ -32 +32 @@
 using HeuristicLab.Persistence.Default.CompositeSerializers.Storable;
 using HeuristicLab.Problems.TestFunctions;
+using Mehroz;
 
 namespace HeuristicLab.Analysis.AlgorithmBehavior.Analyzers {
@@ -93 +94 @@
         var boundsMatrix = run.Parameters["Bounds"] as DoubleMatrix;
         var bounds = new double[] { boundsMatrix[0, 0], boundsMatrix[0, 1] };
-        double hyperCubeVolume = ConvexHullMeasures.CalculateHypercubeVolume(bounds, BestSolutionParameter.Value.Length);
+        Fraction hyperCubeVolumeBig = ConvexHullMeasures.CalculateHypercubeVolumeBig(bounds, BestSolutionParameter.Value.Length);
+        double hyperCubeVolume = double.NaN;
+
+        if (hyperCubeVolumeBig < 10000000) {
+          hyperCubeVolume = ConvexHullMeasures.CalculateHypercubeVolume(bounds, BestSolutionParameter.Value.Length);
+        }
 
         List<double[]> completeHull = null;
@@ -118 +124 @@
 
           if (completeHull != null && completeHull.Any() && completeHull.First().Length < completeHull.Count) {
-            double completeHullVolume = ConvexHullMeasures.CalculateVolume(completeHull);
-            double ratio = completeHullVolume / hyperCubeVolume;
-            dtIncrementalHullVolume.Values.Add(completeHullVolume);
-            dtIncrementalHullVolumeRatio.Values.Add(ratio);
+            if (double.IsNaN(hyperCubeVolume)) {
+              Fraction completeHullVolume = ConvexHullMeasures.CalculateVolumeBig(completeHull);
+              Fraction ratio = completeHullVolume / hyperCubeVolumeBig;
+              dtIncrementalHullVolume.Values.Add(completeHullVolume.ToDouble());
+              dtIncrementalHullVolumeRatio.Values.Add(ratio.ToDouble());
+            } else {
+              double completeHullVolume = ConvexHullMeasures.CalculateVolume(completeHull);
+              double ratio = completeHullVolume / hyperCubeVolume;
+              dtIncrementalHullVolume.Values.Add(completeHullVolume);
+              dtIncrementalHullVolumeRatio.Values.Add(ratio);
+            }
           }
 
@@ -128 +141 @@
 
         if (completeHull != null && completeHull.Any() && completeHull.First().Length < completeHull.Count) {
-          overallVolume.Value = ConvexHullMeasures.CalculateVolume(completeHull);
-          overallVolumeRatio.Value = overallVolume.Value / hyperCubeVolume;
+          if (!double.IsNaN(hyperCubeVolume)) {
+            overallVolume.Value = ConvexHullMeasures.CalculateVolume(completeHull);
+            overallVolumeRatio.Value = overallVolume.Value / hyperCubeVolume;
+          } else {
+            Fraction v = ConvexHullMeasures.CalculateVolumeBig(completeHull);
+            overallVolume.Value = v.ToDouble();
+            overallVolumeRatio.Value = (double)(v / hyperCubeVolumeBig);
+          }
         } else {
           overallVolume.Value = double.NaN;