// File: chhylp123[P]hifiasm/inter.cpp
// Function: retrieve_cigar_err_dir(bit_extz_t*, long, long, long*, long*, long)
int64_t retrieve_cigar_err_dir(bit_extz_t *ez, int64_t s, int64_t e, int64_t *xk, int64_t *ck, int64_t is_back)
{	///[ez->ts, ez->te]/[ez->qs, ez->qe]/[s, e)
    if(!ez->cigar.n) return 0;
    int64_t cn = ez->cigar.n, op, err = 0; int64_t ws, we, os, oe, ovlp;
	if(!is_back) {
		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));
			}
			(*ck)++;
			if((!ovlp) || (!op)) continue;
			err += ovlp;
		}
	} else {
		if(((*ck) < 0) || ((*ck) >= cn)) {
			(*ck) = cn-1; (*xk) = ez->te + 1; ///[ez->ts, ez->te]
		}
		while (((*ck) < cn) && ((*xk) < e)) {
			++(*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) >= 0 && (*xk) > s) {//[s, e)
			we = (*xk); 
			op = ez->cigar.a[(*ck)]>>14;
			if(op!=2) (*xk) -= (ez->cigar.a[(*ck)]&(0x3fff));
			ws = (*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));
			}
			(*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 $0x38, %rsp
movq %rdi, -0x10(%rbp)
movq %rsi, -0x18(%rbp)
movq %rdx, -0x20(%rbp)
movq %rcx, -0x28(%rbp)
movq %r8, -0x30(%rbp)
movq %r9, -0x38(%rbp)
movq -0x10(%rbp), %rax
cmpq $0x0, 0xa8(%rax)
jne 0x2c2a2b
movq $0x0, -0x8(%rbp)
jmp 0x2c2f13
movq -0x10(%rbp), %rax
movq 0xa8(%rax), %rax
movq %rax, -0x40(%rbp)
movq $0x0, -0x50(%rbp)
cmpq $0x0, -0x38(%rbp)
jne 0x2c2ca0
movq -0x30(%rbp), %rax
cmpq $0x0, (%rax)
jl 0x2c2a64
movq -0x30(%rbp), %rax
movq (%rax), %rax
cmpq -0x40(%rbp), %rax
jle 0x2c2a7e
movq -0x30(%rbp), %rax
movq $0x0, (%rax)
movq -0x10(%rbp), %rax
movslq 0x1c(%rax), %rcx
movq -0x28(%rbp), %rax
movq %rcx, (%rax)
jmp 0x2c2a80
movq -0x30(%rbp), %rcx
xorl %eax, %eax
cmpq $0x0, (%rcx)
movb %al, -0x79(%rbp)
jle 0x2c2aa0
movq -0x28(%rbp), %rax
movq (%rax), %rax
cmpq -0x18(%rbp), %rax
setg %al
movb %al, -0x79(%rbp)
movb -0x79(%rbp), %al
testb $0x1, %al
jne 0x2c2aa9
jmp 0x2c2b0d
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, -0x48(%rbp)
cmpq $0x2, -0x48(%rbp)
je 0x2c2b08
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 0x2c2a80
jmp 0x2c2b0f
movq -0x30(%rbp), %rax
movq (%rax), %rcx
xorl %eax, %eax
cmpq -0x40(%rbp), %rcx
movb %al, -0x7a(%rbp)
jge 0x2c2b32
movq -0x28(%rbp), %rax
movq (%rax), %rax
cmpq -0x20(%rbp), %rax
setl %al
movb %al, -0x7a(%rbp)
movb -0x7a(%rbp), %al
testb $0x1, %al
jne 0x2c2b3e
jmp 0x2c2c9b
movq -0x28(%rbp), %rax
movq (%rax), %rax
movq %rax, -0x58(%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, -0x48(%rbp)
cmpq $0x2, -0x48(%rbp)
je 0x2c2b97
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, -0x60(%rbp)
movq -0x18(%rbp), %rax
cmpq -0x58(%rbp), %rax
jl 0x2c2bb9
movq -0x18(%rbp), %rax
movq %rax, -0x88(%rbp)
jmp 0x2c2bc4
movq -0x58(%rbp), %rax
movq %rax, -0x88(%rbp)
movq -0x88(%rbp), %rax
movq %rax, -0x68(%rbp)
movq -0x20(%rbp), %rax
cmpq -0x60(%rbp), %rax
jg 0x2c2be6
movq -0x20(%rbp), %rax
movq %rax, -0x90(%rbp)
jmp 0x2c2bf1
movq -0x60(%rbp), %rax
movq %rax, -0x90(%rbp)
movq -0x90(%rbp), %rax
movq %rax, -0x70(%rbp)
movq -0x70(%rbp), %rax
cmpq -0x68(%rbp), %rax
jle 0x2c2c17
movq -0x70(%rbp), %rax
subq -0x68(%rbp), %rax
movq %rax, -0x98(%rbp)
jmp 0x2c2c22
xorl %eax, %eax
movq %rax, -0x98(%rbp)
jmp 0x2c2c22
movq -0x98(%rbp), %rax
movq %rax, -0x78(%rbp)
cmpq $0x2, -0x48(%rbp)
jne 0x2c2c69
movq -0x58(%rbp), %rax
cmpq -0x18(%rbp), %rax
jl 0x2c2c69
movq -0x58(%rbp), %rax
cmpq -0x20(%rbp), %rax
jge 0x2c2c69
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, -0x78(%rbp)
movq -0x30(%rbp), %rax
movq (%rax), %rcx
addq $0x1, %rcx
movq %rcx, (%rax)
cmpq $0x0, -0x78(%rbp)
je 0x2c2c85
cmpq $0x0, -0x48(%rbp)
jne 0x2c2c8a
jmp 0x2c2b0f
movq -0x78(%rbp), %rax
addq -0x50(%rbp), %rax
movq %rax, -0x50(%rbp)
jmp 0x2c2b0f
jmp 0x2c2f0b
movq -0x30(%rbp), %rax
cmpq $0x0, (%rax)
jl 0x2c2cb7
movq -0x30(%rbp), %rax
movq (%rax), %rax
cmpq -0x40(%rbp), %rax
jl 0x2c2cda
movq -0x40(%rbp), %rcx
subq $0x1, %rcx
movq -0x30(%rbp), %rax
movq %rcx, (%rax)
movq -0x10(%rbp), %rax
movl 0x20(%rax), %eax
addl $0x1, %eax
movslq %eax, %rcx
movq -0x28(%rbp), %rax
movq %rcx, (%rax)
jmp 0x2c2cdc
movq -0x30(%rbp), %rax
movq (%rax), %rcx
xorl %eax, %eax
cmpq -0x40(%rbp), %rcx
movb %al, -0x99(%rbp)
jge 0x2c2d05
movq -0x28(%rbp), %rax
movq (%rax), %rax
cmpq -0x20(%rbp), %rax
setl %al
movb %al, -0x99(%rbp)
movb -0x99(%rbp), %al
testb $0x1, %al
jne 0x2c2d11
jmp 0x2c2d72
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, -0x48(%rbp)
cmpq $0x2, -0x48(%rbp)
je 0x2c2d6d
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)
jmp 0x2c2cdc
jmp 0x2c2d74
movq -0x30(%rbp), %rcx
xorl %eax, %eax
cmpq $0x0, (%rcx)
movb %al, -0x9a(%rbp)
jl 0x2c2d9a
movq -0x28(%rbp), %rax
movq (%rax), %rax
cmpq -0x18(%rbp), %rax
setg %al
movb %al, -0x9a(%rbp)
movb -0x9a(%rbp), %al
testb $0x1, %al
jne 0x2c2da9
jmp 0x2c2f09
movq -0x28(%rbp), %rax
movq (%rax), %rax
movq %rax, -0x60(%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, -0x48(%rbp)
cmpq $0x2, -0x48(%rbp)
je 0x2c2e05
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)
movq -0x28(%rbp), %rax
movq (%rax), %rax
movq %rax, -0x58(%rbp)
movq -0x18(%rbp), %rax
cmpq -0x58(%rbp), %rax
jl 0x2c2e27
movq -0x18(%rbp), %rax
movq %rax, -0xa8(%rbp)
jmp 0x2c2e32
movq -0x58(%rbp), %rax
movq %rax, -0xa8(%rbp)
movq -0xa8(%rbp), %rax
movq %rax, -0x68(%rbp)
movq -0x20(%rbp), %rax
cmpq -0x60(%rbp), %rax
jg 0x2c2e54
movq -0x20(%rbp), %rax
movq %rax, -0xb0(%rbp)
jmp 0x2c2e5f
movq -0x60(%rbp), %rax
movq %rax, -0xb0(%rbp)
movq -0xb0(%rbp), %rax
movq %rax, -0x70(%rbp)
movq -0x70(%rbp), %rax
cmpq -0x68(%rbp), %rax
jle 0x2c2e85
movq -0x70(%rbp), %rax
subq -0x68(%rbp), %rax
movq %rax, -0xb8(%rbp)
jmp 0x2c2e90
xorl %eax, %eax
movq %rax, -0xb8(%rbp)
jmp 0x2c2e90
movq -0xb8(%rbp), %rax
movq %rax, -0x78(%rbp)
cmpq $0x2, -0x48(%rbp)
jne 0x2c2ed7
movq -0x58(%rbp), %rax
cmpq -0x18(%rbp), %rax
jl 0x2c2ed7
movq -0x58(%rbp), %rax
cmpq -0x20(%rbp), %rax
jge 0x2c2ed7
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, -0x78(%rbp)
movq -0x30(%rbp), %rax
movq (%rax), %rcx
addq $-0x1, %rcx
movq %rcx, (%rax)
cmpq $0x0, -0x78(%rbp)
je 0x2c2ef3
cmpq $0x0, -0x48(%rbp)
jne 0x2c2ef8
jmp 0x2c2d74
movq -0x78(%rbp), %rax
addq -0x50(%rbp), %rax
movq %rax, -0x50(%rbp)
jmp 0x2c2d74
jmp 0x2c2f0b
movq -0x50(%rbp), %rax
movq %rax, -0x8(%rbp)
movq -0x8(%rbp), %rax
addq $0x38, %rsp
popq %rbp
retq
nopl (%rax)