+ 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]);
+ }