Skip to content

Commit

Permalink
commented out the debug prints.
Browse files Browse the repository at this point in the history
  • Loading branch information
Robin-Rounthwaite committed Feb 16, 2024
1 parent e2c3bdb commit 5a8bb78
Show file tree
Hide file tree
Showing 2 changed files with 132 additions and 132 deletions.
6 changes: 3 additions & 3 deletions src/algorithms/0_oo_normalize_snarls.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -783,7 +783,7 @@ void SnarlNormalizer::extract_haplotypes(const SubHandleGraph& snarl, const pair
// //todo: end debug_print

// }
cerr << region.first << " " << region.second << " in extract_haplotypes 1" << endl;
// cerr << region.first << " " << region.second << " in extract_haplotypes 1" << endl;
// extract threads.
// haplotypes is of format:
// 0: a set of all the haplotypes which stretch from source to sink, in string format.
Expand All @@ -792,13 +792,13 @@ void SnarlNormalizer::extract_haplotypes(const SubHandleGraph& snarl, const pair
// 2: a vector of all the handles ever touched by the SnarlSequenceFinder.
SnarlSequenceFinder sequence_finder = SnarlSequenceFinder(_graph, snarl, _gbwt_graph, region.first, region.second, false);

cerr << region.first << " " << region.second << " start find gbwt haps" << endl;
// cerr << region.first << " " << region.second << " start find gbwt haps" << endl;

//todo: here is where the exhaustive path finder would be used, if it was working.
tuple<vector<vector<handle_t>>, vector<vector<handle_t>>, unordered_set<id_t>>
gbwt_haplotypes = sequence_finder.find_gbwt_haps();

cerr << region.first << " " << region.second << " end find gbwt haps" << endl;
// cerr << region.first << " " << region.second << " end find gbwt haps" << endl;

unordered_set<id_t> nodes_in_snarl;
snarl.for_each_handle([&](handle_t handle){
Expand Down
Loading

0 comments on commit 5a8bb78

Please sign in to comment.