// File: /workspace/github_cpp_makefile_star2_O1/chhylp123[P]hifiasm/Correct.cpp
// Function: partition_overlaps_advance(overlap_region_alloc*, All_reads*, UC_Read*, UC_Read*, Correct_dumy*, haplotype_evdience_alloc*, int)
void partition_overlaps_advance(overlap_region_alloc* overlap_list, All_reads* R_INF, 
                        UC_Read* g_read, UC_Read* overlap_read, Correct_dumy* dumy, 
                        haplotype_evdience_alloc* hap, int force_repeat)
{
    ResizeInitHaplotypeEvdience(hap);

    uint64_t k, l, m;
    long long window_start, window_end;
    long long num_availiable_win = 0;
    
    Window_Pool w_inf;
    init_Window_Pool(&w_inf, g_read->length, WINDOW, (int)(1.0/asm_opt.max_ov_diff_ec));

    int flag = 0;
    while(get_Window(&w_inf, &window_start, &window_end) && flag != -2)
    {
        dumy->length = 0;
        dumy->lengthNT = 0;
        ///return overlaps that is overlaped with [window_start, window_end]
        flag = get_available_interval(window_start, window_end, overlap_list, dumy);
        switch (flag)
        {
            case 1:    ///found matched overlaps
                break;
            case 0:    ///do not find any matched overlaps
                break;
            case -2: ///do not find any matched overlaps, and the next window also cannot match
                break;
        }
        
        num_availiable_win = num_availiable_win + dumy->length;
        
        ///need to deal with 
        cluster_advance(g_read->seq, window_start, window_end, overlap_list, dumy, R_INF, hap, overlap_read, 1);
    }

    ///very time-consuming
    ///Fix-attention ---> able to be sorted locally
    // qsort(hap->list, hap->length, sizeof(haplotype_evdience), cmp_haplotype_evdience);
    SetSnpMatrix(hap, &(hap->nn_snp), &(overlap_list->length), 0, NULL);
    for (k = 1, l = m = 0; k <= hap->length; ++k) {  
        if (k == hap->length || hap->list[k].site != hap->list[l].site) {
            m += insert_snp_ee(hap, hap->list+l, k-l, hap->list+m, g_read, NULL);
            l = k;
        }
    }
    hap->length = m;

    // generate_haplotypes_naive_advance(hap, overlap_list, NULL);
    generate_haplotypes_naive_HiFi(hap, overlap_list, 0.04, g_read, 1, 0, ((uint64_t)-1));
    // generate_haplotypes_DP(hap, overlap_list, R_INF, g_read->length, force_repeat);
    // generate_haplotypes_naive(hap, overlap_list, R_INF, g_read->length, force_repeat);

    lable_large_indels(overlap_list, g_read->length, dumy, asm_opt.max_ov_diff_ec);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x38, %rsp
movq %r8, %rbx
movq %rcx, 0x30(%rsp)
movq %rdx, %r14
movq %rsi, 0x28(%rsp)
movq %rdi, %r15
xorl %r13d, %r13d
movl %r13d, 0x230(%r9)
movq %r13, 0x220(%r9)
movl %r13d, 0x10(%r9)
movq %r13, 0x8(%r9)
movq %r9, 0x10(%rsp)
leaq 0x18(%r9), %rdi
movl $0x1ac, %edx            # imm = 0x1AC
xorl %esi, %esi
callq 0x20e0
movq %r14, 0x8(%rsp)
movq 0x8(%r14), %rax
movl $0x177, %ecx            # imm = 0x177
cmpq %rcx, %rax
movl $0x177, %r12d           # imm = 0x177
cmovlq %rax, %r12
decq %r12
movq %rax, 0x18(%rsp)
decq %rax
movq %rax, 0x20(%rsp)
leaq 0x30(%rbx), %r14
xorps %xmm0, %xmm0
movl $0x0, 0x4(%rsp)
xorl %ebp, %ebp
movq %r12, %rdx
movq %r13, %rsi
cmpq 0x20(%rsp), %r12
jne 0x581eb
cmpq $0x0, 0x18(%rsp)
sete %al
cmpl $0x1, 0x4(%rsp)
sete %cl
orb %al, %cl
movl $0x1, %ecx
je 0x581ef
xorl %eax, %eax
movq %rsi, %r13
movq %rdx, %r12
jmp 0x5820d
movl 0x4(%rsp), %ecx
movl $0x177, %eax            # imm = 0x177
leaq (%rsi,%rax), %r13
leaq (%rdx,%rax), %r12
cmpq 0x18(%rsp), %r12
cmovgeq 0x20(%rsp), %r12
movb $0x1, %al
movl %ecx, 0x4(%rsp)
cmpl $-0x2, %ebp
je 0x58321
testb %al, %al
je 0x58321
movaps %xmm0, (%r14)
testq %rsi, %rsi
jne 0x5822f
movq $0x0, 0x48(%rbx)
movq 0x48(%rbx), %rax
movq 0x10(%r15), %rcx
cmpq %rcx, %rax
jae 0x5825f
movq (%r15), %rdi
imulq $0xa8, %rax, %rcx
movl 0x4(%rdi,%rcx), %edi
cmpq %rdi, %rdx
jge 0x5829e
movq $0x0, 0x48(%rbx)
movaps %xmm0, (%r14)
xorl %ebp, %ebp
jmp 0x5826c
movq %rcx, 0x48(%rbx)
movaps %xmm0, (%r14)
movl $0xfffffffe, %ebp       # imm = 0xFFFFFFFE
movq 0x8(%rsp), %rax
movq (%rax), %rdi
subq $0x8, %rsp
movq %r15, %rcx
movq %rbx, %r8
movq 0x30(%rsp), %r9
pushq $0x1
pushq 0x40(%rsp)
pushq 0x28(%rsp)
callq 0x50367
xorps %xmm0, %xmm0
addq $0x20, %rsp
jmp 0x581ba
movaps %xmm0, (%r14)
xorl %ebp, %ebp
cmpq 0x10(%r15), %rax
jae 0x5826c
movq (%r15), %rdi
addq %rdi, %rcx
addq $0x4, %rcx
xorl %edi, %edi
xorl %r8d, %r8d
movl (%rcx), %r9d
movl 0x4(%rcx), %r10d
cmpq %rdx, %r10
cmovgeq %rdx, %r10
cmpq %rsi, %r9
movq %rsi, %r11
cmovgq %r9, %r11
cmpq %r11, %r10
jl 0x582fe
cmpb $0x1, 0x28(%rcx)
jne 0x582ec
movq 0x28(%rbx), %r10
movq 0x30(%rbx), %r11
movq %rax, (%r10,%r11,8)
incq 0x30(%rbx)
incq %rdi
testq %r8, %r8
jne 0x582f8
movq %rax, 0x48(%rbx)
movl $0x1, %r8d
cmpq %r9, %rdx
jl 0x58313
incq %rax
addq $0xa8, %rcx
cmpq 0x10(%r15), %rax
jb 0x582b9
xorl %ebp, %ebp
testq %rdi, %rdi
setne %bpl
jmp 0x5826c
movq 0x10(%rsp), %rbp
movl 0x230(%rbp), %eax
cmpq %rax, 0x228(%rbp)
jae 0x58393
decq %rax
movq %rax, %rcx
shrq %rcx
orq %rax, %rcx
movq %rcx, %rax
shrq $0x2, %rax
orq %rcx, %rax
movq %rax, %rcx
shrq $0x4, %rcx
orq %rax, %rcx
movq %rcx, %rax
shrq $0x8, %rax
orq %rcx, %rax
movq %rax, %rcx
shrq $0x10, %rcx
orq %rax, %rcx
leaq 0x1(%rcx), %rax
movq %rax, 0x228(%rbp)
movq 0x218(%rbp), %rdi
leaq 0x8(,%rcx,8), %rax
leaq (%rax,%rax,4), %rsi
callq 0x2370
movq %rax, 0x218(%rbp)
movq $0x0, 0x220(%rbp)
movl 0x10(%r15), %eax
movl %eax, 0x1c8(%rbp)
movl $0x0, 0x1c4(%rbp)
movl 0x10(%rbp), %ecx
testq %rcx, %rcx
je 0x58428
movl $0x1, %r13d
movl $0x18, %r12d
xorl %r14d, %r14d
xorl %eax, %eax
cmpq %rcx, %r13
je 0x583e1
movq (%rbp), %rcx
movl (%rcx,%r12), %edx
leaq (%rax,%rax,2), %rsi
cmpl (%rcx,%rsi,8), %edx
je 0x58413
movq (%rbp), %rcx
leaq (%rax,%rax,2), %rdx
leaq (%rcx,%rdx,8), %rsi
movq %r13, %rdx
subq %rax, %rdx
leaq (%r14,%r14,2), %rax
leaq (%rcx,%rax,8), %rcx
movq %rbp, %rdi
movq 0x8(%rsp), %r8
xorl %r9d, %r9d
callq 0x571e7
cltq
addq %rax, %r14
movq %r13, %rax
leaq 0x1(%r13), %rdx
movl 0x10(%rbp), %ecx
addq $0x18, %r12
cmpq %rcx, %r13
movq %rdx, %r13
jb 0x583cb
jmp 0x5842b
xorl %r14d, %r14d
movl %r14d, 0x10(%rbp)
movsd 0x1cd0b1(%rip), %xmm0   # 0x2254e8
movl $0x1, %ecx
xorps %xmm1, %xmm1
movq %rbp, %rdi
movq %r15, %rsi
movq $-0x1, %r8
callq 0x52bf1
movq 0x8(%rsp), %rax
movq 0x8(%rax), %rsi
movq %r15, %rdi
movq %rbx, %rdx
addq $0x38, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x5198e