void Roe_HLL_solver(double *V_mk, double *F, double gamma, double P_L,
                    double RHO_L, double U_L, double V_L, double P_R,
                    double RHO_R, double U_R, double V_R, double *lambda_max,
                    double delta) {
  double C_L, C_R;
  C_L = sqrt(gamma * P_L / RHO_L);
  C_R = sqrt(gamma * P_R / RHO_R);
  double z = 0.5 * (gamma - 1.0) / gamma;
  double H_L, H_R;
  H_L = gamma / (gamma - 1.0) * P_L / RHO_L + 0.5 * (U_L * U_L);
  H_R = gamma / (gamma - 1.0) * P_R / RHO_R + 0.5 * (U_R * U_R);
  double E_L, E_R;
  E_L = 1.0 / (gamma - 1.0) * P_L / RHO_L + 0.5 * (U_L * U_L);
  E_R = 1.0 / (gamma - 1.0) * P_R / RHO_R + 0.5 * (U_R * U_R);
  double U[3];
  U[0] = RHO_L;
  U[1] = RHO_L * U_L;
  U[2] = RHO_L * E_L;
  double RHO_S, U_S, H_S, C_S;
  RHO_S = sqrt(RHO_L * RHO_R);
  U_S = (U_L * sqrt(RHO_L) + U_R * sqrt(RHO_R)) / (sqrt(RHO_L) + sqrt(RHO_R));
  H_S = (H_L * sqrt(RHO_L) + H_R * sqrt(RHO_R)) / (sqrt(RHO_L) + sqrt(RHO_R));
  C_S = sqrt((gamma - 1.0) * (H_S - 0.5 * U_S * U_S));
  double R[3][3];
  double lambda[3], W[3];
  R[0][0] = 1.0;
  R[0][1] = 1.0;
  R[0][2] = 1.0;
  R[1][0] = U_S - C_S;
  R[1][1] = U_S;
  R[1][2] = U_S + C_S;
  R[2][0] = H_S - U_S * C_S;
  R[2][1] = 0.5 * (U_S * U_S);
  R[2][2] = H_S + U_S * C_S;
  int i, j;
  W[0] = 0.5 * ((P_R - P_L) - RHO_S * C_S * (U_R - U_L)) / (C_S * C_S);
  W[1] = (RHO_R - RHO_L) - (P_R - P_L) / (C_S * C_S);
  W[2] = 0.5 * ((P_R - P_L) + RHO_S * C_S * (U_R - U_L)) / (C_S * C_S);
  lambda[0] = U_S - C_S;
  lambda[1] = U_S;
  lambda[2] = U_S + C_S;
  *lambda_max = fabs(U_S) + C_S;
  for (j = 0; j < 3; j++) {
    if (lambda[j] <= 0) {
      for (i = 0; i < 3; i++) {
        U[i] += W[j] * R[i][j];
      }
    }
  }
  F[0] = U[0];
  F[1] = U[1] / U[0];
  F[2] = (U[2] / U[0] - 0.5 * F[1] * F[1]) * (gamma - 1.0) * F[0];
  double S_L, S_R;
  S_L = (U_S - C_S < U_L - C_L ? U_S - C_S : U_L - C_L);
  S_R = (U_S + C_S > U_R + C_R ? U_S + C_S : U_R + C_R);
  S_L = (0 < S_L ? 0 : S_L);
  S_R = (0 > S_R ? 0 : S_R);
  *V_mk = (RHO_R * (S_R - U_R) * V_R - RHO_L * (S_L - U_L) * V_L) /
          (RHO_R * (S_R - U_R) - RHO_L * (S_L - U_L));
}


x86 assembly (clang, -O0):
<Roe_HLL_solver>:
push   %rbp
mov    %rsp,%rbp
sub    $0x1f0,%rsp
movsd  0x18(%rbp),%xmm8
movsd  0x10(%rbp),%xmm8
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
movsd  %xmm2,-0x28(%rbp)
movsd  %xmm3,-0x30(%rbp)
movsd  %xmm4,-0x38(%rbp)
movsd  %xmm5,-0x40(%rbp)
movsd  %xmm6,-0x48(%rbp)
movsd  %xmm7,-0x50(%rbp)
mov    %rdx,-0x58(%rbp)
movsd  -0x18(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
divsd  -0x28(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x60(%rbp)
movsd  -0x18(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
divsd  -0x48(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x68(%rbp)
movsd  -0x18(%rbp),%xmm1
movsd  0xe7e(%rip),%xmm0
subsd  %xmm0,%xmm1
movsd  0xe6a(%rip),%xmm0
mulsd  %xmm1,%xmm0
divsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x18(%rbp),%xmm1
movsd  -0x18(%rbp),%xmm0
movsd  0xe52(%rip),%xmm2
subsd  %xmm2,%xmm0
divsd  %xmm0,%xmm1
mulsd  -0x20(%rbp),%xmm1
divsd  -0x28(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
movsd  0xe26(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x18(%rbp),%xmm1
movsd  -0x18(%rbp),%xmm0
movsd  0xe0f(%rip),%xmm2
subsd  %xmm2,%xmm0
divsd  %xmm0,%xmm1
mulsd  -0x40(%rbp),%xmm1
divsd  -0x48(%rbp),%xmm1
movsd  -0x50(%rbp),%xmm0
mulsd  -0x50(%rbp),%xmm0
movsd  0xde3(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  0xdd1(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  0xdc5(%rip),%xmm1
divsd  %xmm0,%xmm1
mulsd  -0x20(%rbp),%xmm1
divsd  -0x28(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
movsd  0xd9d(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
movsd  0xd88(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  0xd7c(%rip),%xmm1
divsd  %xmm0,%xmm1
mulsd  -0x40(%rbp),%xmm1
divsd  -0x48(%rbp),%xmm1
movsd  -0x50(%rbp),%xmm0
mulsd  -0x50(%rbp),%xmm0
movsd  0xd54(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
mulsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x1b8(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x1c0(%rbp)
ff
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x1c8(%rbp)
ff
movsd  -0x48(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  -0x1c8(%rbp),%xmm1
ff
movsd  -0x1c0(%rbp),%xmm2
ff
movaps %xmm0,%xmm3
movsd  -0x1b8(%rbp),%xmm0
ff
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x1a8(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x1b0(%rbp)
ff
movsd  -0x48(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  -0x1b0(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x1a8(%rbp),%xmm0
ff
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x190(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x198(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x1a0(%rbp)
ff
movsd  -0x48(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  -0x1a0(%rbp),%xmm1
ff
movsd  -0x198(%rbp),%xmm2
movaps %xmm0,%xmm3
movsd  -0x190(%rbp),%xmm0
ff
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x180(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x188(%rbp)
ff
movsd  -0x48(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  -0x188(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x180(%rbp),%xmm0
ff
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
movsd  0xbac(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  -0xc8(%rbp),%xmm2
ff
movsd  0xb90(%rip),%xmm1
mulsd  -0xc0(%rbp),%xmm1
ff
movsd  -0xc0(%rbp),%xmm3
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0xd0(%rbp)
ff
movsd  0xb50(%rip),%xmm0
movsd  %xmm0,-0x120(%rbp)
ff
movsd  0xb40(%rip),%xmm0
movsd  %xmm0,-0x118(%rbp)
ff
movsd  0xb30(%rip),%xmm0
movsd  %xmm0,-0x110(%rbp)
ff
movsd  -0xc0(%rbp),%xmm0
ff
subsd  -0xd0(%rbp),%xmm0
ff
movsd  %xmm0,-0x108(%rbp)
ff
movsd  -0xc0(%rbp),%xmm0
ff
movsd  %xmm0,-0x100(%rbp)
ff
movsd  -0xc0(%rbp),%xmm0
ff
addsd  -0xd0(%rbp),%xmm0
ff
movsd  %xmm0,-0xf8(%rbp)
ff
movsd  -0xc8(%rbp),%xmm1
ff
movsd  -0xc0(%rbp),%xmm0
ff
movsd  -0xd0(%rbp),%xmm2
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xf0(%rbp)
ff
movsd  -0xc0(%rbp),%xmm1
ff
mulsd  -0xc0(%rbp),%xmm1
ff
movsd  0xa89(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xe8(%rbp)
ff
movsd  -0xc8(%rbp),%xmm1
ff
movsd  -0xc0(%rbp),%xmm0
ff
movsd  -0xd0(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
movsd  -0x40(%rbp),%xmm0
subsd  -0x20(%rbp),%xmm0
movsd  -0xb8(%rbp),%xmm1
ff
mulsd  -0xd0(%rbp),%xmm1
ff
movsd  -0x50(%rbp),%xmm2
subsd  -0x30(%rbp),%xmm2
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  0xa0a(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  -0xd0(%rbp),%xmm1
ff
mulsd  -0xd0(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x160(%rbp)
ff
movsd  -0x48(%rbp),%xmm0
subsd  -0x28(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm1
subsd  -0x20(%rbp),%xmm1
movsd  -0xd0(%rbp),%xmm2
ff
mulsd  -0xd0(%rbp),%xmm2
ff
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x158(%rbp)
ff
movsd  -0x40(%rbp),%xmm0
subsd  -0x20(%rbp),%xmm0
movsd  -0xb8(%rbp),%xmm1
ff
mulsd  -0xd0(%rbp),%xmm1
ff
movsd  -0x50(%rbp),%xmm2
subsd  -0x30(%rbp),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  0x982(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  -0xd0(%rbp),%xmm1
ff
mulsd  -0xd0(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x150(%rbp)
ff
movsd  -0xc0(%rbp),%xmm0
ff
subsd  -0xd0(%rbp),%xmm0
ff
movsd  %xmm0,-0x140(%rbp)
ff
movsd  -0xc0(%rbp),%xmm0
ff
movsd  %xmm0,-0x138(%rbp)
ff
movsd  -0xc0(%rbp),%xmm0
ff
addsd  -0xd0(%rbp),%xmm0
ff
movsd  %xmm0,-0x130(%rbp)
ff
movsd  -0xc0(%rbp),%xmm0
ff
movaps 0x903(%rip),%xmm1
pand   %xmm1,%xmm0
addsd  -0xd0(%rbp),%xmm0
ff
mov    -0x58(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x0,-0x168(%rbp)
cmpl   $0x3,-0x168(%rbp)
jge    17dc <Roe_HLL_solver+0x6cc>
movslq -0x168(%rbp),%rax
xorps  %xmm0,%xmm0
ucomisd -0x140(%rbp,%rax,8),%xmm0
ff ff
jb     17c3 <Roe_HLL_solver+0x6b3>
movl   $0x0,-0x164(%rbp)
cmpl   $0x3,-0x164(%rbp)
jge    17be <Roe_HLL_solver+0x6ae>
movslq -0x168(%rbp),%rax
movsd  -0x160(%rbp,%rax,8),%xmm0
ff ff
movslq -0x164(%rbp),%rcx
lea    -0x120(%rbp),%rax
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x168(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movslq -0x164(%rbp),%rax
movsd  -0xb0(%rbp,%rax,8),%xmm1
ff ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xb0(%rbp,%rax,8)
ff ff
mov    -0x164(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x164(%rbp)
jmp    174b <Roe_HLL_solver+0x63b>
jmp    17c3 <Roe_HLL_solver+0x6b3>
jmp    17c8 <Roe_HLL_solver+0x6b8>
mov    -0x168(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x168(%rbp)
jmp    171b <Roe_HLL_solver+0x60b>
movsd  -0xb0(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0xa8(%rbp),%xmm0
ff
divsd  -0xb0(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x8(%rax)
movsd  -0xa0(%rbp),%xmm1
ff
divsd  -0xb0(%rbp),%xmm1
ff
mov    -0x10(%rbp),%rax
movsd  0x7ef(%rip),%xmm0
mulsd  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0x18(%rbp),%xmm1
movsd  0x7bd(%rip),%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mulsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x10(%rax)
movsd  -0xc0(%rbp),%xmm1
ff
subsd  -0xd0(%rbp),%xmm1
ff
movsd  -0x30(%rbp),%xmm0
subsd  -0x60(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    18b5 <Roe_HLL_solver+0x7a5>
movsd  -0xc0(%rbp),%xmm0
ff
subsd  -0xd0(%rbp),%xmm0
ff
movsd  %xmm0,-0x1d0(%rbp)
ff
jmp    18c7 <Roe_HLL_solver+0x7b7>
movsd  -0x30(%rbp),%xmm0
subsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x1d0(%rbp)
ff
movsd  -0x1d0(%rbp),%xmm0
ff
movsd  %xmm0,-0x170(%rbp)
ff
movsd  -0xc0(%rbp),%xmm0
ff
addsd  -0xd0(%rbp),%xmm0
ff
movsd  -0x50(%rbp),%xmm1
addsd  -0x68(%rbp),%xmm1
ucomisd %xmm1,%xmm0
jbe    1918 <Roe_HLL_solver+0x808>
movsd  -0xc0(%rbp),%xmm0
ff
addsd  -0xd0(%rbp),%xmm0
ff
movsd  %xmm0,-0x1d8(%rbp)
ff
jmp    192a <Roe_HLL_solver+0x81a>
movsd  -0x50(%rbp),%xmm0
addsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x1d8(%rbp)
ff
movsd  -0x1d8(%rbp),%xmm0
ff
movsd  %xmm0,-0x178(%rbp)
ff
movsd  -0x170(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    195f <Roe_HLL_solver+0x84f>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x1e0(%rbp)
ff
jmp    196f <Roe_HLL_solver+0x85f>
movsd  -0x170(%rbp),%xmm0
ff
movsd  %xmm0,-0x1e0(%rbp)
ff
movsd  -0x1e0(%rbp),%xmm0
ff
movsd  %xmm0,-0x170(%rbp)
ff
xorps  %xmm0,%xmm0
ucomisd -0x178(%rbp),%xmm0
ff
jbe    19a0 <Roe_HLL_solver+0x890>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x1e8(%rbp)
ff
jmp    19b0 <Roe_HLL_solver+0x8a0>
movsd  -0x178(%rbp),%xmm0
ff
movsd  %xmm0,-0x1e8(%rbp)
ff
movsd  -0x1e8(%rbp),%xmm0
ff
movsd  %xmm0,-0x178(%rbp)
ff
movsd  -0x48(%rbp),%xmm0
movsd  -0x178(%rbp),%xmm1
ff
subsd  -0x50(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0x10(%rbp),%xmm2
movsd  -0x28(%rbp),%xmm1
movsd  -0x170(%rbp),%xmm3
ff
subsd  -0x30(%rbp),%xmm3
mulsd  %xmm3,%xmm1
mulsd  -0x38(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0x48(%rbp),%xmm1
movsd  -0x178(%rbp),%xmm3
ff
subsd  -0x50(%rbp),%xmm3
movsd  -0x28(%rbp),%xmm2
movsd  -0x170(%rbp),%xmm4
ff
subsd  -0x30(%rbp),%xmm4
mulsd  %xmm4,%xmm2
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,(%rax)
add    $0x1f0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<Roe_HLL_solver>:
push   %r15
push   %r14
push   %rbx
sub    $0xe0,%rsp
mov    %rdx,%r15
movsd  %xmm7,0x90(%rsp)
movsd  %xmm4,0xd8(%rsp)
movaps %xmm3,0xa0(%rsp)
movapd %xmm2,0x30(%rsp)
movapd %xmm0,%xmm3
mov    %rsi,%rbx
mov    %rdi,%r14
movsd  %xmm1,0x98(%rsp)
mulsd  %xmm1,%xmm0
divsd  %xmm2,%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
movsd  %xmm6,(%rsp)
movsd  %xmm5,0x58(%rsp)
movsd  %xmm3,0x10(%rsp)
jb     117b <Roe_HLL_solver+0x6b>
sqrtsd %xmm0,%xmm0
jmp    1195 <Roe_HLL_solver+0x85>
call   1030 <sqrt@plt>
xorpd  %xmm1,%xmm1
movsd  0x10(%rsp),%xmm3
movsd  0x58(%rsp),%xmm5
movsd  (%rsp),%xmm6
movsd  %xmm0,0xd0(%rsp)
mulsd  %xmm5,%xmm3
divsd  %xmm6,%xmm3
ucomisd %xmm1,%xmm3
jb     11b5 <Roe_HLL_solver+0xa5>
xorps  %xmm0,%xmm0
sqrtsd %xmm3,%xmm0
jmp    11c3 <Roe_HLL_solver+0xb3>
movapd %xmm3,%xmm0
call   1030 <sqrt@plt>
movsd  (%rsp),%xmm6
movsd  %xmm0,0xc8(%rsp)
movapd 0x30(%rsp),%xmm1
movapd %xmm1,%xmm0
mulsd  %xmm6,%xmm0
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm0
jb     11ea <Roe_HLL_solver+0xda>
sqrtsd %xmm0,%xmm0
jmp    11f9 <Roe_HLL_solver+0xe9>
call   1030 <sqrt@plt>
xorpd  %xmm2,%xmm2
movapd 0x30(%rsp),%xmm1
movsd  %xmm0,0xc0(%rsp)
ucomisd %xmm2,%xmm1
jb     1211 <Roe_HLL_solver+0x101>
xorps  %xmm2,%xmm2
sqrtsd %xmm1,%xmm2
jmp    121e <Roe_HLL_solver+0x10e>
movapd %xmm1,%xmm0
call   1030 <sqrt@plt>
movapd %xmm0,%xmm2
movsd  (%rsp),%xmm0
movsd  0x10(%rsp),%xmm4
movsd  0xdcf(%rip),%xmm1
addsd  %xmm4,%xmm1
movsd  %xmm1,0x88(%rsp)
xorpd  %xmm3,%xmm3
ucomisd %xmm3,%xmm0
movapd %xmm2,0x70(%rsp)
jb     1257 <Roe_HLL_solver+0x147>
xorps  %xmm1,%xmm1
sqrtsd %xmm0,%xmm1
jmp    126a <Roe_HLL_solver+0x15a>
call   1030 <sqrt@plt>
xorpd  %xmm3,%xmm3
movsd  0x10(%rsp),%xmm4
movapd %xmm0,%xmm1
movapd 0x30(%rsp),%xmm0
divsd  0x88(%rsp),%xmm4
ucomisd %xmm3,%xmm0
movsd  %xmm1,0x80(%rsp)
jb     128e <Roe_HLL_solver+0x17e>
sqrtsd %xmm0,%xmm0
jmp    129f <Roe_HLL_solver+0x18f>
movsd  %xmm4,0x10(%rsp)
call   1030 <sqrt@plt>
movsd  0x10(%rsp),%xmm4
movsd  %xmm0,0x50(%rsp)
movsd  (%rsp),%xmm1
movsd  0x58(%rsp),%xmm0
movsd  0xd58(%rip),%xmm7
movapd 0xa0(%rsp),%xmm5
mulsd  %xmm5,%xmm5
movsd  0xd3b(%rip),%xmm2
movapd %xmm4,%xmm15
mulsd  0x98(%rsp),%xmm15
mulsd  %xmm0,%xmm4
movsd  0x90(%rsp),%xmm6
mulsd  %xmm6,%xmm6
xorpd  %xmm3,%xmm3
ucomisd %xmm3,%xmm1
jb     1300 <Roe_HLL_solver+0x1f0>
xorps  %xmm0,%xmm0
sqrtsd %xmm1,%xmm0
jmp    1354 <Roe_HLL_solver+0x244>
movapd %xmm1,%xmm0
movsd  %xmm15,0x8(%rsp)
movsd  %xmm4,0x10(%rsp)
movapd %xmm5,0x60(%rsp)
movsd  %xmm6,0x28(%rsp)
call   1030 <sqrt@plt>
xorpd  %xmm3,%xmm3
movsd  0x28(%rsp),%xmm6
movsd  0xcdc(%rip),%xmm7
movapd 0x60(%rsp),%xmm5
movsd  0x10(%rsp),%xmm4
movsd  0xcc0(%rip),%xmm2
movsd  0x8(%rsp),%xmm15
movsd  (%rsp),%xmm1
movsd  %xmm0,0xb8(%rsp)
movapd 0x30(%rsp),%xmm0
divsd  0x88(%rsp),%xmm7
mulsd  %xmm2,%xmm5
divsd  %xmm0,%xmm15
divsd  %xmm1,%xmm4
mulsd  %xmm2,%xmm6
ucomisd %xmm3,%xmm0
movapd %xmm5,0x60(%rsp)
jb     1390 <Roe_HLL_solver+0x280>
sqrtsd %xmm0,%xmm9
jmp    13d7 <Roe_HLL_solver+0x2c7>
movsd  %xmm15,0x8(%rsp)
movsd  %xmm4,0x10(%rsp)
movsd  %xmm7,0x48(%rsp)
movsd  %xmm6,0x28(%rsp)
call   1030 <sqrt@plt>
movsd  0x28(%rsp),%xmm6
movsd  0x48(%rsp),%xmm7
movapd 0x60(%rsp),%xmm5
movsd  0x10(%rsp),%xmm4
movsd  0x8(%rsp),%xmm15
movsd  (%rsp),%xmm1
movapd %xmm0,%xmm9
movapd 0x70(%rsp),%xmm3
movsd  0x80(%rsp),%xmm8
mulsd  0x98(%rsp),%xmm7
mulsd  0x90(%rsp),%xmm8
mulsd  0xa0(%rsp),%xmm3
addsd  %xmm5,%xmm15
addsd  %xmm4,%xmm6
xorpd  %xmm4,%xmm4
ucomisd %xmm4,%xmm1
jb     141f <Roe_HLL_solver+0x30f>
xorps  %xmm0,%xmm0
sqrtsd %xmm1,%xmm0
jmp    1480 <Roe_HLL_solver+0x370>
movapd %xmm1,%xmm0
movapd %xmm3,0x70(%rsp)
movsd  %xmm15,0x8(%rsp)
movsd  %xmm8,0x80(%rsp)
movsd  %xmm7,0x48(%rsp)
movsd  %xmm6,0x28(%rsp)
movsd  %xmm9,0x10(%rsp)
call   1030 <sqrt@plt>
xorpd  %xmm4,%xmm4
movsd  0x10(%rsp),%xmm9
movsd  0x28(%rsp),%xmm6
movsd  0x48(%rsp),%xmm7
movsd  0x80(%rsp),%xmm8
movsd  0x8(%rsp),%xmm15
movapd 0x70(%rsp),%xmm3
movapd 0x30(%rsp),%xmm1
movsd  0x50(%rsp),%xmm2
mulsd  %xmm0,%xmm6
divsd  %xmm1,%xmm7
addsd  %xmm8,%xmm3
addsd  0xb8(%rsp),%xmm2
mulsd  %xmm9,%xmm15
ucomisd %xmm4,%xmm1
jb     14b6 <Roe_HLL_solver+0x3a6>
xorps  %xmm4,%xmm4
sqrtsd %xmm1,%xmm4
jmp    1501 <Roe_HLL_solver+0x3f1>
movapd %xmm1,%xmm0
movapd %xmm3,0x70(%rsp)
movsd  %xmm15,0x8(%rsp)
movsd  %xmm2,0x50(%rsp)
movsd  %xmm7,0x48(%rsp)
movsd  %xmm6,0x28(%rsp)
call   1030 <sqrt@plt>
movsd  0x28(%rsp),%xmm6
movsd  0x48(%rsp),%xmm7
movsd  0x50(%rsp),%xmm2
movsd  0x8(%rsp),%xmm15
movapd 0x70(%rsp),%xmm3
movapd %xmm0,%xmm4
movsd  (%rsp),%xmm0
movapd 0x60(%rsp),%xmm1
addsd  %xmm7,%xmm1
divsd  %xmm2,%xmm3
addsd  %xmm6,%xmm15
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm0
movapd %xmm3,0x70(%rsp)
jb     152f <Roe_HLL_solver+0x41f>
sqrtsd %xmm0,%xmm0
jmp    1564 <Roe_HLL_solver+0x454>
movsd  %xmm15,0x8(%rsp)
movapd %xmm1,0x60(%rsp)
movsd  %xmm4,0x10(%rsp)
call   1030 <sqrt@plt>
movsd  0x10(%rsp),%xmm4
xorpd  %xmm2,%xmm2
movapd 0x60(%rsp),%xmm1
movsd  0x8(%rsp),%xmm15
movapd 0x70(%rsp),%xmm3
movapd 0x30(%rsp),%xmm10
addsd  %xmm0,%xmm4
divsd  %xmm4,%xmm15
movapd %xmm10,%xmm5
unpcklpd %xmm10,%xmm5
movapd 0xa0(%rsp),%xmm4
unpcklpd %xmm1,%xmm4
movsd  0xa85(%rip),%xmm0
mulsd  %xmm3,%xmm0
mulsd  %xmm3,%xmm0
addsd  %xmm15,%xmm0
mulsd  0x88(%rsp),%xmm0
ucomisd %xmm2,%xmm0
jb     15b5 <Roe_HLL_solver+0x4a5>
sqrtsd %xmm0,%xmm0
jmp    15ed <Roe_HLL_solver+0x4dd>
movsd  %xmm15,0x8(%rsp)
movapd %xmm5,0x10(%rsp)
movapd %xmm4,0x60(%rsp)
call   1030 <sqrt@plt>
movapd 0x60(%rsp),%xmm4
movapd 0x10(%rsp),%xmm5
movsd  0x8(%rsp),%xmm15
movapd 0x70(%rsp),%xmm3
movapd 0x30(%rsp),%xmm10
movsd  0x58(%rsp),%xmm11
movsd  0xa0b(%rip),%xmm14
mulpd  %xmm4,%xmm5
movapd %xmm3,%xmm13
subsd  %xmm0,%xmm13
movapd %xmm3,%xmm9
mulsd  %xmm0,%xmm9
subsd  0x98(%rsp),%xmm11
movsd  0xc0(%rsp),%xmm1
mulsd  %xmm0,%xmm1
movsd  0x90(%rsp),%xmm4
subsd  0xa0(%rsp),%xmm4
mulsd  %xmm1,%xmm4
movapd %xmm0,%xmm12
mulsd  %xmm0,%xmm12
movapd 0x9cc(%rip),%xmm2
andpd  %xmm3,%xmm2
addsd  %xmm0,%xmm2
movsd  %xmm2,(%r15)
xorpd  %xmm8,%xmm8
ucomisd %xmm13,%xmm8
jae    1674 <Roe_HLL_solver+0x564>
movapd %xmm10,%xmm2
jmp    16af <Roe_HLL_solver+0x59f>
movapd %xmm11,%xmm7
subsd  %xmm4,%xmm7
mulsd  %xmm14,%xmm7
divsd  %xmm12,%xmm7
movapd %xmm15,%xmm6
subsd  %xmm9,%xmm6
movapd %xmm7,%xmm2
addsd  %xmm10,%xmm2
unpcklpd %xmm7,%xmm7
movapd %xmm13,%xmm1
unpcklpd %xmm6,%xmm1
mulpd  %xmm7,%xmm1
addpd  %xmm1,%xmm5
addsd  %xmm3,%xmm0
ucomisd %xmm3,%xmm8
jae    16c2 <Roe_HLL_solver+0x5b2>
movsd  (%rsp),%xmm8
jmp    1701 <Roe_HLL_solver+0x5f1>
movapd %xmm3,%xmm1
mulsd  %xmm3,%xmm1
movsd  (%rsp),%xmm8
movapd %xmm8,%xmm6
subsd  %xmm10,%xmm6
movapd %xmm11,%xmm7
divsd  %xmm12,%xmm7
mulsd  %xmm14,%xmm1
subsd  %xmm7,%xmm6
addsd  %xmm6,%xmm2
unpcklpd %xmm6,%xmm6
unpcklpd %xmm1,%xmm3
mulpd  %xmm6,%xmm3
addpd  %xmm3,%xmm5
movsd  0x100(%rsp),%xmm6
xorpd  %xmm7,%xmm7
ucomisd %xmm0,%xmm7
jb     1741 <Roe_HLL_solver+0x631>
addsd  %xmm11,%xmm4
mulsd  %xmm14,%xmm4
addsd  %xmm15,%xmm9
divsd  %xmm12,%xmm4
addsd  %xmm4,%xmm2
unpcklpd %xmm4,%xmm4
movapd %xmm0,%xmm1
unpcklpd %xmm9,%xmm1
mulpd  %xmm4,%xmm1
addpd  %xmm1,%xmm5
movsd  %xmm2,(%rbx)
movapd %xmm5,%xmm1
divsd  %xmm2,%xmm1
movsd  %xmm1,0x8(%rbx)
unpckhpd %xmm5,%xmm5
divsd  %xmm2,%xmm5
movsd  0x8b6(%rip),%xmm4
mulsd  %xmm1,%xmm4
mulsd  %xmm1,%xmm4
addsd  %xmm5,%xmm4
movsd  0x88(%rsp),%xmm1
mulsd  %xmm4,%xmm1
mulsd  %xmm2,%xmm1
movsd  %xmm1,0x10(%rbx)
movapd 0xa0(%rsp),%xmm2
movapd %xmm2,%xmm1
subsd  0xd0(%rsp),%xmm1
minsd  %xmm1,%xmm13
movsd  0x90(%rsp),%xmm4
movsd  0xc8(%rsp),%xmm1
addsd  %xmm4,%xmm1
maxsd  %xmm1,%xmm0
xorpd  %xmm1,%xmm1
minsd  %xmm13,%xmm1
maxsd  %xmm0,%xmm7
subsd  %xmm4,%xmm7
subsd  %xmm2,%xmm1
xorpd  0x85a(%rip),%xmm1
mulsd  %xmm8,%xmm7
mulsd  %xmm10,%xmm1
movsd  0xd8(%rsp),%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm7,%xmm6
addsd  %xmm0,%xmm6
addsd  %xmm7,%xmm1
divsd  %xmm1,%xmm6
movsd  %xmm6,(%r14)
add    $0xe0,%rsp
pop    %rbx
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @Roe_HLL_solver(double* %0, double* %1, double %2, double %3, double %4, double %5, double %6, double %7, double %8, double %9, double %10, double* %11, double %12) {
%14 = alloca double*
%15 = alloca double*
%16 = alloca double
%17 = alloca double
%18 = alloca double
%19 = alloca double
%20 = alloca double
%21 = alloca double
%22 = alloca double
%23 = alloca double
%24 = alloca double
%25 = alloca double*
%26 = alloca double
%27 = alloca double
%28 = alloca double
%29 = alloca double
%30 = alloca double
%31 = alloca double
%32 = alloca double
%33 = alloca double
%34 = alloca [3 x double]
%35 = alloca double
%36 = alloca double
%37 = alloca double
%38 = alloca double
%39 = alloca [3 x [3 x double]]
%40 = alloca [3 x double]
%41 = alloca [3 x double]
%42 = alloca i32
%43 = alloca i32
%44 = alloca double
%45 = alloca double
store double* %0, double** %14
store double* %1, double** %15
store double %2, double* %16
store double %3, double* %17
store double %4, double* %18
store double %5, double* %19
store double %6, double* %20
store double %7, double* %21
store double %8, double* %22
store double %9, double* %23
store double %10, double* %24
store double* %11, double** %25
store double %12, double* %26
%46 = load double, double* %16
%47 = load double, double* %17
%48 = fmul double %46, %47
%49 = load double, double* %18
%50 = fdiv double %48, %49
%51 = call double @sqrt(double %50)
store double %51, double* %27
%52 = load double, double* %16
%53 = load double, double* %21
%54 = fmul double %52, %53
%55 = load double, double* %22
%56 = fdiv double %54, %55
%57 = call double @sqrt(double %56)
store double %57, double* %28
%58 = load double, double* %16
%59 = fsub double %58, 1.000000e+00
%60 = fmul double 5.000000e-01, %59
%61 = load double, double* %16
%62 = fdiv double %60, %61
store double %62, double* %29
%63 = load double, double* %16
%64 = load double, double* %16
%65 = fsub double %64, 1.000000e+00
%66 = fdiv double %63, %65
%67 = load double, double* %17
%68 = fmul double %66, %67
%69 = load double, double* %18
%70 = fdiv double %68, %69
%71 = load double, double* %19
%72 = load double, double* %19
%73 = fmul double %71, %72
%74 = call double @llvm.fmuladd.f64(double 5.000000e-01, double %73, double %70)
store double %74, double* %30
%75 = load double, double* %16
%76 = load double, double* %16
%77 = fsub double %76, 1.000000e+00
%78 = fdiv double %75, %77
%79 = load double, double* %21
%80 = fmul double %78, %79
%81 = load double, double* %22
%82 = fdiv double %80, %81
%83 = load double, double* %23
%84 = load double, double* %23
%85 = fmul double %83, %84
%86 = call double @llvm.fmuladd.f64(double 5.000000e-01, double %85, double %82)
store double %86, double* %31
%87 = load double, double* %16
%88 = fsub double %87, 1.000000e+00
%89 = fdiv double 1.000000e+00, %88
%90 = load double, double* %17
%91 = fmul double %89, %90
%92 = load double, double* %18
%93 = fdiv double %91, %92
%94 = load double, double* %19
%95 = load double, double* %19
%96 = fmul double %94, %95
%97 = call double @llvm.fmuladd.f64(double 5.000000e-01, double %96, double %93)
store double %97, double* %32
%98 = load double, double* %16
%99 = fsub double %98, 1.000000e+00
%100 = fdiv double 1.000000e+00, %99
%101 = load double, double* %21
%102 = fmul double %100, %101
%103 = load double, double* %22
%104 = fdiv double %102, %103
%105 = load double, double* %23
%106 = load double, double* %23
%107 = fmul double %105, %106
%108 = call double @llvm.fmuladd.f64(double 5.000000e-01, double %107, double %104)
store double %108, double* %33
%109 = load double, double* %18
%110 = getelementptr inbounds [3 x double], [3 x double]* %34, i64 0, i64 0
store double %109, double* %110
%111 = load double, double* %18
%112 = load double, double* %19
%113 = fmul double %111, %112
%114 = getelementptr inbounds [3 x double], [3 x double]* %34, i64 0, i64 1
store double %113, double* %114
%115 = load double, double* %18
%116 = load double, double* %32
%117 = fmul double %115, %116
%118 = getelementptr inbounds [3 x double], [3 x double]* %34, i64 0, i64 2
store double %117, double* %118
%119 = load double, double* %18
%120 = load double, double* %22
%121 = fmul double %119, %120
%122 = call double @sqrt(double %121)
store double %122, double* %35
%123 = load double, double* %19
%124 = load double, double* %18
%125 = call double @sqrt(double %124)
%126 = load double, double* %23
%127 = load double, double* %22
%128 = call double @sqrt(double %127)
%129 = fmul double %126, %128
%130 = call double @llvm.fmuladd.f64(double %123, double %125, double %129)
%131 = load double, double* %18
%132 = call double @sqrt(double %131)
%133 = load double, double* %22
%134 = call double @sqrt(double %133)
%135 = fadd double %132, %134
%136 = fdiv double %130, %135
store double %136, double* %36
%137 = load double, double* %30
%138 = load double, double* %18
%139 = call double @sqrt(double %138)
%140 = load double, double* %31
%141 = load double, double* %22
%142 = call double @sqrt(double %141)
%143 = fmul double %140, %142
%144 = call double @llvm.fmuladd.f64(double %137, double %139, double %143)
%145 = load double, double* %18
%146 = call double @sqrt(double %145)
%147 = load double, double* %22
%148 = call double @sqrt(double %147)
%149 = fadd double %146, %148
%150 = fdiv double %144, %149
store double %150, double* %37
%151 = load double, double* %16
%152 = fsub double %151, 1.000000e+00
%153 = load double, double* %37
%154 = load double, double* %36
%155 = fmul double 5.000000e-01, %154
%156 = load double, double* %36
%157 = fneg double %155
%158 = call double @llvm.fmuladd.f64(double %157, double %156, double %153)
%159 = fmul double %152, %158
%160 = call double @sqrt(double %159)
store double %160, double* %38
%161 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %39, i64 0, i64 0
%162 = getelementptr inbounds [3 x double], [3 x double]* %161, i64 0, i64 0
store double 1.000000e+00, double* %162
%163 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %39, i64 0, i64 0
%164 = getelementptr inbounds [3 x double], [3 x double]* %163, i64 0, i64 1
store double 1.000000e+00, double* %164
%165 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %39, i64 0, i64 0
%166 = getelementptr inbounds [3 x double], [3 x double]* %165, i64 0, i64 2
store double 1.000000e+00, double* %166
%167 = load double, double* %36
%168 = load double, double* %38
%169 = fsub double %167, %168
%170 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %39, i64 0, i64 1
%171 = getelementptr inbounds [3 x double], [3 x double]* %170, i64 0, i64 0
store double %169, double* %171
%172 = load double, double* %36
%173 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %39, i64 0, i64 1
%174 = getelementptr inbounds [3 x double], [3 x double]* %173, i64 0, i64 1
store double %172, double* %174
%175 = load double, double* %36
%176 = load double, double* %38
%177 = fadd double %175, %176
%178 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %39, i64 0, i64 1
%179 = getelementptr inbounds [3 x double], [3 x double]* %178, i64 0, i64 2
store double %177, double* %179
%180 = load double, double* %37
%181 = load double, double* %36
%182 = load double, double* %38
%183 = fneg double %181
%184 = call double @llvm.fmuladd.f64(double %183, double %182, double %180)
%185 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %39, i64 0, i64 2
%186 = getelementptr inbounds [3 x double], [3 x double]* %185, i64 0, i64 0
store double %184, double* %186
%187 = load double, double* %36
%188 = load double, double* %36
%189 = fmul double %187, %188
%190 = fmul double 5.000000e-01, %189
%191 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %39, i64 0, i64 2
%192 = getelementptr inbounds [3 x double], [3 x double]* %191, i64 0, i64 1
store double %190, double* %192
%193 = load double, double* %37
%194 = load double, double* %36
%195 = load double, double* %38
%196 = call double @llvm.fmuladd.f64(double %194, double %195, double %193)
%197 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %39, i64 0, i64 2
%198 = getelementptr inbounds [3 x double], [3 x double]* %197, i64 0, i64 2
store double %196, double* %198
%199 = load double, double* %21
%200 = load double, double* %17
%201 = fsub double %199, %200
%202 = load double, double* %35
%203 = load double, double* %38
%204 = fmul double %202, %203
%205 = load double, double* %23
%206 = load double, double* %19
%207 = fsub double %205, %206
%208 = fneg double %204
%209 = call double @llvm.fmuladd.f64(double %208, double %207, double %201)
%210 = fmul double 5.000000e-01, %209
%211 = load double, double* %38
%212 = load double, double* %38
%213 = fmul double %211, %212
%214 = fdiv double %210, %213
%215 = getelementptr inbounds [3 x double], [3 x double]* %41, i64 0, i64 0
store double %214, double* %215
%216 = load double, double* %22
%217 = load double, double* %18
%218 = fsub double %216, %217
%219 = load double, double* %21
%220 = load double, double* %17
%221 = fsub double %219, %220
%222 = load double, double* %38
%223 = load double, double* %38
%224 = fmul double %222, %223
%225 = fdiv double %221, %224
%226 = fsub double %218, %225
%227 = getelementptr inbounds [3 x double], [3 x double]* %41, i64 0, i64 1
store double %226, double* %227
%228 = load double, double* %21
%229 = load double, double* %17
%230 = fsub double %228, %229
%231 = load double, double* %35
%232 = load double, double* %38
%233 = fmul double %231, %232
%234 = load double, double* %23
%235 = load double, double* %19
%236 = fsub double %234, %235
%237 = call double @llvm.fmuladd.f64(double %233, double %236, double %230)
%238 = fmul double 5.000000e-01, %237
%239 = load double, double* %38
%240 = load double, double* %38
%241 = fmul double %239, %240
%242 = fdiv double %238, %241
%243 = getelementptr inbounds [3 x double], [3 x double]* %41, i64 0, i64 2
store double %242, double* %243
%244 = load double, double* %36
%245 = load double, double* %38
%246 = fsub double %244, %245
%247 = getelementptr inbounds [3 x double], [3 x double]* %40, i64 0, i64 0
store double %246, double* %247
%248 = load double, double* %36
%249 = getelementptr inbounds [3 x double], [3 x double]* %40, i64 0, i64 1
store double %248, double* %249
%250 = load double, double* %36
%251 = load double, double* %38
%252 = fadd double %250, %251
%253 = getelementptr inbounds [3 x double], [3 x double]* %40, i64 0, i64 2
store double %252, double* %253
%254 = load double, double* %36
%255 = call double @llvm.fabs.f64(double %254)
%256 = load double, double* %38
%257 = fadd double %255, %256
%258 = load double*, double** %25
store double %257, double* %258
store i32 0, i32* %43
br label %259
259: ; preds = %294, %13
%260 = load i32, i32* %43
%261 = icmp slt i32 %260, 3
br i1 %261, label %262, label %297
262: ; preds = %259
%263 = load i32, i32* %43
%264 = sext i32 %263 to i64
%265 = getelementptr inbounds [3 x double], [3 x double]* %40, i64 0, i64 %264
%266 = load double, double* %265
%267 = fcmp ole double %266, 0.000000e+00
br i1 %267, label %268, label %293
268: ; preds = %262
store i32 0, i32* %42
br label %269
269: ; preds = %289, %268
%270 = load i32, i32* %42
%271 = icmp slt i32 %270, 3
br i1 %271, label %272, label %292
272: ; preds = %269
%273 = load i32, i32* %43
%274 = sext i32 %273 to i64
%275 = getelementptr inbounds [3 x double], [3 x double]* %41, i64 0, i64 %274
%276 = load double, double* %275
%277 = load i32, i32* %42
%278 = sext i32 %277 to i64
%279 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %39, i64 0, i64 %278
%280 = load i32, i32* %43
%281 = sext i32 %280 to i64
%282 = getelementptr inbounds [3 x double], [3 x double]* %279, i64 0, i64 %281
%283 = load double, double* %282
%284 = load i32, i32* %42
%285 = sext i32 %284 to i64
%286 = getelementptr inbounds [3 x double], [3 x double]* %34, i64 0, i64 %285
%287 = load double, double* %286
%288 = call double @llvm.fmuladd.f64(double %276, double %283, double %287)
store double %288, double* %286
br label %289
289: ; preds = %272
%290 = load i32, i32* %42
%291 = add nsw i32 %290, 1
store i32 %291, i32* %42
br label %269
292: ; preds = %269
br label %293
293: ; preds = %292, %262
br label %294
294: ; preds = %293
%295 = load i32, i32* %43
%296 = add nsw i32 %295, 1
store i32 %296, i32* %43
br label %259
297: ; preds = %259
%298 = getelementptr inbounds [3 x double], [3 x double]* %34, i64 0, i64 0
%299 = load double, double* %298
%300 = load double*, double** %15
%301 = getelementptr inbounds double, double* %300, i64 0
store double %299, double* %301
%302 = getelementptr inbounds [3 x double], [3 x double]* %34, i64 0, i64 1
%303 = load double, double* %302
%304 = getelementptr inbounds [3 x double], [3 x double]* %34, i64 0, i64 0
%305 = load double, double* %304
%306 = fdiv double %303, %305
%307 = load double*, double** %15
%308 = getelementptr inbounds double, double* %307, i64 1
store double %306, double* %308
%309 = getelementptr inbounds [3 x double], [3 x double]* %34, i64 0, i64 2
%310 = load double, double* %309
%311 = getelementptr inbounds [3 x double], [3 x double]* %34, i64 0, i64 0
%312 = load double, double* %311
%313 = fdiv double %310, %312
%314 = load double*, double** %15
%315 = getelementptr inbounds double, double* %314, i64 1
%316 = load double, double* %315
%317 = fmul double 5.000000e-01, %316
%318 = load double*, double** %15
%319 = getelementptr inbounds double, double* %318, i64 1
%320 = load double, double* %319
%321 = fneg double %317
%322 = call double @llvm.fmuladd.f64(double %321, double %320, double %313)
%323 = load double, double* %16
%324 = fsub double %323, 1.000000e+00
%325 = fmul double %322, %324
%326 = load double*, double** %15
%327 = getelementptr inbounds double, double* %326, i64 0
%328 = load double, double* %327
%329 = fmul double %325, %328
%330 = load double*, double** %15
%331 = getelementptr inbounds double, double* %330, i64 2
store double %329, double* %331
%332 = load double, double* %36
%333 = load double, double* %38
%334 = fsub double %332, %333
%335 = load double, double* %19
%336 = load double, double* %27
%337 = fsub double %335, %336
%338 = fcmp olt double %334, %337
br i1 %338, label %339, label %343
339: ; preds = %297
%340 = load double, double* %36
%341 = load double, double* %38
%342 = fsub double %340, %341
br label %347
343: ; preds = %297
%344 = load double, double* %19
%345 = load double, double* %27
%346 = fsub double %344, %345
br label %347
347: ; preds = %343, %339
%348 = phi double [ %342, %339 ], [ %346, %343 ]
store double %348, double* %44
%349 = load double, double* %36
%350 = load double, double* %38
%351 = fadd double %349, %350
%352 = load double, double* %23
%353 = load double, double* %28
%354 = fadd double %352, %353
%355 = fcmp ogt double %351, %354
br i1 %355, label %356, label %360
356: ; preds = %347
%357 = load double, double* %36
%358 = load double, double* %38
%359 = fadd double %357, %358
br label %364
360: ; preds = %347
%361 = load double, double* %23
%362 = load double, double* %28
%363 = fadd double %361, %362
br label %364
364: ; preds = %360, %356
%365 = phi double [ %359, %356 ], [ %363, %360 ]
store double %365, double* %45
%366 = load double, double* %44
%367 = fcmp olt double 0.000000e+00, %366
br i1 %367, label %368, label %369
368: ; preds = %364
br label %371
369: ; preds = %364
%370 = load double, double* %44
br label %371
371: ; preds = %369, %368
%372 = phi double [ 0.000000e+00, %368 ], [ %370, %369 ]
store double %372, double* %44
%373 = load double, double* %45
%374 = fcmp ogt double 0.000000e+00, %373
br i1 %374, label %375, label %376
375: ; preds = %371
br label %378
376: ; preds = %371
%377 = load double, double* %45
br label %378
378: ; preds = %376, %375
%379 = phi double [ 0.000000e+00, %375 ], [ %377, %376 ]
store double %379, double* %45
%380 = load double, double* %22
%381 = load double, double* %45
%382 = load double, double* %23
%383 = fsub double %381, %382
%384 = fmul double %380, %383
%385 = load double, double* %24
%386 = load double, double* %18
%387 = load double, double* %44
%388 = load double, double* %19
%389 = fsub double %387, %388
%390 = fmul double %386, %389
%391 = load double, double* %20
%392 = fmul double %390, %391
%393 = fneg double %392
%394 = call double @llvm.fmuladd.f64(double %384, double %385, double %393)
%395 = load double, double* %22
%396 = load double, double* %45
%397 = load double, double* %23
%398 = fsub double %396, %397
%399 = load double, double* %18
%400 = load double, double* %44
%401 = load double, double* %19
%402 = fsub double %400, %401
%403 = fmul double %399, %402
%404 = fneg double %403
%405 = call double @llvm.fmuladd.f64(double %395, double %398, double %404)
%406 = fdiv double %394, %405
%407 = load double*, double** %14
store double %406, double* %407
ret void
}

LLVM IR (-O2):
define void @Roe_HLL_solver(double* nocapture writeonly %0, double* nocapture %1, double %2, double %3, double %4, double %5, double %6, double %7, double %8, double %9, double %10, double* nocapture writeonly %11, double %12) {
%14 = fmul double %2, %3
%15 = fdiv double %14, %4
%16 = tail call double @sqrt(double %15)
%17 = fmul double %2, %7
%18 = fdiv double %17, %8
%19 = tail call double @sqrt(double %18)
%20 = fadd double %2, -1.000000e+00
%21 = fdiv double %2, %20
%22 = fmul double %21, %3
%23 = fdiv double %22, %4
%24 = fmul double %5, %5
%25 = tail call double @llvm.fmuladd.f64(double %24, double 5.000000e-01, double %23)
%26 = fmul double %21, %7
%27 = fdiv double %26, %8
%28 = fmul double %9, %9
%29 = tail call double @llvm.fmuladd.f64(double %28, double 5.000000e-01, double %27)
%30 = fdiv double 1.000000e+00, %20
%31 = fmul double %30, %3
%32 = fdiv double %31, %4
%33 = tail call double @llvm.fmuladd.f64(double %24, double 5.000000e-01, double %32)
%34 = insertelement <2 x double> poison, double %4, i64 0
%35 = shufflevector <2 x double> %34, <2 x double> poison, <2 x i32> zeroinitializer
%36 = insertelement <2 x double> poison, double %5, i64 0
%37 = insertelement <2 x double> %36, double %33, i64 1
%38 = fmul <2 x double> %35, %37
%39 = fmul double %4, %8
%40 = tail call double @sqrt(double %39)
%41 = tail call double @sqrt(double %4)
%42 = tail call double @sqrt(double %8)
%43 = fmul double %42, %9
%44 = tail call double @llvm.fmuladd.f64(double %5, double %41, double %43)
%45 = tail call double @sqrt(double %4)
%46 = tail call double @sqrt(double %8)
%47 = fadd double %45, %46
%48 = fdiv double %44, %47
%49 = tail call double @sqrt(double %4)
%50 = tail call double @sqrt(double %8)
%51 = fmul double %29, %50
%52 = tail call double @llvm.fmuladd.f64(double %25, double %49, double %51)
%53 = tail call double @sqrt(double %4)
%54 = tail call double @sqrt(double %8)
%55 = fadd double %53, %54
%56 = fdiv double %52, %55
%57 = fmul double %48, -5.000000e-01
%58 = tail call double @llvm.fmuladd.f64(double %57, double %48, double %56)
%59 = fmul double %20, %58
%60 = tail call double @sqrt(double %59)
%61 = fsub double %48, %60
%62 = fadd double %48, %60
%63 = fmul double %48, %48
%64 = fmul double %63, 5.000000e-01
%65 = tail call double @llvm.fmuladd.f64(double %48, double %60, double %56)
%66 = fsub double %7, %3
%67 = fmul double %40, %60
%68 = fsub double %9, %5
%69 = fneg double %67
%70 = tail call double @llvm.fmuladd.f64(double %69, double %68, double %66)
%71 = fmul double %70, 5.000000e-01
%72 = fmul double %60, %60
%73 = fdiv double %71, %72
%74 = fsub double %8, %4
%75 = fdiv double %66, %72
%76 = fsub double %74, %75
%77 = tail call double @llvm.fmuladd.f64(double %67, double %68, double %66)
%78 = fmul double %77, 5.000000e-01
%79 = fdiv double %78, %72
%80 = tail call double @llvm.fabs.f64(double %48)
%81 = fadd double %80, %60
store double %81, double* %11
%82 = fcmp ugt double %61, 0.000000e+00
br i1 %82, label %92, label %83
83: ; preds = %13
%84 = fneg double %48
%85 = tail call double @llvm.fmuladd.f64(double %84, double %60, double %56)
%86 = fadd double %73, %4
%87 = insertelement <2 x double> poison, double %73, i64 0
%88 = shufflevector <2 x double> %87, <2 x double> poison, <2 x i32> zeroinitializer
%89 = insertelement <2 x double> poison, double %61, i64 0
%90 = insertelement <2 x double> %89, double %85, i64 1
%91 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %88, <2 x double> %90, <2 x double> %38)
br label %92
92: ; preds = %83, %13
%93 = phi double [ %4, %13 ], [ %86, %83 ]
%94 = phi <2 x double> [ %38, %13 ], [ %91, %83 ]
%95 = fcmp ugt double %48, 0.000000e+00
br i1 %95, label %103, label %96
96: ; preds = %92
%97 = fadd double %76, %93
%98 = insertelement <2 x double> poison, double %76, i64 0
%99 = shufflevector <2 x double> %98, <2 x double> poison, <2 x i32> zeroinitializer
%100 = insertelement <2 x double> poison, double %48, i64 0
%101 = insertelement <2 x double> %100, double %64, i64 1
%102 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %99, <2 x double> %101, <2 x double> %94)
br label %103
103: ; preds = %96, %92
%104 = phi double [ %93, %92 ], [ %97, %96 ]
%105 = phi <2 x double> [ %94, %92 ], [ %102, %96 ]
%106 = fcmp ugt double %62, 0.000000e+00
br i1 %106, label %114, label %107
107: ; preds = %103
%108 = fadd double %79, %104
%109 = insertelement <2 x double> poison, double %79, i64 0
%110 = shufflevector <2 x double> %109, <2 x double> poison, <2 x i32> zeroinitializer
%111 = insertelement <2 x double> poison, double %62, i64 0
%112 = insertelement <2 x double> %111, double %65, i64 1
%113 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %110, <2 x double> %112, <2 x double> %105)
br label %114
114: ; preds = %107, %103
%115 = phi double [ %104, %103 ], [ %108, %107 ]
%116 = phi <2 x double> [ %105, %103 ], [ %113, %107 ]
store double %115, double* %1
%117 = extractelement <2 x double> %116, i64 0
%118 = fdiv double %117, %115
%119 = getelementptr inbounds double, double* %1, i64 1
store double %118, double* %119
%120 = extractelement <2 x double> %116, i64 1
%121 = fdiv double %120, %115
%122 = fmul double %118, -5.000000e-01
%123 = tail call double @llvm.fmuladd.f64(double %122, double %118, double %121)
%124 = fmul double %20, %123
%125 = fmul double %115, %124
%126 = getelementptr inbounds double, double* %1, i64 2
store double %125, double* %126
%127 = fsub double %5, %16
%128 = fcmp olt double %61, %127
%129 = select i1 %128, double %61, double %127
%130 = fadd double %19, %9
%131 = fcmp ogt double %62, %130
%132 = select i1 %131, double %62, double %130
%133 = fcmp ogt double %129, 0.000000e+00
%134 = select i1 %133, double 0.000000e+00, double %129
%135 = fcmp olt double %132, 0.000000e+00
%136 = select i1 %135, double 0.000000e+00, double %132
%137 = fsub double %136, %9
%138 = fmul double %137, %8
%139 = fsub double %134, %5
%140 = fneg double %139
%141 = fmul double %140, %4
%142 = fmul double %141, %6
%143 = tail call double @llvm.fmuladd.f64(double %138, double %10, double %142)
%144 = tail call double @llvm.fmuladd.f64(double %8, double %137, double %141)
%145 = fdiv double %143, %144
store double %145, double* %0
ret void
}