INNER CODE UNIT · C

st2

lh3/bedtk · bedtk.c:381

		st2 = st1 - win, en2 = en1 + win;
		if (st2 < 0) st2 = 0;
		n_b = cr_overlap(cr, ctg, st2, en2, &b, &m_b);
		if (test_con) {
			for (i = 0, sat = 0; i < n_b; ++i) {
				cr_intv_t *r = &cr->r[b[i]];
				if (cr_st(r) <= st2 && en2 <= cr_en(r)) {
					sat = 1;
					break;
				}
			}
		} else if (min_frac > 0.0) {
			int64_t cov, cnt, depth;
			cov = cal_cov(cr, n_b, b, st1, en1, &cnt, &depth);
			sat = (cov >= (en1 - st1) * min_frac);
		} else sat = (n_b > 0);
		out.l = 0;
		if ((sat && !non_sat) || (!sat && non_sat)) {

View source record →