// File: /workspace/github_cpp_makefile_star2_O1/mouda[P]LTE_random_access-simulate/specialfunctions.cpp
// Function: alglib_impl::invincompletebeta(double, double, double, alglib_impl::ae_state*)
double invincompletebeta(double a, double b, double y, ae_state *_state)
{
    double aaa;
    double bbb;
    double y0;
    double d;
    double yyy;
    double x;
    double x0;
    double x1;
    double lgm;
    double yp;
    double di;
    double dithresh;
    double yl;
    double yh;
    double xt;
    ae_int_t i;
    ae_int_t rflg;
    ae_int_t dir;
    ae_int_t nflg;
    double s;
    ae_int_t mainlooppos;
    ae_int_t ihalve;
    ae_int_t ihalvecycle;
    ae_int_t newt;
    ae_int_t newtcycle;
    ae_int_t breaknewtcycle;
    ae_int_t breakihalvecycle;
    double result;


    i = 0;
    ae_assert(ae_fp_greater_eq(y,0)&&ae_fp_less_eq(y,1), "Domain error in InvIncompleteBeta", _state);
    
    /*
     * special cases
     */
    if( ae_fp_eq(y,0) )
    {
        result = 0;
        return result;
    }
    if( ae_fp_eq(y,1.0) )
    {
        result = 1;
        return result;
    }
    
    /*
     * these initializations are not really necessary,
     * but without them compiler complains about 'possibly uninitialized variables'.
     */
    dithresh = 0;
    rflg = 0;
    aaa = 0;
    bbb = 0;
    y0 = 0;
    x = 0;
    yyy = 0;
    lgm = 0;
    dir = 0;
    di = 0;
    
    /*
     * normal initializations
     */
    x0 = 0.0;
    yl = 0.0;
    x1 = 1.0;
    yh = 1.0;
    nflg = 0;
    mainlooppos = 0;
    ihalve = 1;
    ihalvecycle = 2;
    newt = 3;
    newtcycle = 4;
    breaknewtcycle = 5;
    breakihalvecycle = 6;
    
    /*
     * main loop
     */
    for(;;)
    {
        
        /*
         * start
         */
        if( mainlooppos==0 )
        {
            if( ae_fp_less_eq(a,1.0)||ae_fp_less_eq(b,1.0) )
            {
                dithresh = 1.0e-6;
                rflg = 0;
                aaa = a;
                bbb = b;
                y0 = y;
                x = aaa/(aaa+bbb);
                yyy = incompletebeta(aaa, bbb, x, _state);
                mainlooppos = ihalve;
                continue;
            }
            else
            {
                dithresh = 1.0e-4;
            }
            yp = -invnormaldistribution(y, _state);
            if( ae_fp_greater(y,0.5) )
            {
                rflg = 1;
                aaa = b;
                bbb = a;
                y0 = 1.0-y;
                yp = -yp;
            }
            else
            {
                rflg = 0;
                aaa = a;
                bbb = b;
                y0 = y;
            }
            lgm = (yp*yp-3.0)/6.0;
            x = 2.0/(1.0/(2.0*aaa-1.0)+1.0/(2.0*bbb-1.0));
            d = yp*ae_sqrt(x+lgm, _state)/x-(1.0/(2.0*bbb-1.0)-1.0/(2.0*aaa-1.0))*(lgm+5.0/6.0-2.0/(3.0*x));
            d = 2.0*d;
            if( ae_fp_less(d,ae_log(ae_minrealnumber, _state)) )
            {
                x = 0;
                break;
            }
            x = aaa/(aaa+bbb*ae_exp(d, _state));
            yyy = incompletebeta(aaa, bbb, x, _state);
            yp = (yyy-y0)/y0;
            if( ae_fp_less(ae_fabs(yp, _state),0.2) )
            {
                mainlooppos = newt;
                continue;
            }
            mainlooppos = ihalve;
            continue;
        }
        
        /*
         * ihalve
         */
        if( mainlooppos==ihalve )
        {
            dir = 0;
            di = 0.5;
            i = 0;
            mainlooppos = ihalvecycle;
            continue;
        }
        
        /*
         * ihalvecycle
         */
        if( mainlooppos==ihalvecycle )
        {
            if( i<=99 )
            {
                if( i!=0 )
                {
                    x = x0+di*(x1-x0);
                    if( ae_fp_eq(x,1.0) )
                    {
                        x = 1.0-ae_machineepsilon;
                    }
                    if( ae_fp_eq(x,0.0) )
                    {
                        di = 0.5;
                        x = x0+di*(x1-x0);
                        if( ae_fp_eq(x,0.0) )
                        {
                            break;
                        }
                    }
                    yyy = incompletebeta(aaa, bbb, x, _state);
                    yp = (x1-x0)/(x1+x0);
                    if( ae_fp_less(ae_fabs(yp, _state),dithresh) )
                    {
                        mainlooppos = newt;
                        continue;
                    }
                    yp = (yyy-y0)/y0;
                    if( ae_fp_less(ae_fabs(yp, _state),dithresh) )
                    {
                        mainlooppos = newt;
                        continue;
                    }
                }
                if( ae_fp_less(yyy,y0) )
                {
                    x0 = x;
                    yl = yyy;
                    if( dir<0 )
                    {
                        dir = 0;
                        di = 0.5;
                    }
                    else
                    {
                        if( dir>3 )
                        {
                            di = 1.0-(1.0-di)*(1.0-di);
                        }
                        else
                        {
                            if( dir>1 )
                            {
                                di = 0.5*di+0.5;
                            }
                            else
                            {
                                di = (y0-yyy)/(yh-yl);
                            }
                        }
                    }
                    dir = dir+1;
                    if( ae_fp_greater(x0,0.75) )
                    {
                        if( rflg==1 )
                        {
                            rflg = 0;
                            aaa = a;
                            bbb = b;
                            y0 = y;
                        }
                        else
                        {
                            rflg = 1;
                            aaa = b;
                            bbb = a;
                            y0 = 1.0-y;
                        }
                        x = 1.0-x;
                        yyy = incompletebeta(aaa, bbb, x, _state);
                        x0 = 0.0;
                        yl = 0.0;
                        x1 = 1.0;
                        yh = 1.0;
                        mainlooppos = ihalve;
                        continue;
                    }
                }
                else
                {
                    x1 = x;
                    if( rflg==1&&ae_fp_less(x1,ae_machineepsilon) )
                    {
                        x = 0.0;
                        break;
                    }
                    yh = yyy;
                    if( dir>0 )
                    {
                        dir = 0;
                        di = 0.5;
                    }
                    else
                    {
                        if( dir<-3 )
                        {
                            di = di*di;
                        }
                        else
                        {
                            if( dir<-1 )
                            {
                                di = 0.5*di;
                            }
                            else
                            {
                                di = (yyy-y0)/(yh-yl);
                            }
                        }
                    }
                    dir = dir-1;
                }
                i = i+1;
                mainlooppos = ihalvecycle;
                continue;
            }
            else
            {
                mainlooppos = breakihalvecycle;
                continue;
            }
        }
        
        /*
         * breakihalvecycle
         */
        if( mainlooppos==breakihalvecycle )
        {
            if( ae_fp_greater_eq(x0,1.0) )
            {
                x = 1.0-ae_machineepsilon;
                break;
            }
            if( ae_fp_less_eq(x,0.0) )
            {
                x = 0.0;
                break;
            }
            mainlooppos = newt;
            continue;
        }
        
        /*
         * newt
         */
        if( mainlooppos==newt )
        {
            if( nflg!=0 )
            {
                break;
            }
            nflg = 1;
            lgm = lngamma(aaa+bbb, &s, _state)-lngamma(aaa, &s, _state)-lngamma(bbb, &s, _state);
            i = 0;
            mainlooppos = newtcycle;
            continue;
        }
        
        /*
         * newtcycle
         */
        if( mainlooppos==newtcycle )
        {
            if( i<=7 )
            {
                if( i!=0 )
                {
                    yyy = incompletebeta(aaa, bbb, x, _state);
                }
                if( ae_fp_less(yyy,yl) )
                {
                    x = x0;
                    yyy = yl;
                }
                else
                {
                    if( ae_fp_greater(yyy,yh) )
                    {
                        x = x1;
                        yyy = yh;
                    }
                    else
                    {
                        if( ae_fp_less(yyy,y0) )
                        {
                            x0 = x;
                            yl = yyy;
                        }
                        else
                        {
                            x1 = x;
                            yh = yyy;
                        }
                    }
                }
                if( ae_fp_eq(x,1.0)||ae_fp_eq(x,0.0) )
                {
                    mainlooppos = breaknewtcycle;
                    continue;
                }
                d = (aaa-1.0)*ae_log(x, _state)+(bbb-1.0)*ae_log(1.0-x, _state)+lgm;
                if( ae_fp_less(d,ae_log(ae_minrealnumber, _state)) )
                {
                    break;
                }
                if( ae_fp_greater(d,ae_log(ae_maxrealnumber, _state)) )
                {
                    mainlooppos = breaknewtcycle;
                    continue;
                }
                d = ae_exp(d, _state);
                d = (yyy-y0)/d;
                xt = x-d;
                if( ae_fp_less_eq(xt,x0) )
                {
                    yyy = (x-x0)/(x1-x0);
                    xt = x0+0.5*yyy*(x-x0);
                    if( ae_fp_less_eq(xt,0.0) )
                    {
                        mainlooppos = breaknewtcycle;
                        continue;
                    }
                }
                if( ae_fp_greater_eq(xt,x1) )
                {
                    yyy = (x1-x)/(x1-x0);
                    xt = x1-0.5*yyy*(x1-x);
                    if( ae_fp_greater_eq(xt,1.0) )
                    {
                        mainlooppos = breaknewtcycle;
                        continue;
                    }
                }
                x = xt;
                if( ae_fp_less(ae_fabs(d/x, _state),128.0*ae_machineepsilon) )
                {
                    break;
                }
                i = i+1;
                mainlooppos = newtcycle;
                continue;
            }
            else
            {
                mainlooppos = breaknewtcycle;
                continue;
            }
        }
        
        /*
         * breaknewtcycle
         */
        if( mainlooppos==breaknewtcycle )
        {
            dithresh = 256.0*ae_machineepsilon;
            mainlooppos = ihalve;
            continue;
        }
    }
    
    /*
     * done
     */
    if( rflg!=0 )
    {
        if( ae_fp_less_eq(x,ae_machineepsilon) )
        {
            x = 1.0-ae_machineepsilon;
        }
        else
        {
            x = 1.0-x;
        }
    }
    result = x;
    return result;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0xc8, %rsp
movq %rdi, %rbx
movsd %xmm1, 0x80(%rsp)
movsd %xmm0, 0x88(%rsp)
xorps %xmm0, %xmm0
movsd %xmm0, 0x8(%rsp)
xorpd %xmm1, %xmm1
movsd %xmm2, 0x68(%rsp)
movapd %xmm2, %xmm0
callq 0x2635a
testb %al, %al
je 0xb34c
movsd 0x28e61(%rip), %xmm1    # 0x341a0
movsd 0x68(%rsp), %xmm0
callq 0x2631e
jmp 0xb34e
xorl %eax, %eax
movzbl %al, %edi
leaq 0x2b2aa(%rip), %rsi     # 0x36602
movq %rbx, %rdx
callq 0x2581c
xorpd %xmm1, %xmm1
movsd 0x68(%rsp), %xmm0
callq 0x262ba
testb %al, %al
jne 0xc369
movsd 0x28e21(%rip), %xmm1    # 0x341a0
movsd 0x68(%rsp), %xmm0
movsd %xmm1, 0x8(%rsp)
callq 0x262ba
testb %al, %al
jne 0xc369
movsd 0x28e00(%rip), %xmm2    # 0x341a0
movapd %xmm2, %xmm0
subsd 0x68(%rsp), %xmm0
movsd %xmm0, 0xa8(%rsp)
movsd 0x88(%rsp), %xmm1
movapd %xmm1, %xmm0
addsd 0x80(%rsp), %xmm0
divsd %xmm0, %xmm1
movsd %xmm1, 0xb0(%rsp)
movb $0x1, %cl
xorpd %xmm0, %xmm0
movsd %xmm0, 0x30(%rsp)
xorl %edx, %edx
movl $0x0, %eax
movq %rax, 0x78(%rsp)
movl $0x0, %r14d
movq $0x0, 0xa0(%rsp)
xorl %esi, %esi
xorl %edi, %edi
xorl %r8d, %r8d
xorl %ebp, %ebp
xorl %r9d, %r9d
movsd %xmm0, 0x48(%rsp)
movsd %xmm2, 0x58(%rsp)
movsd %xmm0, 0x70(%rsp)
xorpd %xmm6, %xmm6
xorpd %xmm1, %xmm1
movsd %xmm1, 0x90(%rsp)
movapd %xmm2, %xmm3
xorpd %xmm4, %xmm4
movsd %xmm1, 0x50(%rsp)
xorpd %xmm2, %xmm2
movsd %xmm2, 0x60(%rsp)
movq %r14, %r12
movq %rdx, %r13
movapd %xmm6, %xmm2
testb $0x1, %cl
je 0xb52f
movsd %xmm2, 0x10(%rsp)
movsd %xmm3, 0x20(%rsp)
movsd 0x88(%rsp), %xmm0
movsd 0x28d29(%rip), %xmm1    # 0x341a0
callq 0x2631e
testb %al, %al
jne 0xb49e
movsd 0x80(%rsp), %xmm0
movsd 0x28d0f(%rip), %xmm1    # 0x341a0
callq 0x2631e
testb %al, %al
je 0xb7be
movsd 0x88(%rsp), %xmm0
movsd 0x80(%rsp), %xmm1
movsd 0xb0(%rsp), %xmm2
movq %rbx, %rdi
callq 0xa5e4
movapd %xmm0, %xmm4
movb $0x1, %sil
movq $0x0, 0x78(%rsp)
movsd 0xb0(%rsp), %xmm0
movsd %xmm0, 0x30(%rsp)
movsd 0x10(%rsp), %xmm6
movsd 0x29fea(%rip), %xmm0    # 0x354d8
movq %r13, %rdx
movq %r12, %r14
xorl %ecx, %ecx
xorl %edi, %edi
xorl %r8d, %r8d
xorl %ebp, %ebp
xorl %r9d, %r9d
movsd 0x68(%rsp), %xmm1
movsd %xmm1, 0x50(%rsp)
movsd 0x80(%rsp), %xmm1
movsd 0x88(%rsp), %xmm2
movsd %xmm2, 0x60(%rsp)
movsd 0x20(%rsp), %xmm3
jmp 0xb447
movl %r9d, %eax
movl %ebp, %r10d
movl %r8d, %r11d
movl %edi, %r15d
movapd %xmm0, %xmm5
movb $0x1, %dil
xorl %edx, %edx
movsd 0x28f1c(%rip), %xmm6    # 0x34468
movl $0x0, %r14d
movl $0x0, %ecx
movl $0x0, %r8d
movl $0x0, %ebp
movl $0x0, %r9d
testb $0x1, %sil
movl $0x0, %esi
jne 0xb447
testb $0x1, %r15b
je 0xb81d
movb $0x1, %r8b
xorl %ecx, %ecx
movapd %xmm2, %xmm6
movapd %xmm5, %xmm0
movq %r13, %rdx
movq %r12, %r14
movl $0x0, %edi
cmpq $0x63, %r13
jg 0xb447
testq %r13, %r13
movsd %xmm1, 0x18(%rsp)
movsd %xmm3, 0x20(%rsp)
movsd %xmm5, 0x28(%rsp)
je 0xb759
movsd 0x58(%rsp), %xmm1
movsd 0x48(%rsp), %xmm0
subsd %xmm0, %xmm1
movsd %xmm2, 0x10(%rsp)
movsd %xmm1, 0x8(%rsp)
mulsd %xmm1, %xmm2
addsd %xmm0, %xmm2
movsd %xmm2, 0x30(%rsp)
movapd %xmm2, %xmm0
movsd 0x28bac(%rip), %xmm1    # 0x341a0
callq 0x262ba
movsd 0x29e8f(%rip), %xmm0    # 0x35490
testb %al, %al
jne 0xb60b
movsd 0x30(%rsp), %xmm0
movsd %xmm0, 0x30(%rsp)
xorpd %xmm1, %xmm1
callq 0x262ba
testb %al, %al
je 0xb65f
movsd 0x8(%rsp), %xmm2
movsd 0x28e3c(%rip), %xmm0    # 0x34468
mulsd %xmm0, %xmm2
addsd 0x48(%rsp), %xmm2
xorpd %xmm1, %xmm1
movsd %xmm2, 0x30(%rsp)
movapd %xmm2, %xmm0
callq 0x262ba
movsd 0x28e17(%rip), %xmm0    # 0x34468
movsd %xmm0, 0x10(%rsp)
testb %al, %al
jne 0xc306
movsd 0x60(%rsp), %xmm0
movsd 0x18(%rsp), %xmm1
movsd 0x30(%rsp), %xmm2
movq %rbx, %rdi
callq 0xa5e4
movsd %xmm0, (%rsp)
movsd 0x48(%rsp), %xmm1
addsd 0x58(%rsp), %xmm1
movsd 0x8(%rsp), %xmm0
divsd %xmm1, %xmm0
movq %rbx, %rdi
callq 0x26575
movsd 0x28(%rsp), %xmm1
callq 0x262fe
movsd (%rsp), %xmm4
movsd 0x18(%rsp), %xmm1
movb $0x1, %bpl
xorl %ecx, %ecx
movsd 0x10(%rsp), %xmm6
movsd 0x28(%rsp), %xmm0
movq %r13, %rdx
movq %r12, %r14
movl $0x0, %esi
movl $0x0, %edi
movl $0x0, %r8d
movl $0x0, %r9d
testb %al, %al
movsd 0x20(%rsp), %xmm3
jne 0xb447
movsd 0x50(%rsp), %xmm1
subsd %xmm1, %xmm4
divsd %xmm1, %xmm4
movapd %xmm4, %xmm0
movq %rbx, %rdi
callq 0x26575
movsd 0x28(%rsp), %xmm1
callq 0x262fe
movsd 0x28(%rsp), %xmm0
movsd 0x10(%rsp), %xmm2
movsd (%rsp), %xmm4
movsd 0x20(%rsp), %xmm3
movsd 0x18(%rsp), %xmm1
xorl %ecx, %ecx
movapd %xmm2, %xmm6
movq %r13, %rdx
movq %r12, %r14
movl $0x0, %esi
movl $0x0, %edi
movl $0x0, %r8d
movl $0x0, %r9d
testb %al, %al
jne 0xb447
movsd %xmm2, 0x10(%rsp)
movsd %xmm4, (%rsp)
movapd %xmm4, %xmm0
movsd 0x50(%rsp), %xmm1
callq 0x262fe
testb %al, %al
je 0xbc46
testq %r12, %r12
movsd 0x30(%rsp), %xmm2
js 0xbc99
cmpq $0x4, %r12
movsd 0x20(%rsp), %xmm0
movsd (%rsp), %xmm1
jl 0xbcdb
movsd 0x289f9(%rip), %xmm1    # 0x341a0
movapd %xmm1, %xmm0
subsd 0x10(%rsp), %xmm0
mulsd %xmm0, %xmm0
subsd %xmm0, %xmm1
jmp 0xbcf7
movsd 0x68(%rsp), %xmm0
movq %rbx, %rdi
callq 0x5a37
movaps %xmm0, 0x30(%rsp)
movsd 0x68(%rsp), %xmm0
movsd 0x28c89(%rip), %xmm1    # 0x34468
callq 0x2633e
movl $0x1, %ecx
testb %al, %al
je 0xb8c2
movsd 0xa8(%rsp), %xmm0
movsd %xmm0, 0x50(%rsp)
movsd 0x88(%rsp), %xmm3
movsd 0x80(%rsp), %xmm1
movapd 0x30(%rsp), %xmm0
jmp 0xb8f0
testb $0x1, %r11b
je 0xbab8
movsd %xmm5, 0x28(%rsp)
movsd %xmm2, 0x10(%rsp)
movsd %xmm4, (%rsp)
movsd %xmm3, 0x20(%rsp)
movsd %xmm1, 0x18(%rsp)
movsd 0x48(%rsp), %xmm0
movsd 0x2894e(%rip), %xmm1    # 0x341a0
callq 0x2635a
testb %al, %al
jne 0xc30e
xorps %xmm0, %xmm0
movsd %xmm0, 0x8(%rsp)
xorpd %xmm1, %xmm1
movsd 0x30(%rsp), %xmm0
callq 0x2631e
movb $0x1, %bpl
xorl %ecx, %ecx
movsd 0x10(%rsp), %xmm6
movsd 0x28(%rsp), %xmm0
movq %r13, %rdx
movq %r12, %r14
movl $0x0, %esi
movl $0x0, %edi
movl $0x0, %r8d
movl $0x0, %r9d
testb %al, %al
movsd 0x18(%rsp), %xmm1
movsd 0x20(%rsp), %xmm3
movsd (%rsp), %xmm4
je 0xb447
jmp 0xc32c
movapd 0x30(%rsp), %xmm0
xorpd 0x28910(%rip), %xmm0    # 0x341e0
xorl %ecx, %ecx
movsd 0x68(%rsp), %xmm1
movsd %xmm1, 0x50(%rsp)
movsd 0x80(%rsp), %xmm3
movsd 0x88(%rsp), %xmm1
movapd %xmm0, 0x30(%rsp)
movsd %xmm1, 0x60(%rsp)
movsd %xmm3, 0x18(%rsp)
movq %rcx, 0x78(%rsp)
movapd %xmm0, %xmm2
mulsd %xmm0, %xmm2
movsd 0x29b99(%rip), %xmm0    # 0x354b0
addsd %xmm0, %xmm2
divsd 0x29b95(%rip), %xmm2    # 0x354b8
movapd %xmm2, %xmm4
movsd %xmm2, 0x70(%rsp)
movapd %xmm1, %xmm0
addsd %xmm1, %xmm0
movsd 0x28893(%rip), %xmm2    # 0x341d0
addsd %xmm2, %xmm0
movsd 0x28857(%rip), %xmm1    # 0x341a0
movapd %xmm1, %xmm5
divsd %xmm0, %xmm5
movsd %xmm5, 0x28(%rsp)
movapd %xmm3, %xmm0
addsd %xmm3, %xmm0
addsd %xmm2, %xmm0
divsd %xmm0, %xmm1
movsd %xmm1, (%rsp)
movapd %xmm1, %xmm0
addsd %xmm5, %xmm0
movsd 0x28a64(%rip), %xmm1    # 0x343e0
divsd %xmm0, %xmm1
movsd %xmm1, 0x8(%rsp)
movapd %xmm4, %xmm0
addsd %xmm1, %xmm0
movq %rbx, %rdi
callq 0x2658d
movapd %xmm0, %xmm2
mulsd 0x30(%rsp), %xmm2
movsd 0x8(%rsp), %xmm1
divsd %xmm1, %xmm2
movsd (%rsp), %xmm3
subsd 0x28(%rsp), %xmm3
movsd 0x70(%rsp), %xmm0
addsd 0x29afd(%rip), %xmm0    # 0x354c0
mulsd 0x29ae5(%rip), %xmm1    # 0x354b0
movapd %xmm1, %xmm4
movsd 0x28a09(%rip), %xmm1    # 0x343e0
divsd %xmm4, %xmm1
addsd %xmm0, %xmm1
mulsd %xmm3, %xmm1
subsd %xmm1, %xmm2
addsd %xmm2, %xmm2
movsd %xmm2, 0x30(%rsp)
movsd 0x29a7f(%rip), %xmm0    # 0x35478
movq %rbx, %rdi
callq 0x266d6
movapd %xmm0, %xmm1
movsd 0x30(%rsp), %xmm0
callq 0x262fe
testb %al, %al
jne 0xc322
movsd 0x30(%rsp), %xmm0
movq %rbx, %rdi
callq 0x266e0
movsd 0x18(%rsp), %xmm1
mulsd %xmm1, %xmm0
movsd 0x60(%rsp), %xmm3
addsd %xmm3, %xmm0
movapd %xmm3, %xmm2
divsd %xmm0, %xmm2
movapd %xmm3, %xmm0
movsd %xmm2, 0x30(%rsp)
movq %rbx, %rdi
callq 0xa5e4
movsd %xmm0, (%rsp)
movsd 0x50(%rsp), %xmm1
subsd %xmm1, %xmm0
divsd %xmm1, %xmm0
movq %rbx, %rdi
callq 0x26575
movsd 0x29a51(%rip), %xmm1    # 0x354c8
callq 0x262fe
movsd (%rsp), %xmm4
movl %eax, %ebp
movl %eax, %esi
xorb $0x1, %sil
xorl %ecx, %ecx
movsd 0x10(%rsp), %xmm6
movsd 0x29a37(%rip), %xmm0    # 0x354d0
movq %r13, %rdx
movq %r12, %r14
xorl %edi, %edi
xorl %r8d, %r8d
xorl %r9d, %r9d
movsd 0x18(%rsp), %xmm1
movsd 0x20(%rsp), %xmm3
jmp 0xb447
testb $0x1, %r10b
je 0xbb95
cmpq $0x0, 0xa0(%rsp)
jne 0xc306
movapd %xmm1, %xmm0
addsd 0x60(%rsp), %xmm0
leaq 0xc0(%rsp), %r14
movq %r14, %rdi
movq %rbx, %rsi
movsd %xmm1, 0x18(%rsp)
movsd %xmm3, 0x20(%rsp)
movsd %xmm4, (%rsp)
movsd %xmm2, 0x10(%rsp)
movsd %xmm5, 0x28(%rsp)
callq 0x4e78
movsd %xmm0, 0x70(%rsp)
movsd 0x60(%rsp), %xmm0
movq %r14, %rdi
movq %rbx, %rsi
callq 0x4e78
movsd 0x70(%rsp), %xmm1
subsd %xmm0, %xmm1
movsd %xmm1, 0x70(%rsp)
movsd 0x18(%rsp), %xmm0
movq %r14, %rdi
movq %rbx, %rsi
callq 0x4e78
movsd (%rsp), %xmm4
movsd 0x20(%rsp), %xmm3
movsd 0x18(%rsp), %xmm1
movsd 0x70(%rsp), %xmm2
subsd %xmm0, %xmm2
movsd %xmm2, 0x70(%rsp)
movb $0x1, %r9b
movl $0x1, %eax
movq %rax, 0xa0(%rsp)
xorl %edx, %edx
movsd 0x10(%rsp), %xmm6
movsd 0x28(%rsp), %xmm0
movq %r12, %r14
xorl %ecx, %ecx
xorl %esi, %esi
xorl %edi, %edi
xorl %r8d, %r8d
xorl %ebp, %ebp
jmp 0xb447
movb $0x1, %sil
xorl %ecx, %ecx
movapd %xmm2, %xmm6
movsd 0x298f2(%rip), %xmm0    # 0x35498
movq %r13, %rdx
movq %r12, %r14
movl $0x0, %edi
testb $0x1, %al
je 0xb447
xorl %ecx, %ecx
movapd %xmm2, %xmm6
movapd %xmm5, %xmm0
movq %r13, %rdx
movq %r12, %r14
movl $0x0, %esi
cmpq $0x7, %r13
jg 0xb447
movsd %xmm5, 0x28(%rsp)
movsd %xmm2, 0x10(%rsp)
movsd %xmm3, 0x20(%rsp)
movsd %xmm1, 0x18(%rsp)
testq %r13, %r13
je 0xbc13
movsd 0x60(%rsp), %xmm0
movsd 0x18(%rsp), %xmm1
movsd 0x30(%rsp), %xmm2
movq %rbx, %rdi
callq 0xa5e4
movapd %xmm0, %xmm4
movsd %xmm4, (%rsp)
movapd %xmm4, %xmm0
movsd 0x90(%rsp), %xmm1
callq 0x262fe
testb %al, %al
je 0xbcff
movsd 0x48(%rsp), %xmm0
movsd 0x90(%rsp), %xmm1
jmp 0xbd23
cmpq $0x1, 0x78(%rsp)
movsd 0x30(%rsp), %xmm3
jne 0xbc73
movapd %xmm3, %xmm0
movsd 0x28cc0(%rip), %xmm1    # 0x34920
callq 0x262fe
movsd 0x30(%rsp), %xmm3
testb %al, %al
jne 0xc318
testq %r12, %r12
jle 0xbcaf
movsd 0x287e8(%rip), %xmm6    # 0x34468
xorl %r12d, %r12d
movsd 0x18(%rsp), %xmm1
movsd (%rsp), %xmm4
movsd 0x28(%rsp), %xmm5
jmp 0xc2d8
movsd 0x287c7(%rip), %xmm0    # 0x34468
movsd %xmm0, 0x10(%rsp)
xorl %r12d, %r12d
jmp 0xbd66
cmpq $-0x4, %r12
movsd 0x18(%rsp), %xmm1
movsd 0x90(%rsp), %xmm2
movsd (%rsp), %xmm4
movsd 0x28(%rsp), %xmm5
jg 0xbd2d
movsd 0x10(%rsp), %xmm0
mulsd %xmm0, %xmm0
jmp 0xbd45
cmpq $0x2, %r12
jl 0xbd4e
movsd 0x2877f(%rip), %xmm0    # 0x34468
movsd 0x10(%rsp), %xmm1
mulsd %xmm0, %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, 0x10(%rsp)
jmp 0xbd66
movsd (%rsp), %xmm0
movsd 0x20(%rsp), %xmm1
callq 0x2633e
testb %al, %al
je 0xbea7
movsd 0x58(%rsp), %xmm0
movsd 0x20(%rsp), %xmm1
movsd %xmm1, (%rsp)
jmp 0xbef8
cmpq $-0x2, %r12
jg 0xc2c2
movsd 0x10(%rsp), %xmm0
mulsd 0x28723(%rip), %xmm0    # 0x34468
movapd %xmm0, %xmm6
jmp 0xc2d8
movsd 0x50(%rsp), %xmm3
subsd %xmm1, %xmm3
subsd %xmm1, %xmm0
divsd %xmm0, %xmm3
movsd %xmm3, 0x10(%rsp)
movq %r12, %r14
incq %r14
movapd %xmm2, %xmm0
movsd 0x29730(%rip), %xmm1    # 0x354a8
movsd %xmm2, 0x30(%rsp)
callq 0x2633e
testb %al, %al
je 0xbdb8
cmpq $0x1, 0x78(%rsp)
movsd 0x68(%rsp), %xmm0
jne 0xbdef
movq $0x0, 0x78(%rsp)
movsd %xmm0, 0x50(%rsp)
movsd 0x80(%rsp), %xmm1
movsd 0x88(%rsp), %xmm4
jmp 0xbe1a
movsd 0x30(%rsp), %xmm0
movsd %xmm0, 0x48(%rsp)
movsd (%rsp), %xmm4
movsd %xmm4, 0x90(%rsp)
movsd 0x18(%rsp), %xmm1
movsd 0x20(%rsp), %xmm3
movsd 0x28(%rsp), %xmm5
movsd 0x10(%rsp), %xmm6
jmp 0xc2e8
movl $0x1, %eax
movq %rax, 0x78(%rsp)
movsd 0xa8(%rsp), %xmm0
movsd %xmm0, 0x50(%rsp)
movsd 0x88(%rsp), %xmm1
movsd 0x80(%rsp), %xmm4
movsd 0x30(%rsp), %xmm3
movsd %xmm1, 0x18(%rsp)
movsd 0x28372(%rip), %xmm2    # 0x341a0
subsd %xmm3, %xmm2
movsd %xmm2, 0x20(%rsp)
movsd %xmm4, 0x60(%rsp)
movapd %xmm4, %xmm0
movq %rbx, %rdi
callq 0xa5e4
movsd 0x18(%rsp), %xmm1
movapd %xmm0, %xmm4
movb $0x1, %sil
xorl %ecx, %ecx
movsd 0x20(%rsp), %xmm0
movsd %xmm0, 0x30(%rsp)
xorps %xmm0, %xmm0
movsd %xmm0, 0x48(%rsp)
movsd 0x2832a(%rip), %xmm3    # 0x341a0
movsd %xmm3, 0x58(%rsp)
movsd 0x28(%rsp), %xmm0
xorpd %xmm2, %xmm2
movsd %xmm2, 0x90(%rsp)
movq %r13, %rdx
xorl %edi, %edi
xorl %r8d, %r8d
xorl %ebp, %ebp
xorl %r9d, %r9d
movsd 0x10(%rsp), %xmm6
jmp 0xb447
movsd (%rsp), %xmm0
movsd 0x50(%rsp), %xmm1
callq 0x262fe
movsd 0x30(%rsp), %xmm0
movapd %xmm0, %xmm2
movsd (%rsp), %xmm3
movapd %xmm3, %xmm1
testb %al, %al
jne 0xbee9
movsd 0x48(%rsp), %xmm2
movsd %xmm0, 0x58(%rsp)
movsd 0x90(%rsp), %xmm1
movsd %xmm3, 0x20(%rsp)
movsd %xmm2, 0x48(%rsp)
movsd %xmm1, 0x90(%rsp)
movsd %xmm0, 0x30(%rsp)
movsd 0x2829a(%rip), %xmm1    # 0x341a0
callq 0x262ba
xorl %ecx, %ecx
movsd 0x10(%rsp), %xmm6
movsd 0x28(%rsp), %xmm0
movq %r13, %rdx
movq %r12, %r14
movl $0x0, %esi
movl $0x0, %edi
movl $0x0, %r8d
movl $0x0, %ebp
movl $0x0, %r9d
testb %al, %al
movsd 0x18(%rsp), %xmm1
movsd 0x20(%rsp), %xmm3
movsd (%rsp), %xmm4
jne 0xb447
xorpd %xmm1, %xmm1
movsd 0x30(%rsp), %xmm0
callq 0x262ba
movsd (%rsp), %xmm4
movsd 0x20(%rsp), %xmm3
movsd 0x18(%rsp), %xmm1
xorl %ecx, %ecx
movsd 0x10(%rsp), %xmm6
movsd 0x28(%rsp), %xmm0
movq %r13, %rdx
movq %r12, %r14
movl $0x0, %esi
movl $0x0, %edi
movl $0x0, %r8d
movl $0x0, %ebp
movl $0x0, %r9d
testb %al, %al
jne 0xb447
movsd 0x60(%rsp), %xmm2
movsd 0x28218(%rip), %xmm0    # 0x341d0
addsd %xmm0, %xmm2
movsd %xmm2, 0x98(%rsp)
movsd 0x30(%rsp), %xmm0
movq %rbx, %rdi
callq 0x266d6
movsd %xmm0, 0x8(%rsp)
movsd 0x18(%rsp), %xmm0
addsd 0x281e9(%rip), %xmm0    # 0x341d0
movsd %xmm0, 0xb8(%rsp)
movsd 0x281a8(%rip), %xmm0    # 0x341a0
subsd 0x30(%rsp), %xmm0
movq %rbx, %rdi
callq 0x266d6
mulsd 0xb8(%rsp), %xmm0
movsd 0x8(%rsp), %xmm1
mulsd 0x98(%rsp), %xmm1
addsd %xmm0, %xmm1
addsd 0x70(%rsp), %xmm1
movsd %xmm1, 0x8(%rsp)
movsd 0x29442(%rip), %xmm0    # 0x35478
movq %rbx, %rdi
callq 0x266d6
movapd %xmm0, %xmm1
movsd 0x8(%rsp), %xmm0
callq 0x262fe
testb %al, %al
jne 0xc306
movsd 0x285db(%rip), %xmm0    # 0x34638
movq %rbx, %rdi
callq 0x266d6
movaps %xmm0, %xmm1
movsd 0x8(%rsp), %xmm0
callq 0x2633e
xorl %ecx, %ecx
movsd 0x10(%rsp), %xmm6
movsd 0x28(%rsp), %xmm0
movq %r13, %rdx
movq %r12, %r14
movl $0x0, %esi
movl $0x0, %edi
movl $0x0, %r8d
movl $0x0, %ebp
movl $0x0, %r9d
testb %al, %al
movsd 0x18(%rsp), %xmm1
movsd 0x20(%rsp), %xmm3
movsd (%rsp), %xmm4
jne 0xb447
movsd 0x8(%rsp), %xmm0
movq %rbx, %rdi
callq 0x266e0
movsd (%rsp), %xmm1
subsd 0x50(%rsp), %xmm1
divsd %xmm0, %xmm1
movsd 0x30(%rsp), %xmm0
movsd %xmm1, 0x98(%rsp)
subsd %xmm1, %xmm0
movsd %xmm0, 0x8(%rsp)
movsd 0x48(%rsp), %xmm1
callq 0x2631e
movsd (%rsp), %xmm4
testb %al, %al
je 0xc197
movsd 0x30(%rsp), %xmm3
movsd 0x48(%rsp), %xmm2
subsd %xmm2, %xmm3
movsd 0x58(%rsp), %xmm1
subsd %xmm2, %xmm1
movapd %xmm3, %xmm0
divsd %xmm1, %xmm0
movsd %xmm0, (%rsp)
mulsd 0x28330(%rip), %xmm0    # 0x34468
mulsd %xmm3, %xmm0
addsd %xmm2, %xmm0
xorpd %xmm1, %xmm1
movsd %xmm0, 0x8(%rsp)
callq 0x2631e
movsd (%rsp), %xmm4
movsd 0x20(%rsp), %xmm3
movsd 0x18(%rsp), %xmm1
xorl %ecx, %ecx
movsd 0x10(%rsp), %xmm6
movsd 0x28(%rsp), %xmm0
movq %r13, %rdx
movq %r12, %r14
movl $0x0, %esi
movl $0x0, %edi
movl $0x0, %r8d
movl $0x0, %ebp
movl $0x0, %r9d
testb %al, %al
jne 0xb447
movsd %xmm4, (%rsp)
movsd 0x8(%rsp), %xmm0
movsd 0x58(%rsp), %xmm1
callq 0x2635a
movsd (%rsp), %xmm4
testb %al, %al
je 0xc24c
movsd 0x58(%rsp), %xmm2
movapd %xmm2, %xmm3
subsd 0x30(%rsp), %xmm3
movapd %xmm2, %xmm1
subsd 0x48(%rsp), %xmm1
movapd %xmm3, %xmm0
divsd %xmm1, %xmm0
movsd %xmm0, (%rsp)
mulsd 0x2829f(%rip), %xmm0    # 0x34488
mulsd %xmm3, %xmm0
addsd %xmm2, %xmm0
movsd %xmm0, 0x8(%rsp)
movsd 0x27fa1(%rip), %xmm1    # 0x341a0
callq 0x2635a
movsd (%rsp), %xmm4
movsd 0x20(%rsp), %xmm3
movsd 0x18(%rsp), %xmm1
xorl %ecx, %ecx
movsd 0x10(%rsp), %xmm6
movsd 0x28(%rsp), %xmm0
movq %r13, %rdx
movq %r12, %r14
movl $0x0, %esi
movl $0x0, %edi
movl $0x0, %r8d
movl $0x0, %ebp
movl $0x0, %r9d
testb %al, %al
jne 0xb447
movsd %xmm4, (%rsp)
movsd 0x98(%rsp), %xmm0
divsd 0x8(%rsp), %xmm0
movq %rbx, %rdi
callq 0x26575
movsd 0x29230(%rip), %xmm1    # 0x354a0
callq 0x262fe
testb %al, %al
jne 0xc32c
incq %r13
movb $0x1, %r9b
xorl %ecx, %ecx
movsd 0x8(%rsp), %xmm0
movsd %xmm0, 0x30(%rsp)
movsd 0x10(%rsp), %xmm6
movsd 0x28(%rsp), %xmm0
movq %r13, %rdx
movq %r12, %r14
xorl %esi, %esi
xorl %edi, %edi
xorl %r8d, %r8d
xorl %ebp, %ebp
movsd 0x18(%rsp), %xmm1
movsd 0x20(%rsp), %xmm3
movsd (%rsp), %xmm4
jmp 0xb447
movapd %xmm4, %xmm6
subsd 0x50(%rsp), %xmm6
movapd %xmm4, %xmm0
subsd %xmm2, %xmm0
divsd %xmm0, %xmm6
decq %r12
movsd %xmm3, 0x58(%rsp)
movapd %xmm4, %xmm3
movq %r12, %r14
incq %r13
movb $0x1, %dil
xorl %ecx, %ecx
movapd %xmm5, %xmm0
movq %r13, %rdx
xorl %esi, %esi
xorl %r8d, %r8d
xorl %ebp, %ebp
xorl %r9d, %r9d
jmp 0xb447
movsd 0x30(%rsp), %xmm0
jmp 0xc326
movsd 0x2917a(%rip), %xmm0    # 0x35490
jmp 0xc326
movl $0x1, %eax
movq %rax, 0x78(%rsp)
xorpd %xmm0, %xmm0
movsd %xmm0, 0x8(%rsp)
cmpq $0x0, 0x78(%rsp)
je 0xc369
movsd 0x285e4(%rip), %xmm1    # 0x34920
movsd 0x8(%rsp), %xmm0
callq 0x2631e
testb %al, %al
je 0xc355
movsd 0x2913d(%rip), %xmm0    # 0x35490
jmp 0xc363
movsd 0x27e43(%rip), %xmm0    # 0x341a0
subsd 0x8(%rsp), %xmm0
movsd %xmm0, 0x8(%rsp)
movsd 0x8(%rsp), %xmm0
addq $0xc8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq