INNER CODE UNIT · C
st
lh3/ksw2 · ksw2_extd.c:47
st = i > w? i - w : 0;
en = i + w < qlen - 1? i + w : qlen - 1;
tmp = -(gapoe + gape * i) > -(gapoe2 + gape2 * i)? -(gapoe + gape * i) : -(gapoe2 + gape2 * i);
h1 = st > 0? KSW_NEG_INF : tmp;
f = st > 0? KSW_NEG_INF : tmp - gapoe;
f2 = st > 0? KSW_NEG_INF : tmp - gapoe2;
if (!with_cigar) {
for (j = st; j <= en; ++j) {
eh_t *p = &eh[j];
int32_t h = p->h, h2, e = p->e, e2 = p->e2;
p->h = h1;
h += q[j];
h = h >= e? h : e;
h = h >= f? h : f;
h = h >= e2? h : e2;
h = h >= f2? h : f2;
h1 = h;
max_j = max > h? max_j : j;