Changeset 10208 for branches/HeuristicLab.Analysis.AlgorithmBehavior/PerformanceTests

Timestamp:

12/09/13 16:56:03 (11 years ago)

Author:

ascheibe

Message:

#1886

• added x64 configurations for qhull
• added a wrapper for qhull volume calculation
• removed old volume calculation code

File:

• branches/HeuristicLab.Analysis.AlgorithmBehavior/PerformanceTests/Program.cs (modified) (5 diffs)

branches/HeuristicLab.Analysis.AlgorithmBehavior/PerformanceTests/Program.cs

@@ -1 @@
-using System;
+#region License Information
+/* HeuristicLab
+ * Copyright (C) 2002-2013 Heuristic and Evolutionary Algorithms Laboratory (HEAL)
+ *
+ * This file is part of HeuristicLab.
+ *
+ * HeuristicLab is free software: you can redistribute it and/or modify
+ * it under the terms of the GNU General Public License as published by
+ * the Free Software Foundation, either version 3 of the License, or
+ * (at your option) any later version.
+ *
+ * HeuristicLab is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+ * GNU General Public License for more details.
+ *
+ * You should have received a copy of the GNU General Public License
+ * along with HeuristicLab. If not, see <http://www.gnu.org/licenses/>.
+ */
+#endregion
+
+using System;
+using System.Collections.Generic;
 using System.Diagnostics;
 using System.Linq;
@@ -7 +29 @@
   class Program {
     static void Main(string[] args) {
-      TestVolCalc();
+      // TestSimpleVolumeCalculation();
+      //   TestVolCalc();
+      //TestQHullVolumeCalculation();
+      TestQhull();
     }
 
+    public static void TestQhull() {
+      var points = new List<double[]>();
+      points.Add(new double[] { -0.72045174, 1.44667601, -1.75159125 });
+      points.Add(new double[] { -0.67521503, 0.61127293, -0.33962646 });
+      points.Add(new double[] { 2.02392510, -0.41962562, -1.12211942 });
+      points.Add(new double[] { 0.02647963, -0.24983950, -0.47038916 });
+      points.Add(new double[] { 0.61548803, -0.25039511, 0.97745744 });
+      points.Add(new double[] { -1.65142294, -0.09537136, 1.93154268 });
+      points.Add(new double[] { 1.12618422, 0.49384888, 0.05286691 });
+      points.Add(new double[] { -0.12506782, -2.29161965, -0.09710301 });
+
+      double volume = QhullWrapper.CalculateVolume(points);
+      volume = Math.Round(volume, 6);
+      Console.WriteLine("Volume is: " + volume);
+    }
+
+
     public static void TestVolCalc() {
-      int dim = 4;
-      int size = 20;
+      int dim = 3;
+      int size = 8;
       var data = CreateRandomData(size, dim);
       var dataList = data.ToList();
@@ -27 +69 @@
       //calculate volume of convex hull
       watch.Start();
-      double volume = ConvexHullMeasures.CalculateVolumeNew(convexHull);
+      double volume = ConvexHullMeasures.CalculateVolume(convexHull);
       watch.Stop();
       Console.WriteLine("Volume calculation of convex hull took : " + watch.Elapsed.TotalSeconds + " seconds.");
@@ -35 +77 @@
       //do volume calculation from the whole data set
       watch.Start();
-      volume = ConvexHullMeasures.CalculateVolumeNew(dataList);
+      volume = ConvexHullMeasures.CalculateVolume(dataList);
       watch.Stop();
       Console.WriteLine("Volume calculation of whole data set took : " + watch.Elapsed.TotalSeconds + " seconds.");
@@ -47 +89 @@
       Console.WriteLine("Volume calculation of whole data set with lrs took : " + watch.Elapsed.TotalSeconds + " seconds.");
       Console.WriteLine("Volume of whole data with lrs is: " + volume);
+    }
+
+    private static void TestQHullVolumeCalculation() {
+      const int numPoints = 8;
+      const int dimension = 3;
+
+      var points = new List<double[]>();
+
+      points.Add(new double[] { -0.72045174, 1.44667601, -1.75159125 });
+      points.Add(new double[] { -0.67521503, 0.61127293, -0.33962646 });
+      points.Add(new double[] { 2.02392510, -0.41962562, -1.12211942 });
+      points.Add(new double[] { 0.02647963, -0.24983950, -0.47038916 });
+      points.Add(new double[] { 0.61548803, -0.25039511, 0.97745744 });
+      points.Add(new double[] { -1.65142294, -0.09537136, 1.93154268 });
+      points.Add(new double[] { 1.12618422, 0.49384888, 0.05286691 });
+      points.Add(new double[] { -0.12506782, -2.29161965, -0.09710301 });
+
+      Stopwatch watch = new Stopwatch();
+      //calculate convex hull
+      watch.Start();
+      var convexHull = LPHull.Calculate(points.ToArray());
+      watch.Stop();
+      Console.WriteLine("Convex hull calculation took : " + watch.Elapsed.TotalSeconds + " seconds.");
+
+
+      //calculate volume of convex hull
+      watch.Start();
+      double volume = ConvexHullMeasures.CalculateVolume(convexHull);
+      watch.Stop();
+      Console.WriteLine("Volume calculation of convex hull took : " + watch.Elapsed.TotalSeconds + " seconds.");
+      Console.WriteLine("Volume of convex hull is: " + volume);  //7.960698
+
+
+      //do volume calculation from the whole data set
+      watch.Start();
+      volume = ConvexHullMeasures.CalculateVolume(points);
+      watch.Stop();
+      Console.WriteLine("Volume calculation of whole data set took : " + watch.Elapsed.TotalSeconds + " seconds.");
+      Console.WriteLine("Volume of whole data is: " + volume);
+
+
+      //do volume calculation with lrs
+      watch = new Stopwatch();
+      watch.Start();
+      volume = LiblrsWrapper.CalculateVolume(points);
+      watch.Stop();
+      Console.WriteLine("Volume calculation of whole data set with lrs took : " + watch.Elapsed.TotalSeconds + " seconds.");
+      Console.WriteLine("Volume of whole data with lrs is: " + volume);
+    }
+
+    private static void TestSimpleVolumeCalculation() {
+      const int numPoints = 4;
+      const int dimension = 3;
+
+      double[] num1 = new double[dimension];
+      double[] num2 = new double[dimension];
+      double[] num3 = new double[dimension];
+      double[] num4 = new double[dimension];
+
+      num1[0] = 1;
+      num1[1] = 2;
+      num1[2] = 2;
+
+      num2[0] = 3;
+      num2[1] = 4;
+      num2[2] = 3;
+
+      num3[0] = 3;
+      num3[1] = 6;
+      num3[2] = 3;
+
+      num4[0] = 7;
+      num4[1] = 8;
+      num4[2] = 3;
+
+      var points = new List<double[]>();
+      points.Add(num1);
+      points.Add(num2);
+      points.Add(num3);
+      points.Add(num4);
+
+
+      double volume = ConvexHullMeasures.CalculateSimplexVolume(points);
+
     }