Changeset 1ae9aa for src/Fragmentation/Exporters/unittests/SphericalPointDistributionUnitTest.cpp

Timestamp:

Aug 20, 2014, 1:04:08 PM (11 years ago)

Author:

Frederik Heber <heber@…>

Children:

23c605

Parents:

260540

git-author:

Frederik Heber <heber@…> (06/29/14 21:20:49)

git-committer:

Frederik Heber <heber@…> (08/20/14 13:04:08)

Message:

Extracted joinPoints() function to make it accessible to tests.

• added unit test.
• moved some function definitions around and added documentation.

File:

• src/Fragmentation/Exporters/unittests/SphericalPointDistributionUnitTest.cpp (modified) (1 diff)

src/Fragmentation/Exporters/unittests/SphericalPointDistributionUnitTest.cpp

@@ -245 +245 @@
 }
 
+/** UnitTest for joinPoints()
+ */
+void SphericalPointDistributionTest::joinPointsTest()
+{
+  // test with simple configuration of three points
+  {
+    SphericalPointDistribution::Polygon_t newpolygon;
+    newpolygon += Vector(1.,0.,0.);
+    newpolygon += Vector(0.,1.,0.);
+    newpolygon += Vector(0.,0.,1.);
+    SphericalPointDistribution::Polygon_t expectedpolygon = newpolygon;
+    SphericalPointDistribution::IndexTupleList_t matching;
+    matching += SphericalPointDistribution::IndexList_t(1,0);
+    matching += SphericalPointDistribution::IndexList_t(1,1);
+    matching += SphericalPointDistribution::IndexList_t(1,2);
+    SphericalPointDistribution::IndexList_t IndexList =
+        SphericalPointDistribution::joinPoints(
+            newpolygon,
+            SphericalPointDistribution::VectorArray_t(newpolygon.begin(), newpolygon.end()),
+            matching);
+    SphericalPointDistribution::IndexList_t expected;
+    expected += 0,1,2;
+    CPPUNIT_ASSERT_EQUAL( expected, IndexList );
+    CPPUNIT_ASSERT_EQUAL( expectedpolygon, newpolygon );
+  }
+
+  // test with simple configuration of three points, only two are picked
+  {
+    SphericalPointDistribution::Polygon_t newpolygon;
+    newpolygon += Vector(1.,0.,0.);
+    newpolygon += Vector(0.,1.,0.);
+    newpolygon += Vector(0.,0.,1.);
+    SphericalPointDistribution::Polygon_t expectedpolygon = newpolygon;
+    SphericalPointDistribution::IndexTupleList_t matching;
+    matching += SphericalPointDistribution::IndexList_t(1,1);
+    matching += SphericalPointDistribution::IndexList_t(1,2);
+    SphericalPointDistribution::IndexList_t IndexList =
+        SphericalPointDistribution::joinPoints(
+            newpolygon,
+            SphericalPointDistribution::VectorArray_t(newpolygon.begin(), newpolygon.end()),
+            matching);
+    SphericalPointDistribution::IndexList_t expected;
+    expected += 1,2;
+    CPPUNIT_ASSERT_EQUAL( expected, IndexList );
+    CPPUNIT_ASSERT_EQUAL( expectedpolygon, newpolygon );
+  }
+
+  // test with simple configuration of three points, two are joined
+  {
+    SphericalPointDistribution::Polygon_t newpolygon;
+    newpolygon += Vector(1.,0.,0.);
+    newpolygon += Vector(0.,1.,0.);
+    newpolygon += Vector(0.,0.,1.);
+    SphericalPointDistribution::Polygon_t expectedpolygon;
+    expectedpolygon += Vector(1.,0.,0.);
+    expectedpolygon += Vector(0.,M_SQRT1_2,M_SQRT1_2);
+    SphericalPointDistribution::IndexTupleList_t matching;
+    SphericalPointDistribution::IndexList_t joined;
+    joined += 1,2;
+    matching += SphericalPointDistribution::IndexList_t(1,0);
+    matching += joined;
+    SphericalPointDistribution::IndexList_t IndexList =
+        SphericalPointDistribution::joinPoints(
+            newpolygon,
+            SphericalPointDistribution::VectorArray_t(newpolygon.begin(), newpolygon.end()),
+            matching);
+    SphericalPointDistribution::IndexList_t expected;
+    expected += 0,1;
+    CPPUNIT_ASSERT_EQUAL( expected, IndexList );
+    CPPUNIT_ASSERT_EQUAL( expectedpolygon, newpolygon );
+  }
+
+  // test with simple configuration of six points, two are joined, jumbled indices
+  {
+    SphericalPointDistribution::Polygon_t newpolygon;
+    newpolygon += Vector(1.,0.,1.);
+    newpolygon += Vector(1.,0.,0.);
+    newpolygon += Vector(1.,1.,0.);
+    newpolygon += Vector(0.,1.,0.);
+    newpolygon += Vector(0.,0.,1.);
+    newpolygon += Vector(1.,0.,1.);
+    SphericalPointDistribution::Polygon_t expectedpolygon;
+    expectedpolygon += Vector(1.,0.,1.);
+    expectedpolygon += Vector(1.,0.,0.);
+    expectedpolygon += Vector(1.,1.,0.);
+    expectedpolygon += Vector(1.,0.,1.);
+    expectedpolygon += Vector(0.,M_SQRT1_2,M_SQRT1_2); // new centers go last
+    SphericalPointDistribution::IndexTupleList_t matching;
+    SphericalPointDistribution::IndexList_t joined;
+    joined += 3,4;
+    matching += SphericalPointDistribution::IndexList_t(1,1);
+    matching += joined;
+    SphericalPointDistribution::IndexList_t IndexList =
+        SphericalPointDistribution::joinPoints(
+            newpolygon,
+            SphericalPointDistribution::VectorArray_t(newpolygon.begin(), newpolygon.end()),
+            matching);
+    SphericalPointDistribution::IndexList_t expected;
+    expected += 1,4;
+    CPPUNIT_ASSERT_EQUAL( expected, IndexList );
+    CPPUNIT_ASSERT_EQUAL( expectedpolygon, newpolygon );
+  }
+}
 
 /** UnitTest for matchSphericalPointDistributions() with three points