Changeset 10286

Timestamp:

01/06/14 02:36:31 (11 years ago)

Author:

ascheibe

Message:

#1886 added support for calculating big volumes

Location:

branches/HeuristicLab.Analysis.AlgorithmBehavior

Files:

• AlgorithmBehaviorUnitTests/AlgorithmBehaviorUnitTests.csproj (modified) (1 diff)
• AlgorithmBehaviorUnitTests/BigVolumeTest.cs (added)
• HeuristicLab.Analysis.AlgorithmBehavior.Analyzers/3.3/ConvexHullMeasures.cs (modified) (4 diffs)
• HeuristicLab.Analysis.AlgorithmBehavior.Analyzers/3.3/FractionClass.cs (added)
• HeuristicLab.Analysis.AlgorithmBehavior.Analyzers/3.3/HeuristicLab.Analysis.AlgorithmBehavior.Analyzers-3.3.csproj (modified) (2 diffs)

branches/HeuristicLab.Analysis.AlgorithmBehavior/AlgorithmBehaviorUnitTests/AlgorithmBehaviorUnitTests.csproj

@@ -98 +98 @@
   </Choose>
   <ItemGroup>
+    <Compile Include="BigVolumeTest.cs" />
     <Compile Include="LPConvexHullSortingTest.cs" />
     <Compile Include="LPConvexHullTest.cs" />

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

@@ -23 +23 @@
 using System.Collections.Generic;
 using System.Linq;
+using Mehroz;
 using MIConvexHull;
 
@@ -44 +45 @@
     }
 
+    public static Fraction CalculateSimplexVolumeBig(List<double[]> simplex) {
+      int dim = simplex.First().Length;
+      double[,] diffs = new double[dim, dim];
+
+      for (int i = 1; i < simplex.Count; i++) {
+        for (int j = 0; j < dim; j++) {
+          diffs[j, i - 1] = simplex[i][j] - simplex[0][j];
+        }
+      }
+
+      //calculate determinante from lu decomposition
+      int[] pivots = null;
+      alglib.rmatrixlu(ref diffs, dim, dim, out pivots);
+
+      Fraction det = 1;
+      for (int i = 0; i < dim; i++) {
+        det = det * diffs[i, i];
+      }
+
+      if (det < 0)
+        det = -1 * det;
+
+      Fraction result = det / dim.Fact();
+      return result;
+    }
+
     //calculates the volume of a convex d-polytope
     //decomposition based on delaunay triangulation
@@ -64 +91 @@
       foreach (var simplex in simplices) {
         volume += CalculateSimplexVolume(simplex.ToList());
+      }
+
+      return volume;
+    }
+
+    public static Fraction CalculateVolumeBig(List<double[]> convexHull) {
+      Fraction volume = new Fraction(0.0);
+
+      int dim = convexHull.First().Length;
+      if (dim > convexHull.Count)
+        throw new ArgumentException("Nr. of points for volume calculation must be greater than dimension", "convexHull");
+
+      ITriangulation<DefaultVertex, DefaultTriangulationCell<DefaultVertex>> triangulation = null;
+      try {
+        //Under certain circumstances MIConvexHull is not able to calculate the triangulation
+        triangulation = Triangulation.CreateDelaunay(ConvertToVertex(convexHull));
+      }
+      catch {
+        return Fraction.NaN;
+      }
+      var simplices = ConvertFromTriangulation(triangulation);
+      foreach (var simplex in simplices) {
+        volume += CalculateSimplexVolumeBig(simplex.ToList());
       }
 
@@ -197 +247 @@
 
     public static double CalculateHypercubeVolume(double[] bounds, int dim) {
+      //this is not very smart...
       return CalculateVolume(GenerateHyperCube(bounds, dim));
+    }
+
+    public static Fraction CalculateHypercubeVolumeBig(double[] bounds, int dim) {
+      double diff = bounds[1] - bounds[0];
+
+      Fraction volume = diff;
+      for (int i = 0; i < dim; i++) {
+        volume *= diff;
+      }
+      return volume;
     }
   }

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

@@ -106 +106 @@
     <Reference Include="System.Core" />
     <Reference Include="System.Drawing" />
+    <Reference Include="System.Numerics" />
     <Reference Include="System.Xml.Linq" />
     <Reference Include="System.Data.DataSetExtensions" />
@@ -113 +114 @@
   </ItemGroup>
   <ItemGroup>
+    <Compile Include="FractionClass.cs" />
     <Compile Include="QhullWrapper.cs" />
     <Compile Include="DoubleArrayExtensions.cs" />