From f18dadb1c432598fa8bc05367035e08de8fef1a9 Mon Sep 17 00:00:00 2001 From: Heng Li Date: Tue, 19 Mar 2024 23:47:54 -0400 Subject: [PATCH] r1202: halved RMQ gap cost --- lchain.c | 2 +- minimap.h | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/lchain.c b/lchain.c index 0bb1e73..9c76517 100644 --- a/lchain.c +++ b/lchain.c @@ -281,7 +281,7 @@ mm128_t *mg_lchain_rmq(int max_dist, int max_dist_inner, int bw, int max_chn_ski int64_t j; for (j = i0; j < i; ++j) { q = kmp_alloc_rmq(mp); - q->y = (int32_t)a[j].y, q->i = j, q->pri = -(f[j] + 0.5 * chn_pen_gap * ((int32_t)a[j].x + (int32_t)a[j].y)); + q->y = (int32_t)a[j].y, q->i = j, q->pri = -(f[j] + 0.25 * chn_pen_gap * ((int32_t)a[j].x + (int32_t)a[j].y)); krmq_insert(lc_elem, &root, q, 0); if (max_dist_inner > 0) { r = kmp_alloc_rmq(mp); diff --git a/minimap.h b/minimap.h index 91f664a..5a8f913 100644 --- a/minimap.h +++ b/minimap.h @@ -5,7 +5,7 @@ #include #include -#define MM_VERSION "2.27-r1201-dirty" +#define MM_VERSION "2.27-r1202-dirty" #define MM_F_NO_DIAG (0x001LL) // no exact diagonal hit #define MM_F_NO_DUAL (0x002LL) // skip pairs where query name is lexicographically larger than target name