From a0d62519c1aa556eebb6141be08d844587b91dc7 Mon Sep 17 00:00:00 2001 From: Heng Li Date: Thu, 15 Feb 2018 14:23:42 -0500 Subject: [PATCH] r710: fixed incorrect inversion coordinate (#112) --- align.c | 9 +++++++-- main.c | 2 +- 2 files changed, 8 insertions(+), 3 deletions(-) diff --git a/align.c b/align.c index 264b23d..e7fecc3 100644 --- a/align.c +++ b/align.c @@ -676,8 +676,13 @@ static int mm_align1_inv(void *km, const mm_mapopt_t *opt, const mm_idx_t *mi, i r_inv->rev = !r1->rev; r_inv->rid = r1->rid; r_inv->div = -1.0f; - r_inv->qs = r1->rev? r2->qe + q_off : r1->qe + q_off; - r_inv->qe = r_inv->qs + ez->max_q + 1; + if (r_inv->rev == 0) { + r_inv->qs = r2->qe + q_off; + r_inv->qe = r_inv->qs + ez->max_q + 1; + } else { + r_inv->qe = r2->qs - q_off; + r_inv->qs = r_inv->qe - (ez->max_q + 1); + } r_inv->rs = r1->re + t_off; r_inv->re = r_inv->rs + ez->max_t + 1; mm_update_extra(r_inv, &qseq[q_off], &tseq[t_off], mat, opt->q, opt->e); diff --git a/main.c b/main.c index 6181e34..6a0f333 100644 --- a/main.c +++ b/main.c @@ -6,7 +6,7 @@ #include "mmpriv.h" #include "getopt.h" -#define MM_VERSION "2.8-r708-dirty" +#define MM_VERSION "2.8-r710-dirty" #ifdef __linux__ #include