+ 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);
+ 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(5ul,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(17ul,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);
+ }