// File: /workspace/github_cpp_makefile_star2_O2/chhylp123[P]hifiasm/Correct.cpp
// Function: push_wlst(ul_idx_t const*, hpc_t*, All_reads*, overlap_region*, char*, char*, char*, Correct_dumy*, long, long, long, long, long, long, long, long, long, double, long, unsigned int, double, void*)
uint32_t push_wlst(const ul_idx_t *uref, hpc_t *hpc_g, All_reads *rref, overlap_region* ol, 
                        char* qstr, char *tstr, char *tstr_1, Correct_dumy* dumy,
                        int64_t qs, int64_t qe, int64_t ts, int64_t te, int64_t tl, 
                        int64_t error, int64_t aux_beg, int64_t aux_end, int64_t thres, double e_rate, 
                        int64_t block_s, uint32_t sec_check, double ovlp_cut, void *km)
{

    window_list p, t, *a; int64_t w_e, w_s, ce = qs - 1, cs = ol->x_pos_s, toff, ovl, ualn, aln; 
    uint64_t a_n, k;
    
    p.x_start = qs; p.x_end = qe; p.y_start = ts; p.y_end = te; p.error = error;
    p.extra_begin = aux_beg; p.extra_end = aux_end;
    p.error_threshold = thres; p.cidx = p.clen = 0;
    if(ol->w_list.n > 0) { //utilize the the end pos of pre-window in forward
        w_e = ol->w_list.a[ol->w_list.n-1].x_end; 
        toff = ol->w_list.a[ol->w_list.n-1].y_end + 1 - ol->w_list.a[ol->w_list.n-1].extra_begin;
        while ((w_e < ce) && (toff < tl)) {
            w_s = w_e + 1;
            get_win_id_by_s(ol, w_s, block_s, &w_e);
            // x_start = w_s; x_end = w_e;
            if(aln_wlst_adv(ol, rref, hpc_g, uref, qstr, tstr, tstr_1, dumy,
            ol->y_pos_strand, ol->y_id, w_s, w_e, toff, block_s, e_rate, 0)) {
                toff = ol->w_list.a[ol->w_list.n-1].y_end + 1 - ol->w_list.a[ol->w_list.n-1].extra_begin;
            } else {
                break;
            }
        }
        cs = ol->w_list.a[ol->w_list.n-1].x_end + 1;
    }
    ///utilize the the start pos of next window in backward
    a_n = ol->w_list.n; w_s = qs;
    if(w_s > cs) {
        gen_backtrace_adv(&p, ol, rref, hpc_g, uref, qstr, tstr, tstr_1, dumy, ol->y_pos_strand, ol->y_id);
        p.y_end += p.extra_begin;
        toff = p.y_start - 1;
        while ((w_s > cs) && (toff > 0)) {
            w_e = w_s - 1;
            get_win_id_by_e(ol, w_e, block_s, &w_s);
            // x_start = w_s; x_end = w_e; x_len = x_end + 1 - x_start;
            if(aln_wlst_adv(ol, rref, hpc_g, uref, qstr, tstr, tstr_1, dumy,
                ol->y_pos_strand, ol->y_id, w_s, w_e, toff+1-(w_e+1-w_s), block_s, e_rate, 1)) {
                ///y_start has no shift, but y_end has shift  
                ol->w_list.a[ol->w_list.n-1].y_start -= ol->w_list.a[ol->w_list.n-1].extra_begin; 
                toff = ol->w_list.a[ol->w_list.n-1].y_start - 1;
            } else {
                break;
            }
        }
    }

    ol->align_length += qe + 1 - qs;
    ovl = ol->x_pos_e+1-ol->x_pos_s; ualn = (qe + 1 - ol->x_pos_s) - ol->align_length; aln = ovl-ualn;
    if((!simi_pass(ovl, aln, 0, ovlp_cut, &e_rate)) && (!simi_pass(ovl, aln, sec_check, ovlp_cut, NULL))) {
        kv_push(window_list, ol->w_list, p); 
        return 0;
    }

    if(ol->w_list.n > a_n) {
        a = ol->w_list.a + a_n; a_n = ol->w_list.n - a_n; toff = a_n; a_n >>=1;
        for (k = 0; k < a_n; k++) {
            t = a[k]; a[k] = a[toff-1-k]; a[toff-1-k] = t;
        }
    }
    kv_push(window_list, ol->w_list, p); 
    return 1;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x88, %rsp
movsd %xmm1, 0x58(%rsp)
movsd %xmm0, 0x30(%rsp)
movq %r9, 0x28(%rsp)
movq %r8, 0x20(%rsp)
movq %rcx, %rbx
movq %rdx, 0x18(%rsp)
movq %rsi, 0x10(%rsp)
movq %rdi, 0x8(%rsp)
movq 0x118(%rsp), %r15
movq 0xd8(%rsp), %rsi
movq 0xd0(%rsp), %rdi
movzwl 0x110(%rsp), %eax
movzwl 0x108(%rsp), %ecx
movl 0xe0(%rsp), %edx
movl 0x4(%rbx), %ebp
movl %edi, 0x38(%rsp)
movl %esi, 0x3c(%rsp)
movl %edx, 0x40(%rsp)
movl 0xe8(%rsp), %edx
movl %edx, 0x44(%rsp)
movzwl 0xf8(%rsp), %edx
movw %dx, 0x4c(%rsp)
movzwl 0x100(%rsp), %edx
movw %dx, 0x48(%rsp)
movw %cx, 0x4a(%rsp)
movw %ax, 0x4e(%rsp)
andq $0x0, 0x50(%rsp)
movq 0x48(%rbx), %r13
testq %r13, %r13
je 0x2a3af
leaq -0x1(%rdi), %rbp
movq 0x58(%rbx), %rcx
movq %r13, %rdx
shlq $0x5, %rdx
movslq -0x1c(%rcx,%rdx), %r12
movq %rsp, %rax
movq %r12, (%rax)
movl -0x14(%rcx,%rdx), %eax
movswl -0x10(%rcx,%rdx), %edx
subl %edx, %eax
incl %eax
cmpq %rbp, %r12
jge 0x2a396
movslq %eax, %r14
cmpq 0xf0(%rsp), %r14
jge 0x2a396
incq %r12
movq %rbx, %rdi
movq %r12, %rsi
movq %r15, %rdx
movq %rsp, %rcx
callq 0x51f8b
movl 0x10(%rbx), %eax
movl 0x1c(%rbx), %r10d
subq $0x8, %rsp
movq %rbx, %rdi
movq 0x20(%rsp), %rsi
movq 0x18(%rsp), %rdx
movq 0x10(%rsp), %rcx
movq 0x28(%rsp), %r8
movq 0x30(%rsp), %r9
movsd 0x38(%rsp), %xmm0
pushq $0x0
pushq %r15
pushq %r14
pushq 0x20(%rsp)
pushq %r12
pushq %rax
pushq %r10
pushq 0x108(%rsp)
pushq 0x108(%rsp)
callq 0x56321
addq $0x50, %rsp
movq 0x48(%rbx), %r13
movq 0x58(%rbx), %rcx
testl %eax, %eax
je 0x2a396
movq %r13, %rdx
shlq $0x5, %rdx
movl -0x14(%rcx,%rdx), %eax
movswl -0x10(%rcx,%rdx), %edx
subl %edx, %eax
incl %eax
movq (%rsp), %r12
jmp 0x2a2ec
movq %r13, %rax
shlq $0x5, %rax
movslq -0x1c(%rcx,%rax), %rbp
incq %rbp
movq 0xd0(%rsp), %rdi
jmp 0x2a3b2
xorl %r13d, %r13d
cmpq %rdi, %rbp
jge 0x2a4e2
movl 0x10(%rbx), %eax
movl 0x1c(%rbx), %r10d
subq $0x8, %rsp
leaq 0x40(%rsp), %r12
movq %r12, %rdi
movq %rbx, %rsi
movq 0x20(%rsp), %rdx
movq 0x18(%rsp), %rcx
movq 0x10(%rsp), %r8
movq 0x28(%rsp), %r9
pushq %rax
pushq %r10
pushq 0xe0(%rsp)
pushq 0xe0(%rsp)
pushq 0x50(%rsp)
callq 0x567ea
addq $0x30, %rsp
movswl 0x10(%r12), %eax
addl %eax, 0xc(%r12)
movl 0x8(%r12), %ecx
movq 0xd0(%rsp), %rsi
cmpq %rbp, %rsi
jle 0x2a4e2
cmpl $0x2, %ecx
jl 0x2a4e2
leaq -0x1(%rsi), %r11
movq %r11, (%rsp)
movl 0x10(%rbx), %r10d
movl 0x4(%rbx), %r12d
movq %r12, %rax
xorl %edx, %edx
idivq %r15
subq %r12, %rdx
leaq (%rsi,%rdx), %rax
decq %rax
cqto
idivq %r15
negq %rdx
leaq (%rsi,%rdx), %rax
decq %rax
cmpq %r12, %rax
cmovgq %rax, %r12
movl 0x1c(%rbx), %eax
movl %ecx, %r14d
subq %rsi, %r14
addq %r12, %r14
subq $0x8, %rsp
movq %rbx, %rdi
movq 0x20(%rsp), %rsi
movq 0x18(%rsp), %rdx
movq 0x10(%rsp), %rcx
movq 0x28(%rsp), %r8
movq 0x30(%rsp), %r9
movsd 0x38(%rsp), %xmm0
pushq $0x1
pushq %r15
pushq %r14
pushq %r11
pushq %r12
pushq %r10
pushq %rax
pushq 0x108(%rsp)
pushq 0x108(%rsp)
callq 0x56321
addq $0x50, %rsp
testl %eax, %eax
je 0x2a4e2
movq 0x48(%rbx), %rax
movq 0x58(%rbx), %rdx
shlq $0x5, %rax
movswl -0x10(%rdx,%rax), %esi
movl -0x18(%rdx,%rax), %ecx
subl %esi, %ecx
movl %ecx, -0x18(%rdx,%rax)
movq %r12, %rsi
jmp 0x2a41b
movq 0xd8(%rsp), %rsi
movl %esi, %eax
movq 0xd0(%rsp), %rcx
subl %ecx, %eax
movl 0x28(%rbx), %ecx
addl %ecx, %eax
incl %eax
movl 0x4(%rbx), %ecx
movl 0x8(%rbx), %edx
incl %edx
movl %edx, %edi
subl %ecx, %edi
notq %rsi
addq %rcx, %rsi
addq %rax, %rsi
addq %rdi, %rsi
movl %eax, 0x28(%rbx)
cmpl %ecx, %edx
je 0x2a543
testq %rsi, %rsi
je 0x2a543
cvtsi2sd %rdi, %xmm0
movsd 0x1b8d72(%rip), %xmm1   # 0x1e32a0
subsd 0x30(%rsp), %xmm1
mulsd %xmm0, %xmm1
cvtsi2sd %rsi, %xmm0
ucomisd %xmm1, %xmm0
jae 0x2a595
movl 0x120(%rsp), %edx
movsd 0x58(%rsp), %xmm0
xorl %ecx, %ecx
callq 0x56b28
movq 0x48(%rbx), %rcx
testl %eax, %eax
jne 0x2a599
cmpq 0x50(%rbx), %rcx
jne 0x2a64f
leaq (%rcx,%rcx), %rax
testq %rcx, %rcx
pushq $0x2
popq %rsi
cmovneq %rax, %rsi
movq %rsi, 0x50(%rbx)
movq 0x58(%rbx), %rdi
shlq $0x5, %rsi
callq 0x2370
movq %rax, 0x58(%rbx)
movq 0x48(%rbx), %rcx
jmp 0x2a653
movq 0x48(%rbx), %rcx
movq %rcx, %rax
subq %r13, %rax
jbe 0x2a608
movq 0x58(%rbx), %rdx
shlq $0x5, %r13
addq %rdx, %r13
shrq %rax
shlq $0x5, %rcx
addq %rdx, %rcx
addq $-0x20, %rcx
subq $0x1, %rax
jb 0x2a604
movups (%r13), %xmm0
movups 0x10(%r13), %xmm1
movaps %xmm1, 0x70(%rsp)
movaps %xmm0, 0x60(%rsp)
movups (%rcx), %xmm0
movups 0x10(%rcx), %xmm1
movups %xmm1, 0x10(%r13)
movups %xmm0, (%r13)
movapd 0x60(%rsp), %xmm0
movapd 0x70(%rsp), %xmm1
movupd %xmm1, 0x10(%rcx)
movupd %xmm0, (%rcx)
addq $-0x20, %rcx
addq $0x20, %r13
jmp 0x2a5ba
movq 0x48(%rbx), %rcx
cmpq 0x50(%rbx), %rcx
jne 0x2a637
leaq (%rcx,%rcx), %rax
testq %rcx, %rcx
pushq $0x2
popq %rsi
cmovneq %rax, %rsi
movq %rsi, 0x50(%rbx)
movq 0x58(%rbx), %rdi
shlq $0x5, %rsi
callq 0x2370
movq %rax, 0x58(%rbx)
movq 0x48(%rbx), %rcx
jmp 0x2a63b
movq 0x58(%rbx), %rax
leaq 0x1(%rcx), %rdx
movq %rdx, 0x48(%rbx)
shlq $0x5, %rcx
addq %rax, %rcx
pushq $0x1
popq %rax
jmp 0x2a664
movq 0x58(%rbx), %rax
leaq 0x1(%rcx), %rdx
movq %rdx, 0x48(%rbx)
shlq $0x5, %rcx
addq %rax, %rcx
xorl %eax, %eax
movups 0x38(%rsp), %xmm0
movups 0x48(%rsp), %xmm1
movups %xmm1, 0x10(%rcx)
movups %xmm0, (%rcx)
addq $0x88, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq