xxay/hifiasm
1
0
Fork 0
mirror of https://github.com/chhylp123/hifiasm.git synced 2026-10-07 07:38:11 +08:00
Code Issues Packages Projects Releases Wiki Activity

r230: the first anchor score was not adjusted

Browse Source
This commit is contained in:
Heng Li
2020-04-12 10:39:42 -04:00
parent 75a89c214d
commit 2c4de3a326
2 changed files with 4 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files
CommandLines.h
+1 -1
View File
@@ -3,7 +3,7 @@
#include <pthread.h> #include <pthread.h>
#define HA_VERSION "0.3.0-dirty-r229" #define HA_VERSION "0.3.0-dirty-r230"
#define VERBOSE 0 #define VERBOSE 0
Hash_Table.cpp
+3 -2
View File
@@ -372,7 +372,8 @@ int32_t ha_chain_lis_core(k_mer_hit *a, int32_t n_a, Chain_Data *dp, int32_t min
m = n_a; m = n_a;
} else m = ha_kmer_hit_lis(n_a, a, b, M); } else m = ha_kmer_hit_lis(n_a, a, b, M);
bw_pen = 1.0 / bw_thres; bw_pen = 1.0 / bw_thres;
dp->score[0] = min_sc, dp->pre[0] = -1, dp->indels[0] = 0, dp->self_length[0] = 0; dp->score[0] = a[b[0]].good? min_sc : min_sc>>1;
dp->pre[0] = -1, dp->indels[0] = 0, dp->self_length[0] = 0;
for (i = 1; i < m; ++i) { for (i = 1; i < m; ++i) {
int32_t j0 = b[i-1], j1 = b[i], score, dg; int32_t j0 = b[i-1], j1 = b[i], score, dg;
int32_t dx = (int32_t)a[j1].offset - (int32_t)a[j0].offset; int32_t dx = (int32_t)a[j1].offset - (int32_t)a[j0].offset;
@@ -431,7 +432,7 @@ void chain_DP(k_mer_hit* a, long long a_n, Chain_Data* dp, overlap_region* resul
pos = a[i].offset; pos = a[i].offset;
self_pos = a[i].self_offset; self_pos = a[i].self_offset;
max_j = -1; max_j = -1;
max_score = min_score; max_score = a[i].good? min_score : min_score>>1;
max_indels = 0; max_indels = 0;
max_self_length = 0; max_self_length = 0;