DelauneyMapTest.cs

/*
* This file is part of VECTO.
*
* Copyright © 2012-2016 European Union
*
* Developed by Graz University of Technology,
*              Institute of Internal Combustion Engines and Thermodynamics,
*              Institute of Technical Informatics
*
* VECTO is licensed under the EUPL, Version 1.1 or - as soon they will be approved
* by the European Commission - subsequent versions of the EUPL (the "Licence");
* You may not use VECTO except in compliance with the Licence.
* You may obtain a copy of the Licence at:
*
* https://joinup.ec.europa.eu/community/eupl/og_page/eupl
*
* Unless required by applicable law or agreed to in writing, VECTO
* distributed under the Licence is distributed on an "AS IS" basis,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the Licence for the specific language governing permissions and
* limitations under the Licence.
*
* Authors:
*   Stefan Hausberger, hausberger@ivt.tugraz.at, IVT, Graz University of Technology
*   Christian Kreiner, christian.kreiner@tugraz.at, ITI, Graz University of Technology
*   Michael Krisper, michael.krisper@tugraz.at, ITI, Graz University of Technology
*   Raphael Luz, luz@ivt.tugraz.at, IVT, Graz University of Technology
*   Markus Quaritsch, markus.quaritsch@tugraz.at, IVT, Graz University of Technology
*   Martin Rexeis, rexeis@ivt.tugraz.at, IVT, Graz University of Technology
*/

using System;
using Microsoft.VisualStudio.TestTools.UnitTesting;
using TUGraz.VectoCommon.Exceptions;
using TUGraz.VectoCore.Utils;

namespace TUGraz.VectoCore.Tests.Utils
{
	[TestClass]
	public class DelaunayMapTest
	{
		[TestMethod]
		public void Test_Simple_DelaunayMap()
		{
			var map = new DelaunayMap("TEST");
			map.AddPoint(0, 0, 0);
			map.AddPoint(1, 0, 0);
			map.AddPoint(0, 1, 0);

			map.Triangulate();

			var result = map.Interpolate(0.25, 0.25);

			AssertHelper.AreRelativeEqual(0, result.Value);
		}

		[TestMethod]
		public void Test_DelaunayMapTriangle()
		{
			var map = new DelaunayMap("TEST");
			map.AddPoint(0, 0, 0);
			map.AddPoint(1, 0, 1);
			map.AddPoint(0, 1, 2);

			map.Triangulate();

			// fixed points
			AssertHelper.AreRelativeEqual(0, map.Interpolate(0, 0).Value);
			AssertHelper.AreRelativeEqual(1, map.Interpolate(1, 0).Value);
			AssertHelper.AreRelativeEqual(2, map.Interpolate(0, 1).Value);

			// interpolations
			AssertHelper.AreRelativeEqual(0.5, map.Interpolate(0.5, 0).Value);
			AssertHelper.AreRelativeEqual(1, map.Interpolate(0, 0.5).Value);
			AssertHelper.AreRelativeEqual(1.5, map.Interpolate(0.5, 0.5).Value);

			AssertHelper.AreRelativeEqual(0.25, map.Interpolate(0.25, 0).Value);
			AssertHelper.AreRelativeEqual(0.5, map.Interpolate(0, 0.25).Value);
			AssertHelper.AreRelativeEqual(0.75, map.Interpolate(0.25, 0.25).Value);

			AssertHelper.AreRelativeEqual(0.75, map.Interpolate(0.75, 0).Value);
			AssertHelper.AreRelativeEqual(1.5, map.Interpolate(0, 0.75).Value);

			// extrapolation (should fail)
			Assert.IsNull(map.Interpolate(1, 1));
			Assert.IsNull(map.Interpolate(-1, -1));
			Assert.IsNull(map.Interpolate(1, -1));
			Assert.IsNull(map.Interpolate(-1, 1));
		}

		public void Test_DelaunayMapPlane()
		{
			var map = new DelaunayMap("TEST");
			map.AddPoint(0, 0, 0);
			map.AddPoint(1, 0, 1);
			map.AddPoint(0, 1, 2);
			map.AddPoint(1, 1, 3);

			map.Triangulate();

			// fixed points
			AssertHelper.AreRelativeEqual(0, map.Interpolate(0, 0).Value);
			AssertHelper.AreRelativeEqual(1, map.Interpolate(1, 0).Value);
			AssertHelper.AreRelativeEqual(2, map.Interpolate(0, 1).Value);
			AssertHelper.AreRelativeEqual(3, map.Interpolate(1, 1).Value);

			// interpolations
			AssertHelper.AreRelativeEqual(0.5, map.Interpolate(0.5, 0).Value);
			AssertHelper.AreRelativeEqual(1, map.Interpolate(0, 0.5).Value);
			AssertHelper.AreRelativeEqual(2, map.Interpolate(1, 0.5).Value);
			AssertHelper.AreRelativeEqual(2.5, map.Interpolate(0.5, 1).Value);

			AssertHelper.AreRelativeEqual(1.5, map.Interpolate(0.5, 0.5).Value);

			AssertHelper.AreRelativeEqual(0.75, map.Interpolate(0.25, 0.25).Value);
			AssertHelper.AreRelativeEqual(2.25, map.Interpolate(0.75, 0.75).Value);

			AssertHelper.AreRelativeEqual(1.75, map.Interpolate(0.25, 0.75).Value);
			AssertHelper.AreRelativeEqual(1.25, map.Interpolate(0.75, 0.25).Value);

			// extrapolation (should fail)
			AssertHelper.Exception<VectoException>(() => map.Interpolate(1.5, 0.5), "Interpolation failed.");
			AssertHelper.Exception<VectoException>(() => map.Interpolate(1.5, 1.5), "Interpolation failed.");
			AssertHelper.Exception<VectoException>(() => map.Interpolate(0.5, 1.5), "Interpolation failed.");
			AssertHelper.Exception<VectoException>(() => map.Interpolate(-0.5, 1.5), "Interpolation failed.");
			AssertHelper.Exception<VectoException>(() => map.Interpolate(-0.5, 0.5), "Interpolation failed.");
			AssertHelper.Exception<VectoException>(() => map.Interpolate(-1.5, -1.5), "Interpolation failed.");
			AssertHelper.Exception<VectoException>(() => map.Interpolate(0.5, -0.5), "Interpolation failed.");
			AssertHelper.Exception<VectoException>(() => map.Interpolate(-1.5, -0.5), "Interpolation failed.");
		}

		[TestMethod]
		public void Test_Delaunay_LessThan3Points()
		{
			AssertHelper.Exception<ArgumentException>(() => new DelaunayMap("TEST").Triangulate(),
				"TEST: Triangulation needs at least 3 Points. Got 0 Points.");

			AssertHelper.Exception<ArgumentException>(() => {
				var map1 = new DelaunayMap("TEST");
				map1.AddPoint(1, 0, 0);
				map1.Triangulate();
			}, "TEST: Triangulation needs at least 3 Points. Got 1 Points.");

			AssertHelper.Exception<ArgumentException>(() => {
				var map2 = new DelaunayMap("TEST");
				map2.AddPoint(1, 0, 0);
				map2.AddPoint(0, 1, 0);
				map2.Triangulate();
			}, "TEST: Triangulation needs at least 3 Points. Got 2 Points.");

			var map = new DelaunayMap("TEST");
			map.AddPoint(1, 0, 0);
			map.AddPoint(0, 1, 0);
			map.AddPoint(0, 0, 1);
			map.Triangulate();
		}

		[TestMethod]
		public void Test_Delaunay_DuplicatePoints()
		{
			var map = new DelaunayMap("TEST");
			map.AddPoint(0, 0, 0);
			map.AddPoint(1, 0, 1);
			map.AddPoint(1, 1, 3);
			map.AddPoint(0, 1, 2);
			map.AddPoint(1, 1, 5);

			AssertHelper.Exception<VectoException>(() => { map.Triangulate(); },
				"TEST: Input Data for Delaunay map contains duplicates! \n1 / 1");
		}
	}
}