// File: /workspace/github_c_makefile_star2_O1/mohuangrui[P]ArtraCFD/cfd_commons.c
// Function: SymmetricAverage
void SymmetricAverage(const int averager, const Real gamma,
        const Real UL[restrict], const Real UR[restrict], Real Uo[restrict])
{
    const Real rhoL = UL[0];
    const Real uL = UL[1] / UL[0];
    const Real vL = UL[2] / UL[0];
    const Real wL = UL[3] / UL[0];
    const Real hTL = (UL[4] / UL[0]) * gamma - 0.5 * (uL * uL + vL * vL + wL * wL) * (gamma - 1.0);
    const Real rhoR = UR[0];
    const Real uR = UR[1] / UR[0];
    const Real vR = UR[2] / UR[0];
    const Real wR = UR[3] / UR[0];
    const Real hTR = (UR[4] / UR[0]) * gamma - 0.5 * (uR * uR + vR * vR + wR * wR) * (gamma - 1.0);
    Real D = 0.0;
    switch (averager) {
        case 0: /* arithmetic mean */
            Uo[1] = 0.5 * (uL + uR); /* u average */
            Uo[2] = 0.5 * (vL + vR); /* v average */
            Uo[3] = 0.5 * (wL + wR); /* w average */
            Uo[4] = 0.5 * (hTL + hTR); /* hT average */
            break;
        case 1: /* Roe average */
            D = sqrt(rhoR / rhoL);
            Uo[1] = (uL + D * uR) / (1.0 + D); /* u average */
            Uo[2] = (vL + D * vR) / (1.0 + D); /* v average */
            Uo[3] = (wL + D * wR) / (1.0 + D); /* w average */
            Uo[4] = (hTL + D * hTR) / (1.0 + D); /* hT average */
            break;
        default:
            break;
    }
    Uo[5] = sqrt((gamma - 1.0) * (Uo[4] - 0.5 * (Uo[1] * Uo[1] + Uo[2] * Uo[2] + Uo[3] * Uo[3]))); /* the speed of sound */
    return;
}

// Assembly:
pushq %rbx
subq $0x50, %rsp
movq %rcx, %rbx
movsd (%rsi), %xmm2
movsd 0x8(%rsi), %xmm10
divsd %xmm2, %xmm10
movsd 0x10(%rsi), %xmm9
divsd %xmm2, %xmm9
movsd 0x18(%rsi), %xmm8
divsd %xmm2, %xmm8
movsd 0x20(%rsi), %xmm7
divsd %xmm2, %xmm7
movapd %xmm9, %xmm1
mulsd %xmm9, %xmm1
movapd %xmm10, %xmm3
mulsd %xmm10, %xmm3
addsd %xmm1, %xmm3
movapd %xmm8, %xmm1
mulsd %xmm8, %xmm1
addsd %xmm3, %xmm1
movsd 0x15c90(%rip), %xmm6    # 0x1d5b0
addsd %xmm0, %xmm6
movsd 0x15c8c(%rip), %xmm3    # 0x1d5b8
mulsd %xmm3, %xmm1
mulsd %xmm6, %xmm1
mulsd %xmm0, %xmm7
addsd %xmm1, %xmm7
movsd (%rdx), %xmm1
movsd 0x8(%rdx), %xmm14
divsd %xmm1, %xmm14
movsd 0x10(%rdx), %xmm13
divsd %xmm1, %xmm13
movsd 0x18(%rdx), %xmm12
divsd %xmm1, %xmm12
movsd 0x20(%rdx), %xmm11
divsd %xmm1, %xmm11
movapd %xmm13, %xmm4
mulsd %xmm13, %xmm4
movapd %xmm14, %xmm5
mulsd %xmm14, %xmm5
addsd %xmm4, %xmm5
movapd %xmm12, %xmm4
mulsd %xmm12, %xmm4
addsd %xmm5, %xmm4
mulsd %xmm3, %xmm4
mulsd %xmm6, %xmm4
mulsd %xmm0, %xmm11
addsd %xmm4, %xmm11
cmpl $0x1, %edi
je 0x79fc
testl %edi, %edi
jne 0x7af9
addsd %xmm14, %xmm10
movsd 0x15c02(%rip), %xmm0    # 0x1d5c0
mulsd %xmm0, %xmm10
movsd %xmm10, 0x8(%rbx)
addsd %xmm13, %xmm9
mulsd %xmm0, %xmm9
movsd %xmm9, 0x10(%rbx)
addsd %xmm12, %xmm8
mulsd %xmm0, %xmm8
movsd %xmm8, 0x18(%rbx)
addsd %xmm11, %xmm7
mulsd %xmm0, %xmm7
movapd %xmm7, %xmm11
jmp 0x7af3
divsd %xmm2, %xmm1
xorpd %xmm0, %xmm0
ucomisd %xmm0, %xmm1
jb 0x7a16
xorps %xmm0, %xmm0
sqrtsd %xmm1, %xmm0
jmp 0x7a99
movapd %xmm1, %xmm0
movsd %xmm6, 0x48(%rsp)
movsd %xmm7, 0x40(%rsp)
movsd %xmm8, 0x38(%rsp)
movsd %xmm9, 0x30(%rsp)
movsd %xmm10, 0x28(%rsp)
movsd %xmm11, 0x20(%rsp)
movsd %xmm12, 0x18(%rsp)
movsd %xmm13, 0x10(%rsp)
movsd %xmm14, 0x8(%rsp)
callq 0x2210
movsd 0x8(%rsp), %xmm14
movsd 0x10(%rsp), %xmm13
movsd 0x18(%rsp), %xmm12
movsd 0x20(%rsp), %xmm11
movsd 0x28(%rsp), %xmm10
movsd 0x30(%rsp), %xmm9
movsd 0x38(%rsp), %xmm8
movsd 0x40(%rsp), %xmm7
movsd 0x48(%rsp), %xmm6
mulsd %xmm0, %xmm14
addsd %xmm10, %xmm14
movsd 0x1255d(%rip), %xmm1    # 0x1a008
addsd %xmm0, %xmm1
divsd %xmm1, %xmm14
movsd %xmm14, 0x8(%rbx)
mulsd %xmm0, %xmm13
addsd %xmm9, %xmm13
divsd %xmm1, %xmm13
movsd %xmm13, 0x10(%rbx)
mulsd %xmm0, %xmm12
addsd %xmm8, %xmm12
divsd %xmm1, %xmm12
movsd %xmm12, 0x18(%rbx)
mulsd %xmm0, %xmm11
addsd %xmm7, %xmm11
divsd %xmm1, %xmm11
movsd %xmm11, 0x20(%rbx)
movsd 0x8(%rbx), %xmm0
movsd 0x10(%rbx), %xmm1
mulsd %xmm1, %xmm1
mulsd %xmm0, %xmm0
addsd %xmm1, %xmm0
movsd 0x18(%rbx), %xmm1
mulsd %xmm1, %xmm1
addsd %xmm0, %xmm1
mulsd 0x15a94(%rip), %xmm1    # 0x1d5b8
addsd 0x20(%rbx), %xmm1
mulsd %xmm1, %xmm6
xorpd %xmm0, %xmm0
ucomisd %xmm0, %xmm6
jb 0x7b40
xorps %xmm0, %xmm0
sqrtsd %xmm6, %xmm0
jmp 0x7b49
movapd %xmm6, %xmm0
callq 0x2210
movsd %xmm0, 0x28(%rbx)
addq $0x50, %rsp
popq %rbx
retq