+ void testInversionDistance (){
+ GeneOrder go(_validPerm.begin(),_validPerm.end());
+ size_t d = inversionDistance(go);
+ CPPUNIT_ASSERT_EQUAL((size_t)0u,d);
+
+ GeneOrder go2(_validPerm2.begin(),_validPerm2.end());
+ d = inversionDistance(go2);
+ CPPUNIT_ASSERT_EQUAL((size_t)5u,d);
+
+ GeneOrder go3(_validPerm3.begin(),_validPerm3.end());
+ d = inversionDistance(go3);
+ CPPUNIT_ASSERT_EQUAL((size_t)13u,d);
+ }
+
+ void testFindComponents (){
+ GeneOrder go(_validPerm.begin(),_validPerm.end());
+ vector<Component> v = findComponents(go);
+ CPPUNIT_ASSERT_EQUAL((size_t)4u,v.size());
+ Component go10(0,1,0,0,1);
+ Component go11(1,2,0,1,2);
+ Component go12(2,3,0,2,3);
+ Component go13(3,4,0,3,4);
+ CPPUNIT_ASSERT(go10 == v[0]);
+ CPPUNIT_ASSERT(go11 == v[1]);
+ CPPUNIT_ASSERT(go12 == v[2]);
+ CPPUNIT_ASSERT(go13 == v[3]);
+
+ GeneOrder go2(_validPerm3.begin(),_validPerm3.end());
+ v = findComponents(go2);
+ CPPUNIT_ASSERT_EQUAL((size_t)6u,v.size());
+ Component go20(1,2,0,2,3);
+ Component go21(0,4,0,0,4);
+ Component go22(4,7,1,4,7);
+ Component go23(-15,-12,-1,8,11);
+ Component go24(-12,-9,-1,11,14);
+ Component go25(7,16,0,7,16);
+ CPPUNIT_ASSERT(go20 == v[0]);
+ CPPUNIT_ASSERT(go21 == v[1]);
+ CPPUNIT_ASSERT(go22 == v[2]);
+ CPPUNIT_ASSERT(go23 == v[3]);
+ CPPUNIT_ASSERT(go24 == v[4]);
+ CPPUNIT_ASSERT(go25 == v[5]);
+ }
+