INNER CODE UNIT · C
b
lh3/ksw2 · cli.c:40
b = b > 0? -b : b;
for (i = 0; i < m - 1; ++i) {
for (j = 0; j < m - 1; ++j)
mat[i * m + j] = i == j? a : b;
mat[i * m + m - 1] = 0;
}
for (j = 0; j < m; ++j)
mat[(m - 1) * m + j] = 0;
}
static void global_aln(const char *algo, void *km, const char *qseq_, const char *tseq_, int8_t m, const int8_t *mat, int8_t q, int8_t e, int8_t q2, int8_t e2,
int w, int zdrop, int flag, ksw_extz_t *ez)
{
int i, qlen, tlen;
uint8_t *qseq, *tseq;
ez->max_q = ez->max_t = ez->mqe_t = ez->mte_q = -1;
ez->max = 0, ez->mqe = ez->mte = KSW_NEG_INF;
ez->n_cigar = 0;