INNER CODE UNIT · C
f
lh3/ksw2 · ksw2_extd.c:51
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;
max = max > h? max : h;
h -= gapoe;
e -= gape;
e = e > h? e : h;