// File: /workspace/github_c_makefile_star2_O0/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 %rbp
movq %rsp, %rbp
subq $0x90, %rsp
movl %edi, -0x4(%rbp)
movsd %xmm0, -0x10(%rbp)
movq %rsi, -0x18(%rbp)
movq %rdx, -0x20(%rbp)
movq %rcx, -0x28(%rbp)
movq -0x18(%rbp), %rax
movsd (%rax), %xmm0
movsd %xmm0, -0x30(%rbp)
movq -0x18(%rbp), %rax
movsd (%rax), %xmm1
movsd 0x8(%rax), %xmm0
divsd %xmm1, %xmm0
movsd %xmm0, -0x38(%rbp)
movq -0x18(%rbp), %rax
movsd (%rax), %xmm1
movsd 0x10(%rax), %xmm0
divsd %xmm1, %xmm0
movsd %xmm0, -0x40(%rbp)
movq -0x18(%rbp), %rax
movsd (%rax), %xmm1
movsd 0x18(%rax), %xmm0
divsd %xmm1, %xmm0
movsd %xmm0, -0x48(%rbp)
movq -0x18(%rbp), %rax
movsd (%rax), %xmm1
movsd 0x20(%rax), %xmm0
divsd %xmm1, %xmm0
movsd -0x10(%rbp), %xmm2
movsd -0x38(%rbp), %xmm3
movsd -0x40(%rbp), %xmm1
mulsd %xmm1, %xmm1
mulsd %xmm3, %xmm3
addsd %xmm1, %xmm3
movsd -0x48(%rbp), %xmm1
mulsd %xmm1, %xmm1
addsd %xmm3, %xmm1
movsd 0x1ec49(%rip), %xmm3    # 0x27708
mulsd %xmm3, %xmm1
movsd 0x1ec45(%rip), %xmm4    # 0x27710
movaps %xmm2, %xmm5
subsd %xmm4, %xmm5
mulsd %xmm5, %xmm1
mulsd %xmm2, %xmm0
subsd %xmm1, %xmm0
movsd %xmm0, -0x50(%rbp)
movq -0x20(%rbp), %rax
movsd (%rax), %xmm0
movsd %xmm0, -0x58(%rbp)
movq -0x20(%rbp), %rax
movsd (%rax), %xmm1
movsd 0x8(%rax), %xmm0
divsd %xmm1, %xmm0
movsd %xmm0, -0x60(%rbp)
movq -0x20(%rbp), %rax
movsd (%rax), %xmm1
movsd 0x10(%rax), %xmm0
divsd %xmm1, %xmm0
movsd %xmm0, -0x68(%rbp)
movq -0x20(%rbp), %rax
movsd (%rax), %xmm1
movsd 0x18(%rax), %xmm0
divsd %xmm1, %xmm0
movsd %xmm0, -0x70(%rbp)
movq -0x20(%rbp), %rax
movsd (%rax), %xmm1
movsd 0x20(%rax), %xmm0
divsd %xmm1, %xmm0
movsd -0x10(%rbp), %xmm2
movsd -0x60(%rbp), %xmm5
movsd -0x68(%rbp), %xmm1
mulsd %xmm1, %xmm1
mulsd %xmm5, %xmm5
addsd %xmm1, %xmm5
movsd -0x70(%rbp), %xmm1
mulsd %xmm1, %xmm1
addsd %xmm5, %xmm1
mulsd %xmm3, %xmm1
movaps %xmm2, %xmm3
subsd %xmm4, %xmm3
mulsd %xmm3, %xmm1
mulsd %xmm2, %xmm0
subsd %xmm1, %xmm0
movsd %xmm0, -0x78(%rbp)
movq $0x0, -0x80(%rbp)
movl -0x4(%rbp), %eax
movl %eax, -0x84(%rbp)
testl %eax, %eax
je 0x8bb2
jmp 0x8b9e
movl -0x84(%rbp), %eax
subl $0x1, %eax
je 0x8c33
jmp 0x8d0d
movsd -0x38(%rbp), %xmm1
addsd -0x60(%rbp), %xmm1
movsd 0x1eb44(%rip), %xmm0    # 0x27708
mulsd %xmm1, %xmm0
movq -0x28(%rbp), %rax
movsd %xmm0, 0x8(%rax)
movsd -0x40(%rbp), %xmm1
addsd -0x68(%rbp), %xmm1
movsd 0x1eb25(%rip), %xmm0    # 0x27708
mulsd %xmm1, %xmm0
movq -0x28(%rbp), %rax
movsd %xmm0, 0x10(%rax)
movsd -0x48(%rbp), %xmm1
addsd -0x70(%rbp), %xmm1
movsd 0x1eb06(%rip), %xmm0    # 0x27708
mulsd %xmm1, %xmm0
movq -0x28(%rbp), %rax
movsd %xmm0, 0x18(%rax)
movsd -0x50(%rbp), %xmm1
addsd -0x78(%rbp), %xmm1
movsd 0x1eae7(%rip), %xmm0    # 0x27708
mulsd %xmm1, %xmm0
movq -0x28(%rbp), %rax
movsd %xmm0, 0x20(%rax)
jmp 0x8d0f
movsd -0x58(%rbp), %xmm0
divsd -0x30(%rbp), %xmm0
callq 0x2200
movsd %xmm0, -0x80(%rbp)
movsd -0x38(%rbp), %xmm1
movsd -0x80(%rbp), %xmm0
movsd -0x60(%rbp), %xmm2
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd 0x1eaaa(%rip), %xmm1    # 0x27710
addsd -0x80(%rbp), %xmm1
divsd %xmm1, %xmm0
movq -0x28(%rbp), %rax
movsd %xmm0, 0x8(%rax)
movsd -0x40(%rbp), %xmm1
movsd -0x80(%rbp), %xmm0
movsd -0x68(%rbp), %xmm2
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd 0x1ea79(%rip), %xmm1    # 0x27710
addsd -0x80(%rbp), %xmm1
divsd %xmm1, %xmm0
movq -0x28(%rbp), %rax
movsd %xmm0, 0x10(%rax)
movsd -0x48(%rbp), %xmm1
movsd -0x80(%rbp), %xmm0
movsd -0x70(%rbp), %xmm2
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd 0x1ea48(%rip), %xmm1    # 0x27710
addsd -0x80(%rbp), %xmm1
divsd %xmm1, %xmm0
movq -0x28(%rbp), %rax
movsd %xmm0, 0x18(%rax)
movsd -0x50(%rbp), %xmm1
movsd -0x80(%rbp), %xmm0
movsd -0x78(%rbp), %xmm2
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd 0x1ea17(%rip), %xmm1    # 0x27710
addsd -0x80(%rbp), %xmm1
divsd %xmm1, %xmm0
movq -0x28(%rbp), %rax
movsd %xmm0, 0x20(%rax)
jmp 0x8d0f
jmp 0x8d0f
movsd -0x10(%rbp), %xmm0
movsd 0x1e9f4(%rip), %xmm1    # 0x27710
subsd %xmm1, %xmm0
movq -0x28(%rbp), %rax
movsd 0x20(%rax), %xmm2
movq -0x28(%rbp), %rax
movsd 0x8(%rax), %xmm3
movq -0x28(%rbp), %rax
movsd 0x8(%rax), %xmm4
movq -0x28(%rbp), %rax
movsd 0x10(%rax), %xmm1
movq -0x28(%rbp), %rax
mulsd 0x10(%rax), %xmm1
mulsd %xmm4, %xmm3
addsd %xmm1, %xmm3
movq -0x28(%rbp), %rax
movsd 0x18(%rax), %xmm1
movq -0x28(%rbp), %rax
movsd 0x18(%rax), %xmm4
mulsd %xmm4, %xmm1
addsd %xmm3, %xmm1
movsd 0x1e9a1(%rip), %xmm3    # 0x27718
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
callq 0x2200
movq -0x28(%rbp), %rax
movsd %xmm0, 0x28(%rax)
addq $0x90, %rsp
popq %rbp
retq
nopw (%rax,%rax)