From bed4b3da86b4c2c27ae31ad76904d08404bcfb36 Mon Sep 17 00:00:00 2001
From: Heng Li <lh3@live.co.uk>
Date: Thu, 30 Sep 2010 14:12:02 +0000
Subject: [PATCH] round to the nearest integer

---
 kaln.c | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/kaln.c b/kaln.c
index b52d73a..9fe7d6f 100644
--- a/kaln.c
+++ b/kaln.c
@@ -542,7 +542,7 @@ int ka_prob_glocal(const uint8_t *_ref, int l_ref, const uint8_t *_query, int l_
 		}
 		max /= sum; sum *= s[i]; // if everything works as is expected, sum == 1.0
 		if (state) state[i-1] = max_k;
-		if (q) k = -4.343 * log(1. - max), q[i-1] = k > 100? 99 : k;
+		if (q) k = (int)(-4.343 * log(1. - max) + .499), q[i-1] = k > 100? 99 : k;
 #ifdef _MAIN
 		fprintf(stderr, "(%.10lg,%.10lg) (%d,%d:%d)~%lg\n", pb, sum, i-1, max_k>>2, max_k&3, max); // DEBUG
 #endif
-- 
2.39.5