// File: /workspace/github_cpp_makefile_star2_O3/chhylp123[P]hifiasm/hic.cpp
// Function: get_LCA(unsigned int, unsigned long, unsigned int, unsigned long, unsigned char*, unsigned long, MT*, bubble_type*, unsigned long*)
uint64_t get_LCA(uint32_t x, uint64_t xLen, uint32_t y, uint64_t yLen, uint8_t* dis, uint64_t n, MT* M, bubble_type* bub, uint64_t* min_rev)
{
    if(IF_BUB(x>>1, *bub) && IF_BUB(y>>1, *bub) && bub->index[x>>1] == bub->index[y>>1])
    {
        return get_LCA_bubble(x, xLen, y, yLen, dis, n, M, bub, min_rev);
    }
    else
    {
        memset(dis, -1, sizeof(uint8_t)*n);
    }

    uint32_t j, v;
    uint64_t u, d = (uint64_t)-1, tmp;
    uint8_t rev;
    
    v = x;
    for (j = 0; j < M->matrix.a[v].a.n; j++)
    {
        u = M->matrix.a[v].a.a[j] >> M->uID_shift;
        d = M->matrix.a[v].a.a[j] & M->dis_mode;
        dis[u] = dis[u] >> 4;
    }

    v = y;
    for (j = 0; j < M->matrix.a[v].a.n; j++)
    {
        u = M->matrix.a[v].a.a[j] >> M->uID_shift;
        d = M->matrix.a[v].a.a[j] & M->dis_mode;
        dis[u] = dis[u] >> 4;
    }

    uint64_t x_i = 0, y_i = 0, d_x, d_y, min_d = (uint64_t)-1;
    uint32_t min_j = (uint32_t)-1;
    (*min_rev) = (uint64_t)-1;
    for (j = 0; j < n; j++)
    {
        if(dis[j] != 0)
        {
            dis[j] = (uint8_t)-1;
            continue;
        }

        for (; x_i < M->matrix.a[x].a.n; x_i++)
        {
            u = M->matrix.a[x].a.a[x_i] >> M->uID_shift;
            d = M->matrix.a[x].a.a[x_i] & M->dis_mode;
            if(u == j) break;
        }
        if(x_i == M->matrix.a[x].a.n && M->matrix.a[x].a.n != 0) fprintf(stderr, "ERROR X\n");
        d_x = d;

        for (; y_i < M->matrix.a[y].a.n; y_i++)
        {
            u = M->matrix.a[y].a.a[y_i] >> M->uID_shift;
            d = M->matrix.a[y].a.a[y_i] & M->dis_mode;
            if(u == j) break;
        }
        if(y_i == M->matrix.a[y].a.n && M->matrix.a[y].a.n != 0) fprintf(stderr, "ERROR Y\n");
        d_y = d;


        tmp = LCA_distance(d_x, d_y, xLen, yLen, &rev);
        if(tmp < min_d) min_d = tmp, (*min_rev) = rev, min_j = j;
    }

    if(min_j == x || min_j == y) return (uint64_t)-1;

    return min_d;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x58, %rsp
movl %edx, %ebx
movl %edi, %r14d
movq 0xa0(%rsp), %rax
movq %rax, 0x48(%rsp)
movq 0x98(%rsp), %rdx
movq 0x90(%rsp), %rbp
movq (%rdx), %rax
movq 0x80(%rdx), %rdx
shrl %edi
movl (%rax,%rdi,4), %edi
incq %rdx
cmpq %rdi, %rdx
jbe 0x142a2f
movl %ebx, %edx
shrl %edx
cmpl (%rax,%rdx,4), %edi
jne 0x142a2f
movl %r14d, %edi
movl %ebx, %edx
addq $0x58, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x1426b0
movq %rsi, 0x10(%rsp)
movq %rcx, 0x20(%rsp)
movq %r8, 0x28(%rsp)
movq %r8, %rdi
movl $0xff, %esi
movq %r9, 0x18(%rsp)
movq %r9, %rdx
callq 0x20e0
movl %r14d, 0x40(%rsp)
movl %r14d, %edi
movq 0x10(%rbp), %rax
shlq $0x5, %rdi
cmpq $0x0, (%rax,%rdi)
je 0x142aa4
movl $0x1, %edx
xorl %ecx, %ecx
movq 0x28(%rsp), %r8
movq 0x10(%rax,%rdi), %rax
movq (%rax,%rcx,8), %rsi
movb 0x18(%rbp), %cl
movq %rsi, %rax
shrq %cl, %rax
movq 0x20(%rbp), %r14
shrb $0x4, (%r8,%rax)
movl %edx, %ecx
movq 0x10(%rbp), %rax
incl %edx
cmpq %rcx, (%rax,%rdi)
ja 0x142a76
andq %rsi, %r14
jmp 0x142ab0
movq $-0x1, %r14
movq 0x28(%rsp), %r8
movq 0x20(%rsp), %r9
movq 0x18(%rsp), %r10
movq 0x10(%rsp), %r11
movl %ebx, %r15d
shlq $0x5, %r15
cmpq $0x0, (%rax,%r15)
je 0x142b00
movl $0x1, %edx
xorl %ecx, %ecx
movq 0x10(%rax,%r15), %rax
movq (%rax,%rcx,8), %rsi
movb 0x18(%rbp), %cl
movq %rsi, %rax
shrq %cl, %rax
movq 0x20(%rbp), %r14
shrb $0x4, (%r8,%rax)
movl %edx, %ecx
movq 0x10(%rbp), %rax
incl %edx
cmpq %rcx, (%rax,%r15)
ja 0x142ad4
andq %rsi, %r14
movl %ebx, 0x44(%rsp)
movq 0x48(%rsp), %rax
movq $-0x1, (%rax)
testq %r10, %r10
je 0x142ca3
movq %r15, 0x50(%rsp)
xorl %r12d, %r12d
movq $-0x1, 0x30(%rsp)
movl $0xffffffff, 0xc(%rsp)  # imm = 0xFFFFFFFF
xorl %ebx, %ebx
xorl %r15d, %r15d
movq %rdi, 0x38(%rsp)
cmpb $0x0, (%r8,%r12)
je 0x142b4d
movb $-0x1, (%r8,%r12)
jmp 0x142c30
movq 0x10(%rbp), %rax
movq (%rax,%rdi), %rdx
cmpq %rdx, %r15
jae 0x142b8b
leaq (%rax,%rdi), %rcx
movq 0x10(%rcx), %rsi
movq 0x18(%rbp), %rcx
movq (%rsi,%r15,8), %r13
movq %r13, %rdi
shrq %cl, %rdi
cmpq %r12, %rdi
je 0x142b80
incq %r15
cmpq %r15, %rdx
jne 0x142b66
movq %rdx, %r15
andq 0x20(%rbp), %r13
movq 0x38(%rsp), %rdi
jmp 0x142b8e
movq %r14, %r13
cmpq %rdx, %r15
setne %cl
testq %rdx, %rdx
sete %dl
orb %cl, %dl
je 0x142c3e
movq 0x50(%rsp), %rcx
movq (%rax,%rcx), %rdx
cmpq %rdx, %rbx
jae 0x142bdb
addq %rcx, %rax
movq 0x10(%rax), %rax
movq 0x18(%rbp), %rcx
movq (%rax,%rbx,8), %r14
movq %r14, %rsi
shrq %cl, %rsi
cmpq %r12, %rsi
je 0x142bd5
incq %rbx
cmpq %rbx, %rdx
jne 0x142bbb
movq %rdx, %rbx
andq 0x20(%rbp), %r14
jmp 0x142bde
movq %r13, %r14
cmpq %rdx, %rbx
setne %al
testq %rdx, %rdx
sete %cl
orb %al, %cl
je 0x142c80
movq %r11, %rax
subq %r13, %rax
movq %r9, %rcx
subq %r14, %rcx
cmpq %r14, %r13
cmovlq %rax, %rcx
movq %r14, %rax
cmovgq %r13, %rax
addq %rcx, %rax
cmpq 0x30(%rsp), %rax
jae 0x142c30
xorl %ecx, %ecx
cmpq %r14, %r13
setl %cl
movq 0x48(%rsp), %rdx
movq %rcx, (%rdx)
movq %rax, 0x30(%rsp)
movl %r12d, 0xc(%rsp)
incl %r12d
cmpq %r10, %r12
jb 0x142b3c
jmp 0x142cb4
movq 0x15239b(%rip), %rax    # 0x294fe0
movq (%rax), %rcx
movl $0x8, %esi
movl $0x1, %edx
leaq 0xf8aae(%rip), %rdi     # 0x23b707
callq 0x2410
movq 0x38(%rsp), %rdi
movq 0x10(%rsp), %r11
movq 0x18(%rsp), %r10
movq 0x20(%rsp), %r9
movq 0x28(%rsp), %r8
movq 0x10(%rbp), %rax
jmp 0x142ba2
callq 0x1c3af
movq 0x38(%rsp), %rdi
movq 0x10(%rsp), %r11
movq 0x18(%rsp), %r10
movq 0x20(%rsp), %r9
movq 0x28(%rsp), %r8
jmp 0x142bf2
movl $0xffffffff, 0xc(%rsp)  # imm = 0xFFFFFFFF
movq $-0x1, 0x30(%rsp)
movl 0xc(%rsp), %edx
cmpl 0x44(%rsp), %edx
movq $-0x1, %rcx
movq 0x30(%rsp), %rax
cmoveq %rcx, %rax
cmpl 0x40(%rsp), %edx
cmoveq %rcx, %rax
addq $0x58, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq