// File: /workspace/github_cpp_makefile_star2_O2/chhylp123[P]hifiasm/inter.cpp
// Function: worker_for_ul_recorrect_alignment(void*, long, int)
static void worker_for_ul_recorrect_alignment(void *data, long i, int tid) // callback for kt_for()
{
    utepdat_t *s = (utepdat_t*)data;
    ha_ovec_buf_t *b = s->hab[tid];
    glchain_t *bl = &(s->ll[tid]);
    int64_t /**rid = s->id+i,**/ winLen = MIN((((double)THRESHOLD_MAX_SIZE)/s->opt->diff_ec_ul), WINDOW), is_correct;
    // uint64_t align = 0;
    int fully_cov, abnormal;
	// if(UL_INF.a[s->id+i].rlen != s->len[i]) {
	// 	fprintf(stderr, "[M::%s] rid:%ld, s->len:%lu, UL_INF->rlen:%u\n", __func__, s->id+i, s->len[i], UL_INF.a[s->id+i].rlen);
	// }
	is_correct = ck_ul_alignment(&(UL_INF.a[s->id+i]));
	if(is_correct) {
		assert((UL_INF.a[s->id+i].rlen == s->len[i]) && (!s->seq[i]));
		return;
	}
	assert(UL_INF.a[s->id+i].rlen&((uint32_t)(0x80000000)));
	UL_INF.a[s->id+i].rlen<<=1; UL_INF.a[s->id+i].rlen>>=1;
	assert((UL_INF.a[s->id+i].rlen == s->len[i]));
	
    // void *km = s->buf?(s->buf[tid]?s->buf[tid]->km:NULL):NULL;
    // if(s->id+i!=41927 && s->id+i!=47072 && s->id+i!=67641 && s->id+i!=90305 && s->id+i!=698342 && s->id+i!=329421) {
	// 	return;
	// }
	// if(s->id+i!=41927) return;

    // fprintf(stderr, "\n[M::%s] rid:%ld, len:%lu\n", __func__, s->id+i, s->len[i]);
    // if (memcmp(UL_INF.nid.a[s->id+i].a, "d0aab024-b3a7-40fb-83cc-22c3d6d951f8", UL_INF.nid.a[s->id+i].n-1)) return;
    // fprintf(stderr, "[M::%s::] ==> len: %lu\n", __func__, s->len[i]);
    ha_get_ul_candidates_interface(b->abl, i, s->seq[i], s->len[i], s->opt->w, s->opt->k, s->uu, &b->olist, &b->olist_hp, &b->clist, s->opt->bw_thres, 
        s->opt->max_n_chain, 1, NULL/**&(b->k_flag)**/, &b->r_buf, &(b->tmp_region), NULL, &(b->sp), 1, NULL);
    
    clear_Cigar_record(&b->cigar1);
    clear_Round2_alignment(&b->round2);
    // return;
    // b->num_correct_base += overlap_statistics(&b->olist, NULL, 0);   

    b->self_read.seq = s->seq[i]; b->self_read.length = s->len[i]; b->self_read.size = 0;
    correct_ul_overlap(&b->olist, s->uu, &b->self_read, &b->correct, &b->ovlp_read, &b->POA_Graph, &b->DAGCon,
            &b->cigar1, &b->hap, &b->round2, &b->r_buf, &(b->tmp_region.w_list), 0, 1, &fully_cov, &abnormal, s->opt->diff_ec_ul, winLen, NULL);

    // uint64_t k;
    // for (k = 0; k < b->olist.length; k++) {
    //  if(b->olist.list[k].is_match == 1) b->num_correct_base += b->olist.list[k].x_pos_e+1-b->olist.list[k].x_pos_s;
    //  if(b->olist.list[k].is_match == 2) b->num_recorrect_base += b->olist.list[k].x_pos_e+1-b->olist.list[k].x_pos_s;
    // }


    // gl_chain_refine(&b->olist, &b->correct, &b->hap, bl, s->uu, s->opt->diff_ec_ul, winLen, s->len[i], km);
    gl_chain_refine_advance_combine(s->buf[tid], &(UL_INF.a[s->id+i]), &b->olist, &b->correct, &b->hap, &(s->sps[tid]), bl, &(s->gdp[tid]), s->uu, s->opt->diff_ec_ul, winLen, s->len[i], s->uopt, s->id+i, tid, NULL);
    // return;
    // b->num_read_base += b->self_read.length;
    // b->num_correct_base += b->correct.corrected_base;
    // b->num_recorrect_base += b->round2.dumy.corrected_base;
    memset(&b->self_read, 0, sizeof(b->self_read));



	is_correct = ck_ul_alignment(&(UL_INF.a[s->id+i]));
	if(is_correct) b->num_correct_base++;
    s->hab[tid]->num_read_base++;

	// fprintf(stderr, "[M::%s] rid:%ld, dd:%u\n", __func__, s->id+i, UL_INF.a[s->id+i].dd);
    // int64_t mem[6], mem_hab[6];
    // if(get_utepdat_t_mem_tid(s, tid, mem, mem_hab)>((int64_t)5*(int64_t)1073741824)) {
    //     fprintf(stderr, "[M::%s::tid->%d::rid->%ld] buffer[0]: %.3fGB(%.3fGB::%.3fGB::%.3fGB::%.3fGB::%.3fGB), buffer[1]: %.3fGB, buffer[2]: %.3fGB, buffer[3]: %.3fGB, buffer[4]: %.3fGB, buffer[5]: %.3fGB\n", 
    //             __func__, tid, i, mem[0]/1073741824.0, 
	// 			mem_hab[0]/1073741824.0, mem_hab[1]/1073741824.0, mem_hab[2]/1073741824.0, 
	// 			mem_hab[3]/1073741824.0, mem_hab[4]/1073741824.0,
	// 			mem[1]/1073741824.0, mem[2]/1073741824.0,
    //             mem[3]/1073741824.0, mem[4]/1073741824.0, mem[5]/1073741824.0);
    // }

    // align = kv_ul_ov_t_statistics(&(bl->tk), i, &(b->num_recorrect_base));
    // if(align == s->len[i]) {
    //  free(s->seq[i]); s->seq[i] = NULL;
    // }
    // b->num_correct_base += align;

    // uint64_t k;
    // b->num_read_base += overlap_statistics(&b->olist, NULL, NULL, 1);
    // for (k = 0; k < bl->tk.n; k++) {
    //  if(bl->tk.a[k].sec == 0) b->num_correct_base += bl->tk.a[k].qe - bl->tk.a[k].qs;
    //  if(bl->tk.a[k].sec > 0) b->num_recorrect_base += bl->tk.a[k].qe - bl->tk.a[k].qs;
    // }
    // for (k = 0; k < bl->lo.n; k++) {
    //  b->num_read_base += bl->lo.a[k].qe - bl->lo.a[k].qs;
    // }

    // uint32_t l1 = overlap_statistics(&b->olist, s->uu->ug, 1), l2 = overlap_statistics(&b->olist, s->uu->ug, 2);
    // 
    // if(l1 == 0 && l2 > 0) fprintf(stderr, "[M::%s::%lu::no_match]\n", UL_INF.nid.a[s->id+i].a, s->len[i]);
    // fprintf(stderr, "[M::%s::%lu::] l1->%u; l2->%u\n", UL_INF.nid.a[s->id+i].a, s->len[i], l1, l2);
    // fprintf(stderr, "[M::%s::rid->%ld] done\n", __func__, s->id+i);
	// exit(1);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x58, %rsp
movq %rsi, %r13
movq %rdi, %rbx
movq (%rdi), %rcx
movslq %edx, %r14
imulq $0x68, %r14, %r15
addq 0x88(%rdi), %r15
movq 0x80(%rdi), %rax
movsd 0x409ad(%rip), %xmm0    # 0x1e3570
movq %rcx, 0x10(%rsp)
divsd 0x68(%rcx), %xmm0
movsd %xmm0, (%rsp)
movq (%rax,%r14,8), %rax
movq %rax, 0x48(%rsp)
leaq 0xb1736(%rip), %rax     # 0x254318
imulq $0x58, 0x50(%rdi), %rbp
addq 0x48(%rax), %rbp
imulq $0x58, %rsi, %r12
leaq (%r12,%rbp), %rdi
callq 0x17853e
movl 0x18(%r12,%rbp), %ecx
testl %eax, %eax
je 0x1a2c39
movq 0x48(%rbx), %rax
cmpq %rcx, (%rax,%r13,8)
jne 0x1a2c1a
movq 0x58(%rbx), %rax
cmpq $0x0, (%rax,%r13,8)
je 0x1a2eb3
leaq 0x5437f(%rip), %rdi     # 0x1f6fa0
leaq 0x4eefb(%rip), %rsi     # 0x1f1b23
leaq 0x543a7(%rip), %rcx     # 0x1f6fd6
movl $0x315d, %edx           # imm = 0x315D
callq 0x2110
movsd (%rsp), %xmm0
movq 0x10(%rsp), %r11
movq %r15, 0x40(%rsp)
movq 0x48(%rsp), %rsi
movq %r14, 0x38(%rsp)
testl %ecx, %ecx
jns 0x1a2ec2
leaq (%r12,%rbp), %rax
addq $0x18, %rax
andl $0x7fffffff, %ecx       # imm = 0x7FFFFFFF
movl %ecx, (%rax)
movq 0x48(%rbx), %rax
cmpq %rcx, (%rax,%r13,8)
movq %r12, 0x30(%rsp)
jne 0x1a2ee1
movq %r13, %r14
minsd 0x408f0(%rip), %xmm0    # 0x1e3578
cvttsd2si %xmm0, %rax
movq %rax, 0x28(%rsp)
movq %rsi, %r13
movq 0xe8(%rsi), %rdi
movq 0x58(%rbx), %rax
movq (%rax,%r14,8), %rdx
movslq (%r11), %r8
movslq 0x4(%r11), %r9
leaq 0xa0(%rsi), %rax
movq %rax, (%rsp)
leaq 0xc0(%rsi), %rax
movq %rax, 0x8(%rsp)
leaq 0x48(%rsi), %r10
movsd 0x60(%r11), %xmm0
movl 0x80(%r11), %r11d
leaq 0x12160(%rsi), %r12
movq %r12, 0x10(%rsp)
leaq 0x121a0(%rsi), %r15
leaq 0x12250(%rsi), %rbp
subq $0x8, %rsp
pushq $0x1
popq %rax
movq %r14, %rsi
pushq $0x0
pushq %rax
pushq %rbp
pushq $0x0
pushq %r15
pushq %r12
pushq $0x0
pushq %rax
pushq %r11
pushq %r10
pushq 0x60(%rsp)
movq 0x60(%rsp), %rbp
pushq %rbp
pushq 0x30(%rbx)
callq 0xe57e0
addq $0x70, %rsp
leaq 0x108(%r13), %rdi
movq %rdi, 0x20(%rsp)
callq 0x1b25d
leaq 0x9330(%r13), %r15
movq %r15, %rdi
callq 0x1b1e8
movq 0x58(%rbx), %rax
movq (%rax,%r14,8), %rax
leaq 0x8(%r13), %rdx
movq %rax, 0x8(%r13)
movq 0x48(%rbx), %rax
movq (%rax,%r14,8), %rax
movq %rax, 0x10(%r13)
andq $0x0, 0x18(%r13)
movq (%rbx), %rdi
movq 0x30(%rbx), %rsi
leaq 0x2c0(%r13), %rcx
movq %rcx, 0x8(%rsp)
leaq 0x28(%r13), %r8
leaq 0x148(%r13), %r9
leaq 0x200(%r13), %rax
leaq 0x9040(%r13), %r12
movq %r12, 0x18(%rsp)
leaq 0x121e8(%r13), %r10
movsd 0x68(%rdi), %xmm0
leaq 0x50(%rsp), %r11
movq %rbp, %rdi
pushq $0x0
movq 0x30(%rsp), %rbp
pushq %rbp
pushq %r11
leaq 0x6c(%rsp), %r11
pushq %r11
pushq $0x1
popq %r11
pushq %r11
pushq $0x0
pushq %r10
pushq 0x48(%rsp)
pushq %r15
pushq %r12
pushq 0x70(%rsp)
pushq %rax
callq 0x29452
addq $0x60, %rsp
movq 0x78(%rbx), %rax
movq 0x38(%rsp), %r12
movq (%rax,%r12,8), %rdi
movq 0x50(%rbx), %rax
addq %r14, %rax
imulq $0x58, %rax, %rsi
leaq 0xb1524(%rip), %r15     # 0x254318
addq 0x48(%r15), %rsi
imulq $0x18, %r12, %r9
addq 0x68(%rbx), %r9
imulq $0xc0, %r12, %rcx
addq 0x90(%rbx), %rcx
movq 0x30(%rbx), %rdx
movq (%rbx), %r8
movq 0x28(%rbx), %r10
movsd 0x68(%r8), %xmm0
movq 0x48(%rbx), %r8
movq (%r8,%r14,8), %r8
subq $0x50, %rsp
movq %rax, 0x30(%rsp)
movq %r10, 0x28(%rsp)
movq %r8, 0x20(%rsp)
movq %rbp, 0x18(%rsp)
movq %rdx, 0x10(%rsp)
movq %rcx, 0x8(%rsp)
andq $0x0, 0x40(%rsp)
movq 0x90(%rsp), %rax
movq %rax, (%rsp)
movq 0x50(%rsp), %rdx
movq 0x58(%rsp), %rcx
movq 0x68(%rsp), %r8
callq 0x16a8b5
addq $0x50, %rsp
xorpd %xmm0, %xmm0
movupd %xmm0, 0x8(%r13)
movupd %xmm0, 0x18(%r13)
imulq $0x58, 0x50(%rbx), %rdi
addq 0x48(%r15), %rdi
addq 0x30(%rsp), %rdi
callq 0x17853e
testl %eax, %eax
je 0x1a2ea1
incq 0xf8(%r13)
movq 0x80(%rbx), %rax
movq (%rax,%r12,8), %rax
incq 0xf0(%rax)
addq $0x58, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
leaq 0x54147(%rip), %rdi     # 0x1f7010
leaq 0x4ec53(%rip), %rsi     # 0x1f1b23
leaq 0x540ff(%rip), %rcx     # 0x1f6fd6
movl $0x3160, %edx           # imm = 0x3160
callq 0x2110
leaq 0x54158(%rip), %rdi     # 0x1f7040
leaq 0x4ec34(%rip), %rsi     # 0x1f1b23
leaq 0x540e0(%rip), %rcx     # 0x1f6fd6
movl $0x3162, %edx           # imm = 0x3162
callq 0x2110