From 2763fa6d4cc3a7befb4d7f513ccca81088755eff Mon Sep 17 00:00:00 2001 From: Heng Li Date: Wed, 20 Feb 2013 20:14:04 -0500 Subject: [PATCH] bugfix: point address changes --- ksw.c | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/ksw.c b/ksw.c index 8d741a6..742fec9 100644 --- a/ksw.c +++ b/ksw.c @@ -447,7 +447,7 @@ static inline uint32_t *push_cigar(int *n_cigar, int *m_cigar, uint32_t *cigar, if (*n_cigar == 0 || op != (cigar[(*n_cigar) - 1]&0xf)) { if (*n_cigar == *m_cigar) { *m_cigar = *m_cigar? (*m_cigar)<<1 : 4; - cigar = realloc(cigar, (*m_cigar) << 4); + cigar = realloc(cigar, (*m_cigar) << 2); } cigar[(*n_cigar)++] = len<<4 | op; } else cigar[(*n_cigar)-1] += len<<4; @@ -520,8 +520,8 @@ int ksw_global(int qlen, const uint8_t *query, int tlen, const uint8_t *target, else if (which == 1) cigar = push_cigar(&n_cigar, &m_cigar, cigar, 2, 1), --i; else cigar = push_cigar(&n_cigar, &m_cigar, cigar, 1, 1), --k; } - if (i >= 0) push_cigar(&n_cigar, &m_cigar, cigar, 2, i + 1); - if (k >= 0) push_cigar(&n_cigar, &m_cigar, cigar, 1, k + 1); + if (i >= 0) cigar = push_cigar(&n_cigar, &m_cigar, cigar, 2, i + 1); + if (k >= 0) cigar = push_cigar(&n_cigar, &m_cigar, cigar, 1, k + 1); for (i = 0; i < n_cigar>>1; ++i) // reverse CIGAR tmp = cigar[i], cigar[i] = cigar[n_cigar-1-i], cigar[n_cigar-1-i] = tmp; *n_cigar_ = n_cigar, *cigar_ = cigar;