Commit 3c138657 authored by Mathieu Giraud's avatar Mathieu Giraud

tests: test best_overlap_split

Currently, this test is not passing.
parent 5e6bc91f
......@@ -11,6 +11,29 @@
using namespace std;
void testDynprog()
{
// aaaaaaaTaa
// cGcccccccc
string refA = "AAAAAAAAAA";
string refC = "CCCCCCCCCC";
// here seqA and seqC do not actullay overlap, but the goal is to find the good split point
string seqA = "AAAAAAATAA";
string seqC = "CGCCCCCCCC";
int trimA, trimC ;
best_overlap_split(5, seqA, seqC, refA, refC, &trimA, &trimC, VDJ);
cout << trimA << " / " << trimC << endl ;
TAP_TEST(trimA == 3, TEST_FINE_SEGMENT_OVERLAP, "number of trim nucleotides");
TAP_TEST(trimC == 2, TEST_FINE_SEGMENT_OVERLAP, "number of trim nucleotides");
}
void testFineSegment()
{
Fasta seqV("../../germline/IGHV.fa", 2);
......@@ -322,5 +345,6 @@ void testSegment() {
testSegmentationCause();
testExtractor();
testProbability();
testDynprog();
testFineSegment();
}
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment