// File: /workspace/github_cpp_makefile_star2_O3/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
movl $0xffffffff, %eax       # imm = 0xFFFFFFFF
jne 0xed7dc
movsd (%rdx), %xmm14
movsd 0x3f3a1(%rip), %xmm0    # 0x12c4a8
mulsd %xmm14, %xmm0
addsd 0x3f39c(%rip), %xmm0    # 0x12c4b0
movsd 0x3f39c(%rip), %xmm2    # 0x12c4b8
mulsd %xmm14, %xmm2
mulsd %xmm14, %xmm2
addsd %xmm0, %xmm2
movapd %xmm14, %xmm1
movsd 0x3f399(%rip), %xmm4    # 0x12c4d0
mulsd %xmm14, %xmm4
movsd 0x3f394(%rip), %xmm3    # 0x12c4d8
mulsd %xmm14, %xmm3
mulsd %xmm14, %xmm3
movsd 0x3f392(%rip), %xmm0    # 0x12c4e8
mulsd %xmm14, %xmm0
mulsd %xmm14, %xmm0
movsd %xmm0, -0x18(%rbp)
movsd 0x3f383(%rip), %xmm6    # 0x12c4f0
mulsd %xmm14, %xmm6
movsd 0x3f38e(%rip), %xmm7    # 0x12c508
mulsd %xmm14, %xmm7
movsd 0x3f390(%rip), %xmm8    # 0x12c518
mulsd %xmm14, %xmm8
mulsd %xmm14, %xmm8
movsd 0x3f395(%rip), %xmm9    # 0x12c530
mulsd %xmm14, %xmm9
movsd 0x3f397(%rip), %xmm10   # 0x12c540
mulsd %xmm14, %xmm10
mulsd %xmm14, %xmm10
movsd 0x3f39c(%rip), %xmm11   # 0x12c558
mulsd %xmm14, %xmm11
movsd 0x3f39e(%rip), %xmm12   # 0x12c568
mulsd %xmm14, %xmm12
mulsd %xmm14, %xmm12
movsd 0x3f3a3(%rip), %xmm13   # 0x12c580
mulsd %xmm14, %xmm13
movapd %xmm14, %xmm0
mulsd %xmm14, %xmm0
mulsd %xmm0, %xmm1
movsd 0x3f2c7(%rip), %xmm15   # 0x12c4c0
mulsd %xmm1, %xmm15
addsd %xmm2, %xmm15
movsd (%rsi), %xmm2
mulsd 0x31329(%rip), %xmm2    # 0x11e538
addsd 0x3f2b1(%rip), %xmm2    # 0x12c4c8
addsd %xmm2, %xmm4
movsd 0x3f36d(%rip), %xmm2    # 0x12c590
mulsd %xmm14, %xmm2
mulsd %xmm14, %xmm2
addsd %xmm4, %xmm3
movapd 0x1bf77(%rip), %xmm4    # 0x1091b0
xorpd %xmm4, %xmm14
mulsd %xmm0, %xmm14
movsd 0x3f295(%rip), %xmm5    # 0x12c4e0
mulsd %xmm14, %xmm5
addsd %xmm3, %xmm5
subsd %xmm15, %xmm5
movsd %xmm5, 0x28e3ff(%rip)   # 0x37b660
addsd -0x18(%rbp), %xmm6
movsd 0x3f28a(%rip), %xmm3    # 0x12c4f8
mulsd %xmm1, %xmm3
addsd %xmm6, %xmm3
xorpd %xmm0, %xmm4
mulsd %xmm0, %xmm4
movapd %xmm0, %xmm5
mulsd %xmm0, %xmm5
movsd 0x3f272(%rip), %xmm0    # 0x12c500
mulsd %xmm5, %xmm0
addsd %xmm3, %xmm0
addsd %xmm15, %xmm0
addsd 0x3f26d(%rip), %xmm7    # 0x12c510
mulsd 0x3f275(%rip), %xmm1    # 0x12c520
addsd %xmm7, %xmm8
addsd %xmm8, %xmm1
movsd 0x3f26b(%rip), %xmm3    # 0x12c528
mulsd %xmm4, %xmm3
addsd %xmm1, %xmm3
movsd %xmm3, 0x28e32b(%rip)   # 0x37b5f8
addsd 0x3f262(%rip), %xmm9    # 0x12c538
addsd %xmm9, %xmm10
movsd 0x3f265(%rip), %xmm1    # 0x12c548
mulsd %xmm14, %xmm1
addsd %xmm10, %xmm1
movsd 0x3f25b(%rip), %xmm3    # 0x12c550
mulsd %xmm5, %xmm3
addsd %xmm1, %xmm3
movsd %xmm3, 0x28e303(%rip)   # 0x37b608
addsd 0x3f252(%rip), %xmm11   # 0x12c560
addsd %xmm11, %xmm12
movsd 0x3f255(%rip), %xmm1    # 0x12c570
mulsd %xmm14, %xmm1
addsd %xmm12, %xmm1
mulsd 0x3f24b(%rip), %xmm5    # 0x12c578
addsd %xmm1, %xmm5
movsd %xmm5, 0x28e32f(%rip)   # 0x37b668
addsd 0x3f246(%rip), %xmm13   # 0x12c588
addsd %xmm13, %xmm2
mulsd 0x3f248(%rip), %xmm14   # 0x12c598
addsd %xmm2, %xmm14
mulsd 0x3f243(%rip), %xmm4    # 0x12c5a0
addsd %xmm14, %xmm4
movsd %xmm4, 0x28e2ce(%rip)   # 0x37b638
movsd 0x30dc6(%rip), %xmm1    # 0x11e138
callq 0x94d0
movsd %xmm0, 0x28e291(%rip)   # 0x37b610
movsd 0x28e2d9(%rip), %xmm0   # 0x37b660
movsd 0x30da9(%rip), %xmm1    # 0x11e138
callq 0x94d0
movsd %xmm0, 0x28e2c4(%rip)   # 0x37b660
movsd 0x28e254(%rip), %xmm0   # 0x37b5f8
movsd 0x30d8c(%rip), %xmm1    # 0x11e138
callq 0x94d0
movsd %xmm0, 0x28e23f(%rip)   # 0x37b5f8
movsd 0x28e247(%rip), %xmm0   # 0x37b608
movsd 0x30d6f(%rip), %xmm1    # 0x11e138
callq 0x94d0
movsd %xmm0, 0x28e232(%rip)   # 0x37b608
movsd 0x28e28a(%rip), %xmm0   # 0x37b668
movsd 0x30d52(%rip), %xmm1    # 0x11e138
callq 0x94d0
movsd %xmm0, 0x28e275(%rip)   # 0x37b668
movsd 0x28e23d(%rip), %xmm0   # 0x37b638
movsd 0x30d35(%rip), %xmm1    # 0x11e138
callq 0x94d0
movsd %xmm0, 0x28e228(%rip)   # 0x37b638
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, 0x28e1ff(%rip)   # 0x37b650
xorps %xmm2, %xmm2
sqrtsd %xmm4, %xmm2
divsd %xmm3, %xmm2
movsd %xmm2, 0x28e1e4(%rip)   # 0x37b648
movapd %xmm1, %xmm4
divsd %xmm2, %xmm4
divsd %xmm3, %xmm4
movsd %xmm4, 0x28e1a0(%rip)   # 0x37b618
movapd %xmm0, %xmm4
divsd %xmm2, %xmm4
divsd %xmm3, %xmm4
movsd %xmm4, 0x28e194(%rip)   # 0x37b620
callq 0x9050
movsd %xmm0, 0x28e1bf(%rip)   # 0x37b658
xorpd %xmm0, %xmm0
movupd %xmm0, (%rbx)
movq $0x0, 0x10(%rbx)
movl $0x1, 0x28e14d(%rip)    # 0x37b600
movsd 0x28e195(%rip), %xmm2   # 0x37b650
movl $0x1, %eax
movsd 0x28e180(%rip), %xmm1   # 0x37b648
leaq 0x825e1(%rip), %r14     # 0x16fab0
leal (%rax,%rax,8), %eax
cltq
movsd -0x48(%r14,%rax,8), %xmm0
mulsd 0x28e12d(%rip), %xmm0   # 0x37b610
addsd 0x28e175(%rip), %xmm0   # 0x37b660
movsd -0x40(%r14,%rax,8), %xmm3
mulsd 0x28e0fe(%rip), %xmm3   # 0x37b5f8
addsd %xmm0, %xmm3
movsd -0x38(%r14,%rax,8), %xmm0
mulsd 0x28e0fb(%rip), %xmm0   # 0x37b608
addsd %xmm3, %xmm0
movsd -0x30(%r14,%rax,8), %xmm3
mulsd 0x28e148(%rip), %xmm3   # 0x37b668
addsd %xmm0, %xmm3
movsd -0x28(%r14,%rax,8), %xmm0
mulsd 0x28e105(%rip), %xmm0   # 0x37b638
addsd %xmm3, %xmm0
mulsd 0x1e3d1(%rip), %xmm0    # 0x10b910
movsd %xmm0, 0x28e0f9(%rip)   # 0x37b640
movsd -0x20(%r14,%rax,8), %xmm3
addsd %xmm3, %xmm3
mulsd %xmm2, %xmm3
mulsd %xmm1, %xmm3
movsd %xmm3, -0x18(%rbp)
addsd 0x28e0f1(%rip), %xmm0   # 0x37b658
callq 0x9080
movslq 0x28e08d(%rip), %rax    # 0x37b600
leaq (%rax,%rax,8), %rax
movsd -0x18(%r14,%rax,8), %xmm1
addsd %xmm1, %xmm1
mulsd 0x28e0c6(%rip), %xmm1   # 0x37b650
mulsd 0x28e0b6(%rip), %xmm1   # 0x37b648
movsd %xmm1, -0x20(%rbp)
movsd %xmm0, -0x28(%rbp)
movsd 0x28e09c(%rip), %xmm0   # 0x37b640
addsd 0x28e0ac(%rip), %xmm0   # 0x37b658
callq 0x9710
mulsd -0x20(%rbp), %xmm0
movsd -0x28(%rbp), %xmm1
mulsd -0x18(%rbp), %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, 0x28e064(%rip)   # 0x37b630
movsd 0x28e074(%rip), %xmm1   # 0x37b648
movsd 0x28e074(%rip), %xmm0   # 0x37b650
movslq 0x28e01d(%rip), %rax    # 0x37b600
leaq (%rax,%rax,8), %rax
mulsd %xmm0, %xmm0
mulsd %xmm1, %xmm1
subsd %xmm0, %xmm1
movsd %xmm1, -0x18(%rbp)
movsd -0x10(%r14,%rax,8), %xmm0
mulsd %xmm1, %xmm0
movsd %xmm0, -0x28(%rbp)
movsd 0x28e030(%rip), %xmm0   # 0x37b640
addsd 0x28e040(%rip), %xmm0   # 0x37b658
callq 0x9080
movslq 0x28dfdc(%rip), %rax    # 0x37b600
leaq (%rax,%rax,8), %rax
movsd -0x18(%rbp), %xmm1
mulsd -0x8(%r14,%rax,8), %xmm1
movsd %xmm1, -0x18(%rbp)
movsd %xmm0, -0x20(%rbp)
movsd 0x28dffa(%rip), %xmm0   # 0x37b640
addsd 0x28e00a(%rip), %xmm0   # 0x37b658
callq 0x9710
mulsd -0x18(%rbp), %xmm0
movsd -0x20(%rbp), %xmm1
mulsd -0x28(%rbp), %xmm1
addsd %xmm0, %xmm1
movslq 0x28df93(%rip), %rax    # 0x37b600
leaq (%rax,%rax,8), %rax
movsd -0x10(%r14,%rax,8), %xmm0
mulsd 0x28dfd0(%rip), %xmm0   # 0x37b650
movsd %xmm0, -0x18(%rbp)
movsd %xmm1, 0x28df9b(%rip)   # 0x37b628
movsd 0x28dfab(%rip), %xmm0   # 0x37b640
addsd 0x28dfbb(%rip), %xmm0   # 0x37b658
callq 0x9710
movsd %xmm0, -0x28(%rbp)
movslq 0x28df52(%rip), %rax    # 0x37b600
leaq (%rax,%rax,8), %rax
movsd 0x28df96(%rip), %xmm1   # 0x37b650
movsd 0x28df7e(%rip), %xmm0   # 0x37b640
addsd 0x28df8e(%rip), %xmm0   # 0x37b658
mulsd -0x8(%r14,%rax,8), %xmm1
movsd %xmm1, -0x20(%rbp)
callq 0x9080
mulsd -0x20(%rbp), %xmm0
movsd -0x28(%rbp), %xmm9
mulsd -0x18(%rbp), %xmm9
subsd %xmm0, %xmm9
movsd 0x28df37(%rip), %xmm0   # 0x37b630
movsd 0x28df17(%rip), %xmm4   # 0x37b618
movapd %xmm0, %xmm5
mulsd %xmm4, %xmm5
movsd 0x28df37(%rip), %xmm1   # 0x37b648
mulsd %xmm1, %xmm5
addsd (%rbx), %xmm5
movsd 0x28deff(%rip), %xmm6   # 0x37b620
movapd %xmm9, %xmm2
mulsd %xmm6, %xmm2
subsd %xmm2, %xmm5
movsd 0x28def2(%rip), %xmm3   # 0x37b628
movsd 0x28df12(%rip), %xmm2   # 0x37b650
movapd %xmm3, %xmm7
xorpd 0x1ba66(%rip), %xmm7    # 0x1091b0
mulsd %xmm2, %xmm7
movapd %xmm7, %xmm8
mulsd %xmm4, %xmm8
addsd %xmm5, %xmm8
movsd %xmm8, (%rbx)
movapd %xmm0, %xmm5
mulsd %xmm6, %xmm5
mulsd %xmm1, %xmm5
addsd 0x8(%rbx), %xmm5
mulsd %xmm4, %xmm9
addsd %xmm5, %xmm9
mulsd %xmm6, %xmm7
addsd %xmm9, %xmm7
movsd %xmm7, 0x8(%rbx)
mulsd %xmm2, %xmm0
addsd 0x10(%rbx), %xmm0
mulsd %xmm1, %xmm3
addsd %xmm0, %xmm3
movsd %xmm3, 0x10(%rbx)
movl 0x28de59(%rip), %ecx    # 0x37b600
leal 0x1(%rcx), %eax
movl %eax, 0x28de50(%rip)    # 0x37b600
cmpl $0x1f, %ecx
jl 0xed4cf
xorl %eax, %eax
movsd 0x1b875(%rip), %xmm0    # 0x109038
xorl %ecx, %ecx
movsd (%rbx,%rcx,8), %xmm1
divsd %xmm0, %xmm1
movsd %xmm1, (%rbx,%rcx,8)
incq %rcx
cmpq $0x3, %rcx
jne 0xed7c5
addq $0x20, %rsp
popq %rbx
popq %r14
popq %rbp
retq