r333: also subsample high-freq mz during indexing

This commit is contained in:
Heng Li
2021-01-22 23:17:25 -05:00
parent 3b8f73fbab
commit 61befd357a
2 changed files with 2 additions and 2 deletions

View File

@@ -3,7 +3,7 @@
#include <pthread.h>
#define HA_VERSION "0.13-r332-dirty"
#define HA_VERSION "0.13-r333-dirty"
#define VERBOSE 0

View File

@@ -556,7 +556,7 @@ static void worker_for_mz(void *data, long i, int tid)
ha_mz1_v *b = &s->mz_buf[tid];
s->mz_buf[tid].n = 0;
///s->p->opt->w = 51, s->p->opt->k
ha_sketch(s->seq[i], s->len[i], s->p->opt->w, s->p->opt->k, s->n_seq0 + i, s->p->opt->is_HPC, b, s->p->flt_tab, 0, 0, 0);
ha_sketch(s->seq[i], s->len[i], s->p->opt->w, s->p->opt->k, s->n_seq0 + i, s->p->opt->is_HPC, b, s->p->flt_tab, asm_opt.mz_sample_dist, 0, 0);
s->mz[i].n = s->mz[i].m = b->n;
MALLOC(s->mz[i].a, b->n);
memcpy(s->mz[i].a, b->a, b->n * sizeof(ha_mz1_t));