r1202: halved RMQ gap cost

lh3 · Mar 20, 2024 · f18dadb · f18dadb
1 parent a83b8fe
commit f18dadb
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git 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
@@ -5,7 +5,7 @@
 #include <stdio.h>
 #include <sys/types.h>
 
-#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