changed use of IntervalTree to accomodate changes in method names

andersen-lab · Jan 4, 2021 · 9c69970 · 9c69970
1 parent 272544d
commit 9c69970
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/src/trim_primer_quality.cpp b/src/trim_primer_quality.cpp
@@ -331,7 +331,7 @@ int get_bigger_primer(std::vector<primer> primers){
   return max_primer_len;
 }
 
-// check if read overlaps with any of the amplicons
+// check if read is enveloped by any of the amplicons
 bool amplicon_filter(IntervalTree amplicons, bam1_t* r){
   Interval fragment_coords = Interval(0, 1);
   if(r->core.isize > 0){
@@ -342,7 +342,7 @@ bool amplicon_filter(IntervalTree amplicons, bam1_t* r){
     fragment_coords.high = bam_endpos(r);
   }
   // debugging
-  bool amplicon_flag = amplicons.overlapSearch(fragment_coords);
+  bool amplicon_flag = amplicons.envelopSearch(fragment_coords);
   return amplicon_flag;
 }
 

diff --git a/tests/test_interval_tree.cpp b/tests/test_interval_tree.cpp
@@ -7,7 +7,7 @@ int test_itree_overlap(IntervalTree tree, Interval queries[], int num_tests, boo
     int result = 0;
     for (int i = 0; i < num_tests; i++)
     {
-        result = tree.overlapSearch(queries[i]);
+        result = tree.envelopSearch(queries[i]);
         if (result != expected[i])
         {
             std::cout << "Interval Tree overlap behavior incorrect for interval " << queries[i].low << ":" << queries[i].high