// File: /workspace/github_cpp_makefile_star2_O2/mouda[P]LTE_random_access-simulate/alglibinternal.cpp
// Function: alglib_impl::hsschur_internalauxschur(bool, bool, long, long, long, alglib_impl::ae_matrix*, alglib_impl::ae_vector*, alglib_impl::ae_vector*, long, long, alglib_impl::ae_matrix*, alglib_impl::ae_vector*, alglib_impl::ae_vector*, alglib_impl::ae_vector*, alglib_impl::ae_vector*, long*, alglib_impl::ae_state*)
static void hsschur_internalauxschur(ae_bool wantt,
     ae_bool wantz,
     ae_int_t n,
     ae_int_t ilo,
     ae_int_t ihi,
     /* Real    */ ae_matrix* h,
     /* Real    */ ae_vector* wr,
     /* Real    */ ae_vector* wi,
     ae_int_t iloz,
     ae_int_t ihiz,
     /* Real    */ ae_matrix* z,
     /* Real    */ ae_vector* work,
     /* Real    */ ae_vector* workv3,
     /* Real    */ ae_vector* workc1,
     /* Real    */ ae_vector* works1,
     ae_int_t* info,
     ae_state *_state)
{
    ae_int_t i;
    ae_int_t i1;
    ae_int_t i2;
    ae_int_t itn;
    ae_int_t its;
    ae_int_t j;
    ae_int_t k;
    ae_int_t l;
    ae_int_t m;
    ae_int_t nh;
    ae_int_t nr;
    ae_int_t nz;
    double ave;
    double cs;
    double disc;
    double h00;
    double h10;
    double h11;
    double h12;
    double h21;
    double h22;
    double h33;
    double h33s;
    double h43h34;
    double h44;
    double h44s;
    double ovfl;
    double s;
    double smlnum;
    double sn;
    double sum;
    double t1;
    double t2;
    double t3;
    double tst1;
    double unfl;
    double v1;
    double v2;
    double v3;
    ae_bool failflag;
    double dat1;
    double dat2;
    ae_int_t p1;
    double him1im1;
    double him1i;
    double hiim1;
    double hii;
    double wrim1;
    double wri;
    double wiim1;
    double wii;
    double ulp;

    *info = 0;

    *info = 0;
    dat1 = 0.75;
    dat2 = -0.4375;
    ulp = ae_machineepsilon;
    
    /*
     * Quick return if possible
     */
    if( n==0 )
    {
        return;
    }
    if( ilo==ihi )
    {
        wr->ptr.p_double[ilo] = h->ptr.pp_double[ilo][ilo];
        wi->ptr.p_double[ilo] = 0;
        return;
    }
    nh = ihi-ilo+1;
    nz = ihiz-iloz+1;
    
    /*
     * Set machine-dependent constants for the stopping criterion.
     * If norm(H) <= sqrt(OVFL), overflow should not occur.
     */
    unfl = ae_minrealnumber;
    ovfl = 1/unfl;
    smlnum = unfl*(nh/ulp);
    
    /*
     * I1 and I2 are the indices of the first row and last column of H
     * to which transformations must be applied. If eigenvalues only are
     * being computed, I1 and I2 are set inside the main loop.
     */
    i1 = 1;
    i2 = n;
    
    /*
     * ITN is the total number of QR iterations allowed.
     */
    itn = 30*nh;
    
    /*
     * The main loop begins here. I is the loop index and decreases from
     * IHI to ILO in steps of 1 or 2. Each iteration of the loop works
     * with the active submatrix in rows and columns L to I.
     * Eigenvalues I+1 to IHI have already converged. Either L = ILO or
     * H(L,L-1) is negligible so that the matrix splits.
     */
    i = ihi;
    for(;;)
    {
        l = ilo;
        if( i<ilo )
        {
            return;
        }
        
        /*
         * Perform QR iterations on rows and columns ILO to I until a
         * submatrix of order 1 or 2 splits off at the bottom because a
         * subdiagonal element has become negligible.
         */
        failflag = ae_true;
        for(its=0; its<=itn; its++)
        {
            
            /*
             * Look for a single small subdiagonal element.
             */
            for(k=i; k>=l+1; k--)
            {
                tst1 = ae_fabs(h->ptr.pp_double[k-1][k-1], _state)+ae_fabs(h->ptr.pp_double[k][k], _state);
                if( ae_fp_eq(tst1,0) )
                {
                    tst1 = upperhessenberg1norm(h, l, i, l, i, work, _state);
                }
                if( ae_fp_less_eq(ae_fabs(h->ptr.pp_double[k][k-1], _state),ae_maxreal(ulp*tst1, smlnum, _state)) )
                {
                    break;
                }
            }
            l = k;
            if( l>ilo )
            {
                
                /*
                 * H(L,L-1) is negligible
                 */
                h->ptr.pp_double[l][l-1] = 0;
            }
            
            /*
             * Exit from loop if a submatrix of order 1 or 2 has split off.
             */
            if( l>=i-1 )
            {
                failflag = ae_false;
                break;
            }
            
            /*
             * Now the active submatrix is in rows and columns L to I. If
             * eigenvalues only are being computed, only the active submatrix
             * need be transformed.
             */
            if( its==10||its==20 )
            {
                
                /*
                 * Exceptional shift.
                 */
                s = ae_fabs(h->ptr.pp_double[i][i-1], _state)+ae_fabs(h->ptr.pp_double[i-1][i-2], _state);
                h44 = dat1*s+h->ptr.pp_double[i][i];
                h33 = h44;
                h43h34 = dat2*s*s;
            }
            else
            {
                
                /*
                 * Prepare to use Francis' double shift
                 * (i.e. 2nd degree generalized Rayleigh quotient)
                 */
                h44 = h->ptr.pp_double[i][i];
                h33 = h->ptr.pp_double[i-1][i-1];
                h43h34 = h->ptr.pp_double[i][i-1]*h->ptr.pp_double[i-1][i];
                s = h->ptr.pp_double[i-1][i-2]*h->ptr.pp_double[i-1][i-2];
                disc = (h33-h44)*0.5;
                disc = disc*disc+h43h34;
                if( ae_fp_greater(disc,0) )
                {
                    
                    /*
                     * Real roots: use Wilkinson's shift twice
                     */
                    disc = ae_sqrt(disc, _state);
                    ave = 0.5*(h33+h44);
                    if( ae_fp_greater(ae_fabs(h33, _state)-ae_fabs(h44, _state),0) )
                    {
                        h33 = h33*h44-h43h34;
                        h44 = h33/(hsschur_extschursign(disc, ave, _state)+ave);
                    }
                    else
                    {
                        h44 = hsschur_extschursign(disc, ave, _state)+ave;
                    }
                    h33 = h44;
                    h43h34 = 0;
                }
            }
            
            /*
             * Look for two consecutive small subdiagonal elements.
             */
            for(m=i-2; m>=l; m--)
            {
                
                /*
                 * Determine the effect of starting the double-shift QR
                 * iteration at row M, and see if this would make H(M,M-1)
                 * negligible.
                 */
                h11 = h->ptr.pp_double[m][m];
                h22 = h->ptr.pp_double[m+1][m+1];
                h21 = h->ptr.pp_double[m+1][m];
                h12 = h->ptr.pp_double[m][m+1];
                h44s = h44-h11;
                h33s = h33-h11;
                v1 = (h33s*h44s-h43h34)/h21+h12;
                v2 = h22-h11-h33s-h44s;
                v3 = h->ptr.pp_double[m+2][m+1];
                s = ae_fabs(v1, _state)+ae_fabs(v2, _state)+ae_fabs(v3, _state);
                v1 = v1/s;
                v2 = v2/s;
                v3 = v3/s;
                workv3->ptr.p_double[1] = v1;
                workv3->ptr.p_double[2] = v2;
                workv3->ptr.p_double[3] = v3;
                if( m==l )
                {
                    break;
                }
                h00 = h->ptr.pp_double[m-1][m-1];
                h10 = h->ptr.pp_double[m][m-1];
                tst1 = ae_fabs(v1, _state)*(ae_fabs(h00, _state)+ae_fabs(h11, _state)+ae_fabs(h22, _state));
                if( ae_fp_less_eq(ae_fabs(h10, _state)*(ae_fabs(v2, _state)+ae_fabs(v3, _state)),ulp*tst1) )
                {
                    break;
                }
            }
            
            /*
             * Double-shift QR step
             */
            for(k=m; k<=i-1; k++)
            {
                
                /*
                 * The first iteration of this loop determines a reflection G
                 * from the vector V and applies it from left and right to H,
                 * thus creating a nonzero bulge below the subdiagonal.
                 *
                 * Each subsequent iteration determines a reflection G to
                 * restore the Hessenberg form in the (K-1)th column, and thus
                 * chases the bulge one step toward the bottom of the active
                 * submatrix. NR is the order of G.
                 */
                nr = ae_minint(3, i-k+1, _state);
                if( k>m )
                {
                    for(p1=1; p1<=nr; p1++)
                    {
                        workv3->ptr.p_double[p1] = h->ptr.pp_double[k+p1-1][k-1];
                    }
                }
                generatereflection(workv3, nr, &t1, _state);
                if( k>m )
                {
                    h->ptr.pp_double[k][k-1] = workv3->ptr.p_double[1];
                    h->ptr.pp_double[k+1][k-1] = 0;
                    if( k<i-1 )
                    {
                        h->ptr.pp_double[k+2][k-1] = 0;
                    }
                }
                else
                {
                    if( m>l )
                    {
                        h->ptr.pp_double[k][k-1] = -h->ptr.pp_double[k][k-1];
                    }
                }
                v2 = workv3->ptr.p_double[2];
                t2 = t1*v2;
                if( nr==3 )
                {
                    v3 = workv3->ptr.p_double[3];
                    t3 = t1*v3;
                    
                    /*
                     * Apply G from the left to transform the rows of the matrix
                     * in columns K to I2.
                     */
                    for(j=k; j<=i2; j++)
                    {
                        sum = h->ptr.pp_double[k][j]+v2*h->ptr.pp_double[k+1][j]+v3*h->ptr.pp_double[k+2][j];
                        h->ptr.pp_double[k][j] = h->ptr.pp_double[k][j]-sum*t1;
                        h->ptr.pp_double[k+1][j] = h->ptr.pp_double[k+1][j]-sum*t2;
                        h->ptr.pp_double[k+2][j] = h->ptr.pp_double[k+2][j]-sum*t3;
                    }
                    
                    /*
                     * Apply G from the right to transform the columns of the
                     * matrix in rows I1 to min(K+3,I).
                     */
                    for(j=i1; j<=ae_minint(k+3, i, _state); j++)
                    {
                        sum = h->ptr.pp_double[j][k]+v2*h->ptr.pp_double[j][k+1]+v3*h->ptr.pp_double[j][k+2];
                        h->ptr.pp_double[j][k] = h->ptr.pp_double[j][k]-sum*t1;
                        h->ptr.pp_double[j][k+1] = h->ptr.pp_double[j][k+1]-sum*t2;
                        h->ptr.pp_double[j][k+2] = h->ptr.pp_double[j][k+2]-sum*t3;
                    }
                    if( wantz )
                    {
                        
                        /*
                         * Accumulate transformations in the matrix Z
                         */
                        for(j=iloz; j<=ihiz; j++)
                        {
                            sum = z->ptr.pp_double[j][k]+v2*z->ptr.pp_double[j][k+1]+v3*z->ptr.pp_double[j][k+2];
                            z->ptr.pp_double[j][k] = z->ptr.pp_double[j][k]-sum*t1;
                            z->ptr.pp_double[j][k+1] = z->ptr.pp_double[j][k+1]-sum*t2;
                            z->ptr.pp_double[j][k+2] = z->ptr.pp_double[j][k+2]-sum*t3;
                        }
                    }
                }
                else
                {
                    if( nr==2 )
                    {
                        
                        /*
                         * Apply G from the left to transform the rows of the matrix
                         * in columns K to I2.
                         */
                        for(j=k; j<=i2; j++)
                        {
                            sum = h->ptr.pp_double[k][j]+v2*h->ptr.pp_double[k+1][j];
                            h->ptr.pp_double[k][j] = h->ptr.pp_double[k][j]-sum*t1;
                            h->ptr.pp_double[k+1][j] = h->ptr.pp_double[k+1][j]-sum*t2;
                        }
                        
                        /*
                         * Apply G from the right to transform the columns of the
                         * matrix in rows I1 to min(K+3,I).
                         */
                        for(j=i1; j<=i; j++)
                        {
                            sum = h->ptr.pp_double[j][k]+v2*h->ptr.pp_double[j][k+1];
                            h->ptr.pp_double[j][k] = h->ptr.pp_double[j][k]-sum*t1;
                            h->ptr.pp_double[j][k+1] = h->ptr.pp_double[j][k+1]-sum*t2;
                        }
                        if( wantz )
                        {
                            
                            /*
                             * Accumulate transformations in the matrix Z
                             */
                            for(j=iloz; j<=ihiz; j++)
                            {
                                sum = z->ptr.pp_double[j][k]+v2*z->ptr.pp_double[j][k+1];
                                z->ptr.pp_double[j][k] = z->ptr.pp_double[j][k]-sum*t1;
                                z->ptr.pp_double[j][k+1] = z->ptr.pp_double[j][k+1]-sum*t2;
                            }
                        }
                    }
                }
            }
        }
        if( failflag )
        {
            
            /*
             * Failure to converge in remaining number of iterations
             */
            *info = i;
            return;
        }
        if( l==i )
        {
            
            /*
             * H(I,I-1) is negligible: one eigenvalue has converged.
             */
            wr->ptr.p_double[i] = h->ptr.pp_double[i][i];
            wi->ptr.p_double[i] = 0;
        }
        else
        {
            if( l==i-1 )
            {
                
                /*
                 * H(I-1,I-2) is negligible: a pair of eigenvalues have converged.
                 *
                 *        Transform the 2-by-2 submatrix to standard Schur form,
                 *        and compute and store the eigenvalues.
                 */
                him1im1 = h->ptr.pp_double[i-1][i-1];
                him1i = h->ptr.pp_double[i-1][i];
                hiim1 = h->ptr.pp_double[i][i-1];
                hii = h->ptr.pp_double[i][i];
                hsschur_aux2x2schur(&him1im1, &him1i, &hiim1, &hii, &wrim1, &wiim1, &wri, &wii, &cs, &sn, _state);
                wr->ptr.p_double[i-1] = wrim1;
                wi->ptr.p_double[i-1] = wiim1;
                wr->ptr.p_double[i] = wri;
                wi->ptr.p_double[i] = wii;
                h->ptr.pp_double[i-1][i-1] = him1im1;
                h->ptr.pp_double[i-1][i] = him1i;
                h->ptr.pp_double[i][i-1] = hiim1;
                h->ptr.pp_double[i][i] = hii;
                if( wantt )
                {
                    
                    /*
                     * Apply the transformation to the rest of H.
                     */
                    if( i2>i )
                    {
                        workc1->ptr.p_double[1] = cs;
                        works1->ptr.p_double[1] = sn;
                        applyrotationsfromtheleft(ae_true, i-1, i, i+1, i2, workc1, works1, h, work, _state);
                    }
                    workc1->ptr.p_double[1] = cs;
                    works1->ptr.p_double[1] = sn;
                    applyrotationsfromtheright(ae_true, i1, i-2, i-1, i, workc1, works1, h, work, _state);
                }
                if( wantz )
                {
                    
                    /*
                     * Apply the transformation to Z.
                     */
                    workc1->ptr.p_double[1] = cs;
                    works1->ptr.p_double[1] = sn;
                    applyrotationsfromtheright(ae_true, iloz, iloz+nz-1, i-1, i, workc1, works1, z, work, _state);
                }
            }
        }
        
        /*
         * Decrement number of remaining iterations, and return to start of
         * the main loop with new value of I.
         */
        itn = itn-its;
        i = l-1;
    }
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x168, %rsp            # imm = 0x168
movq %r9, %r14
movq %rdx, %r12
movl %esi, 0x94(%rsp)
movl %edi, 0xdc(%rsp)
movq 0x1e0(%rsp), %rax
andq $0x0, (%rax)
movq 0x1a8(%rsp), %rdx
movq 0x1a0(%rsp), %rsi
movq %r8, %rax
subq %rcx, %rax
jne 0x1a0f2
movq 0x38(%r14), %rax
movq (%rax,%rcx,8), %rax
movsd (%rax,%rcx,8), %xmm0
movq 0x28(%rsi), %rax
movsd %xmm0, (%rax,%rcx,8)
movq 0x28(%rdx), %rax
andq $0x0, (%rax,%rcx,8)
jmp 0x1b681
movq %r8, %r13
movq 0x1b8(%rsp), %rbx
incq %rax
cvtsi2sd %rax, %xmm0
divsd 0x166e3(%rip), %xmm0    # 0x307f0
movq 0x1b0(%rsp), %rdx
mulsd 0x17153(%rip), %xmm0    # 0x31270
movsd %xmm0, 0x148(%rsp)
imulq $0x1e, %rax, %rax
xorl %esi, %esi
testq %rdx, %rdx
cmovgq %rdx, %rsi
xorpd %xmm0, %xmm0
movhpd 0x160b9(%rip), %xmm0    # xmm0 = xmm0[0],mem[0]
movapd %xmm0, 0x150(%rsp)
movq %rcx, 0xc0(%rsp)
movq %rsi, 0x110(%rsp)
movq %r12, 0xd0(%rsp)
movq %rax, 0x138(%rsp)
cmpq %rcx, %r13
jl 0x1b681
leaq -0x1(%r13), %rax
movq %rax, 0xb0(%rsp)
leaq -0x2(%r13), %rax
movq %rax, 0x100(%rsp)
testq %r13, %r13
movl $0x0, %edi
cmovgq %r13, %rdi
leaq 0x1(%r13), %rax
movq %rax, 0x108(%rsp)
leaq (,%r13,8), %rbp
movq %rcx, 0x20(%rsp)
xorl %eax, %eax
movq %rbp, 0x70(%rsp)
movq %r13, 0xb8(%rsp)
movq %rdi, 0xa0(%rsp)
movq 0x138(%rsp), %rdx
subq %rax, %rdx
jl 0x1b676
movq %rdx, 0x48(%rsp)
movq %rax, 0xf8(%rsp)
xorl %ecx, %ecx
movq %rcx, %r15
leaq (%rcx,%r13), %rax
movq %rax, 0x60(%rsp)
cmpq 0x20(%rsp), %rax
jle 0x1a2fb
movq 0x38(%r14), %rax
addq %rbp, %rax
movq -0x8(%rax,%r15,8), %rax
addq %rbp, %rax
movsd -0x8(%rax,%r15,8), %xmm0
movq %r15, %r13
movq 0x1e8(%rsp), %r15
movq %r15, %rdi
callq 0x2440c
movsd %xmm0, (%rsp)
movq 0x38(%r14), %rax
addq %rbp, %rax
movq (%rax,%r13,8), %rax
addq %rbp, %rax
movq %r13, %rbp
movsd (%rax,%r13,8), %xmm0
movq 0xb8(%rsp), %r12
movq %r15, %rdi
callq 0x2440c
addsd (%rsp), %xmm0
xorps %xmm1, %xmm1
movsd %xmm0, (%rsp)
callq 0x241f5
testb %al, %al
je 0x1a292
subq $0x8, %rsp
movq %r14, %rdi
movq 0x28(%rsp), %rcx
movq %rcx, %rsi
movq %r12, %rdx
movq %r12, %r8
movq 0x1c8(%rsp), %r9
pushq %r15
callq 0x15682
popq %rax
popq %rcx
movsd %xmm0, (%rsp)
movq 0x38(%r14), %rax
movq 0x70(%rsp), %r13
addq %r13, %rax
movq (%rax,%rbp,8), %rax
addq %r13, %rax
movsd -0x8(%rax,%rbp,8), %xmm0
movq %r15, %rdi
callq 0x2440c
movsd %xmm0, 0x10(%rsp)
movsd (%rsp), %xmm0
mulsd 0x1652a(%rip), %xmm0    # 0x307f0
movsd 0x148(%rsp), %xmm1
movq %r15, %rdi
movq %rbp, %r15
movq %r13, %rbp
movq %r12, %r13
callq 0x244b6
movapd %xmm0, %xmm1
movsd 0x10(%rsp), %xmm0
callq 0x24259
leaq -0x1(%r15), %rcx
testb %al, %al
je 0x1a1e5
movq 0xc0(%rsp), %rcx
cmpq %rcx, 0x60(%rsp)
jle 0x1a31e
movq 0x38(%r14), %rax
addq %rbp, %rax
movq (%rax,%r15,8), %rax
addq %rbp, %rax
andq $0x0, -0x8(%rax,%r15,8)
movq 0x60(%rsp), %rax
cmpq 0xb0(%rsp), %rax
movq 0x1a0(%rsp), %rbp
movsd 0x1614d(%rip), %xmm2    # 0x30488
movq 0xd0(%rsp), %r12
jge 0x1ad2f
movq 0xf8(%rsp), %rax
cmpq $0x14, %rax
je 0x1a361
cmpq $0xa, %rax
jne 0x1a654
movq 0x38(%r14), %rax
movq (%rax,%r13,8), %rax
movsd -0x8(%rax,%r13,8), %xmm0
movq 0x1e8(%rsp), %rbp
movq %rbp, %rdi
callq 0x2440c
movsd %xmm0, (%rsp)
movq 0x38(%r14), %rax
movq -0x8(%rax,%r13,8), %rax
movsd -0x10(%rax,%r13,8), %xmm0
movq %rbp, %rdi
callq 0x2440c
addsd (%rsp), %xmm0
movq 0x38(%r14), %rax
movq (%rax,%r13,8), %rax
movapd %xmm0, %xmm2
mulsd 0x186c2(%rip), %xmm2    # 0x32a78
mulsd %xmm0, %xmm2
mulsd 0x16ee6(%rip), %xmm0    # 0x312a8
addsd (%rax,%r13,8), %xmm0
movapd %xmm0, %xmm1
xorpd 0x15ddc(%rip), %xmm2    # 0x301b0
addq $0x2, %r15
movq %r15, 0x140(%rsp)
xorl %r15d, %r15d
movapd %xmm2, 0x120(%rsp)
movsd %xmm1, 0x88(%rsp)
movsd %xmm0, 0x98(%rsp)
movq 0x100(%rsp), %rax
addq %r15, %rax
movq %rax, 0x50(%rsp)
cmpq 0x60(%rsp), %rax
jl 0x1a64a
movq 0x38(%r14), %rax
movq 0x70(%rsp), %r13
addq %r13, %rax
movq -0x10(%rax,%r15,8), %rcx
addq %r13, %rcx
movsd -0x10(%rcx,%r15,8), %xmm4
movq -0x8(%rax,%r15,8), %rdx
addq %r13, %rdx
movsd -0x8(%rdx,%r15,8), %xmm3
subsd %xmm4, %xmm0
subsd %xmm4, %xmm1
movsd %xmm3, 0x30(%rsp)
movsd %xmm4, 0xe0(%rsp)
subsd %xmm4, %xmm3
subsd %xmm1, %xmm3
mulsd %xmm0, %xmm1
addsd %xmm2, %xmm1
divsd -0x10(%rdx,%r15,8), %xmm1
addsd -0x8(%rcx,%r15,8), %xmm1
movsd %xmm1, (%rsp)
subsd %xmm0, %xmm3
movsd %xmm3, 0x10(%rsp)
movq (%rax,%r15,8), %rax
addq %r13, %rax
movsd -0x8(%rax,%r15,8), %xmm0
movsd %xmm0, 0x20(%rsp)
movapd %xmm1, %xmm0
movq %rbp, %rdi
callq 0x2440c
movsd %xmm0, 0x48(%rsp)
movsd 0x10(%rsp), %xmm0
movq %rbp, %rdi
callq 0x2440c
addsd 0x48(%rsp), %xmm0
movsd %xmm0, 0x48(%rsp)
movsd 0x20(%rsp), %xmm0
movq %rbp, %rdi
callq 0x2440c
addsd 0x48(%rsp), %xmm0
movsd (%rsp), %xmm3
divsd %xmm0, %xmm3
movsd 0x10(%rsp), %xmm1
divsd %xmm0, %xmm1
movsd 0x20(%rsp), %xmm2
divsd %xmm0, %xmm2
movq 0x1c8(%rsp), %rcx
movq 0x28(%rcx), %rax
movsd %xmm3, 0x8(%rax)
movq 0x28(%rcx), %rax
movsd %xmm1, 0x10(%rsp)
movsd %xmm1, 0x10(%rax)
movq 0x28(%rcx), %rax
movsd %xmm2, 0x20(%rsp)
movsd %xmm2, 0x18(%rax)
cmpq %r15, 0x140(%rsp)
je 0x1a792
movapd %xmm3, %xmm0
movq 0x38(%r14), %rax
addq %r13, %rax
movq -0x18(%rax,%r15,8), %rcx
addq %r13, %rcx
movsd -0x18(%rcx,%r15,8), %xmm1
movsd %xmm1, (%rsp)
movq -0x10(%rax,%r15,8), %rax
addq %r13, %rax
movsd -0x18(%rax,%r15,8), %xmm1
movsd %xmm1, 0x48(%rsp)
movq %rbp, %rdi
callq 0x2440c
movsd %xmm0, 0xc8(%rsp)
movsd (%rsp), %xmm0
movq %rbp, %rdi
callq 0x2440c
movsd %xmm0, (%rsp)
movsd 0xe0(%rsp), %xmm0
movq %rbp, %rdi
callq 0x2440c
movsd (%rsp), %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, (%rsp)
movsd 0x30(%rsp), %xmm0
movq %rbp, %rdi
callq 0x2440c
movsd (%rsp), %xmm1
addsd %xmm0, %xmm1
mulsd 0xc8(%rsp), %xmm1
movsd %xmm1, (%rsp)
movsd 0x48(%rsp), %xmm0
movq %rbp, %rdi
callq 0x2440c
movsd %xmm0, 0x48(%rsp)
movsd 0x10(%rsp), %xmm0
movq %rbp, %rdi
callq 0x2440c
movsd %xmm0, 0x10(%rsp)
movsd 0x20(%rsp), %xmm0
movq %rbp, %rdi
callq 0x2440c
addsd 0x10(%rsp), %xmm0
mulsd 0x48(%rsp), %xmm0
movsd (%rsp), %xmm1
mulsd 0x161d1(%rip), %xmm1    # 0x307f0
callq 0x24259
decq %r15
testb %al, %al
movapd 0x120(%rsp), %xmm2
movsd 0x88(%rsp), %xmm1
movsd 0x98(%rsp), %xmm0
je 0x1a3fe
movq 0x50(%rsp), %rax
jmp 0x1a797
movq 0x38(%r14), %rax
movq -0x8(%rax,%r13,8), %rcx
movq (%rax,%r13,8), %rax
movsd -0x8(%rax,%r13,8), %xmm4
movsd (%rax,%r13,8), %xmm3
movsd -0x8(%rcx,%r13,8), %xmm1
mulsd (%rcx,%r13,8), %xmm4
movsd %xmm1, 0x88(%rsp)
movapd %xmm1, %xmm0
movsd %xmm3, 0x98(%rsp)
subsd %xmm3, %xmm0
mulsd %xmm2, %xmm0
mulsd %xmm0, %xmm0
movapd %xmm4, 0x120(%rsp)
addsd %xmm4, %xmm0
xorpd %xmm1, %xmm1
movsd %xmm0, (%rsp)
callq 0x24279
testb %al, %al
movq 0x1e8(%rsp), %rbp
je 0x1acee
movsd (%rsp), %xmm0
movq %rbp, %rdi
callq 0x24424
movsd %xmm0, 0x10(%rsp)
movsd 0x98(%rsp), %xmm1
movsd 0x88(%rsp), %xmm0
addsd %xmm0, %xmm1
mulsd 0x15d8f(%rip), %xmm1    # 0x30488
movsd %xmm1, (%rsp)
movq %rbp, %rdi
callq 0x2440c
movsd %xmm0, 0x20(%rsp)
movsd 0x98(%rsp), %xmm0
movq %rbp, %rdi
callq 0x2440c
movsd 0x20(%rsp), %xmm2
subsd %xmm0, %xmm2
xorpd %xmm1, %xmm1
movapd %xmm2, %xmm0
callq 0x24279
testb %al, %al
je 0x1ad0e
movsd 0x88(%rsp), %xmm0
mulsd 0x98(%rsp), %xmm0
subsd 0x120(%rsp), %xmm0
movsd %xmm0, 0x88(%rsp)
movsd 0x10(%rsp), %xmm0
movsd (%rsp), %xmm1
movq %rbp, %rdi
callq 0x230b0
addsd (%rsp), %xmm0
movsd 0x88(%rsp), %xmm1
divsd %xmm0, %xmm1
xorpd %xmm2, %xmm2
movapd %xmm1, %xmm0
jmp 0x1a3cc
movq 0x60(%rsp), %rax
movq %rax, 0x50(%rsp)
movq %rax, %r8
movq 0xb8(%rsp), %r13
cmpq %r13, %r8
jge 0x1accf
movq %r8, %r15
movq %r13, %r12
movq 0x108(%rsp), %rsi
subq %r8, %rsi
pushq $0x3
popq %rdi
movq %rbp, %rdx
callq 0x244ab
movq %rbp, %rcx
movq %rax, %rbp
movq 0x50(%rsp), %r13
cmpq %r13, %r15
jle 0x1a82b
testq %rbp, %rbp
movl $0x0, %eax
cmovgq %rbp, %rax
leaq (,%r15,8), %rcx
xorl %edx, %edx
movq 0x1c8(%rsp), %rdi
movq %r12, %r13
movq 0xd0(%rsp), %r12
cmpq %rdx, %rax
je 0x1a87c
movq 0x38(%r14), %rsi
addq %rcx, %rsi
movq (%rsi,%rdx,8), %rsi
movsd -0x8(%rsi,%r15,8), %xmm0
movq 0x28(%rdi), %rsi
movsd %xmm0, 0x8(%rsi,%rdx,8)
incq %rdx
jmp 0x1a805
movq 0x1c8(%rsp), %rbx
movq %rbx, %rdi
movq %rbp, %rsi
leaq 0x118(%rsp), %rdx
callq 0x16e77
cmpq 0x60(%rsp), %r13
movq %r12, %r13
movq 0xd0(%rsp), %r12
jle 0x1a8e4
movq 0x38(%r14), %rax
movq (%rax,%r15,8), %rax
movsd -0x8(%rax,%r15,8), %xmm0
xorpd 0x1593d(%rip), %xmm0    # 0x301b0
movlpd %xmm0, -0x8(%rax,%r15,8)
jmp 0x1a8e4
movq 0x1c8(%rsp), %rdi
movq %rbp, %rsi
leaq 0x118(%rsp), %rdx
movq 0x1e8(%rsp), %rcx
callq 0x16e77
movq 0x1c8(%rsp), %rbx
movq 0x28(%rbx), %rax
movsd 0x8(%rax), %xmm0
movq 0x38(%r14), %rax
movq (%rax,%r15,8), %rax
movsd %xmm0, -0x8(%rax,%r15,8)
movq 0x38(%r14), %rax
movq 0x8(%rax,%r15,8), %rax
andq $0x0, -0x8(%rax,%r15,8)
cmpq 0xb0(%rsp), %r15
jge 0x1a8e4
movq 0x38(%r14), %rax
movq 0x10(%rax,%r15,8), %rax
andq $0x0, -0x8(%rax,%r15,8)
movq 0x28(%rbx), %rax
movsd 0x10(%rax), %xmm3
movsd 0x118(%rsp), %xmm4
movapd %xmm3, %xmm5
mulsd %xmm4, %xmm5
cmpq $0x2, %rbp
movq 0x110(%rsp), %rsi
movq 0xa0(%rsp), %rdi
je 0x1ab88
cmpq $0x3, %rbp
movq 0x1b8(%rsp), %rbx
jne 0x1acbc
movsd 0x18(%rax), %xmm6
movapd %xmm4, %xmm7
mulsd %xmm6, %xmm7
leaq 0x1(%r15), %rdi
movq %r15, %rax
cmpq %r12, %rax
jg 0x1a9cb
movq 0x38(%r14), %rcx
movq (%rcx,%r15,8), %rdx
movq 0x8(%rcx,%r15,8), %rsi
movsd (%rdx,%rax,8), %xmm0
movsd (%rsi,%rax,8), %xmm1
mulsd %xmm3, %xmm1
addsd %xmm0, %xmm1
movq 0x10(%rcx,%r15,8), %rcx
movsd (%rcx,%rax,8), %xmm2
mulsd %xmm6, %xmm2
addsd %xmm1, %xmm2
movapd %xmm2, %xmm1
mulsd %xmm4, %xmm1
subsd %xmm1, %xmm0
movsd %xmm0, (%rdx,%rax,8)
movq 0x38(%r14), %rcx
movq 0x8(%rcx,%r15,8), %rcx
movsd (%rcx,%rax,8), %xmm0
movapd %xmm2, %xmm1
mulsd %xmm5, %xmm1
subsd %xmm1, %xmm0
movsd %xmm0, (%rcx,%rax,8)
movq 0x38(%r14), %rcx
movq 0x10(%rcx,%r15,8), %rcx
movsd (%rcx,%rax,8), %xmm0
mulsd %xmm7, %xmm2
subsd %xmm2, %xmm0
movsd %xmm0, (%rcx,%rax,8)
incq %rax
jmp 0x1a93e
movsd %xmm7, 0x30(%rsp)
movsd %xmm6, 0x48(%rsp)
movq %rdi, 0xe0(%rsp)
movsd %xmm5, 0x20(%rsp)
movsd %xmm4, 0x10(%rsp)
movsd %xmm3, (%rsp)
movq %r13, %rbx
leaq 0x3(%r15), %rbp
pushq $0x1
popq %r13
movq %rbp, %rdi
movq %rbx, %rsi
movq 0x1e8(%rsp), %rdx
callq 0x244ab
cmpq %rax, %r13
jg 0x1aaa6
movq 0x38(%r14), %rax
movq (%rax,%r13,8), %rax
movsd (%rax,%r15,8), %xmm0
movsd 0x8(%rax,%r15,8), %xmm1
mulsd (%rsp), %xmm1
addsd %xmm0, %xmm1
movsd 0x10(%rax,%r15,8), %xmm2
mulsd 0x48(%rsp), %xmm2
addsd %xmm1, %xmm2
movapd %xmm2, %xmm1
mulsd 0x10(%rsp), %xmm1
subsd %xmm1, %xmm0
movsd %xmm0, (%rax,%r15,8)
movq 0x38(%r14), %rax
movq (%rax,%r13,8), %rax
movsd 0x8(%rax,%r15,8), %xmm0
movapd %xmm2, %xmm1
mulsd 0x20(%rsp), %xmm1
subsd %xmm1, %xmm0
movsd %xmm0, 0x8(%rax,%r15,8)
movq 0x38(%r14), %rax
movq (%rax,%r13,8), %rax
movsd 0x10(%rax,%r15,8), %xmm0
mulsd 0x30(%rsp), %xmm2
subsd %xmm2, %xmm0
movsd %xmm0, 0x10(%rax,%r15,8)
incq %r13
jmp 0x1a9fb
cmpb $0x0, 0x94(%rsp)
movq %rbx, %r13
movq 0x1e8(%rsp), %rbp
movq 0x1b8(%rsp), %rbx
movq 0x110(%rsp), %rsi
movsd (%rsp), %xmm3
movsd 0x10(%rsp), %xmm4
movsd 0x20(%rsp), %xmm5
movq 0xe0(%rsp), %r8
movsd 0x48(%rsp), %xmm6
movsd 0x30(%rsp), %xmm7
je 0x1a7a7
xorl %eax, %eax
cmpq %rax, %rsi
je 0x1a7a7
movq 0x38(%rbx), %rcx
movq 0x8(%rcx,%rax,8), %rcx
movsd (%rcx,%r15,8), %xmm0
movsd 0x8(%rcx,%r15,8), %xmm1
mulsd %xmm3, %xmm1
addsd %xmm0, %xmm1
movsd 0x10(%rcx,%r15,8), %xmm2
mulsd %xmm6, %xmm2
addsd %xmm1, %xmm2
movapd %xmm2, %xmm1
mulsd %xmm4, %xmm1
subsd %xmm1, %xmm0
movsd %xmm0, (%rcx,%r15,8)
movq 0x38(%rbx), %rcx
movq 0x8(%rcx,%rax,8), %rcx
movsd 0x8(%rcx,%r15,8), %xmm0
movapd %xmm2, %xmm1
mulsd %xmm5, %xmm1
subsd %xmm1, %xmm0
movsd %xmm0, 0x8(%rcx,%r15,8)
movq 0x38(%rbx), %rcx
movq 0x8(%rcx,%rax,8), %rcx
movsd 0x10(%rcx,%r15,8), %xmm0
mulsd %xmm7, %xmm2
subsd %xmm2, %xmm0
movsd %xmm0, 0x10(%rcx,%r15,8)
incq %rax
jmp 0x1aaf6
leaq 0x1(%r15), %r8
movq %r15, %rax
movq 0x1e8(%rsp), %rbp
movq 0x1b8(%rsp), %rbx
cmpq %r12, %rax
jg 0x1abf4
movq 0x38(%r14), %rcx
movq (%rcx,%r15,8), %rdx
movq 0x8(%rcx,%r15,8), %rcx
movsd (%rdx,%rax,8), %xmm0
movsd (%rcx,%rax,8), %xmm1
mulsd %xmm3, %xmm1
addsd %xmm0, %xmm1
movapd %xmm1, %xmm2
mulsd %xmm4, %xmm2
subsd %xmm2, %xmm0
movsd %xmm0, (%rdx,%rax,8)
movq 0x38(%r14), %rcx
movq 0x8(%rcx,%r15,8), %rcx
movsd (%rcx,%rax,8), %xmm0
mulsd %xmm5, %xmm1
subsd %xmm1, %xmm0
movsd %xmm0, (%rcx,%rax,8)
incq %rax
jmp 0x1ab9f
xorl %eax, %eax
cmpq %rax, %rdi
je 0x1ac4f
movq 0x38(%r14), %rcx
movq 0x8(%rcx,%rax,8), %rcx
movsd (%rcx,%r15,8), %xmm0
movsd 0x8(%rcx,%r15,8), %xmm1
mulsd %xmm3, %xmm1
addsd %xmm0, %xmm1
movapd %xmm1, %xmm2
mulsd %xmm4, %xmm2
subsd %xmm2, %xmm0
movsd %xmm0, (%rcx,%r15,8)
movq 0x38(%r14), %rcx
movq 0x8(%rcx,%rax,8), %rcx
movsd 0x8(%rcx,%r15,8), %xmm0
mulsd %xmm5, %xmm1
subsd %xmm1, %xmm0
movsd %xmm0, 0x8(%rcx,%r15,8)
incq %rax
jmp 0x1abf6
cmpb $0x0, 0x94(%rsp)
je 0x1a7a7
xorl %eax, %eax
cmpq %rax, %rsi
je 0x1a7a7
movq 0x38(%rbx), %rcx
movq 0x8(%rcx,%rax,8), %rcx
movsd (%rcx,%r15,8), %xmm0
movsd 0x8(%rcx,%r15,8), %xmm1
mulsd %xmm3, %xmm1
addsd %xmm0, %xmm1
movapd %xmm1, %xmm2
mulsd %xmm4, %xmm2
subsd %xmm2, %xmm0
movsd %xmm0, (%rcx,%r15,8)
movq 0x38(%rbx), %rcx
movq 0x8(%rcx,%rax,8), %rcx
movsd 0x8(%rcx,%r15,8), %xmm0
mulsd %xmm5, %xmm1
subsd %xmm1, %xmm0
movsd %xmm0, 0x8(%rcx,%r15,8)
incq %rax
jmp 0x1ac5f
incq %r15
movq %r15, %r8
movq 0x1e8(%rsp), %rbp
jmp 0x1a7a7
movq 0xf8(%rsp), %rax
incq %rax
movq 0x60(%rsp), %rcx
movq %rcx, 0x20(%rsp)
movq 0x70(%rsp), %rbp
jmp 0x1a1c5
movapd 0x120(%rsp), %xmm2
movsd 0x88(%rsp), %xmm1
movsd 0x98(%rsp), %xmm0
jmp 0x1a3cc
movsd 0x10(%rsp), %xmm0
movsd (%rsp), %xmm1
movq %rbp, %rdi
callq 0x230b0
addsd (%rsp), %xmm0
xorpd %xmm2, %xmm2
jmp 0x1a3c8
testq %r15, %r15
je 0x1adb2
cmpq $-0x1, %r15
jne 0x1b2d6
movq 0x38(%r14), %rax
movq -0x8(%rax,%r13,8), %rcx
movq (%rax,%r13,8), %rax
movsd -0x8(%rcx,%r13,8), %xmm0
movaps %xmm0, 0x20(%rsp)
movsd (%rcx,%r13,8), %xmm0
movaps %xmm0, 0x10(%rsp)
movsd -0x8(%rax,%r13,8), %xmm0
movsd (%rax,%r13,8), %xmm1
movapd %xmm1, 0x60(%rsp)
xorpd %xmm1, %xmm1
movapd %xmm0, (%rsp)
callq 0x241f5
testb %al, %al
je 0x1ade0
movaps 0x20(%rsp), %xmm0
movaps %xmm0, 0x50(%rsp)
movaps 0x60(%rsp), %xmm0
movaps %xmm0, 0x20(%rsp)
movsd 0x15455(%rip), %xmm0    # 0x301f8
movaps %xmm0, 0x30(%rsp)
movapd (%rsp), %xmm0
jmp 0x1b0a9
movq 0x38(%r14), %rax
movq (%rax,%r13,8), %rax
movsd (%rax,%r13,8), %xmm0
movq 0x28(%rbp), %rax
movsd %xmm0, (%rax,%r13,8)
movq 0x1a8(%rsp), %rax
movq 0x28(%rax), %rax
andq $0x0, (%rax,%r13,8)
jmp 0x1b2d6
xorpd %xmm1, %xmm1
movaps 0x10(%rsp), %xmm0
callq 0x241f5
testb %al, %al
je 0x1ae27
movapd (%rsp), %xmm0
xorpd 0x153b1(%rip), %xmm0    # 0x301b0
movaps 0x60(%rsp), %xmm1
movaps %xmm1, 0x50(%rsp)
movapd %xmm0, 0x10(%rsp)
xorpd %xmm0, %xmm0
movapd 0x150(%rsp), %xmm1
movapd %xmm1, 0x30(%rsp)
jmp 0x1b0a9
movapd 0x20(%rsp), %xmm0
subsd 0x60(%rsp), %xmm0
xorpd %xmm1, %xmm1
movsd %xmm0, 0x50(%rsp)
callq 0x241f5
testb %al, %al
je 0x1ae6e
movaps 0x10(%rsp), %xmm0
callq 0x230eb
movq %rax, %rbp
movapd (%rsp), %xmm0
callq 0x230eb
cmpq %rax, %rbp
movq 0x1a0(%rsp), %rbp
jne 0x1ad87
movsd 0x50(%rsp), %xmm1
movsd 0x1560c(%rip), %xmm0    # 0x30488
mulsd %xmm0, %xmm1
movapd %xmm1, 0x30(%rsp)
movaps 0x10(%rsp), %xmm0
movq 0x1e8(%rsp), %rdi
callq 0x2440c
movsd %xmm0, 0xa0(%rsp)
movaps (%rsp), %xmm0
movq 0x1e8(%rsp), %rdi
callq 0x2440c
movaps %xmm0, %xmm1
movsd 0xa0(%rsp), %xmm0
movq 0x1e8(%rsp), %rdi
callq 0x244b6
movsd %xmm0, 0xa0(%rsp)
movaps 0x10(%rsp), %xmm0
movq 0x1e8(%rsp), %rdi
callq 0x2440c
movsd %xmm0, 0x70(%rsp)
movaps (%rsp), %xmm0
movq 0x1e8(%rsp), %rdi
callq 0x2440c
movaps %xmm0, %xmm1
movsd 0x70(%rsp), %xmm0
movq 0x1e8(%rsp), %rdi
callq 0x244bb
movsd %xmm0, 0x70(%rsp)
movaps 0x10(%rsp), %xmm0
callq 0x230eb
cvtsi2sd %rax, %xmm0
mulsd 0x70(%rsp), %xmm0
movsd %xmm0, 0x70(%rsp)
movaps (%rsp), %xmm0
callq 0x230eb
cvtsi2sd %rax, %xmm0
mulsd 0x70(%rsp), %xmm0
movsd %xmm0, 0x70(%rsp)
movapd 0x30(%rsp), %xmm0
movq 0x1e8(%rsp), %rdi
callq 0x2440c
movsd 0xa0(%rsp), %xmm1
movq 0x1e8(%rsp), %rdi
callq 0x244b6
movapd %xmm0, %xmm3
movapd 0x30(%rsp), %xmm1
movapd %xmm1, %xmm0
divsd %xmm3, %xmm0
movsd 0xa0(%rsp), %xmm2
movsd %xmm3, 0xc8(%rsp)
divsd %xmm3, %xmm2
mulsd 0x70(%rsp), %xmm2
mulsd %xmm1, %xmm0
addsd %xmm2, %xmm0
movapd %xmm0, 0xe0(%rsp)
movsd 0x17ab2(%rip), %xmm1    # 0x32a70
callq 0x24295
testb %al, %al
je 0x1b2e6
movsd 0xc8(%rsp), %xmm0
movq 0x1e8(%rsp), %rdi
callq 0x24424
movsd %xmm0, 0x20(%rsp)
movapd 0xe0(%rsp), %xmm0
movq 0x1e8(%rsp), %rdi
callq 0x24424
mulsd 0x20(%rsp), %xmm0
movaps 0x30(%rsp), %xmm1
movq 0x1e8(%rsp), %rdi
callq 0x230b0
movapd %xmm0, %xmm2
addsd 0x30(%rsp), %xmm2
movapd %xmm2, 0x30(%rsp)
movapd 0x60(%rsp), %xmm1
movapd %xmm1, %xmm0
addsd %xmm2, %xmm0
movapd %xmm0, 0x50(%rsp)
movsd 0xa0(%rsp), %xmm0
divsd %xmm2, %xmm0
mulsd 0x70(%rsp), %xmm0
subsd %xmm0, %xmm1
movapd %xmm1, 0x60(%rsp)
movapd (%rsp), %xmm0
movapd %xmm2, %xmm1
movq 0x1e8(%rsp), %rdi
callq 0x15e46
movapd (%rsp), %xmm1
movapd 0x30(%rsp), %xmm2
unpcklpd %xmm1, %xmm2            # xmm2 = xmm2[0],xmm1[0]
unpcklpd %xmm0, %xmm0            # xmm0 = xmm0[0,0]
divpd %xmm0, %xmm2
movapd %xmm2, 0x30(%rsp)
movapd 0x10(%rsp), %xmm0
subsd %xmm1, %xmm0
movapd %xmm0, 0x10(%rsp)
xorpd %xmm0, %xmm0
movapd 0x60(%rsp), %xmm1
movapd %xmm1, 0x20(%rsp)
movapd %xmm0, (%rsp)
xorpd %xmm1, %xmm1
callq 0x241f5
xorpd %xmm0, %xmm0
xorpd %xmm3, %xmm3
testb %al, %al
jne 0x1b11b
movaps 0x10(%rsp), %xmm0
movq %r15, %r13
movq 0x1e8(%rsp), %r15
movq %r15, %rdi
callq 0x2440c
movq %r15, %rdi
callq 0x24424
movsd %xmm0, 0x60(%rsp)
movapd (%rsp), %xmm0
movq %r15, %rdi
callq 0x2440c
movq %r15, %rdi
movq %r13, %r15
movq 0xb8(%rsp), %r13
callq 0x24424
mulsd 0x60(%rsp), %xmm0
movapd %xmm0, %xmm3
xorpd 0x15095(%rip), %xmm3    # 0x301b0
movq 0x28(%rbp), %rax
movapd 0x50(%rsp), %xmm2
movsd %xmm2, -0x8(%rax,%r13,8)
movq 0x1a8(%rsp), %rcx
movq 0x28(%rcx), %rax
movsd %xmm0, -0x8(%rax,%r13,8)
movq 0x28(%rbp), %rax
movapd 0x20(%rsp), %xmm1
movsd %xmm1, (%rax,%r13,8)
movq 0x28(%rcx), %rax
movsd %xmm3, (%rax,%r13,8)
movq 0x38(%r14), %rax
movq -0x8(%rax,%r13,8), %rax
movsd %xmm2, -0x8(%rax,%r13,8)
movq 0x38(%r14), %rax
movq -0x8(%rax,%r13,8), %rax
movaps 0x10(%rsp), %xmm0
movsd %xmm0, (%rax,%r13,8)
movq 0x38(%r14), %rax
movq (%rax,%r13,8), %rax
movapd (%rsp), %xmm0
movsd %xmm0, -0x8(%rax,%r13,8)
movq 0x38(%r14), %rax
movq (%rax,%r13,8), %rax
movsd %xmm1, (%rax,%r13,8)
cmpb $0x0, 0xdc(%rsp)
je 0x1b264
cmpq %r12, %r13
jge 0x1b20b
movq 0x1d0(%rsp), %r9
movq 0x28(%r9), %rax
movapd 0x30(%rsp), %xmm0
movlpd %xmm0, 0x8(%rax)
movq 0x1d8(%rsp), %r10
movq 0x28(%r10), %rax
movhpd %xmm0, 0x8(%rax)
subq $0x8, %rsp
pushq $0x1
popq %rdi
movq 0xb8(%rsp), %rsi
movq %r13, %rdx
movq 0x110(%rsp), %rcx
movq %r12, %r8
pushq 0x1c8(%rsp)
pushq %r14
pushq %r10
callq 0x18253
addq $0x20, %rsp
movq 0x1d0(%rsp), %r9
movq 0x28(%r9), %rax
movapd 0x30(%rsp), %xmm0
movlpd %xmm0, 0x8(%rax)
movq 0x1d8(%rsp), %r10
movq 0x28(%r10), %rax
movhpd %xmm0, 0x8(%rax)
subq $0x8, %rsp
pushq $0x1
popq %rdi
pushq $0x1
popq %rsi
movq 0x108(%rsp), %rdx
movq 0xb8(%rsp), %rcx
movq %r13, %r8
pushq 0x1c8(%rsp)
pushq %r14
pushq %r10
callq 0x187e0
addq $0x20, %rsp
cmpb $0x0, 0x94(%rsp)
movq 0xc0(%rsp), %rcx
je 0x1b2d6
movq 0x1d0(%rsp), %r9
movq 0x28(%r9), %rax
movapd 0x30(%rsp), %xmm0
movlpd %xmm0, 0x8(%rax)
movq 0x1d8(%rsp), %r10
movq 0x28(%r10), %rax
movhpd %xmm0, 0x8(%rax)
subq $0x8, %rsp
pushq $0x1
popq %rdi
pushq $0x1
popq %rsi
movq 0x1b8(%rsp), %rdx
movq 0xb8(%rsp), %rcx
movq %r13, %r8
pushq 0x1c8(%rsp)
pushq %rbx
pushq %r10
callq 0x187e0
movq 0xe0(%rsp), %rcx
addq $0x20, %rsp
addq %r15, %r13
decq %r13
movq 0x48(%rsp), %rax
jmp 0x1a160
movapd 0x10(%rsp), %xmm0
addsd (%rsp), %xmm0
movsd %xmm0, 0xa0(%rsp)
movsd 0x50(%rsp), %xmm1
movq 0x1e8(%rsp), %rdi
callq 0x15e46
movsd %xmm0, 0x50(%rsp)
movsd 0xa0(%rsp), %xmm0
movq 0x1e8(%rsp), %rdi
callq 0x2440c
divsd 0x50(%rsp), %xmm0
movsd 0x14ec1(%rip), %xmm1    # 0x301f8
addsd %xmm1, %xmm0
movsd 0x15145(%rip), %xmm1    # 0x30488
mulsd %xmm1, %xmm0
movq 0x1e8(%rsp), %rdi
callq 0x24424
movapd %xmm0, %xmm1
movapd %xmm0, 0x70(%rsp)
movapd 0x14e4a(%rip), %xmm0    # 0x301b0
movapd 0x30(%rsp), %xmm2
xorpd %xmm0, %xmm2
movsd 0x50(%rsp), %xmm0
mulsd %xmm1, %xmm0
divsd %xmm0, %xmm2
movapd %xmm2, 0x30(%rsp)
movsd 0x14e6c(%rip), %xmm0    # 0x301f8
movsd 0xa0(%rsp), %xmm1
movq 0x1e8(%rsp), %rdi
callq 0x230b0
mulsd 0x30(%rsp), %xmm0
movapd 0x10(%rsp), %xmm3
movapd %xmm3, %xmm1
mulsd %xmm0, %xmm1
movapd 0x20(%rsp), %xmm2
movapd %xmm2, %xmm5
movapd 0x70(%rsp), %xmm8
mulsd %xmm8, %xmm5
addsd %xmm1, %xmm5
mulsd %xmm8, %xmm3
mulsd %xmm0, %xmm2
subsd %xmm2, %xmm3
movapd (%rsp), %xmm7
movapd %xmm7, %xmm2
movapd 0x14dc2(%rip), %xmm4    # 0x301b0
xorpd %xmm4, %xmm2
movapd %xmm3, %xmm1
movapd %xmm3, %xmm6
xorpd %xmm4, %xmm1
movapd 0x60(%rsp), %xmm3
unpcklpd %xmm3, %xmm3            # xmm3 = xmm3[0,0]
movapd %xmm8, %xmm4
unpcklpd %xmm0, %xmm4            # xmm4 = xmm4[0],xmm0[0]
mulpd %xmm4, %xmm3
unpcklpd %xmm7, %xmm2            # xmm2 = xmm2[0],xmm7[0]
movapd %xmm0, %xmm7
unpcklpd %xmm8, %xmm7            # xmm7 = xmm7[0],xmm8[0]
mulpd %xmm7, %xmm2
addpd %xmm3, %xmm2
movapd %xmm0, %xmm3
mulsd %xmm2, %xmm3
mulsd %xmm8, %xmm6
addsd %xmm3, %xmm6
movapd %xmm6, 0x10(%rsp)
movapd %xmm4, 0x30(%rsp)
movapd %xmm4, %xmm3
mulpd %xmm2, %xmm3
unpckhpd %xmm2, %xmm2            # xmm2 = xmm2[1,1]
mulsd %xmm8, %xmm2
unpcklpd %xmm5, %xmm1            # xmm1 = xmm1[0],xmm5[0]
movapd %xmm0, 0x60(%rsp)
mulsd %xmm0, %xmm5
subsd %xmm5, %xmm2
movapd %xmm7, 0x20(%rsp)
mulpd %xmm7, %xmm1
addpd %xmm3, %xmm1
movapd %xmm1, %xmm0
unpckhpd %xmm1, %xmm0            # xmm0 = xmm0[1],xmm1[1]
addsd %xmm1, %xmm0
mulsd 0x14ffc(%rip), %xmm0    # 0x30488
movapd %xmm0, 0x50(%rsp)
xorpd %xmm1, %xmm1
movapd %xmm2, (%rsp)
movapd %xmm2, %xmm0
callq 0x24217
testb %al, %al
je 0x1b60a
xorpd %xmm1, %xmm1
movaps 0x10(%rsp), %xmm0
callq 0x24217
testb %al, %al
je 0x1b619
movaps 0x10(%rsp), %xmm0
callq 0x230eb
movq %rax, %rbp
movaps (%rsp), %xmm0
callq 0x230eb
cmpq %rax, %rbp
jne 0x1b65f
movaps 0x10(%rsp), %xmm0
movq %r15, %r13
movq 0x1e8(%rsp), %r15
movq %r15, %rdi
callq 0x2440c
movq %r15, %rdi
callq 0x24424
movaps %xmm0, 0x30(%rsp)
movapd (%rsp), %xmm0
movq %r15, %rdi
callq 0x2440c
movq %r15, %rdi
callq 0x24424
movapd %xmm0, %xmm1
movapd %xmm0, 0xa0(%rsp)
movapd 0x30(%rsp), %xmm0
mulsd %xmm1, %xmm0
movaps (%rsp), %xmm1
movq %r15, %rdi
callq 0x230b0
movsd %xmm0, 0x20(%rsp)
movapd 0x10(%rsp), %xmm0
addsd (%rsp), %xmm0
movq %r15, %rdi
callq 0x2440c
movq %r15, %rdi
movq %r13, %r15
movq 0xb8(%rsp), %r13
callq 0x24424
movsd 0x14c86(%rip), %xmm1    # 0x301f8
divsd %xmm0, %xmm1
movapd 0x50(%rsp), %xmm3
movapd %xmm3, %xmm0
movsd 0x20(%rsp), %xmm2
addsd %xmm2, %xmm3
subsd %xmm2, %xmm0
movapd %xmm0, 0x20(%rsp)
movapd 0x10(%rsp), %xmm0
subsd (%rsp), %xmm0
movapd %xmm0, 0x10(%rsp)
movapd 0x60(%rsp), %xmm0
movapd %xmm0, %xmm2
xorpd 0x14bf9(%rip), %xmm2    # 0x301b0
movapd 0x30(%rsp), %xmm4
unpcklpd 0xa0(%rsp), %xmm4       # xmm4 = xmm4[0],mem[0]
unpcklpd %xmm1, %xmm1            # xmm1 = xmm1[0,0]
mulpd %xmm4, %xmm1
movapd 0x70(%rsp), %xmm4
unpcklpd %xmm4, %xmm4            # xmm4 = xmm4[0,0]
mulpd %xmm1, %xmm4
shufpd $0x1, %xmm1, %xmm1      # xmm1 = xmm1[1,0]
unpcklpd %xmm0, %xmm2            # xmm2 = xmm2[0],xmm0[0]
mulpd %xmm1, %xmm2
addpd %xmm2, %xmm4
xorpd %xmm0, %xmm0
movapd %xmm3, 0x50(%rsp)
movapd %xmm4, 0x30(%rsp)
movq 0x1a0(%rsp), %rbp
jmp 0x1b0a9
movaps 0x50(%rsp), %xmm0
movaps %xmm0, 0x20(%rsp)
jmp 0x1ada8
movapd 0x14b8f(%rip), %xmm0    # 0x301b0
movapd (%rsp), %xmm1
xorpd %xmm0, %xmm1
movapd 0x60(%rsp), %xmm2
xorpd %xmm0, %xmm2
movapd 0x20(%rsp), %xmm3
movsd %xmm2, %xmm3            # xmm3 = xmm2[0],xmm3[1]
movapd %xmm1, 0x10(%rsp)
xorpd %xmm0, %xmm0
movapd 0x50(%rsp), %xmm1
movapd %xmm1, 0x20(%rsp)
movapd %xmm3, 0x30(%rsp)
jmp 0x1b0a9
movaps 0x50(%rsp), %xmm0
movaps %xmm0, 0x20(%rsp)
movq 0x1a0(%rsp), %rbp
jmp 0x1ada8
movq 0x1e0(%rsp), %rax
movq %r13, (%rax)
addq $0x168, %rsp            # imm = 0x168
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq