@@ -632,9 +632,13 @@ bool RGFACover::add_interval(vector<pair<step_handle_t, step_handle_t>>& thread_
632632#pragma omp critical(cerr)
633633 cerr << " adding interval " << graph->get_path_name (graph->get_path_handle_of_step (new_interval.first ))
634634 << (graph->get_is_reverse (graph->get_handle_of_step (new_interval.first )) ? " <" : " >" )
635- << " :" << graph->get_id (graph->get_handle_of_step (new_interval.first ))
636- << " -" << (graph->get_is_reverse (graph->get_handle_of_step (new_interval.second )) ? " <" : " >" )
637- << graph->get_id (graph->get_handle_of_step (new_interval.second )) << endl;
635+ << " :" << graph->get_id (graph->get_handle_of_step (new_interval.first ));
636+ if (new_interval.second == graph->path_end (graph->get_path_handle_of_step (new_interval.first ))) {
637+ cerr << " PATH_END" << endl;
638+ } else {
639+ cerr << " -" << (graph->get_is_reverse (graph->get_handle_of_step (new_interval.second )) ? " <" : " >" )
640+ << graph->get_id (graph->get_handle_of_step (new_interval.second )) << endl;
641+ }
638642#endif
639643 bool merged = false ;
640644 int64_t merged_interval_idx = -1 ;
@@ -655,9 +659,13 @@ bool RGFACover::add_interval(vector<pair<step_handle_t, step_handle_t>>& thread_
655659#pragma omp critical(cerr)
656660 cerr << " prev interval found" << graph->get_path_name (graph->get_path_handle_of_step (prev_interval.first ))
657661 << " :" << (graph->get_is_reverse (graph->get_handle_of_step (prev_interval.first )) ? " <" : " >" )
658- << graph->get_id (graph->get_handle_of_step (prev_interval.first ))
659- << " -" << (graph->get_is_reverse (graph->get_handle_of_step (prev_interval.second )) ? " <" : " >" )
660- << graph->get_id (graph->get_handle_of_step (prev_interval.second )) << endl;
662+ << graph->get_id (graph->get_handle_of_step (prev_interval.first ));
663+ if (prev_interval.second == graph->path_end (graph->get_path_handle_of_step (prev_interval.first ))) {
664+ cerr << " PATH_END" << endl;
665+ } else {
666+ cerr << " -" << (graph->get_is_reverse (graph->get_handle_of_step (prev_interval.second )) ? " <" : " >" )
667+ << graph->get_id (graph->get_handle_of_step (prev_interval.second )) << endl;
668+ }
661669#endif
662670 prev_interval.second = new_interval.second ;
663671 merged = true ;
0 commit comments