// File: chhylp123[P]hifiasm/Correct.cpp
// Function: retrieve_cigar_err(bit_extz_t*, long, long, long*, long*)
int64_t retrieve_cigar_err(bit_extz_t *ez, int64_t s, int64_t e, int64_t *xk, int64_t *ck)
{
    if(!ez->cigar.n) return 0;
    int64_t cn = ez->cigar.n, op, err = 0; int64_t ws, we, os, oe, ovlp;
    if(((*ck) < 0) || ((*ck) > cn)) {//(*ck) == cn is allowed
        (*ck) = 0; (*xk) = ez->ts; 
    }

    while ((*ck) > 0 && (*xk) > s) {
        --(*ck);
        op = ez->cigar.a[(*ck)]>>14;
        if(op!=2) (*xk) -= (ez->cigar.a[(*ck)]&(0x3fff));
    }

    //some cigar will span s or e
    while ((*ck) < cn && (*xk) < e) {//[s, e)
        ws = (*xk); 
        op = ez->cigar.a[(*ck)]>>14;
        if(op!=2) (*xk) += (ez->cigar.a[(*ck)]&(0x3fff));
        we = (*xk);
        os = MAX(s, ws); oe = MIN(e, we);
        ovlp = ((oe>os)? (oe-os):0);
        if((op==2) && (ws>=s) && (ws<e)) {
            ovlp = (ez->cigar.a[(*ck)]&(0x3fff));
        }

        // if(s == 22694 && e == 38018) {
        //     fprintf(stderr, "[M::%s]\tw::[%ld,\t%ld)\tc::%ld\top::%ld\terr::%ld\n", __func__, ws, we, (*ck), op, err + (op?ovlp:0));
        // } 
        (*ck)++;
        if((!ovlp) || (!op)) continue;
        err += ovlp;
    }

    // int64_t debug_err = retrieve_cigar_err_debug(ez, s, e);
    // if(!(err == debug_err)) {
    //     fprintf(stderr, "[M::%s::] err::%ld, debug_err::%ld, s::%ld, e::%ld, ez->ts::%u, ez->te::%u\n", 
    //     __func__, err, debug_err, s, e, ez->ts, ez->te);
    // }
    // assert(err == debug_err);
    return err;
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x10, %rsp
movq %rdi, -0x10(%rbp)
movq %rsi, -0x18(%rbp)
movq %rdx, -0x20(%rbp)
movq %rcx, -0x28(%rbp)
movq %r8, -0x30(%rbp)
movq -0x10(%rbp), %rax
cmpq $0x0, 0xa8(%rax)
jne 0x6b4a7
movq $0x0, -0x8(%rbp)
jmp 0x6b70b
movq -0x10(%rbp), %rax
movq 0xa8(%rax), %rax
movq %rax, -0x38(%rbp)
movq $0x0, -0x48(%rbp)
movq -0x30(%rbp), %rax
cmpq $0x0, (%rax)
jl 0x6b4d5
movq -0x30(%rbp), %rax
movq (%rax), %rax
cmpq -0x38(%rbp), %rax
jle 0x6b4ef
movq -0x30(%rbp), %rax
movq $0x0, (%rax)
movq -0x10(%rbp), %rax
movslq 0x1c(%rax), %rcx
movq -0x28(%rbp), %rax
movq %rcx, (%rax)
jmp 0x6b4f1
movq -0x30(%rbp), %rcx
xorl %eax, %eax
cmpq $0x0, (%rcx)
movb %al, -0x71(%rbp)
jle 0x6b511
movq -0x28(%rbp), %rax
movq (%rax), %rax
cmpq -0x18(%rbp), %rax
setg %al
movb %al, -0x71(%rbp)
movb -0x71(%rbp), %al
testb $0x1, %al
jne 0x6b51a
jmp 0x6b57e
movq -0x30(%rbp), %rax
movq (%rax), %rcx
addq $-0x1, %rcx
movq %rcx, (%rax)
movq -0x10(%rbp), %rax
movq 0xb8(%rax), %rax
movq -0x30(%rbp), %rcx
movq (%rcx), %rcx
movzwl (%rax,%rcx,2), %eax
sarl $0xe, %eax
cltq
movq %rax, -0x40(%rbp)
cmpq $0x2, -0x40(%rbp)
je 0x6b579
movq -0x10(%rbp), %rax
movq 0xb8(%rax), %rax
movq -0x30(%rbp), %rcx
movq (%rcx), %rcx
movzwl (%rax,%rcx,2), %eax
andl $0x3fff, %eax           # imm = 0x3FFF
movslq %eax, %rdx
movq -0x28(%rbp), %rax
movq (%rax), %rcx
subq %rdx, %rcx
movq %rcx, (%rax)
jmp 0x6b4f1
jmp 0x6b580
movq -0x30(%rbp), %rax
movq (%rax), %rcx
xorl %eax, %eax
cmpq -0x38(%rbp), %rcx
movb %al, -0x72(%rbp)
jge 0x6b5a3
movq -0x28(%rbp), %rax
movq (%rax), %rax
cmpq -0x20(%rbp), %rax
setl %al
movb %al, -0x72(%rbp)
movb -0x72(%rbp), %al
testb $0x1, %al
jne 0x6b5af
jmp 0x6b703
movq -0x28(%rbp), %rax
movq (%rax), %rax
movq %rax, -0x50(%rbp)
movq -0x10(%rbp), %rax
movq 0xb8(%rax), %rax
movq -0x30(%rbp), %rcx
movq (%rcx), %rcx
movzwl (%rax,%rcx,2), %eax
sarl $0xe, %eax
cltq
movq %rax, -0x40(%rbp)
cmpq $0x2, -0x40(%rbp)
je 0x6b608
movq -0x10(%rbp), %rax
movq 0xb8(%rax), %rax
movq -0x30(%rbp), %rcx
movq (%rcx), %rcx
movzwl (%rax,%rcx,2), %eax
andl $0x3fff, %eax           # imm = 0x3FFF
movslq %eax, %rcx
movq -0x28(%rbp), %rax
addq (%rax), %rcx
movq %rcx, (%rax)
movq -0x28(%rbp), %rax
movq (%rax), %rax
movq %rax, -0x58(%rbp)
movq -0x18(%rbp), %rax
cmpq -0x50(%rbp), %rax
jl 0x6b627
movq -0x18(%rbp), %rax
movq %rax, -0x80(%rbp)
jmp 0x6b62f
movq -0x50(%rbp), %rax
movq %rax, -0x80(%rbp)
movq -0x80(%rbp), %rax
movq %rax, -0x60(%rbp)
movq -0x20(%rbp), %rax
cmpq -0x58(%rbp), %rax
jg 0x6b64e
movq -0x20(%rbp), %rax
movq %rax, -0x88(%rbp)
jmp 0x6b659
movq -0x58(%rbp), %rax
movq %rax, -0x88(%rbp)
movq -0x88(%rbp), %rax
movq %rax, -0x68(%rbp)
movq -0x68(%rbp), %rax
cmpq -0x60(%rbp), %rax
jle 0x6b67f
movq -0x68(%rbp), %rax
subq -0x60(%rbp), %rax
movq %rax, -0x90(%rbp)
jmp 0x6b68a
xorl %eax, %eax
movq %rax, -0x90(%rbp)
jmp 0x6b68a
movq -0x90(%rbp), %rax
movq %rax, -0x70(%rbp)
cmpq $0x2, -0x40(%rbp)
jne 0x6b6d1
movq -0x50(%rbp), %rax
cmpq -0x18(%rbp), %rax
jl 0x6b6d1
movq -0x50(%rbp), %rax
cmpq -0x20(%rbp), %rax
jge 0x6b6d1
movq -0x10(%rbp), %rax
movq 0xb8(%rax), %rax
movq -0x30(%rbp), %rcx
movq (%rcx), %rcx
movzwl (%rax,%rcx,2), %eax
andl $0x3fff, %eax           # imm = 0x3FFF
cltq
movq %rax, -0x70(%rbp)
movq -0x30(%rbp), %rax
movq (%rax), %rcx
addq $0x1, %rcx
movq %rcx, (%rax)
cmpq $0x0, -0x70(%rbp)
je 0x6b6ed
cmpq $0x0, -0x40(%rbp)
jne 0x6b6f2
jmp 0x6b580
movq -0x70(%rbp), %rax
addq -0x48(%rbp), %rax
movq %rax, -0x48(%rbp)
jmp 0x6b580
movq -0x48(%rbp), %rax
movq %rax, -0x8(%rbp)
movq -0x8(%rbp), %rax
addq $0x10, %rsp
popq %rbp
retq
nopw %cs:(%rax,%rax)