// File: /workspace/github_cpp_makefile_star2_O2/jkingan[P]ladyheather/heathmsc.cpp
// Function: step2diu_(double*, double*, double*, double*)
int step2diu_(double *xsta, double *fhr, double *t, double *xcorsta)
{
    /* Initialized data */

    static double deg2rad = .017453292519943295769;

    /* F2C generated locals */
    double d__1, d__2, d__3, d__4;

    /* Local variables */
    static double rsta, h;
    static int i, j;
    static double p, s, cosla, sinla, de, 
            dn, dr, pr, ps, thetaf, cosphi, sinphi, zla, tau, zns;

   if(xsta == 0) return (-1);
   if(fhr == 0) return (-1);
   if(t == 0) return (-1);
   if(xcorsta == 0) return (-1);

/* last change:  vd   17 may 00   1:20 pm */
/* these are the subroutines for the step2 of the tidal corrections. */

/* they are called to account for the frequency dependence */
/* of the love numbers. */

    /* Function Body */
/*** note, following table is derived from dehanttideinelMJD.f (2000oct30 16:10) */
/*** has minor differences from that of dehanttideinel.f (2000apr17 14:10) */
/*** D.M. edited to strictly follow published table 7.5a (2006aug08 13:46) */
/* cf. table 7.5a of IERS conventions 2003 (TN.32, pg.82) */
/* columns are s,h,p,N',ps, dR(ip),dR(op),dT(ip),dT(op) */
/* units of mm */
/*   data datdi/ */
/*  * -3., 0., 2., 0., 0.,-0.01,-0.01, 0.0 , 0.0, */
/*  * -3., 2., 0., 0., 0.,-0.01,-0.01, 0.0 , 0.0, */
/*  * -2., 0., 1.,-1., 0.,-0.02,-0.01, 0.0 , 0.0, */
/*  * -2., 0., 1., 0., 0.,-0.08, 0.00, 0.01, 0.01,      !*** table 7.5a
 */
/*  * -2., 2.,-1., 0., 0.,-0.02,-0.01, 0.0 , 0.0, */
/*  * -1., 0., 0.,-1., 0.,-0.10, 0.00, 0.00, 0.00,      !*** table 7.5a
 */
/*  * -1., 0., 0., 0., 0.,-0.51, 0.00,-0.02, 0.03,      !*** table 7.5a
 */
/*  * -1., 2., 0., 0., 0., 0.01, 0.0 , 0.0 , 0.0, */
/*  *  0.,-2., 1., 0., 0., 0.01, 0.0 , 0.0 , 0.0, */
/*  *  0., 0.,-1., 0., 0., 0.02, 0.01, 0.0 , 0.0, */
/*  *  0., 0., 1., 0., 0., 0.06, 0.00, 0.00, 0.00,      !*** table 7.5a
 */
/*  *  0., 0., 1., 1., 0., 0.01, 0.0 , 0.0 , 0.0, */
/*  *  0., 2.,-1., 0., 0., 0.01, 0.0 , 0.0 , 0.0, */
/*  *  1.,-3., 0., 0., 1.,-0.06, 0.00, 0.00, 0.00,      !*** table 7.5a
 */
/*  *  1.,-2., 0., 1., 0., 0.01, 0.0 , 0.0 , 0.0, */
/*  *  1.,-2., 0., 0., 0.,-1.23,-0.07, 0.06, 0.01,      !*** table 7.5a
 */
/*  *  1.,-1., 0., 0.,-1., 0.02, 0.0 , 0.0 , 0.0, */
/*  *  1.,-1., 0., 0., 1., 0.04, 0.0 , 0.0 , 0.0, */
/*  *  1., 0., 0.,-1., 0.,-0.22, 0.01, 0.01, 0.00,      !*** table 7.5a
 */
/***  *  1., 0., 0., 0., 0.,12.00,-0.78,-0.67,-0.03,      !**** table 7.5a */
/*  *  1., 0., 0., 1., 0., 1.73,-0.12,-0.10, 0.00,      !*** table 7.5a */
/*  *  1., 0., 0., 2., 0.,-0.04, 0.0 , 0.0 , 0.0, */
/*  *  1., 1., 0., 0.,-1.,-0.50,-0.01, 0.03, 0.00,      !*** table 7.5a */
/*  *  1., 1., 0., 0., 1., 0.01, 0.0 , 0.0 , 0.0, */
/***  *  1., 1., 0., 1.,-1.,-0.01, 0.0 , 0.0 , 0.0,       !*** v.dehant 2007 */
/*  *  1., 2.,-2., 0., 0.,-0.01, 0.0 , 0.0 , 0.0, */
/*  *  1., 2., 0., 0., 0.,-0.11, 0.01, 0.01, 0.00,      !*** table 7.5a */
/*  *  2.,-2., 1., 0., 0.,-0.01, 0.0 , 0.0 , 0.0, */
/*  *  2., 0.,-1., 0., 0.,-0.02, 0.02, 0.0 , 0.01, */
/*  *  3., 0., 0., 0., 0., 0.0 , 0.01, 0.0 , 0.01, */
/*  *  3., 0., 0., 1., 0., 0.0 , 0.01, 0.0 , 0.0/ */
/* Computing 3rd power */
    d__1 = *t, d__2 = d__1;
    s = *t * 481267.88194 + 218.31664563 - *t * .0014663889 * *t + d__2 * (
            d__1 * d__1) * 1.85139e-6;
/* Computing 3rd power */
    d__1 = *t, d__2 = d__1;
    tau = *fhr * 15. + 280.4606184 + *t * 36000.7700536 + *t * 3.8793e-4 * *t 
            - d__2 * (d__1 * d__1) * 2.58e-8 - s;
/* Computing 3rd power */
    d__1 = *t, d__2 = d__1;
/* Computing 4th power */
    d__3 = *t, d__3 *= d__3;
    pr = *t * 1.396971278f + *t * 3.08889e-4f * *t + d__2 * (d__1 * d__1) * 
            2.1e-8f + d__3 * d__3 * 7e-9f;
    s += pr;
/* Computing 3rd power */
    d__1 = *t, d__2 = d__1;
/* Computing 4th power */
    d__3 = *t, d__3 *= d__3;
    h = *t * 36000.7697489 + 280.46645 + *t * 3.0322222e-4 * *t + d__2 * (
            d__1 * d__1) * 2e-8f - d__3 * d__3 * 6.54e-9f;
/* Computing 3rd power */
    d__1 = *t, d__2 = d__1;
/* Computing 4th power */
    d__3 = *t, d__3 *= d__3;
    p = *t * 4069.01363525 + 83.35324312 - *t * .01032172222 * *t - d__2 * (
            d__1 * d__1) * 1.24991e-5 + d__3 * d__3 * 5.263e-8;
/* Computing 3rd power */
    d__1 = *t, d__2 = d__1;
/* Computing 4th power */
    d__3 = *t, d__3 *= d__3;
    zns = *t * 1934.13626197 + 234.95544499 - *t * .00207561111 * *t - d__2 * 
            (d__1 * d__1) * 2.13944e-6 + d__3 * d__3 * 1.65e-8;
/* Computing 3rd power */
    d__1 = *t, d__2 = d__1;
/* Computing 4th power */
    d__3 = *t, d__3 *= d__3;
    ps = *t * 1.71945766667 + 282.93734098 + *t * 4.5688889e-4 * *t - d__2 * (
            d__1 * d__1) * 1.778e-8 - d__3 * d__3 * 3.34e-9;
/* reduce angles to between 0 and 360 */
    s = d_mod(&s, &c_b56);
    tau = d_mod(&tau, &c_b56);
    h = d_mod(&h, &c_b56);
    p = d_mod(&p, &c_b56);
    zns = d_mod(&zns, &c_b56);
    ps = d_mod(&ps, &c_b56);
/* Computing 2nd power */
    d__1 = xsta[0];
/* Computing 2nd power */
    d__2 = xsta[1];
/* Computing 2nd power */
    d__3 = xsta[2];
    rsta = sqrt(d__1 * d__1 + d__2 * d__2 + d__3 * d__3);
    sinphi = xsta[2] / rsta;
/* Computing 2nd power */
    d__1 = xsta[0];
/* Computing 2nd power */
    d__2 = xsta[1];
    cosphi = sqrt(d__1 * d__1 + d__2 * d__2) / rsta;
    cosla = xsta[0] / cosphi / rsta;
    sinla = xsta[1] / cosphi / rsta;
    zla = atan2(xsta[1], xsta[0]);
    for(i=0; i<=2; ++i) {
        xcorsta[i] = 0.;
    }
    for(j=1; j<=31; ++j) {
        thetaf = (tau + datdi[j * 9 - 9] * s + datdi[j * 9 - 8] * h + datdi[j 
                * 9 - 7] * p + datdi[j * 9 - 6] * zns + datdi[j * 9 - 5] * ps)
                 * deg2rad;
        dr = datdi[j * 9 - 4] * 2. * sinphi * cosphi * sin(thetaf + zla) + 
                datdi[j * 9 - 3] * 2. * sinphi * cosphi * cos(thetaf + zla);
/* Computing 2nd power */
        d__1 = cosphi;
/* Computing 2nd power */
        d__2 = sinphi;
/* Computing 2nd power */
        d__3 = cosphi;
/* Computing 2nd power */
        d__4 = sinphi;
        dn = datdi[j * 9 - 2] * (d__1 * d__1 - d__2 * d__2) * sin(thetaf + 
                zla) + datdi[j * 9 - 1] * (d__3 * d__3 - d__4 * d__4) * cos(
                thetaf + zla);
/***** following correction by V.Dehant to match eq.16b, p.81, 2003 Co
nventions*/
/* ****   de=datdi(8,j)*sinphi*cos(thetaf+zla)+ */
        de = datdi[j * 9 - 2] * sinphi * cos(thetaf + zla) - datdi[j * 9 - 1] 
                * sinphi * sin(thetaf + zla);
        xcorsta[0] = xcorsta[0] + dr * cosla * cosphi - de * sinla - dn * 
                sinphi * cosla;
        xcorsta[1] = xcorsta[1] + dr * sinla * cosphi + de * cosla - dn * 
                sinphi * sinla;
        xcorsta[2] = xcorsta[2] + dr * sinphi + dn * cosphi;
    }
    for(i=0; i<= 2; ++i) {
        xcorsta[i] /= 1e3;
    }
    return 0;
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r14
pushq %rbx
subq $0x20, %rsp
movq %rcx, %rbx
movq %rdi, %r14
testq %rdi, %rdi
sete %al
testq %rsi, %rsi
sete %cl
orb %al, %cl
testq %rdx, %rdx
sete %al
testq %rbx, %rbx
sete %dil
orb %al, %dil
orb %cl, %dil
je 0xd8a21
pushq $-0x1
popq %rax
addq $0x20, %rsp
popq %rbx
popq %r14
popq %rbp
retq
movsd (%rdx), %xmm14
movsd 0x3d9aa(%rip), %xmm0    # 0x1163d8
mulsd %xmm14, %xmm0
addsd 0x3d9a5(%rip), %xmm0    # 0x1163e0
movsd 0x3d9a5(%rip), %xmm2    # 0x1163e8
mulsd %xmm14, %xmm2
mulsd %xmm14, %xmm2
addsd %xmm0, %xmm2
movapd %xmm14, %xmm1
movsd 0x3d9a2(%rip), %xmm4    # 0x116400
mulsd %xmm14, %xmm4
movsd 0x3d99d(%rip), %xmm3    # 0x116408
mulsd %xmm14, %xmm3
mulsd %xmm14, %xmm3
movsd 0x3d99b(%rip), %xmm0    # 0x116418
mulsd %xmm14, %xmm0
mulsd %xmm14, %xmm0
movsd %xmm0, -0x18(%rbp)
movsd 0x3d98c(%rip), %xmm6    # 0x116420
mulsd %xmm14, %xmm6
movsd 0x3d997(%rip), %xmm7    # 0x116438
mulsd %xmm14, %xmm7
movsd 0x3d999(%rip), %xmm8    # 0x116448
mulsd %xmm14, %xmm8
mulsd %xmm14, %xmm8
movsd 0x3d99e(%rip), %xmm9    # 0x116460
mulsd %xmm14, %xmm9
movsd 0x3d9a0(%rip), %xmm10   # 0x116470
mulsd %xmm14, %xmm10
mulsd %xmm14, %xmm10
movsd 0x3d9a5(%rip), %xmm11   # 0x116488
mulsd %xmm14, %xmm11
movsd 0x3d9a7(%rip), %xmm12   # 0x116498
mulsd %xmm14, %xmm12
mulsd %xmm14, %xmm12
movsd 0x3d9ac(%rip), %xmm13   # 0x1164b0
mulsd %xmm14, %xmm13
movapd %xmm14, %xmm0
mulsd %xmm14, %xmm0
mulsd %xmm0, %xmm1
movsd 0x3d8d0(%rip), %xmm15   # 0x1163f0
mulsd %xmm1, %xmm15
addsd %xmm2, %xmm15
movsd (%rsi), %xmm2
mulsd 0x2f9f2(%rip), %xmm2    # 0x108528
addsd 0x3d8ba(%rip), %xmm2    # 0x1163f8
addsd %xmm2, %xmm4
movsd 0x3d976(%rip), %xmm2    # 0x1164c0
mulsd %xmm14, %xmm2
mulsd %xmm14, %xmm2
addsd %xmm4, %xmm3
movapd 0x1a650(%rip), %xmm4    # 0xf31b0
xorpd %xmm4, %xmm14
mulsd %xmm0, %xmm14
movsd 0x3d89e(%rip), %xmm5    # 0x116410
mulsd %xmm14, %xmm5
addsd %xmm3, %xmm5
subsd %xmm15, %xmm5
movsd %xmm5, 0x28cad8(%rip)   # 0x365660
addsd -0x18(%rbp), %xmm6
movsd 0x3d893(%rip), %xmm3    # 0x116428
mulsd %xmm1, %xmm3
addsd %xmm6, %xmm3
xorpd %xmm0, %xmm4
mulsd %xmm0, %xmm4
movapd %xmm0, %xmm5
mulsd %xmm0, %xmm5
movsd 0x3d87b(%rip), %xmm0    # 0x116430
mulsd %xmm5, %xmm0
addsd %xmm3, %xmm0
addsd %xmm15, %xmm0
addsd 0x3d876(%rip), %xmm7    # 0x116440
mulsd 0x3d87e(%rip), %xmm1    # 0x116450
addsd %xmm7, %xmm8
addsd %xmm8, %xmm1
movsd 0x3d874(%rip), %xmm3    # 0x116458
mulsd %xmm4, %xmm3
addsd %xmm1, %xmm3
movsd %xmm3, 0x28ca04(%rip)   # 0x3655f8
addsd 0x3d86b(%rip), %xmm9    # 0x116468
addsd %xmm9, %xmm10
movsd 0x3d86e(%rip), %xmm1    # 0x116478
mulsd %xmm14, %xmm1
addsd %xmm10, %xmm1
movsd 0x3d864(%rip), %xmm3    # 0x116480
mulsd %xmm5, %xmm3
addsd %xmm1, %xmm3
movsd %xmm3, 0x28c9dc(%rip)   # 0x365608
addsd 0x3d85b(%rip), %xmm11   # 0x116490
addsd %xmm11, %xmm12
movsd 0x3d85e(%rip), %xmm1    # 0x1164a0
mulsd %xmm14, %xmm1
addsd %xmm12, %xmm1
mulsd 0x3d854(%rip), %xmm5    # 0x1164a8
addsd %xmm1, %xmm5
movsd %xmm5, 0x28ca08(%rip)   # 0x365668
addsd 0x3d84f(%rip), %xmm13   # 0x1164b8
addsd %xmm13, %xmm2
mulsd 0x3d851(%rip), %xmm14   # 0x1164c8
addsd %xmm2, %xmm14
mulsd 0x3d84c(%rip), %xmm4    # 0x1164d0
addsd %xmm14, %xmm4
movsd %xmm4, 0x28c9a7(%rip)   # 0x365638
movsd 0x2f517(%rip), %xmm1    # 0x1081b0
callq 0x94d0
movsd %xmm0, 0x28c96a(%rip)   # 0x365610
movsd 0x28c9b2(%rip), %xmm0   # 0x365660
movsd 0x2f4fa(%rip), %xmm1    # 0x1081b0
callq 0x94d0
movsd %xmm0, 0x28c99d(%rip)   # 0x365660
movsd 0x28c92d(%rip), %xmm0   # 0x3655f8
movsd 0x2f4dd(%rip), %xmm1    # 0x1081b0
callq 0x94d0
movsd %xmm0, 0x28c918(%rip)   # 0x3655f8
movsd 0x28c920(%rip), %xmm0   # 0x365608
movsd 0x2f4c0(%rip), %xmm1    # 0x1081b0
callq 0x94d0
movsd %xmm0, 0x28c90b(%rip)   # 0x365608
movsd 0x28c963(%rip), %xmm0   # 0x365668
movsd 0x2f4a3(%rip), %xmm1    # 0x1081b0
callq 0x94d0
movsd %xmm0, 0x28c94e(%rip)   # 0x365668
movsd 0x28c916(%rip), %xmm0   # 0x365638
movsd 0x2f486(%rip), %xmm1    # 0x1081b0
callq 0x94d0
movsd %xmm0, 0x28c901(%rip)   # 0x365638
movsd (%r14), %xmm1
movsd 0x8(%r14), %xmm0
movsd 0x10(%r14), %xmm2
movapd %xmm0, %xmm3
mulsd %xmm0, %xmm3
movapd %xmm1, %xmm4
mulsd %xmm1, %xmm4
addsd %xmm3, %xmm4
movapd %xmm2, %xmm3
mulsd %xmm2, %xmm3
addsd %xmm4, %xmm3
sqrtsd %xmm3, %xmm3
divsd %xmm3, %xmm2
movsd %xmm2, 0x28c8d8(%rip)   # 0x365650
sqrtsd %xmm4, %xmm2
divsd %xmm3, %xmm2
movsd %xmm2, 0x28c8c0(%rip)   # 0x365648
movapd %xmm1, %xmm4
divsd %xmm2, %xmm4
divsd %xmm3, %xmm4
movsd %xmm4, 0x28c87c(%rip)   # 0x365618
movapd %xmm0, %xmm4
divsd %xmm2, %xmm4
divsd %xmm3, %xmm4
movsd %xmm4, 0x28c870(%rip)   # 0x365620
callq 0x9050
movsd %xmm0, 0x28c89b(%rip)   # 0x365658
xorl %eax, %eax
cmpq $0x3, %rax
je 0xd8dcf
andq $0x0, (%rbx,%rax,8)
incq %rax
jmp 0xd8dbf
pushq $0x1
popq %rax
leaq 0x80cd7(%rip), %r14     # 0x159ab0
movl %eax, 0x28c821(%rip)    # 0x365600
cmpl $0x20, %eax
jge 0xd90d5
leal (%rax,%rax,8), %eax
cltq
movsd -0x48(%r14,%rax,8), %xmm0
movsd -0x40(%r14,%rax,8), %xmm1
mulsd 0x28c80d(%rip), %xmm0   # 0x365610
addsd 0x28c855(%rip), %xmm0   # 0x365660
mulsd 0x28c7e5(%rip), %xmm1   # 0x3655f8
movsd -0x38(%r14,%rax,8), %xmm2
mulsd 0x28c7e6(%rip), %xmm2   # 0x365608
addsd %xmm0, %xmm1
movsd -0x30(%r14,%rax,8), %xmm3
mulsd 0x28c833(%rip), %xmm3   # 0x365668
addsd %xmm1, %xmm2
movsd -0x28(%r14,%rax,8), %xmm0
mulsd 0x28c7f0(%rip), %xmm0   # 0x365638
addsd %xmm2, %xmm3
addsd %xmm3, %xmm0
mulsd 0x1cb30(%rip), %xmm0    # 0xf5988
movsd %xmm0, 0x28c7e0(%rip)   # 0x365640
movsd -0x20(%r14,%rax,8), %xmm1
addsd %xmm1, %xmm1
mulsd 0x28c7dd(%rip), %xmm1   # 0x365650
mulsd 0x28c7cd(%rip), %xmm1   # 0x365648
movsd %xmm1, -0x18(%rbp)
addsd 0x28c7d0(%rip), %xmm0   # 0x365658
callq 0x9080
movslq 0x28c76c(%rip), %rax    # 0x365600
imulq $0x48, %rax, %rax
movsd -0x18(%rax,%r14), %xmm1
addsd %xmm1, %xmm1
mulsd 0x28c7a5(%rip), %xmm1   # 0x365650
mulsd 0x28c795(%rip), %xmm1   # 0x365648
movsd %xmm1, -0x20(%rbp)
movsd %xmm0, -0x28(%rbp)
movsd 0x28c77b(%rip), %xmm0   # 0x365640
addsd 0x28c78b(%rip), %xmm0   # 0x365658
callq 0x9710
mulsd -0x20(%rbp), %xmm0
movsd -0x28(%rbp), %xmm1
mulsd -0x18(%rbp), %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, 0x28c743(%rip)   # 0x365630
movsd 0x28c753(%rip), %xmm1   # 0x365648
movsd 0x28c753(%rip), %xmm0   # 0x365650
movslq 0x28c6fc(%rip), %rax    # 0x365600
imulq $0x48, %rax, %rax
mulsd %xmm0, %xmm0
mulsd %xmm1, %xmm1
subsd %xmm0, %xmm1
movsd %xmm1, -0x18(%rbp)
movsd -0x10(%rax,%r14), %xmm0
mulsd %xmm1, %xmm0
movsd %xmm0, -0x28(%rbp)
movsd 0x28c70f(%rip), %xmm0   # 0x365640
addsd 0x28c71f(%rip), %xmm0   # 0x365658
callq 0x9080
movslq 0x28c6bb(%rip), %rax    # 0x365600
imulq $0x48, %rax, %rax
movsd -0x18(%rbp), %xmm1
mulsd -0x8(%rax,%r14), %xmm1
movsd %xmm1, -0x18(%rbp)
movsd %xmm0, -0x20(%rbp)
movsd 0x28c6d9(%rip), %xmm0   # 0x365640
addsd 0x28c6e9(%rip), %xmm0   # 0x365658
callq 0x9710
mulsd -0x18(%rbp), %xmm0
movsd -0x20(%rbp), %xmm1
mulsd -0x28(%rbp), %xmm1
addsd %xmm0, %xmm1
movslq 0x28c672(%rip), %rax    # 0x365600
imulq $0x48, %rax, %rax
movsd -0x10(%rax,%r14), %xmm0
mulsd 0x28c6af(%rip), %xmm0   # 0x365650
movsd %xmm0, -0x18(%rbp)
movsd %xmm1, 0x28c67a(%rip)   # 0x365628
movsd 0x28c68a(%rip), %xmm0   # 0x365640
addsd 0x28c69a(%rip), %xmm0   # 0x365658
callq 0x9710
movsd %xmm0, -0x28(%rbp)
movslq 0x28c631(%rip), %rax    # 0x365600
imulq $0x48, %rax, %rax
movsd 0x28c675(%rip), %xmm1   # 0x365650
movsd 0x28c65d(%rip), %xmm0   # 0x365640
addsd 0x28c66d(%rip), %xmm0   # 0x365658
mulsd -0x8(%rax,%r14), %xmm1
movsd %xmm1, -0x20(%rbp)
callq 0x9080
mulsd -0x20(%rbp), %xmm0
movsd -0x28(%rbp), %xmm6
mulsd -0x18(%rbp), %xmm6
movsd 0x28c61d(%rip), %xmm1   # 0x365630
movsd 0x28c5fd(%rip), %xmm3   # 0x365618
movapd %xmm1, %xmm4
mulsd %xmm3, %xmm4
movsd 0x28c61d(%rip), %xmm2   # 0x365648
mulsd %xmm2, %xmm4
addsd (%rbx), %xmm4
subsd %xmm0, %xmm6
movsd 0x28c5e1(%rip), %xmm5   # 0x365620
movapd %xmm6, %xmm0
movapd %xmm6, %xmm9
mulsd %xmm5, %xmm0
subsd %xmm0, %xmm4
movsd 0x28c5d0(%rip), %xmm0   # 0x365628
movsd 0x28c5f0(%rip), %xmm6   # 0x365650
movapd %xmm0, %xmm7
xorpd 0x1a144(%rip), %xmm7    # 0xf31b0
mulsd %xmm6, %xmm7
movapd %xmm7, %xmm8
mulsd %xmm3, %xmm8
addsd %xmm4, %xmm8
movsd %xmm8, (%rbx)
movapd %xmm1, %xmm4
mulsd %xmm5, %xmm4
mulsd %xmm2, %xmm4
addsd 0x8(%rbx), %xmm4
movapd %xmm9, %xmm8
mulsd %xmm3, %xmm8
addsd %xmm4, %xmm8
mulsd %xmm5, %xmm7
addsd %xmm8, %xmm7
movsd %xmm7, 0x8(%rbx)
mulsd %xmm6, %xmm1
addsd 0x10(%rbx), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, 0x10(%rbx)
movl 0x28c532(%rip), %eax    # 0x365600
incl %eax
jmp 0xd8dd9
xorl %eax, %eax
movsd 0x19f59(%rip), %xmm0    # 0xf3038
xorl %ecx, %ecx
cmpq $0x3, %rcx
je 0xd8a18
movsd (%rbx,%rcx,8), %xmm1
divsd %xmm0, %xmm1
movsd %xmm1, (%rbx,%rcx,8)
incq %rcx
jmp 0xd90e1