// File: /workspace/github_c_makefile_star2_O3/mohuangrui[P]ArtraCFD/solid_dynamics.c
// Function: IntegrateSurfaceForce
void IntegrateSurfaceForce(Space *space, const Model *model)
{
    const Partition *const part = &(space->part);
    const Node *const node = space->node;
    Geometry *const geo = &(space->geo);
    const IntVec nMin = {part->ns[PIN][X][MIN], part->ns[PIN][Y][MIN], part->ns[PIN][Z][MIN]};
    const IntVec nMax = {part->ns[PIN][X][MAX], part->ns[PIN][Y][MAX], part->ns[PIN][Z][MAX]};
    const RealVec sMin = {part->domain[X][MIN], part->domain[Y][MIN], part->domain[Z][MIN]};
    const RealVec d = {part->d[X], part->d[Y], part->d[Z]};
    const RealVec dd = {part->dd[X], part->dd[Y], part->dd[Z]};
    const IntVec ng = {part->ng[X], part->ng[Y], part->ng[Z]};
    const Real zero = 0.0;
    const Real percent = FLT_EPSILON * FLT_EPSILON;
    Polyhedron *poly = NULL;
    int idx = 0; /* linear array index math variable */
    int box[DIMS][LIMIT] = {{0}}; /* bounding box in node space */
    int lidN = 0; /* count total number of interfacial nodes */
    int gstN = 0; /* count total number of ghost nodes */
    RealVec pG = {zero}; /* ghost point */
    RealVec pO = {zero}; /* boundary point */
    RealVec pI = {zero}; /* image point */
    RealVec N = {zero}; /* normal */
    Real Uo[DIMUo] = {zero};
    RealVec V = {zero}; /* velocity vector */
    RealVec r = {zero}; /* position vector */
    RealVec Fp = {zero}; /* pressure force */
    RealVec Fv = {zero}; /* viscous force */
    RealVec Fs = {zero}; /* surface force */
    RealVec Tt = {zero}; /* torque */
    RealVec fvar = {zero}; /* force offset, mean, variance */
    Real Vn = zero; /* velocity projection */
    Real mu = zero; /* viscosity */
    Real ds = zero; /* infinitesimal area for integration */
    for (int n = 0; n < geo->totN; ++n) {
        poly = geo->poly + n;
        if (0 < poly->state) { /* surface force negligible */
            continue;
        }
        /* reset some non accumulative information to zero */
        memset(poly->Fp, 0, DIMS * sizeof(*poly->Fp));
        memset(poly->Fv, 0, DIMS * sizeof(*poly->Fv));
        memset(poly->Tt, 0, DIMS * sizeof(*poly->Tt));
        memset(fvar, 0, DIMS * sizeof(*fvar));
        lidN = 0;
        gstN = 0;
        /* determine search range according to bounding box of polyhedron and valid node space */
        for (int s = 0; s < DIMS; ++s) {
            box[s][MIN] = ConfineSpace(MapNode(poly->box[s][MIN], sMin[s], dd[s], ng[s]), nMin[s], nMax[s]);
            box[s][MAX] = ConfineSpace(MapNode(poly->box[s][MAX], sMin[s], dd[s], ng[s]), nMin[s], nMax[s]) + 1;
        }
        for (int k = box[Z][MIN]; k < box[Z][MAX]; ++k) {
            for (int j = box[Y][MIN]; j < box[Y][MAX]; ++j) {
                for (int i = box[X][MIN]; i < box[X][MAX]; ++i) {
                    idx = IndexNode(k, j, i, part->n[Y], part->n[X]);
                    if ((2 == node[idx].lid) && (n + 1 == node[idx].did)) {
                        ++lidN; /* an interfacial node of current geometry */
                    }
                    if ((2 != node[idx].gst) || (n + 1 != node[idx].did)) {
                        continue;
                    }
                    ++gstN; /* a ghost node of current geometry */
                    /* surface force exerted by fluid (pressure + shear force) */
                    pG[X] = MapPoint(i, sMin[X], d[X], ng[X]);
                    pG[Y] = MapPoint(j, sMin[Y], d[Y], ng[Y]);
                    pG[Z] = MapPoint(k, sMin[Z], d[Z], ng[Z]);
                    ComputeGeometricData(pG, node[idx].fid, poly, pO, pI, N);
                    r[X] = pO[X] - poly->O[X];
                    r[Y] = pO[Y] - poly->O[Y];
                    r[Z] = pO[Z] - poly->O[Z];
                    MapPrimitive(model->gamma, model->gasR, node[idx].U[TO], Uo);
                    Fp[X] = Uo[4] * N[X];
                    Fp[Y] = Uo[4] * N[Y];
                    Fp[Z] = Uo[4] * N[Z];
                    if (1 == gstN) {
                        fvar[0] = Uo[4];
                    }
                    fvar[1] = fvar[1] + Uo[4] - fvar[0];
                    fvar[2] = fvar[2] + (Uo[4] - fvar[0]) * (Uo[4] - fvar[0]);
                    if ((zero < model->refMu) && (zero < poly->cf)) {
                        mu = model->refMu * Viscosity(Uo[5] * model->refT);
                        Cross(poly->W[TO], r, V);
                        V[X] = Uo[1] - (poly->V[TO][X] + V[X]);
                        V[Y] = Uo[2] - (poly->V[TO][Y] + V[Y]);
                        V[Z] = Uo[3] - (poly->V[TO][Z] + V[Z]);
                        Vn = Dot(V, N);
                        Fv[X] = mu * (V[X] - Vn * N[X]) / Dist(pG, pO);
                        Fv[Y] = mu * (V[Y] - Vn * N[Y]) / Dist(pG, pO);
                        Fv[Z] = mu * (V[Z] - Vn * N[Z]) / Dist(pG, pO);
                    } else {
                        memset(Fv, 0, DIMS * sizeof(*Fv));
                    }
                    Fs[X] = Fp[X] + Fv[X];
                    Fs[Y] = Fp[Y] + Fv[Y];
                    Fs[Z] = Fp[Z] + Fv[Z];
                    Cross(r, Fs, Tt);
                    /* integration sum */
                    for (int s = 0; s < DIMS; ++s) {
                        poly->Fp[s] = poly->Fp[s] + Fp[s];
                        poly->Fv[s] = poly->Fv[s] + Fv[s];
                        poly->Tt[s] = poly->Tt[s] + Tt[s];
                    }
                }
            }
        }
        /* calibrate the sum of discrete forces into integration */
        if ((0 == lidN) || (0 == gstN)) { /* no surface force exerted */
            continue;
        }
        ds = poly->area / lidN;
        fvar[2] = (fvar[2] - fvar[1] * fvar[1] / gstN) / gstN; /* variance */
        fvar[1] = fvar[1] / gstN + fvar[0]; /* mean */
        if (percent * fvar[1] * fvar[1] > fvar[2]) { /* recover equilibrium state and ignore integration error */
            ds = zero;
        }
        for (int s = 0; s < DIMS; ++s) {
            poly->Fp[s] = -poly->Fp[s] * ds;
            poly->Fv[s] = -poly->Fv[s] * ds;
            poly->Tt[s] = -poly->Tt[s] * ds;
        }
    }
    return;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x338, %rsp            # imm = 0x338
movq %rsi, 0x60(%rsp)
movq (%rdi), %r13
movl 0x90(%rdi), %eax
movl %eax, 0x124(%rsp)
movl 0x98(%rdi), %eax
movl %eax, 0x128(%rsp)
movl 0xa0(%rdi), %eax
movl %eax, 0x12c(%rsp)
movl 0x94(%rdi), %eax
movl %eax, 0x118(%rsp)
movl 0x9c(%rdi), %eax
movl %eax, 0x11c(%rsp)
movl 0xa4(%rdi), %eax
movl %eax, 0x120(%rsp)
movsd 0x3f0(%rdi), %xmm0
movsd %xmm0, 0x1f8(%rsp)
movsd %xmm0, 0x2a0(%rsp)
movsd 0x400(%rdi), %xmm0
movsd %xmm0, 0x1f0(%rsp)
movsd %xmm0, 0x2a8(%rsp)
movsd 0x410(%rdi), %xmm0
movsd %xmm0, 0x1e8(%rsp)
movsd %xmm0, 0x2b0(%rsp)
movsd 0x58(%rdi), %xmm0
movsd %xmm0, 0x1e0(%rsp)
movsd 0x60(%rdi), %xmm0
movsd %xmm0, 0x1d8(%rsp)
movsd 0x68(%rdi), %xmm0
movsd %xmm0, 0x1d0(%rsp)
movups 0x70(%rdi), %xmm0
movaps %xmm0, 0x320(%rsp)
movsd 0x80(%rdi), %xmm0
movsd %xmm0, 0x330(%rsp)
movl 0x40(%rdi), %eax
movl %eax, 0x7c(%rsp)
movl %eax, 0x10c(%rsp)
movl 0x44(%rdi), %eax
movl %eax, 0x78(%rsp)
movl %eax, 0x110(%rsp)
movl 0x48(%rdi), %eax
movl %eax, 0x74(%rsp)
movl %eax, 0x114(%rsp)
xorl %eax, %eax
movq %rax, 0xb0(%rsp)
xorpd %xmm0, %xmm0
movapd %xmm0, 0xa0(%rsp)
movq %rax, 0x1a0(%rsp)
movapd %xmm0, 0x190(%rsp)
movq %rax, 0x180(%rsp)
movapd %xmm0, 0x170(%rsp)
movq %rax, 0x310(%rsp)
movapd %xmm0, 0x300(%rsp)
movq %rax, 0xe0(%rsp)
movapd %xmm0, 0xd0(%rsp)
movapd %xmm0, 0x270(%rsp)
movapd %xmm0, 0x260(%rsp)
movapd %xmm0, 0x250(%rsp)
movq %rax, 0x50(%rsp)
movapd %xmm0, 0x40(%rsp)
movq %rax, 0x160(%rsp)
movapd %xmm0, 0x150(%rsp)
movq %rax, 0x230(%rsp)
movapd %xmm0, 0x220(%rsp)
movq %rax, 0x140(%rsp)
movapd %xmm0, 0x130(%rsp)
movq %rax, 0x210(%rsp)
movapd %xmm0, 0x200(%rsp)
movq %rax, 0x290(%rsp)
movapd %xmm0, 0x280(%rsp)
cmpl $0x0, 0x8(%rdi)
jle 0x17e2f
movq %rdi, %r14
movl $0x1d8, %ecx            # imm = 0x1D8
movl $0x188, %edx            # imm = 0x188
xorl %esi, %esi
movq %r13, 0x100(%rsp)
movq 0x18(%r14), %rbp
imulq $0x238, %rsi, %rax      # imm = 0x238
cmpl $0x0, 0xc(%rbp,%rax)
jg 0x17e14
movq %rsi, 0xe8(%rsp)
movq %rdx, 0xf0(%rsp)
addq %rbp, %rax
xorpd %xmm0, %xmm0
movupd %xmm0, 0x188(%rax)
movupd %xmm0, 0x178(%rax)
movupd %xmm0, 0x168(%rax)
movupd %xmm0, 0x158(%rax)
movq %rax, 0x10(%rsp)
movq $0x0, 0x198(%rax)
movq %rcx, 0xf8(%rsp)
leaq (%rcx,%rbp), %rax
movq %rax, 0x30(%rsp)
xorl %ebx, %ebx
movq 0x30(%rsp), %rax
movsd -0x8(%rax,%rbx,4), %xmm0
movsd 0x2a0(%rsp,%rbx,2), %xmm1
movsd %xmm1, 0x8(%rsp)
movsd 0x320(%rsp,%rbx,2), %xmm2
movsd %xmm2, 0x28(%rsp)
movl 0x10c(%rsp,%rbx), %edi
movl %edi, 0x20(%rsp)
callq 0x7d83
movl 0x124(%rsp,%rbx), %r12d
movl 0x118(%rsp,%rbx), %r15d
movl %eax, %edi
movl %r12d, %esi
movl %r15d, %edx
callq 0x7d9a
movl %eax, 0xa0(%rsp,%rbx,2)
movq 0x30(%rsp), %rax
movsd (%rax,%rbx,4), %xmm0
movsd 0x8(%rsp), %xmm1
movsd 0x28(%rsp), %xmm2
movl 0x20(%rsp), %edi
callq 0x7d83
movl %eax, %edi
movl %r12d, %esi
movl %r15d, %edx
callq 0x7d9a
incl %eax
movl %eax, 0xa4(%rsp,%rbx,2)
addq $0x4, %rbx
cmpq $0xc, %rbx
jne 0x1765f
movl 0xb0(%rsp), %eax
movl 0xb4(%rsp), %ecx
movl %eax, 0x8(%rsp)
movl %ecx, 0x6c(%rsp)
cmpl %ecx, %eax
movq 0xf8(%rsp), %rcx
movq 0xf0(%rsp), %rdx
movq 0xe8(%rsp), %rsi
movq 0x10(%rsp), %rax
jge 0x17e14
movl 0xa8(%rsp), %ecx
movl %ecx, 0x68(%rsp)
movl 0xac(%rsp), %ecx
movl %ecx, 0x18(%rsp)
movl 0xa0(%rsp), %ecx
movl %ecx, 0x70(%rsp)
movl 0xa4(%rsp), %ecx
movl %ecx, 0x28(%rsp)
leaq 0x1(%rsi), %rcx
movq %rcx, 0x20(%rsp)
addq $0xb0, %rax
movq %rax, 0x1c0(%rsp)
addq %rdx, %rbp
movl $0x0, 0x4(%rsp)
xorpd %xmm0, %xmm0
xorpd %xmm1, %xmm1
movapd %xmm1, 0x1b0(%rsp)
movapd %xmm1, 0x90(%rsp)
movl $0x0, 0x1c(%rsp)
movl 0x68(%rsp), %eax
movl %eax, %r15d
cmpl 0x18(%rsp), %eax
jge 0x17d11
movl 0x70(%rsp), %eax
movl %eax, %r12d
cmpl 0x28(%rsp), %eax
jge 0x17d03
movapd %xmm0, 0x30(%rsp)
movl 0x34(%r14), %r8d
movl 0x38(%r14), %ecx
movl 0x8(%rsp), %edi
movl %r15d, %esi
movl %r12d, %edx
callq 0x7d53
cltq
movq %rax, %rcx
shlq $0x7, %rcx
leaq (%rcx,%rax,8), %rax
leaq (%rax,%r13), %rbx
cmpl $0x2, 0x8(%r13,%rax)
jne 0x17812
movl (%rbx), %eax
xorl %ecx, %ecx
cmpq %rax, 0x20(%rsp)
sete %cl
addl %ecx, 0x1c(%rsp)
cmpl $0x2, 0xc(%rbx)
movapd 0x30(%rsp), %xmm0
jne 0x17cf5
movl (%rbx), %eax
cmpq %rax, 0x20(%rsp)
jne 0x17cf5
movl %r12d, %edi
movsd 0x1f8(%rsp), %xmm0
movsd 0x1e0(%rsp), %xmm1
movl 0x7c(%rsp), %esi
callq 0x7db8
movsd %xmm0, 0x190(%rsp)
movl %r15d, %edi
movsd 0x1f0(%rsp), %xmm0
movsd 0x1d8(%rsp), %xmm1
movl 0x78(%rsp), %esi
callq 0x7db8
movsd %xmm0, 0x198(%rsp)
movl 0x8(%rsp), %edi
movsd 0x1e8(%rsp), %xmm0
movsd 0x1d0(%rsp), %xmm1
movl 0x74(%rsp), %esi
callq 0x7db8
movsd %xmm0, 0x1a0(%rsp)
movl 0x4(%rbx), %esi
leaq 0x190(%rsp), %rdi
movq 0x10(%rsp), %r13
movq %r13, %rdx
leaq 0x170(%rsp), %rcx
leaq 0x300(%rsp), %r8
leaq 0xd0(%rsp), %r9
callq 0xb5ce
movsd 0x170(%rsp), %xmm0
subsd 0x20(%r13), %xmm0
movsd 0x178(%rsp), %xmm1
movsd %xmm0, 0x150(%rsp)
subsd 0x28(%r13), %xmm1
movsd %xmm1, 0x158(%rsp)
movsd 0x180(%rsp), %xmm0
subsd 0x30(%r13), %xmm0
movsd %xmm0, 0x160(%rsp)
movq 0x60(%rsp), %rax
movsd 0x40(%rax), %xmm0
movsd 0x48(%rax), %xmm1
addq $0x10, %rbx
movq %rbx, %rdi
leaq 0x250(%rsp), %rsi
callq 0x7c16
movsd 0x270(%rsp), %xmm3
movapd %xmm3, %xmm2
unpcklpd %xmm3, %xmm2            # xmm2 = xmm2[0],xmm3[0]
mulpd 0xd0(%rsp), %xmm2
movapd %xmm2, 0x220(%rsp)
movsd 0xe0(%rsp), %xmm4
mulsd %xmm3, %xmm4
movsd %xmm4, 0x230(%rsp)
movapd %xmm3, %xmm1
cmpl $0x0, 0x4(%rsp)
je 0x1798d
movapd 0x90(%rsp), %xmm1
movq 0x60(%rsp), %rax
movsd 0x38(%rax), %xmm5
ucomisd 0x1671(%rip), %xmm5     # 0x19010
movsd %xmm1, 0x90(%rsp)
jbe 0x17bfe
movq 0x10(%rsp), %rax
movsd 0x1b8(%rax), %xmm0
ucomisd 0x164d(%rip), %xmm0     # 0x19010
jbe 0x17bfe
movsd 0x278(%rsp), %xmm0
movq 0x60(%rsp), %rax
mulsd 0x70(%rax), %xmm0
movapd %xmm2, 0x2f0(%rsp)
movapd %xmm3, 0x2e0(%rsp)
movsd %xmm4, 0x1c8(%rsp)
movapd %xmm5, 0x80(%rsp)
callq 0x7bdf
movapd 0x80(%rsp), %xmm1
mulsd %xmm0, %xmm1
movapd %xmm1, 0x80(%rsp)
movq 0x1c0(%rsp), %rdi
leaq 0x150(%rsp), %rsi
leaq 0x40(%rsp), %rbx
movq %rbx, %rdx
callq 0x7eee
movq 0x10(%rsp), %rax
movsd 0x80(%rax), %xmm0
addsd 0x40(%rsp), %xmm0
movsd 0x258(%rsp), %xmm1
subsd %xmm0, %xmm1
movsd 0x260(%rsp), %xmm0
movsd %xmm1, 0x40(%rsp)
movsd 0x88(%rax), %xmm1
addsd 0x48(%rsp), %xmm1
subsd %xmm1, %xmm0
movsd %xmm0, 0x48(%rsp)
movsd 0x268(%rsp), %xmm0
movsd 0x90(%rax), %xmm1
addsd 0x50(%rsp), %xmm1
subsd %xmm1, %xmm0
movsd %xmm0, 0x50(%rsp)
movq %rbx, %rdi
leaq 0xd0(%rsp), %rsi
callq 0x7e36
movapd %xmm0, 0x240(%rsp)
movaps 0x40(%rsp), %xmm0
movaps %xmm0, 0xc0(%rsp)
movaps 0xd0(%rsp), %xmm0
movaps %xmm0, 0x2c0(%rsp)
leaq 0x190(%rsp), %rbx
movq %rbx, %rdi
leaq 0x170(%rsp), %r13
movq %r13, %rsi
callq 0x7eb6
movaps %xmm0, 0x2d0(%rsp)
movapd 0x240(%rsp), %xmm0
movapd 0x2c0(%rsp), %xmm1
movhpd 0xd8(%rsp), %xmm1       # xmm1 = xmm1[0],mem[0]
unpcklpd %xmm0, %xmm0            # xmm0 = xmm0[0,0]
mulpd %xmm0, %xmm1
movapd 0xc0(%rsp), %xmm0
movhpd 0x48(%rsp), %xmm0       # xmm0 = xmm0[0],mem[0]
subpd %xmm1, %xmm0
movapd 0x80(%rsp), %xmm1
unpcklpd %xmm1, %xmm1            # xmm1 = xmm1[0,0]
mulpd %xmm0, %xmm1
movapd %xmm1, 0xc0(%rsp)
movq %rbx, %rdi
movq %r13, %rsi
callq 0x7eb6
movapd 0x2d0(%rsp), %xmm1
unpcklpd %xmm0, %xmm1            # xmm1 = xmm1[0],xmm0[0]
movapd 0xc0(%rsp), %xmm0
divpd %xmm1, %xmm0
movapd %xmm0, 0xc0(%rsp)
movapd %xmm0, 0x130(%rsp)
movsd 0x50(%rsp), %xmm1
movapd 0x240(%rsp), %xmm0
mulsd 0xe0(%rsp), %xmm0
subsd %xmm0, %xmm1
mulsd 0x80(%rsp), %xmm1
movsd %xmm1, 0x80(%rsp)
movq %rbx, %rdi
movq %r13, %rsi
callq 0x7eb6
movsd 0x80(%rsp), %xmm6
movapd 0xc0(%rsp), %xmm5
movsd 0x1c8(%rsp), %xmm4
movapd 0x2e0(%rsp), %xmm3
movapd 0x2f0(%rsp), %xmm2
movsd 0x90(%rsp), %xmm1
divsd %xmm0, %xmm6
movsd %xmm6, 0x140(%rsp)
jmp 0x17c1f
xorpd %xmm0, %xmm0
movapd %xmm0, 0x130(%rsp)
movq $0x0, 0x140(%rsp)
xorpd %xmm5, %xmm5
xorpd %xmm6, %xmm6
movq 0x100(%rsp), %r13
incl 0x4(%rsp)
movapd 0x1b0(%rsp), %xmm0
addsd %xmm3, %xmm0
subsd %xmm1, %xmm0
movapd %xmm0, 0x1b0(%rsp)
subsd %xmm1, %xmm3
mulsd %xmm3, %xmm3
movapd 0x30(%rsp), %xmm0
addsd %xmm3, %xmm0
movapd %xmm0, 0x30(%rsp)
addpd %xmm5, %xmm2
movapd %xmm2, 0x200(%rsp)
addsd %xmm6, %xmm4
movsd %xmm4, 0x210(%rsp)
leaq 0x150(%rsp), %rdi
leaq 0x200(%rsp), %rsi
leaq 0x280(%rsp), %rdx
callq 0x7eee
xorl %eax, %eax
movsd -0x30(%rbp,%rax,8), %xmm0
addsd 0x220(%rsp,%rax,8), %xmm0
movsd -0x18(%rbp,%rax,8), %xmm1
movsd %xmm0, -0x30(%rbp,%rax,8)
addsd 0x130(%rsp,%rax,8), %xmm1
movsd %xmm1, -0x18(%rbp,%rax,8)
movsd (%rbp,%rax,8), %xmm0
addsd 0x280(%rsp,%rax,8), %xmm0
movsd %xmm0, (%rbp,%rax,8)
incq %rax
cmpq $0x3, %rax
jne 0x17c96
movsd 0x90(%rsp), %xmm0
movaps %xmm0, 0x90(%rsp)
movapd 0x30(%rsp), %xmm0
incl %r12d
cmpl 0x28(%rsp), %r12d
jne 0x177cc
incl %r15d
cmpl 0x18(%rsp), %r15d
jne 0x177bb
movl 0x8(%rsp), %eax
incl %eax
movl %eax, 0x8(%rsp)
cmpl 0x6c(%rsp), %eax
jne 0x177aa
movl 0x1c(%rsp), %r8d
testl %r8d, %r8d
sete %al
movl 0x4(%rsp), %edi
testl %edi, %edi
sete %cl
orb %al, %cl
movq 0xf8(%rsp), %rcx
movq 0xf0(%rsp), %rdx
movq 0xe8(%rsp), %rsi
movq 0x10(%rsp), %rax
jne 0x17e14
movapd %xmm0, %xmm4
movsd 0x1c0(%rax), %xmm0
xorps %xmm1, %xmm1
cvtsi2sd %r8d, %xmm1
divsd %xmm1, %xmm0
movapd 0x1b0(%rsp), %xmm3
movapd %xmm3, %xmm1
xorps %xmm2, %xmm2
cvtsi2sd %edi, %xmm2
mulsd %xmm3, %xmm1
divsd %xmm2, %xmm1
subsd %xmm1, %xmm4
divsd %xmm2, %xmm4
divsd %xmm2, %xmm3
movapd 0x90(%rsp), %xmm2
addsd %xmm3, %xmm2
movapd %xmm2, %xmm1
mulsd 0x81e1(%rip), %xmm1     # 0x1ff98
mulsd %xmm2, %xmm1
cmpltsd %xmm1, %xmm4
andnpd %xmm0, %xmm4
movq $-0x3, %rax
movsd -0x18(%rbp,%rax,8), %xmm0
movsd (%rbp,%rax,8), %xmm1
movsd 0x18(%rbp,%rax,8), %xmm2
movapd 0x483b(%rip), %xmm3     # 0x1c620
xorpd %xmm3, %xmm0
mulsd %xmm4, %xmm0
movsd %xmm0, -0x18(%rbp,%rax,8)
xorpd %xmm3, %xmm1
mulsd %xmm4, %xmm1
movsd %xmm1, (%rbp,%rax,8)
xorpd %xmm3, %xmm2
mulsd %xmm4, %xmm2
movsd %xmm2, 0x18(%rbp,%rax,8)
incq %rax
jne 0x17dcb
incq %rsi
movslq 0x8(%r14), %rax
movl $0x238, %edi            # imm = 0x238
addq %rdi, %rcx
addq %rdi, %rdx
cmpq %rax, %rsi
jl 0x175ef
addq $0x338, %rsp            # imm = 0x338
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
nopl (%rax)