+ void testLongestSequences (){
+ GeneOrder go(_validPerm.begin(),_validPerm.end());
+ pair<int,int> p = longestSequences(go);
+ CPPUNIT_ASSERT_EQUAL(5,p.first);
+ CPPUNIT_ASSERT_EQUAL(1,p.second);
+
+ GeneOrder go2(_validPerm2.begin(),_validPerm2.end());
+ p = longestSequences(go2);
+ CPPUNIT_ASSERT_EQUAL(7,p.first);
+ CPPUNIT_ASSERT_EQUAL(2,p.second);
+ }
+ void testFindIntervals (){
+ GeneOrder go(_validPerm.begin(),_validPerm.end());
+ vector<Interval> v = findIntervals(go);
+ CPPUNIT_ASSERT_EQUAL((size_t)4u,v.size());
+ Interval go10(1,1);
+ Interval go12(3,3);
+ CPPUNIT_ASSERT(go10 == v[0]);
+ CPPUNIT_ASSERT(go12 == v[2]);
+
+ GeneOrder go2(_validPerm2.begin(),_validPerm2.end());
+ v = findIntervals(go2);
+ CPPUNIT_ASSERT_EQUAL((size_t)9u,v.size());
+ Interval go20(1,3,true);
+ Interval go21(2,2);
+ Interval go22(1,4,true);
+ Interval go23(5,3);
+ Interval go25(6,7);
+ Interval go26(8,8);
+ Interval go27(9,7,true);
+ CPPUNIT_ASSERT(go20 == v[0]);
+ CPPUNIT_ASSERT(go21 == v[1]);
+ CPPUNIT_ASSERT(go22 == v[2]);
+ CPPUNIT_ASSERT(go23 == v[3]);
+ CPPUNIT_ASSERT(go25 == v[5]);
+ CPPUNIT_ASSERT(go26 == v[6]);
+ CPPUNIT_ASSERT(go27 == v[7]);
+
+ /*GeneOrder go2(_validPerm3.begin(),_validPerm3.end());
+ v = findIntervals(go2);
+ CPPUNIT_ASSERT_EQUAL(16ul,v.size());
+ Interval go20(1,2,true);
+ Interval go22(4,2);
+ Interval go215(8,16);
+ CPPUNIT_ASSERT(go20 == v[0]);
+ CPPUNIT_ASSERT(go22 == v[2]);
+ CPPUNIT_ASSERT(go215 == v[15]);*/
+ }
+ void testFindIntervalsAtPoints (){
+ GeneOrder go(_validPerm.begin(),_validPerm.end());
+ vector<Interval> v = findIntervals(go);
+ v = findIntervalsAtPoints(v);
+ CPPUNIT_ASSERT_EQUAL((size_t)5u,v.size());
+ Interval go10(0,0);
+ Interval go12(2,2);
+ CPPUNIT_ASSERT(go10 == v[1]);
+ CPPUNIT_ASSERT(go12 == v[3]);
+
+ GeneOrder go2(_validPerm3.begin(),_validPerm3.end());
+ v = findIntervals(go2);
+ v = findIntervalsAtPoints(v);
+ CPPUNIT_ASSERT_EQUAL((size_t)17u,v.size());
+ Interval go20(0,3);
+ Interval go22(1,1);
+ CPPUNIT_ASSERT(go20 == v[1]);
+ CPPUNIT_ASSERT(go22 == v[3]);
+ }
+ void testCountCycles (){
+ GeneOrder go(_validPerm.begin(),_validPerm.end());
+ int c = countCycles(go);
+ CPPUNIT_ASSERT_EQUAL(4,c);
+
+ GeneOrder go2(_validPerm3.begin(),_validPerm3.end());
+ c = countCycles(go2);
+ CPPUNIT_ASSERT_EQUAL(6,c);
+ }
+
+ 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]);
+ }
+