countCycles implemented and passess test

[germs.git] / src / test / genealgorithmstest.cpp

diff --git a/src/test/genealgorithmstest.cpp b/src/test/genealgorithmstest.cpp
index f7bd2fc13bf3d0c7551536f03db2c348c442b127..3f57e0fa3e47a246c02a557e7b0fbeea6996bc98 100644 (file)
--- a/src/test/genealgorithmstest.cpp
+++ b/src/test/genealgorithmstest.cpp
@@ -21,6 +21,10 @@ class TESTNAME : public CPPUNIT_NS::TestFixture
 {
   CPPUNIT_TEST_SUITE( TESTNAME );
   CPPUNIT_TEST( testRobinsonSchensted );
+  CPPUNIT_TEST( testLongestSequences );
+  CPPUNIT_TEST( testFindIntervals );
+  CPPUNIT_TEST( testFindIntervalsAtPoints );
+  CPPUNIT_TEST( testCountCycles );
   CPPUNIT_TEST_SUITE_END();
 
 protected:
@@ -65,6 +69,63 @@ protected:
                CPPUNIT_ASSERT(equal(v[1].begin(),v[1].end(),second));
        }
 
+       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);
+       }
+
 };
 
 CPPUNIT_TEST_SUITE_REGISTRATION( TESTNAME );