+ 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<pair<size_t,size_t> > v = findIntervals(go);
+ CPPUNIT_ASSERT_EQUAL(4ul,v.size());
+ Interval go10(1,1);
+ Interval go12(3,3);
+ CPPUNIT_ASSERT(go10 == v[0]);
+ CPPUNIT_ASSERT(go12 == v[2]);
+
+ GeneOrder go2(_validPerm3.begin(),_validPerm3.end());
+ v = findIntervals(go2);
+ CPPUNIT_ASSERT_EQUAL(16ul,v.size());
+ Interval go20(1,2);
+ Interval go22(4,2);
+ CPPUNIT_ASSERT(go20 == v[0]);
+ CPPUNIT_ASSERT(go22 == v[2]);
+ }
+ void testFindIntervalsAtPoints (){
+ GeneOrder go(_validPerm.begin(),_validPerm.end());
+ vector<Interval> v = findIntervals(go);
+ v = findIntervalsAtPoints(v);
+ CPPUNIT_ASSERT_EQUAL(4ul,v.size());
+ Interval go10(0,0);
+ Interval go12(2,2);
+ CPPUNIT_ASSERT(go10 == v[0]);
+ CPPUNIT_ASSERT(go12 == v[2]);
+
+ GeneOrder go2(_validPerm3.begin(),_validPerm3.end());
+ v = findIntervals(go2);
+ v = findIntervalsAtPoints(v);
+ CPPUNIT_ASSERT_EQUAL(16ul,v.size());
+ Interval go20(0,3);
+ Interval go22(1,1);
+ CPPUNIT_ASSERT(go20 == v[0]);
+ CPPUNIT_ASSERT(go22 == v[2]);
+ }
+ 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);
+ }
+