source: branches/HeuristicLab.Analysis.AlgorithmBehavior/AlgorithmBehaviorUnitTests/VolumeCalculationTest.cs @ 10207

using System;
using System.Collections.Generic;
using HeuristicLab.Analysis.AlgorithmBehavior.Analyzers;
using HeuristicLab.Common;
using Microsoft.VisualStudio.TestTools.UnitTesting;

namespace AlgorithmBehaviorUnitTests {
  [TestClass]
  public class VolumeCalculationTest {
    [TestMethod]
    public void TestVolumeCalculation() {
      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 });

      var convexHull = LPHull.Calculate(points.ToArray());
      double volume = ConvexHullMeasures.CalculateVolumeNew(convexHull);
      volume = Math.Round(volume, 6);
      //volume computed with R and the geometry package which uses qhull
      Assert.IsTrue(volume.IsAlmost(7.960698));
      volume = ConvexHullMeasures.CalculateVolumeNew(points);
      volume = Math.Round(volume, 6);
      Assert.IsTrue(volume.IsAlmost(7.960698));
    }
  }
}