// File: /workspace/github_cpp_makefile_star2_O3/chhylp123[P]hifiasm/Correct.cpp
// Function: partition_overlaps_advance_back(overlap_region_alloc*, All_reads*, UC_Read*, UC_Read*, Correct_dumy*, haplotype_evdience_alloc*, int)
void partition_overlaps_advance_back(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);

    long long i;
    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), 1, NULL);

    uint64_t pre_site = (uint64_t)-1;
    uint64_t num_of_snps = 0;
    long long pre_i = -1;
    long long sub_length;
    haplotype_evdience* sub_list;

    ////split reads
    for (i = 0; i < hap->length; i++)
    {
        if(pre_site != hap->list[i].site)
        {
            if(i != 0)
            {
                sub_list = hap->list + pre_i;
                sub_length = i - pre_i;
                split_sub_list(hap, sub_list, sub_length, overlap_list, R_INF, g_read);
            }
            num_of_snps++;
            pre_site = hap->list[i].site;
            pre_i = i;
        }
    }

    if(pre_i != -1)
    {
        sub_list = hap->list + pre_i;
        sub_length = i - pre_i;
        split_sub_list(hap, sub_list, sub_length, overlap_list, R_INF, g_read);
    }

    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 $0x48, %rsp
movq %r8, %rbx
movq %rcx, 0x40(%rsp)
movq %rdx, %r12
movq %rsi, 0x38(%rsp)
movq %rdi, %r15
xorl %r14d, %r14d
movl %r14d, 0x230(%r9)
movq %r14, 0x220(%r9)
movl %r14d, 0x10(%r9)
movq %r14, 0x8(%r9)
movq %r9, 0x18(%rsp)
leaq 0x18(%r9), %rdi
movl $0x1ac, %edx            # imm = 0x1AC
xorl %esi, %esi
callq 0x20e0
movq %r12, 0x10(%rsp)
movq 0x8(%r12), %rax
movl $0x177, %ecx            # imm = 0x177
cmpq %rcx, %rax
movl $0x177, %edx            # imm = 0x177
cmovlq %rax, %rdx
decq %rdx
movq %rax, 0x20(%rsp)
leaq -0x1(%rax), %rcx
leaq 0x30(%rbx), %r8
xorps %xmm0, %xmm0
movl $0x0, 0xc(%rsp)
xorl %ebp, %ebp
movq %rcx, 0x30(%rsp)
movq %r8, 0x28(%rsp)
cmpq %rcx, %rdx
jne 0x563b3
cmpq $0x0, 0x20(%rsp)
je 0x564eb
cmpl $0x1, 0xc(%rsp)
movl $0x1, 0xc(%rsp)
je 0x564eb
movl $0x177, %eax            # imm = 0x177
leaq (%rdx,%rax), %r12
cmpq 0x20(%rsp), %r12
cmovgeq %rcx, %r12
cmpl $-0x2, %ebp
je 0x564eb
movaps %xmm0, (%r8)
testq %r14, %r14
je 0x563dd
movq 0x48(%rbx), %rax
jmp 0x563e7
movq $0x0, 0x48(%rbx)
xorl %eax, %eax
movq 0x10(%r15), %rcx
cmpq %rcx, %rax
jae 0x56416
movq (%r15), %rcx
imulq $0xa8, %rax, %rsi
movl 0x4(%rcx,%rsi), %edi
cmpq %rdi, %rdx
jge 0x56425
movq $0x0, 0x48(%rbx)
movaps %xmm0, (%r8)
xorl %ebp, %ebp
jmp 0x5649d
movq %rcx, 0x48(%rbx)
movaps %xmm0, (%r8)
movl $0xfffffffe, %ebp       # imm = 0xFFFFFFFE
jmp 0x5649d
movaps %xmm0, (%r8)
addq %rsi, %rcx
addq $0x4, %rcx
xorl %edi, %edi
xorl %esi, %esi
xorl %r8d, %r8d
movl (%rcx), %r9d
movl 0x4(%rcx), %r10d
cmpq %rdx, %r10
cmovaeq %rdx, %r10
cmpq %r14, %r9
movq %r14, %r11
cmovgq %r9, %r11
cmpq %r11, %r10
jb 0x5647f
cmpb $0x1, 0x28(%rcx)
jne 0x5646d
movq 0x28(%rbx), %r10
movq %rax, (%r10,%rdi,8)
movq 0x30(%rbx), %rdi
incq %rdi
movq %rdi, 0x30(%rbx)
incq %rsi
testq %r8, %r8
jne 0x56479
movq %rax, 0x48(%rbx)
movl $0x1, %r8d
cmpq %r9, %rdx
jl 0x56494
incq %rax
addq $0xa8, %rcx
cmpq 0x10(%r15), %rax
jb 0x56437
xorl %ebp, %ebp
testq %rsi, %rsi
setne %bpl
movl $0x177, %eax            # imm = 0x177
leaq (%r14,%rax), %r13
movq 0x10(%rsp), %rax
movq (%rax), %rdi
subq $0x8, %rsp
movq %r14, %rsi
movq %r15, %rcx
movq %rbx, %r8
movq 0x40(%rsp), %r9
pushq $0x1
pushq 0x50(%rsp)
pushq 0x30(%rsp)
callq 0x4f4dd
addq $0x20, %rsp
movq %r13, %r14
movq %r12, %rdx
movq 0x30(%rsp), %rcx
movq 0x28(%rsp), %r8
xorps %xmm0, %xmm0
jmp 0x5638f
movq 0x18(%rsp), %r13
movl 0x230(%r13), %r14d
cmpq %r14, 0x228(%r13)
jae 0x56565
decq %r14
movq %r14, %rax
shrq %rax
orq %r14, %rax
movq %rax, %rcx
shrq $0x2, %rcx
orq %rax, %rcx
movq %rcx, %rax
shrq $0x4, %rax
orq %rcx, %rax
movq %rax, %rcx
shrq $0x8, %rcx
orq %rax, %rcx
movq %rcx, %rax
shrq $0x10, %rax
orq %rcx, %rax
leaq 0x1(%rax), %rcx
movq %rcx, 0x228(%r13)
movq 0x218(%r13), %rdi
leaq 0x8(,%rax,8), %rax
leaq (%rax,%rax,4), %rsi
callq 0x2370
movq %rax, 0x218(%r13)
movl 0x230(%r13), %r14d
movq $0x0, 0x220(%r13)
movq 0x10(%r15), %r12
movl %r12d, 0x1c8(%r13)
movl $0x0, 0x1c4(%r13)
imulq %r12, %r14
movl 0x1d8(%r13), %eax
cmpq %rax, %r14
jbe 0x565b8
movl %r14d, 0x1d8(%r13)
movq 0x1d0(%r13), %rdi
movl %r14d, %esi
callq 0x2370
movq %rax, 0x1d0(%r13)
movq %rax, %rdi
jmp 0x565bf
movq 0x1d0(%r13), %rdi
movl $0xff, %esi
movq %r14, %rdx
callq 0x20e0
movl 0x1e8(%r13), %eax
cmpq %rax, %r12
jbe 0x565f5
movl %r12d, 0x1e8(%r13)
movq 0x1e0(%r13), %rdi
movl %r12d, %esi
callq 0x2370
movq %rax, 0x1e0(%r13)
movl 0x10(%r13), %eax
testl %eax, %eax
je 0x56684
movq (%r13), %rcx
movq $-0x1, %rdi
xorl %r12d, %r12d
movq $-0x1, %rdx
xorl %r14d, %r14d
movl (%rcx,%r12), %esi
cmpq %rsi, %rdx
je 0x56651
testq %r12, %r12
je 0x5664a
leaq (%rdi,%rdi,2), %rax
leaq (%rcx,%rax,8), %rsi
movq %r14, %rdx
subq %rdi, %rdx
movq %r13, %rdi
movq 0x10(%rsp), %r9
callq 0x4fd28
movq (%r13), %rcx
movl 0x10(%r13), %eax
movl (%rcx,%r12), %edx
movq %r14, %rdi
incq %r14
movl %eax, %esi
addq $0x18, %r12
cmpq %rsi, %r14
jb 0x56619
cmpq $-0x1, %rdi
je 0x56684
leaq (%rdi,%rdi,2), %rsi
shlq $0x3, %rsi
addq (%r13), %rsi
subq %rdi, %r14
movq %r13, %rdi
movq %r14, %rdx
movq 0x10(%rsp), %r9
callq 0x4fd28
movq %r13, %rdi
movq %r15, %rsi
callq 0x50e1e
movq 0x10(%rsp), %rax
movq 0x8(%rax), %rsi
movq %r15, %rdi
movq %rbx, %rdx
addq $0x48, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x50add