// File: /workspace/github_cpp_makefile_star2_O3/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), %r15
movq 0xd8(%rsp), %rcx
movq 0xd0(%rsp), %rdi
movq 0xc0(%rsp), %rdx
movl 0xe0(%rsp), %eax
movl 0x4(%rbx), %r12d
movl %edi, 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), %r14
testq %r14, %r14
je 0x5a2b0
movq 0xf0(%rsp), %rsi
leaq -0x1(%rdi), %r12
movq 0x58(%rbx), %rcx
movq %r14, %rdx
shlq $0x5, %rdx
movl -0x1c(%rcx,%rdx), %eax
movslq %eax, %r13
movslq -0x14(%rcx,%rdx), %r10
incq %r10
cmpq %r13, %r12
setg %cl
cmpq %rsi, %r10
setl %dl
andb %cl, %dl
cmpb $0x1, %dl
jne 0x5a2a9
movq %rsi, %rbp
movq %r13, %r11
incq %r11
movl 0x8(%rbx), %r13d
movl 0x4(%rbx), %ecx
movq %rcx, %rax
xorl %edx, %edx
idivq %r15
subq %rdx, %rcx
movq %r11, %rax
subq %rcx, %rax
cqto
idivq %r15
incq %rax
imulq %r15, %rax
addq %rcx, %rax
decq %rax
cmpq %r13, %rax
cmovlq %rax, %r13
movl 0x1c(%rbx), %eax
movl 0x10(%rbx), %r14d
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 %r15
pushq %r10
pushq %r13
pushq %r11
pushq %r14
pushq %rax
movl 0x108(%rsp), %eax
pushq %rax
pushq 0x108(%rsp)
callq 0x8a4b1
addq $0x50, %rsp
movq 0x48(%rbx), %r14
movq 0x58(%rbx), %rcx
testl %eax, %eax
je 0x5a296
cmpq %r12, %r13
jge 0x5a296
movq %r14, %rax
shlq $0x5, %rax
movslq -0x14(%rcx,%rax), %r10
incq %r10
cmpq %rbp, %r10
jl 0x5a1e4
movq %r14, %rax
shlq $0x5, %rax
movl -0x1c(%rcx,%rax), %eax
movq 0xd0(%rsp), %rdi
incl %eax
movslq %eax, %r12
jmp 0x5a2b3
xorl %r14d, %r14d
cmpq %rdi, %r12
jge 0x5a3b4
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 0x8a9ca
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 %r15
movq %rcx, %r13
subq %rdx, %r13
movq %r10, %rsi
subq %r13, %rsi
movq %rsi, %rax
cqto
idivq %r15
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 0x5a3b4
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 %r15
pushq %rax
pushq %r10
pushq %r13
pushq %r11
pushq %rbp
movl 0x108(%rsp), %eax
pushq %rax
pushq 0x108(%rsp)
callq 0x8a4b1
addq $0x50, %rsp
testl %eax, %eax
je 0x5a3b4
movq 0x48(%rbx), %rax
movq 0x58(%rbx), %rcx
shlq $0x5, %rax
movl -0x18(%rcx,%rax), %edi
cmpq %r12, %r13
jg 0x5a305
movq 0xd8(%rsp), %rsi
movl %esi, %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 0x5a4a3
movq 0x48(%rbx), %rcx
movq %rcx, %rax
subq %r14, %rax
jbe 0x5a455
cmpq $0x2, %rax
jb 0x5a455
movq 0x58(%rbx), %rdx
shlq $0x5, %r14
addq %rdx, %r14
shrq %rax
shlq $0x5, %rcx
addq %rdx, %rcx
addq $-0x20, %rcx
movups (%r14), %xmm0
movups 0x10(%r14), %xmm1
movaps %xmm1, 0x70(%rsp)
movaps %xmm0, 0x60(%rsp)
movups (%rcx), %xmm0
movups 0x10(%rcx), %xmm1
movups %xmm1, 0x10(%r14)
movups %xmm0, (%r14)
movapd 0x60(%rsp), %xmm0
movapd 0x70(%rsp), %xmm1
movupd %xmm1, 0x10(%rcx)
movupd %xmm0, (%rcx)
addq $0x20, %r14
addq $-0x20, %rcx
decq %rax
jne 0x5a40c
movq 0x48(%rbx), %rcx
cmpq 0x50(%rbx), %rcx
jne 0x5a486
leaq (%rcx,%rcx), %rax
testq %rcx, %rcx
movl $0x2, %esi
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 0x5a48a
movq 0x58(%rbx), %rax
leaq 0x1(%rcx), %rdx
movq %rdx, 0x48(%rbx)
shlq $0x5, %rcx
addq %rax, %rcx
movl $0x1, %eax
jmp 0x5a561
movl 0x4(%rbx), %edx
movl 0x8(%rbx), %ecx
incl %ecx
cmpl %edx, %ecx
je 0x5a517
subl %edx, %ecx
notq %rsi
addq %rdx, %rsi
addq %rax, %rsi
addq %rcx, %rsi
je 0x5a517
movl %ecx, %eax
xorps %xmm0, %xmm0
cvtsi2sd %rax, %xmm0
movsd 0x1d6e3f(%rip), %xmm2   # 0x231310
subsd 0x8(%rsp), %xmm2
mulsd %xmm0, %xmm2
cvtsi2sd %rsi, %xmm1
ucomisd %xmm2, %xmm1
jae 0x5a3e1
cmpl $0x0, 0x110(%rsp)
je 0x5a507
cmpq $0x1c2, %rsi            # imm = 0x1C2
jl 0x5a517
mulsd 0x1d70eb(%rip), %xmm0   # 0x2315f0
jmp 0x5a50d
mulsd 0x58(%rsp), %xmm0
ucomisd %xmm0, %xmm1
jae 0x5a3e1
movq 0x48(%rbx), %rcx
cmpq 0x50(%rbx), %rcx
jne 0x5a54c
leaq (%rcx,%rcx), %rax
testq %rcx, %rcx
movl $0x2, %esi
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 0x5a550
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