Compare commits

...

2 Commits

Author SHA1 Message Date
zzh b38214edd9 同上 2026-01-06 14:40:04 +08:00
zzh 90f6bc5cc0 测试过了5.1GB的数据,串行和多线程结果一致 2026-01-06 14:39:48 +08:00
3 changed files with 6 additions and 6 deletions

View File

@ -5,6 +5,6 @@ set(CMAKE_CXX_STANDARD_REQUIRED ON)
# set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -pthread")
# set(CMAKE_BUILD_TYPE Debug)
# set(CMAKE_BUILD_TYPE Release)
add_definitions(-DSHOW_PERF)
#add_definitions(-DSHOW_PERF)
# add_definitions(-DSHOW_PERF=1)
add_subdirectory(src)

View File

@ -187,11 +187,11 @@ int SerialApplyBQSR(AuxVar &aux) {
recalibratedQuals[offset] = quantizedQualityScore;
}
fprintf(gf[4], "%s %d %ld ", bam_get_qname(sd.bw->b), sd.bw->b->core.flag, sd.rid);
for (size_t si = 0; si < sd.read_len; ++si) {
fprintf(gf[4], "%d ", (int)recalibratedQuals[si]);
}
fprintf(gf[4], "\n");
// fprintf(gf[4], "%s %d %ld ", bam_get_qname(sd.bw->b), sd.bw->b->core.flag, sd.rid);
// for (size_t si = 0; si < sd.read_len; ++si) {
// fprintf(gf[4], "%d ", (int)recalibratedQuals[si]);
// }
// fprintf(gf[4], "\n");
if (sam_write1(nsgv::gOutBamFp, nsgv::gOutBamHeader, bw->b) < 0) {
spdlog::error("failed writing sam record to \"{}\"", nsgv::gBqsrArg.OUTPUT_FILE.c_str());

Binary file not shown.