// File: chhylp123[P]hifiasm/hic.cpp
// Function: LeastSquare_advance(trans_idx*, ha_ug_index*, unsigned long)
void LeastSquare_advance(trans_idx* dis, ha_ug_index* idx, uint64_t med)
{
    #define SCAL_RATE 1000
    long double t1=0, t2=0, t3=0, t4=0, x, y;
    uint64_t i, m, ava_size;

    for (i = m = 0; i < dis->n; i++)
    {
        x = ((double)(dis->a[i].beg + dis->a[i].end))/2;
        y = ((double)(dis->a[i].cnt_1))/((double)(dis->a[i].cnt_0 + dis->a[i].cnt_1));
        if(dis->a[i].beg >= med) break;

        t1 += x*x;
        t2 += x;
        t3 += x*y;
        t4 += y;
        m++;
    }

    if(i < dis->n)
    {
        uint64_t beg, end, cnt_0, cnt_1;
        for (beg = dis->a[i].beg, end = dis->a[i].end, cnt_0 = cnt_1 = 0; i < dis->n; i++)
        {
            cnt_0 += dis->a[i].cnt_0;
            cnt_1 += dis->a[i].cnt_1;
            beg = MIN(beg, dis->a[i].beg);
            end = MAX(end, dis->a[i].end);
        }

        x = ((double)(beg + end))/2;
        y = ((double)(cnt_1))/((double)(cnt_0 + cnt_1));

        t1 += x*x;
        t2 += x;
        t3 += x*y;
        t4 += y;
        m++;
    }



    if(t2 > t4)
    {
        idx->frac = t2/t4;
        if(idx->frac > SCAL_RATE) idx->frac = idx->frac / SCAL_RATE;
        if(idx->frac < 1) idx->frac = 1;
    }

    ava_size = m;
    t1 /= (idx->frac*idx->frac);
    t2 /= idx->frac;
    t3 /= idx->frac;
    if((t1*ava_size - t2*t2) != 0)
    {
        idx->a = (t3*ava_size - t2*t4) / (t1*ava_size - t2*t2);
    }
    if((t1*ava_size - t2*t2) != 0)
    {
        idx->b = (t1*t4 - t2*t3) / (t1*ava_size - t2*t2);
    }
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x90, %rsp
movq %rdi, -0x8(%rbp)
movq %rsi, -0x10(%rbp)
movq %rdx, -0x18(%rbp)
fldz
fld %st(0)
fstpt -0x30(%rbp)
fld %st(0)
fstpt -0x40(%rbp)
fld %st(0)
fstpt -0x50(%rbp)
fstpt -0x60(%rbp)
movq $0x0, -0x90(%rbp)
movq $0x0, -0x88(%rbp)
movq -0x88(%rbp), %rax
movq -0x8(%rbp), %rcx
cmpq 0x8(%rcx), %rax
jae 0x1dbc29
movq -0x8(%rbp), %rax
movq (%rax), %rcx
movq -0x88(%rbp), %rax
leaq (%rax,%rax,4), %rdx
movq (%rcx,%rdx,8), %rax
movq 0x8(%rcx,%rdx,8), %rcx
addq %rcx, %rax
movq %rax, %xmm0
movaps 0x1a6500(%rip), %xmm3   # 0x382010
punpckldq %xmm3, %xmm0    # xmm0 = xmm0[0],xmm3[0],xmm0[1],xmm3[1]
movapd 0x1a6504(%rip), %xmm2   # 0x382020
subpd %xmm2, %xmm0
movaps %xmm0, %xmm1
unpckhpd %xmm0, %xmm0            # xmm0 = xmm0[1,1]
addsd %xmm1, %xmm0
movsd 0x1a673d(%rip), %xmm1   # 0x382270
mulsd %xmm1, %xmm0
movsd %xmm0, -0xc0(%rbp)
fldl -0xc0(%rbp)
fstpt -0x70(%rbp)
movq -0x8(%rbp), %rax
movq (%rax), %rcx
movq -0x88(%rbp), %rax
leaq (%rax,%rax,4), %rdx
movq 0x18(%rcx,%rdx,8), %rax
movq 0x20(%rcx,%rdx,8), %rcx
movq %rcx, %xmm0
punpckldq %xmm3, %xmm0    # xmm0 = xmm0[0],xmm3[0],xmm0[1],xmm3[1]
subpd %xmm2, %xmm0
movaps %xmm0, %xmm1
unpckhpd %xmm0, %xmm0            # xmm0 = xmm0[1,1]
addsd %xmm1, %xmm0
addq %rcx, %rax
movq %rax, %xmm1
punpckldq %xmm3, %xmm1    # xmm1 = xmm1[0],xmm3[0],xmm1[1],xmm3[1]
subpd %xmm2, %xmm1
movaps %xmm1, %xmm2
unpckhpd %xmm1, %xmm1            # xmm1 = xmm1[1,1]
addsd %xmm2, %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, -0xc8(%rbp)
fldl -0xc8(%rbp)
fstpt -0x80(%rbp)
movq -0x8(%rbp), %rax
movq (%rax), %rax
imulq $0x28, -0x88(%rbp), %rcx
addq %rcx, %rax
movq (%rax), %rax
cmpq -0x18(%rbp), %rax
jb 0x1dbbc9
jmp 0x1dbc29
fldt -0x70(%rbp)
fldt -0x30(%rbp)
fxch %st(1)
fmul %st, %st(0)
faddp %st, %st(1)
fstpt -0x30(%rbp)
fldt -0x70(%rbp)
fldt -0x40(%rbp)
faddp %st, %st(1)
fstpt -0x40(%rbp)
fldt -0x70(%rbp)
fldt -0x80(%rbp)
fldt -0x50(%rbp)
fxch %st(2)
fmulp %st, %st(1)
faddp %st, %st(1)
fstpt -0x50(%rbp)
fldt -0x80(%rbp)
fldt -0x60(%rbp)
faddp %st, %st(1)
fstpt -0x60(%rbp)
movq -0x90(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x90(%rbp)
movq -0x88(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x88(%rbp)
jmp 0x1dbad1
movq -0x88(%rbp), %rax
movq -0x8(%rbp), %rcx
cmpq 0x8(%rcx), %rax
jae 0x1dbea7
movq -0x8(%rbp), %rax
movq (%rax), %rax
imulq $0x28, -0x88(%rbp), %rcx
addq %rcx, %rax
movq (%rax), %rax
movq %rax, -0xa0(%rbp)
movq -0x8(%rbp), %rax
movq (%rax), %rax
imulq $0x28, -0x88(%rbp), %rcx
addq %rcx, %rax
movq 0x8(%rax), %rax
movq %rax, -0xa8(%rbp)
movq $0x0, -0xb8(%rbp)
movq $0x0, -0xb0(%rbp)
movq -0x88(%rbp), %rax
movq -0x8(%rbp), %rcx
cmpq 0x8(%rcx), %rax
jae 0x1dbdb3
movq -0x8(%rbp), %rax
movq (%rax), %rax
imulq $0x28, -0x88(%rbp), %rcx
addq %rcx, %rax
movq 0x18(%rax), %rax
addq -0xb0(%rbp), %rax
movq %rax, -0xb0(%rbp)
movq -0x8(%rbp), %rax
movq (%rax), %rax
imulq $0x28, -0x88(%rbp), %rcx
addq %rcx, %rax
movq 0x20(%rax), %rax
addq -0xb8(%rbp), %rax
movq %rax, -0xb8(%rbp)
movq -0xa0(%rbp), %rax
movq -0x8(%rbp), %rcx
movq (%rcx), %rcx
imulq $0x28, -0x88(%rbp), %rdx
addq %rdx, %rcx
cmpq (%rcx), %rax
ja 0x1dbd18
movq -0xa0(%rbp), %rax
movq %rax, -0x100(%rbp)
jmp 0x1dbd34
movq -0x8(%rbp), %rax
movq (%rax), %rax
imulq $0x28, -0x88(%rbp), %rcx
addq %rcx, %rax
movq (%rax), %rax
movq %rax, -0x100(%rbp)
movq -0x100(%rbp), %rax
movq %rax, -0xa0(%rbp)
movq -0xa8(%rbp), %rax
movq -0x8(%rbp), %rcx
movq (%rcx), %rcx
imulq $0x28, -0x88(%rbp), %rdx
addq %rdx, %rcx
cmpq 0x8(%rcx), %rax
jb 0x1dbd71
movq -0xa8(%rbp), %rax
movq %rax, -0x108(%rbp)
jmp 0x1dbd8e
movq -0x8(%rbp), %rax
movq (%rax), %rax
imulq $0x28, -0x88(%rbp), %rcx
addq %rcx, %rax
movq 0x8(%rax), %rax
movq %rax, -0x108(%rbp)
movq -0x108(%rbp), %rax
movq %rax, -0xa8(%rbp)
movq -0x88(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x88(%rbp)
jmp 0x1dbc8d
movq -0xa0(%rbp), %rax
movq -0xa8(%rbp), %rcx
addq %rcx, %rax
movq %rax, %xmm0
movaps 0x1a6240(%rip), %xmm3   # 0x382010
punpckldq %xmm3, %xmm0    # xmm0 = xmm0[0],xmm3[0],xmm0[1],xmm3[1]
movapd 0x1a6244(%rip), %xmm2   # 0x382020
subpd %xmm2, %xmm0
movaps %xmm0, %xmm1
unpckhpd %xmm0, %xmm0            # xmm0 = xmm0[1,1]
addsd %xmm1, %xmm0
movsd 0x1a647d(%rip), %xmm1   # 0x382270
mulsd %xmm1, %xmm0
movsd %xmm0, -0xd0(%rbp)
fldl -0xd0(%rbp)
fstpt -0x70(%rbp)
movq -0xb8(%rbp), %rcx
movq %rcx, %xmm0
punpckldq %xmm3, %xmm0    # xmm0 = xmm0[0],xmm3[0],xmm0[1],xmm3[1]
subpd %xmm2, %xmm0
movaps %xmm0, %xmm1
unpckhpd %xmm0, %xmm0            # xmm0 = xmm0[1,1]
addsd %xmm1, %xmm0
movq -0xb0(%rbp), %rax
addq %rcx, %rax
movq %rax, %xmm1
punpckldq %xmm3, %xmm1    # xmm1 = xmm1[0],xmm3[0],xmm1[1],xmm3[1]
subpd %xmm2, %xmm1
movaps %xmm1, %xmm2
unpckhpd %xmm1, %xmm1            # xmm1 = xmm1[1,1]
addsd %xmm2, %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, -0xd8(%rbp)
fldl -0xd8(%rbp)
fstpt -0x80(%rbp)
fldt -0x70(%rbp)
fldt -0x30(%rbp)
fxch %st(1)
fmul %st, %st(0)
faddp %st, %st(1)
fstpt -0x30(%rbp)
fldt -0x70(%rbp)
fldt -0x40(%rbp)
faddp %st, %st(1)
fstpt -0x40(%rbp)
fldt -0x70(%rbp)
fldt -0x80(%rbp)
fldt -0x50(%rbp)
fxch %st(2)
fmulp %st, %st(1)
faddp %st, %st(1)
fstpt -0x50(%rbp)
fldt -0x80(%rbp)
fldt -0x60(%rbp)
faddp %st, %st(1)
fstpt -0x60(%rbp)
movq -0x90(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x90(%rbp)
fldt -0x40(%rbp)
fldt -0x60(%rbp)
fxch %st(1)
fucompi %st(1), %st
fstp %st(0)
jbe 0x1dbf1d
jmp 0x1dbeb7
fldt -0x40(%rbp)
fldt -0x60(%rbp)
fdivrp %st, %st(1)
movq -0x10(%rbp), %rax
fstpt 0xa0(%rax)
movq -0x10(%rbp), %rax
fldt 0xa0(%rax)
flds 0x1affdf(%rip)          # 0x38beb8
fxch %st(1)
fucompi %st(1), %st
fstp %st(0)
jbe 0x1dbefb
jmp 0x1dbee3
movq -0x10(%rbp), %rax
fldt 0xa0(%rax)
flds 0x1affc5(%rip)          # 0x38beb8
fdivrp %st, %st(1)
fstpt 0xa0(%rax)
movq -0x10(%rbp), %rax
fldt 0xa0(%rax)
fld1
fucompi %st(1), %st
fstp %st(0)
jbe 0x1dbf1b
jmp 0x1dbf0f
movq -0x10(%rbp), %rax
fld1
fstpt 0xa0(%rax)
jmp 0x1dbf1d
movq -0x90(%rbp), %rax
movq %rax, -0x98(%rbp)
movq -0x10(%rbp), %rax
fldt 0xa0(%rax)
fmul %st, %st(0)
fldt -0x30(%rbp)
fdivp %st, %st(1)
fstpt -0x30(%rbp)
movq -0x10(%rbp), %rax
fldt 0xa0(%rax)
fldt -0x40(%rbp)
fdivp %st, %st(1)
fstpt -0x40(%rbp)
movq -0x10(%rbp), %rax
fldt 0xa0(%rax)
fldt -0x50(%rbp)
fdivp %st, %st(1)
fstpt -0x50(%rbp)
fldt -0x30(%rbp)
movq -0x98(%rbp), %rax
movq %rax, -0xe0(%rbp)
testq %rax, %rax
sets %al
movzbl %al, %eax
movl %eax, %ecx
fildll -0xe0(%rbp)
leaq 0x1afef4(%rip), %rax    # 0x38be80
flds (%rax,%rcx,4)
faddp %st, %st(1)
fldt -0x40(%rbp)
fmul %st, %st(0)
fxch %st(2)
fmulp %st, %st(1)
fsubp %st, %st(1)
fldz
fxch %st(1)
fucompi %st(1), %st
fstp %st(0)
jne 0x1dbfaa
jp 0x1dbfaa
jmp 0x1dbfff
fldt -0x50(%rbp)
movq -0x98(%rbp), %rax
movq %rax, -0xe8(%rbp)
testq %rax, %rax
sets %al
movzbl %al, %eax
movl %eax, %ecx
fildll -0xe8(%rbp)
leaq 0x1afead(%rip), %rax    # 0x38be80
flds (%rax,%rcx,4)
faddp %st, %st(1)
fldt -0x40(%rbp)
fldt -0x60(%rbp)
fmul %st(1), %st
fxch %st(3)
fmul %st(2), %st
fsubp %st, %st(3)
fldt -0x30(%rbp)
fxch %st(1)
fmul %st, %st(0)
fxch %st(1)
fmulp %st, %st(2)
fsubrp %st, %st(1)
fdivrp %st, %st(1)
movq -0x10(%rbp), %rax
fstpt 0x80(%rax)
fldt -0x30(%rbp)
movq -0x98(%rbp), %rax
movq %rax, -0xf0(%rbp)
testq %rax, %rax
sets %al
movzbl %al, %eax
movl %eax, %ecx
fildll -0xf0(%rbp)
leaq 0x1afe58(%rip), %rax    # 0x38be80
flds (%rax,%rcx,4)
faddp %st, %st(1)
fldt -0x40(%rbp)
fmul %st, %st(0)
fxch %st(2)
fmulp %st, %st(1)
fsubp %st, %st(1)
fldz
fxch %st(1)
fucompi %st(1), %st
fstp %st(0)
jne 0x1dc046
jp 0x1dc046
jmp 0x1dc09b
fldt -0x30(%rbp)
fldt -0x60(%rbp)
fldt -0x40(%rbp)
fldt -0x50(%rbp)
fmul %st(1), %st
fxch %st(2)
fmul %st(3), %st
fsubp %st, %st(2)
movq -0x98(%rbp), %rax
movq %rax, -0xf8(%rbp)
testq %rax, %rax
sets %al
movzbl %al, %eax
movl %eax, %ecx
fildll -0xf8(%rbp)
leaq 0x1afe00(%rip), %rax    # 0x38be80
flds (%rax,%rcx,4)
faddp %st, %st(1)
fxch %st(1)
fmul %st, %st(0)
fxch %st(3)
fmulp %st, %st(1)
fsubp %st, %st(2)
fdivp %st, %st(1)
movq -0x10(%rbp), %rax
fstpt 0x90(%rax)
addq $0x90, %rsp
popq %rbp
retq
nopw %cs:(%rax,%rax)