// File: /workspace/github_cpp_makefile_star2_O1/chhylp123[P]hifiasm/Correct.cpp
// Function: push_wlst_exz(ul_idx_t const*, hpc_t*, All_reads*, overlap_region*, char*, char*, bit_extz_t*, unsigned int, long, long, long, long, long, long, long, double, long, unsigned int, double, long, void*)
uint32_t push_wlst_exz(const ul_idx_t *uref, hpc_t *hpc_g, All_reads *rref, overlap_region* ol, 
                        char* qstr, char *tstr, bit_extz_t *exz, uint32_t max_err,
                        int64_t qs, int64_t qe, int64_t ts, int64_t te, int64_t tl, 
                        int64_t aux_beg, int64_t aux_end, double e_rate, int64_t block_s, uint32_t sec_check, double ovlp_cut, int64_t force_aln, 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, ys; 
    uint64_t a_n, k;
    
    p.x_start = qs; p.x_end = qe; p.y_start = ts; p.y_end = te; p.error = exz->err;
    p.extra_begin = aux_beg; p.extra_end = aux_end; p.error_threshold = exz->thre; 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;
        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_exz(ol, rref, hpc_g, uref, qstr, tstr, exz, max_err,
            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;
            } 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_exz(&p, ol, rref, hpc_g, uref, qstr, tstr, exz, ol->y_pos_strand, ol->y_id);
        toff = p.y_start - 1;
        while (w_s > cs) {
            w_e = w_s - 1;
            get_win_id_by_e(ol, w_e, block_s, &w_s); ys = toff+1-(w_e+1-w_s);
            // x_start = w_s; x_end = w_e; x_len = x_end + 1 - x_start;
            if((ys >= 0) && aln_wlst_adv_exz(ol, rref, hpc_g, uref, qstr, tstr, exz, max_err,
                ol->y_pos_strand, ol->y_id, w_s, w_e, ys, block_s, e_rate, 1)) {
                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((!force_aln) && (!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, 0x8(%rsp)
movq %r9, 0x30(%rsp)
movq %r8, 0x28(%rsp)
movq %rcx, %rbx
movq %rdx, 0x20(%rsp)
movq %rsi, 0x18(%rsp)
movq %rdi, 0x10(%rsp)
movq 0x108(%rsp), %r14
movq 0xd8(%rsp), %rcx
movq 0xd0(%rsp), %rsi
movq 0xc0(%rsp), %rdx
movl 0xe0(%rsp), %eax
movl 0x4(%rbx), %r12d
movl %esi, 0x38(%rsp)
movl %ecx, 0x3c(%rsp)
movl %eax, 0x40(%rsp)
movl 0xe8(%rsp), %eax
movl %eax, 0x44(%rsp)
movzwl 0x2c(%rdx), %eax
movw %ax, 0x4c(%rsp)
movzwl 0xf8(%rsp), %eax
movw %ax, 0x48(%rsp)
movzwl 0x100(%rsp), %eax
movw %ax, 0x4a(%rsp)
movzwl 0x28(%rdx), %eax
movw %ax, 0x4e(%rsp)
movq $0x0, 0x50(%rsp)
movq 0x48(%rbx), %rax
testq %rax, %rax
je 0x5ad63
movq 0xf0(%rsp), %rdx
leaq -0x1(%rsi), %r12
movq 0x58(%rbx), %rcx
shlq $0x5, %rax
movslq -0x1c(%rcx,%rax), %r15
movslq -0x14(%rcx,%rax), %r10
incq %r10
cmpq %r15, %r12
setg %al
cmpq %rdx, %r10
setl %cl
andb %al, %cl
cmpb $0x1, %cl
jne 0x5ad47
movq %rdx, %rbp
movq %r15, %r11
incq %r11
movl 0x8(%rbx), %r15d
movl 0x4(%rbx), %ecx
movq %rcx, %rax
xorl %edx, %edx
idivq %r14
subq %rdx, %rcx
movq %r11, %rax
subq %rcx, %rax
cqto
idivq %r14
incq %rax
imulq %r14, %rax
addq %rcx, %rax
decq %rax
cmpq %r15, %rax
cmovlq %rax, %r15
movl 0x1c(%rbx), %eax
movl 0x10(%rbx), %r13d
subq $0x8, %rsp
movq %rbx, %rdi
movq 0x28(%rsp), %rsi
movq 0x20(%rsp), %rdx
movq 0x18(%rsp), %rcx
movq 0x30(%rsp), %r8
movq 0x38(%rsp), %r9
movsd 0x10(%rsp), %xmm0
pushq $0x0
pushq %r14
pushq %r10
pushq %r15
pushq %r11
pushq %r13
pushq %rax
movl 0x108(%rsp), %eax
pushq %rax
pushq 0x108(%rsp)
callq 0x8b0aa
addq $0x50, %rsp
testl %eax, %eax
je 0x5ad47
cmpq %r12, %r15
jge 0x5ad47
movq 0x48(%rbx), %rax
movq 0x58(%rbx), %rcx
shlq $0x5, %rax
movslq -0x14(%rcx,%rax), %r10
incq %r10
cmpq %rbp, %r10
jl 0x5ac98
movq 0x48(%rbx), %rax
movq 0x58(%rbx), %rcx
shlq $0x5, %rax
movslq -0x1c(%rcx,%rax), %r12
incq %r12
movq 0xd0(%rsp), %rsi
movq 0x48(%rbx), %r15
cmpq %rsi, %r12
jge 0x5ae68
movl 0x10(%rbx), %eax
movl 0x1c(%rbx), %r10d
leaq 0x38(%rsp), %r13
movq %r13, %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 0xd0(%rsp)
pushq 0x48(%rsp)
callq 0x8b5b0
addq $0x20, %rsp
movl 0x8(%r13), %edi
movq 0xd0(%rsp), %r13
movq %r13, %r8
leaq -0x1(%r13), %r10
movl 0x4(%rbx), %ecx
movq %rcx, %rax
xorl %edx, %edx
idivq %r14
movq %rcx, %r13
subq %rdx, %r13
movq %r10, %rsi
subq %r13, %rsi
movq %rsi, %rax
cqto
idivq %r14
subq %rdx, %r13
addq %rsi, %r13
cmpq %rcx, %r13
cmovleq %rcx, %r13
movslq %edi, %rcx
movq %r13, %rax
subq %r8, %rax
addq %rcx, %rax
js 0x5ae68
movl 0x10(%rbx), %r11d
movl 0x1c(%rbx), %ebp
subq $0x8, %rsp
movq %rbx, %rdi
movq 0x28(%rsp), %rsi
movq 0x20(%rsp), %rdx
movq 0x18(%rsp), %rcx
movq 0x30(%rsp), %r8
movq 0x38(%rsp), %r9
movsd 0x10(%rsp), %xmm0
pushq $0x1
pushq %r14
pushq %rax
pushq %r10
pushq %r13
pushq %r11
pushq %rbp
movl 0x108(%rsp), %eax
pushq %rax
pushq 0x108(%rsp)
callq 0x8b0aa
addq $0x50, %rsp
testl %eax, %eax
je 0x5ae68
movq 0x48(%rbx), %rax
movq 0x58(%rbx), %rcx
shlq $0x5, %rax
movl -0x18(%rcx,%rax), %edi
cmpq %r12, %r13
jg 0x5adb9
movq 0xd8(%rsp), %r8
movl %r8d, %eax
movq 0xd0(%rsp), %rcx
subl %ecx, %eax
movl 0x28(%rbx), %ecx
addl %ecx, %eax
incl %eax
movl %eax, 0x28(%rbx)
cmpq $0x0, 0x118(%rsp)
je 0x5af19
movq 0x48(%rbx), %rcx
movq %rcx, %rax
subq %r15, %rax
jbe 0x5af0e
cmpq $0x2, %rax
jb 0x5af0e
movq 0x58(%rbx), %rdx
shlq $0x5, %r15
addq %rdx, %r15
shrq %rax
cmpq $0x1, %rax
adcq $0x0, %rax
shlq $0x5, %rcx
addq %rdx, %rcx
addq $-0x20, %rcx
movups (%r15), %xmm0
movups 0x10(%r15), %xmm1
movaps %xmm1, 0x70(%rsp)
movaps %xmm0, 0x60(%rsp)
movups (%rcx), %xmm0
movups 0x10(%rcx), %xmm1
movups %xmm1, 0x10(%r15)
movups %xmm0, (%r15)
movapd 0x60(%rsp), %xmm0
movapd 0x70(%rsp), %xmm1
movupd %xmm1, 0x10(%rcx)
movupd %xmm0, (%rcx)
addq $0x20, %r15
addq $-0x20, %rcx
decq %rax
jne 0x5aec9
movq 0x50(%rbx), %rax
movl $0x1, %ebp
jmp 0x5af45
movl 0x4(%rbx), %edx
movl 0x8(%rbx), %esi
incl %esi
movl %esi, %ecx
subl %edx, %ecx
notq %r8
addq %rdx, %r8
addq %rax, %r8
addq %rcx, %r8
sete %dil
cmpl %edx, %esi
sete %al
orb %dil, %al
je 0x5afab
movq 0x50(%rbx), %rax
xorl %ebp, %ebp
cmpq %rax, 0x48(%rbx)
jne 0x5af70
leaq (%rax,%rax), %rcx
testq %rax, %rax
movl $0x2, %esi
cmovneq %rcx, %rsi
movq %rsi, 0x50(%rbx)
movq 0x58(%rbx), %rdi
shlq $0x5, %rsi
callq 0x2370
movq %rax, 0x58(%rbx)
movq 0x48(%rbx), %rax
movq 0x58(%rbx), %rcx
leaq 0x1(%rax), %rdx
movq %rdx, 0x48(%rbx)
shlq $0x5, %rax
movups 0x38(%rsp), %xmm0
movups 0x48(%rsp), %xmm1
movups %xmm1, 0x10(%rcx,%rax)
movups %xmm0, (%rcx,%rax)
movl %ebp, %eax
addq $0x88, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
movl %ecx, %ecx
xorps %xmm0, %xmm0
cvtsi2sd %rcx, %xmm0
movsd 0x1ca2b3(%rip), %xmm1   # 0x225270
subsd 0x8(%rsp), %xmm1
mulsd %xmm0, %xmm1
xorps %xmm0, %xmm0
cvtsi2sd %r8, %xmm0
ucomisd %xmm1, %xmm0
jae 0x5ae96
testb %al, %al
jne 0x5af3f
cmpl $0x0, 0x110(%rsp)
xorps %xmm1, %xmm1
cvtsi2sd %rcx, %xmm1
je 0x5b00a
cmpq $0x1c2, %r8             # imm = 0x1C2
jl 0x5af3f
mulsd 0x1ca4e8(%rip), %xmm1   # 0x2254f0
jmp 0x5b010
mulsd 0x58(%rsp), %xmm1
ucomisd %xmm1, %xmm0
jae 0x5ae96
jmp 0x5af3f