From 494f4360d91c015bcef9ed060f372b61bcba5656 Mon Sep 17 00:00:00 2001 From: Evelin Date: Thu, 12 Sep 2024 15:20:32 +0200 Subject: [PATCH] Which bucket --- include/valik/search/search_local.hpp | 7 ++++++- 1 file changed, 6 insertions(+), 1 deletion(-) diff --git a/include/valik/search/search_local.hpp b/include/valik/search/search_local.hpp index 4fe8737d..6d926d0c 100644 --- a/include/valik/search/search_local.hpp +++ b/include/valik/search/search_local.hpp @@ -270,6 +270,7 @@ bool search_local(search_arguments & arguments, search_time_statistics & time_st uint64_t correct_bin{0}; uint64_t incorrect_bin{0}; + // what about the last k-mer in the QGramDir? for (uint64_t i{1}; i < seqan2::length(indexDir(index)); i++) { uint64_t bucketBegin = indexDir(index)[i - 1]; @@ -277,7 +278,11 @@ bool search_local(search_arguments & arguments, search_time_statistics & time_st if (bucketBegin > bucketEnd) { incorrect_bin++; - std::cout << "Error: bucketBegin > bucketEnd\n" << bucketBegin << " > " << bucketEnd << "\n"; + std::cout << "Error: bucketBegin > bucketEnd\n"; + std::cout << "h1=" << i << '\n'; + std::cout << "previous bucket begin\t" << indexDir(index)[i-1] << '\n'; + std::cout << "current bucket begin\t" << indexDir(index)[i] << '\n'; + std::cout << "next bucket begin\t" << indexDir(index)[i+1] << '\n'; } else {