#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ bcon_calc_g(double *, double *, double *, double *,
                               double *);
int /*<<< orphan*/ lower_g(double *, double **, double *);
int /*<<< orphan*/ ncov_calc(double **, double *);
double pressure_rho0_u(double, double);
int /*<<< orphan*/ ucon_calc_g(double *, double **, double **, double *);
void primtoU_g(double prim[8], double gcov[4][4], double gcon[4][4],
               double U[8]) {
  int i, j;
  double rho0;
  double ucon[4], ucov[4], bcon[4], bcov[4], ncov[4];
  double gamma, n_dot_b, bsq, u, p, w;
  ucon_calc_g(prim, gcov, gcon, ucon);
  lower_g(ucon, gcov, ucov);
  ncov_calc(gcon, ncov);
  gamma = -ncov[0] * ucon[0];
  bcon_calc_g(prim, ucon, ucov, ncov, bcon);
  lower_g(bcon, gcov, bcov);
  n_dot_b = 0.;
  for (i = 0; i < 4; i++)
    n_dot_b += ncov[i] * bcon[i];
  bsq = 0.;
  for (i = 0; i < 4; i++)
    bsq += bcov[i] * bcon[i];
  rho0 = prim[0];
  u = prim[1];
  p = pressure_rho0_u(rho0, u);
  w = rho0 + u + p;
  U[0] = gamma * rho0;
  for (i = 0; i < 4; i++)
    U[1 + i] = gamma * (w + bsq) * ucov[i] - (p + bsq / 2.) * ncov[i] +
               n_dot_b * bcov[i];
  U[5] = prim[5];
  U[6] = prim[6];
  U[7] = prim[7];
  return;
}


x86 assembly (clang, -O0):
<primtoU_g>:
push   %rbp
mov    %rsp,%rbp
sub    $0x100,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    -0x18(%rbp),%rdx
lea    -0x50(%rbp),%rcx
call   1070 <ucon_calc_g@plt>
lea    -0x50(%rbp),%rdi
mov    -0x10(%rbp),%rsi
lea    -0x70(%rbp),%rdx
call   1050 <lower_g@plt>
mov    -0x18(%rbp),%rdi
lea    -0xd0(%rbp),%rsi
call   1030 <ncov_calc@plt>
movsd  -0xd0(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
mov    -0x8(%rbp),%rdi
lea    -0x50(%rbp),%rsi
lea    -0x70(%rbp),%rdx
lea    -0xd0(%rbp),%rcx
lea    -0x90(%rbp),%r8
call   1060 <bcon_calc_g@plt>
lea    -0x90(%rbp),%rdi
mov    -0x10(%rbp),%rsi
lea    -0xb0(%rbp),%rdx
call   1050 <lower_g@plt>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
movl   $0x0,-0x24(%rbp)
cmpl   $0x4,-0x24(%rbp)
jge    125f <primtoU_g+0x10f>
movslq -0x24(%rbp),%rax
movsd  -0xd0(%rbp,%rax,8),%xmm0
ff ff
movslq -0x24(%rbp),%rax
movsd  -0x90(%rbp,%rax,8),%xmm2
ff ff
movsd  -0xe0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1215 <primtoU_g+0xc5>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xe8(%rbp)
ff
movl   $0x0,-0x24(%rbp)
cmpl   $0x4,-0x24(%rbp)
jge    12bb <primtoU_g+0x16b>
movslq -0x24(%rbp),%rax
movsd  -0xb0(%rbp,%rax,8),%xmm0
ff ff
movslq -0x24(%rbp),%rax
movsd  -0x90(%rbp,%rax,8),%xmm2
ff ff
movsd  -0xe8(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xe8(%rbp)
ff
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1271 <primtoU_g+0x121>
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0xf0(%rbp)
ff
movsd  -0x30(%rbp),%xmm0
movsd  -0xf0(%rbp),%xmm1
ff
call   1040 <pressure_rho0_u@plt>
movsd  %xmm0,-0xf8(%rbp)
ff
movsd  -0x30(%rbp),%xmm0
addsd  -0xf0(%rbp),%xmm0
ff
addsd  -0xf8(%rbp),%xmm0
ff
movsd  %xmm0,-0x100(%rbp)
ff
movsd  -0xd8(%rbp),%xmm0
ff
mulsd  -0x30(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x0,-0x24(%rbp)
cmpl   $0x4,-0x24(%rbp)
jge    13e5 <primtoU_g+0x295>
movsd  -0xd8(%rbp),%xmm1
ff
movsd  -0x100(%rbp),%xmm0
ff
addsd  -0xe8(%rbp),%xmm0
ff
mulsd  %xmm0,%xmm1
movslq -0x24(%rbp),%rax
movsd  -0x70(%rbp,%rax,8),%xmm2
movsd  -0xf8(%rbp),%xmm0
ff
movsd  -0xe8(%rbp),%xmm3
ff
movsd  0xc8c(%rip),%xmm4
divsd  %xmm4,%xmm3
addsd  %xmm3,%xmm0
movslq -0x24(%rbp),%rax
mulsd  -0xd0(%rbp,%rax,8),%xmm0
ff ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0xe0(%rbp),%xmm0
ff
movslq -0x24(%rbp),%rax
movsd  -0xb0(%rbp,%rax,8),%xmm2
ff ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    132c <primtoU_g+0x1dc>
mov    -0x8(%rbp),%rax
movsd  0x28(%rax),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,0x28(%rax)
mov    -0x8(%rbp),%rax
movsd  0x30(%rax),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,0x30(%rax)
mov    -0x8(%rbp),%rax
movsd  0x38(%rax),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,0x38(%rax)
add    $0x100,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<primtoU_g>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xf8,%rsp
mov    %rcx,%r14
mov    %rdx,%rbx
mov    %rsi,%r12
mov    %rdi,%r15
lea    0xd0(%rsp),%r13
mov    %r13,%rcx
call   1070 <ucon_calc_g@plt>
lea    0xb0(%rsp),%rdx
mov    %r13,%rdi
mov    %r12,%rsi
call   1050 <lower_g@plt>
mov    %rsp,%rbp
mov    %rbx,%rdi
mov    %rbp,%rsi
call   1030 <ncov_calc@plt>
movsd  (%rsp),%xmm0
xorpd  0xe55(%rip),%xmm0
mulsd  0xd0(%rsp),%xmm0
movapd %xmm0,0xa0(%rsp)
lea    0x50(%rsp),%rbx
mov    %r15,%rdi
mov    %r13,%rsi
lea    0xb0(%rsp),%rdx
mov    %rbp,%rcx
mov    %rbx,%r8
call   1060 <bcon_calc_g@plt>
lea    0x30(%rsp),%rdx
mov    %rbx,%rdi
mov    %r12,%rsi
call   1050 <lower_g@plt>
movsd  0x50(%rsp),%xmm2
movsd  0x58(%rsp),%xmm1
movsd  0x60(%rsp),%xmm0
movsd  0x68(%rsp),%xmm6
movsd  0x30(%rsp),%xmm3
mulsd  %xmm2,%xmm3
xorpd  %xmm4,%xmm4
addsd  %xmm4,%xmm3
movsd  0x38(%rsp),%xmm5
mulsd  %xmm1,%xmm5
addsd  %xmm3,%xmm5
movsd  0x40(%rsp),%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm5,%xmm3
movsd  0x48(%rsp),%xmm5
mulsd  %xmm6,%xmm5
addsd  %xmm3,%xmm5
movapd %xmm5,0x80(%rsp)
mulsd  (%rsp),%xmm2
addsd  %xmm4,%xmm2
mulsd  0x8(%rsp),%xmm1
addsd  %xmm2,%xmm1
mulsd  0x10(%rsp),%xmm0
mulsd  0x18(%rsp),%xmm6
addsd  %xmm1,%xmm0
addsd  %xmm0,%xmm6
movapd %xmm6,0x90(%rsp)
movsd  (%r15),%xmm0
movsd  %xmm0,0x28(%rsp)
movsd  0x8(%r15),%xmm1
movaps %xmm1,0x70(%rsp)
call   1040 <pressure_rho0_u@plt>
movapd 0x70(%rsp),%xmm4
movsd  0x28(%rsp),%xmm2
addsd  %xmm2,%xmm4
addsd  %xmm0,%xmm4
movapd 0xa0(%rsp),%xmm1
mulsd  %xmm1,%xmm2
movsd  %xmm2,(%r14)
movapd 0x80(%rsp),%xmm2
addsd  %xmm2,%xmm4
mulsd  %xmm1,%xmm4
mulsd  0xd3e(%rip),%xmm2
addsd  %xmm0,%xmm2
unpcklpd %xmm2,%xmm2
movapd (%rsp),%xmm0
mulpd  %xmm2,%xmm0
movapd %xmm2,%xmm3
unpcklpd %xmm4,%xmm4
movapd 0xb0(%rsp),%xmm1
mulpd  %xmm4,%xmm1
subpd  %xmm0,%xmm1
movapd 0x90(%rsp),%xmm2
unpcklpd %xmm2,%xmm2
movapd 0x30(%rsp),%xmm0
mulpd  %xmm2,%xmm0
addpd  %xmm1,%xmm0
movupd %xmm0,0x8(%r14)
mulpd  0x10(%rsp),%xmm3
mulpd  0xc0(%rsp),%xmm4
subpd  %xmm3,%xmm4
mulpd  0x40(%rsp),%xmm2
addpd  %xmm4,%xmm2
movupd %xmm2,0x18(%r14)
movsd  0x28(%r15),%xmm0
movsd  %xmm0,0x28(%r14)
movsd  0x30(%r15),%xmm0
movsd  %xmm0,0x30(%r14)
movsd  0x38(%r15),%xmm0
movsd  %xmm0,0x38(%r14)
add    $0xf8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @primtoU_g(double* %0, [4 x double]* %1, [4 x double]* %2, double* %3) {
%5 = alloca double*
%6 = alloca [4 x double]*
%7 = alloca [4 x double]*
%8 = alloca double*
%9 = alloca i32
%10 = alloca i32
%11 = alloca double
%12 = alloca [4 x double]
%13 = alloca [4 x double]
%14 = alloca [4 x double]
%15 = alloca [4 x double]
%16 = alloca [4 x double]
%17 = alloca double
%18 = alloca double
%19 = alloca double
%20 = alloca double
%21 = alloca double
%22 = alloca double
store double* %0, double** %5
store [4 x double]* %1, [4 x double]** %6
store [4 x double]* %2, [4 x double]** %7
store double* %3, double** %8
%23 = load double*, double** %5
%24 = load [4 x double]*, [4 x double]** %6
%25 = bitcast [4 x double]* %24 to double**
%26 = load [4 x double]*, [4 x double]** %7
%27 = bitcast [4 x double]* %26 to double**
%28 = getelementptr inbounds [4 x double], [4 x double]* %12, i64 0, i64 0
%29 = call i32 @ucon_calc_g(double* %23, double** %25, double** %27, double* %28)
%30 = getelementptr inbounds [4 x double], [4 x double]* %12, i64 0, i64 0
%31 = load [4 x double]*, [4 x double]** %6
%32 = bitcast [4 x double]* %31 to double**
%33 = getelementptr inbounds [4 x double], [4 x double]* %13, i64 0, i64 0
%34 = call i32 @lower_g(double* %30, double** %32, double* %33)
%35 = load [4 x double]*, [4 x double]** %7
%36 = bitcast [4 x double]* %35 to double**
%37 = getelementptr inbounds [4 x double], [4 x double]* %16, i64 0, i64 0
%38 = call i32 @ncov_calc(double** %36, double* %37)
%39 = getelementptr inbounds [4 x double], [4 x double]* %16, i64 0, i64 0
%40 = load double, double* %39
%41 = fneg double %40
%42 = getelementptr inbounds [4 x double], [4 x double]* %12, i64 0, i64 0
%43 = load double, double* %42
%44 = fmul double %41, %43
store double %44, double* %17
%45 = load double*, double** %5
%46 = getelementptr inbounds [4 x double], [4 x double]* %12, i64 0, i64 0
%47 = getelementptr inbounds [4 x double], [4 x double]* %13, i64 0, i64 0
%48 = getelementptr inbounds [4 x double], [4 x double]* %16, i64 0, i64 0
%49 = getelementptr inbounds [4 x double], [4 x double]* %14, i64 0, i64 0
%50 = call i32 @bcon_calc_g(double* %45, double* %46, double* %47, double* %48, double* %49)
%51 = getelementptr inbounds [4 x double], [4 x double]* %14, i64 0, i64 0
%52 = load [4 x double]*, [4 x double]** %6
%53 = bitcast [4 x double]* %52 to double**
%54 = getelementptr inbounds [4 x double], [4 x double]* %15, i64 0, i64 0
%55 = call i32 @lower_g(double* %51, double** %53, double* %54)
store double 0.000000e+00, double* %18
store i32 0, i32* %9
br label %56
56: ; preds = %70, %4
%57 = load i32, i32* %9
%58 = icmp slt i32 %57, 4
br i1 %58, label %59, label %73
59: ; preds = %56
%60 = load i32, i32* %9
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds [4 x double], [4 x double]* %16, i64 0, i64 %61
%63 = load double, double* %62
%64 = load i32, i32* %9
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds [4 x double], [4 x double]* %14, i64 0, i64 %65
%67 = load double, double* %66
%68 = load double, double* %18
%69 = call double @llvm.fmuladd.f64(double %63, double %67, double %68)
store double %69, double* %18
br label %70
70: ; preds = %59
%71 = load i32, i32* %9
%72 = add nsw i32 %71, 1
store i32 %72, i32* %9
br label %56
73: ; preds = %56
store double 0.000000e+00, double* %19
store i32 0, i32* %9
br label %74
74: ; preds = %88, %73
%75 = load i32, i32* %9
%76 = icmp slt i32 %75, 4
br i1 %76, label %77, label %91
77: ; preds = %74
%78 = load i32, i32* %9
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds [4 x double], [4 x double]* %15, i64 0, i64 %79
%81 = load double, double* %80
%82 = load i32, i32* %9
%83 = sext i32 %82 to i64
%84 = getelementptr inbounds [4 x double], [4 x double]* %14, i64 0, i64 %83
%85 = load double, double* %84
%86 = load double, double* %19
%87 = call double @llvm.fmuladd.f64(double %81, double %85, double %86)
store double %87, double* %19
br label %88
88: ; preds = %77
%89 = load i32, i32* %9
%90 = add nsw i32 %89, 1
store i32 %90, i32* %9
br label %74
91: ; preds = %74
%92 = load double*, double** %5
%93 = getelementptr inbounds double, double* %92, i64 0
%94 = load double, double* %93
store double %94, double* %11
%95 = load double*, double** %5
%96 = getelementptr inbounds double, double* %95, i64 1
%97 = load double, double* %96
store double %97, double* %20
%98 = load double, double* %11
%99 = load double, double* %20
%100 = call double @pressure_rho0_u(double %98, double %99)
store double %100, double* %21
%101 = load double, double* %11
%102 = load double, double* %20
%103 = fadd double %101, %102
%104 = load double, double* %21
%105 = fadd double %103, %104
store double %105, double* %22
%106 = load double, double* %17
%107 = load double, double* %11
%108 = fmul double %106, %107
%109 = load double*, double** %8
%110 = getelementptr inbounds double, double* %109, i64 0
store double %108, double* %110
store i32 0, i32* %9
br label %111
111: ; preds = %146, %91
%112 = load i32, i32* %9
%113 = icmp slt i32 %112, 4
br i1 %113, label %114, label %149
114: ; preds = %111
%115 = load double, double* %17
%116 = load double, double* %22
%117 = load double, double* %19
%118 = fadd double %116, %117
%119 = fmul double %115, %118
%120 = load i32, i32* %9
%121 = sext i32 %120 to i64
%122 = getelementptr inbounds [4 x double], [4 x double]* %13, i64 0, i64 %121
%123 = load double, double* %122
%124 = load double, double* %21
%125 = load double, double* %19
%126 = fdiv double %125, 2.000000e+00
%127 = fadd double %124, %126
%128 = load i32, i32* %9
%129 = sext i32 %128 to i64
%130 = getelementptr inbounds [4 x double], [4 x double]* %16, i64 0, i64 %129
%131 = load double, double* %130
%132 = fmul double %127, %131
%133 = fneg double %132
%134 = call double @llvm.fmuladd.f64(double %119, double %123, double %133)
%135 = load double, double* %18
%136 = load i32, i32* %9
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds [4 x double], [4 x double]* %15, i64 0, i64 %137
%139 = load double, double* %138
%140 = call double @llvm.fmuladd.f64(double %135, double %139, double %134)
%141 = load double*, double** %8
%142 = load i32, i32* %9
%143 = add nsw i32 1, %142
%144 = sext i32 %143 to i64
%145 = getelementptr inbounds double, double* %141, i64 %144
store double %140, double* %145
br label %146
146: ; preds = %114
%147 = load i32, i32* %9
%148 = add nsw i32 %147, 1
store i32 %148, i32* %9
br label %111
149: ; preds = %111
%150 = load double*, double** %5
%151 = getelementptr inbounds double, double* %150, i64 5
%152 = load double, double* %151
%153 = load double*, double** %8
%154 = getelementptr inbounds double, double* %153, i64 5
store double %152, double* %154
%155 = load double*, double** %5
%156 = getelementptr inbounds double, double* %155, i64 6
%157 = load double, double* %156
%158 = load double*, double** %8
%159 = getelementptr inbounds double, double* %158, i64 6
store double %157, double* %159
%160 = load double*, double** %5
%161 = getelementptr inbounds double, double* %160, i64 7
%162 = load double, double* %161
%163 = load double*, double** %8
%164 = getelementptr inbounds double, double* %163, i64 7
store double %162, double* %164
ret void
}

LLVM IR (-O2):
define void @primtoU_g(double* %0, [4 x double]* %1, [4 x double]* %2, double* nocapture writeonly %3) {
%5 = alloca [4 x double]
%6 = alloca [4 x double]
%7 = alloca [4 x double]
%8 = alloca [4 x double]
%9 = alloca [4 x double]
%10 = bitcast [4 x double]* %5 to i8*
%11 = bitcast [4 x double]* %6 to i8*
%12 = bitcast [4 x double]* %7 to i8*
%13 = bitcast [4 x double]* %8 to i8*
%14 = bitcast [4 x double]* %9 to i8*
%15 = bitcast [4 x double]* %1 to double**
%16 = bitcast [4 x double]* %2 to double**
%17 = getelementptr inbounds [4 x double], [4 x double]* %5, i64 0, i64 0
%18 = call i32 @ucon_calc_g(double* %0, double** %15, double** %16, double* %17)
%19 = getelementptr inbounds [4 x double], [4 x double]* %6, i64 0, i64 0
%20 = call i32 @lower_g(double* %17, double** %15, double* %19)
%21 = getelementptr inbounds [4 x double], [4 x double]* %9, i64 0, i64 0
%22 = call i32 @ncov_calc(double** %16, double* %21)
%23 = load double, double* %21
%24 = load double, double* %17
%25 = getelementptr inbounds [4 x double], [4 x double]* %7, i64 0, i64 0
%26 = call i32 @bcon_calc_g(double* %0, double* %17, double* %19, double* %21, double* %25)
%27 = getelementptr inbounds [4 x double], [4 x double]* %8, i64 0, i64 0
%28 = call i32 @lower_g(double* %25, double** %15, double* %27)
%29 = load double, double* %21
%30 = load double, double* %25
%31 = getelementptr inbounds [4 x double], [4 x double]* %9, i64 0, i64 1
%32 = load double, double* %31
%33 = getelementptr inbounds [4 x double], [4 x double]* %7, i64 0, i64 1
%34 = load double, double* %33
%35 = getelementptr inbounds [4 x double], [4 x double]* %9, i64 0, i64 2
%36 = load double, double* %35
%37 = getelementptr inbounds [4 x double], [4 x double]* %7, i64 0, i64 2
%38 = load double, double* %37
%39 = getelementptr inbounds [4 x double], [4 x double]* %9, i64 0, i64 3
%40 = load double, double* %39
%41 = getelementptr inbounds [4 x double], [4 x double]* %7, i64 0, i64 3
%42 = load double, double* %41
%43 = fneg double %23
%44 = load double, double* %27
%45 = call double @llvm.fmuladd.f64(double %44, double %30, double 0.000000e+00)
%46 = getelementptr inbounds [4 x double], [4 x double]* %8, i64 0, i64 1
%47 = load double, double* %46
%48 = call double @llvm.fmuladd.f64(double %47, double %34, double %45)
%49 = getelementptr inbounds [4 x double], [4 x double]* %8, i64 0, i64 2
%50 = load double, double* %49
%51 = call double @llvm.fmuladd.f64(double %50, double %38, double %48)
%52 = getelementptr inbounds [4 x double], [4 x double]* %8, i64 0, i64 3
%53 = load double, double* %52
%54 = call double @llvm.fmuladd.f64(double %53, double %42, double %51)
%55 = call double @llvm.fmuladd.f64(double %29, double %30, double 0.000000e+00)
%56 = call double @llvm.fmuladd.f64(double %32, double %34, double %55)
%57 = call double @llvm.fmuladd.f64(double %36, double %38, double %56)
%58 = call double @llvm.fmuladd.f64(double %40, double %42, double %57)
%59 = fmul double %24, %43
%60 = load double, double* %0
%61 = getelementptr inbounds double, double* %0, i64 1
%62 = load double, double* %61
%63 = call double @pressure_rho0_u(double %60, double %62)
%64 = fadd double %60, %62
%65 = fadd double %64, %63
%66 = fmul double %59, %60
store double %66, double* %3
%67 = fadd double %54, %65
%68 = fmul double %59, %67
%69 = fmul double %54, 5.000000e-01
%70 = fadd double %69, %63
%71 = fneg double %70
%72 = getelementptr inbounds double, double* %3, i64 1
%73 = bitcast [4 x double]* %6 to <2 x double>*
%74 = load <2 x double>, <2 x double>* %73
%75 = bitcast [4 x double]* %9 to <2 x double>*
%76 = load <2 x double>, <2 x double>* %75
%77 = insertelement <2 x double> poison, double %71, i64 0
%78 = shufflevector <2 x double> %77, <2 x double> poison, <2 x i32> zeroinitializer
%79 = fmul <2 x double> %76, %78
%80 = insertelement <2 x double> poison, double %68, i64 0
%81 = shufflevector <2 x double> %80, <2 x double> poison, <2 x i32> zeroinitializer
%82 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %81, <2 x double> %74, <2 x double> %79)
%83 = bitcast [4 x double]* %8 to <2 x double>*
%84 = load <2 x double>, <2 x double>* %83
%85 = insertelement <2 x double> poison, double %58, i64 0
%86 = shufflevector <2 x double> %85, <2 x double> poison, <2 x i32> zeroinitializer
%87 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %86, <2 x double> %84, <2 x double> %82)
%88 = bitcast double* %72 to <2 x double>*
store <2 x double> %87, <2 x double>* %88
%89 = getelementptr inbounds [4 x double], [4 x double]* %6, i64 0, i64 2
%90 = getelementptr inbounds double, double* %3, i64 3
%91 = bitcast double* %89 to <2 x double>*
%92 = load <2 x double>, <2 x double>* %91
%93 = bitcast double* %35 to <2 x double>*
%94 = load <2 x double>, <2 x double>* %93
%95 = fmul <2 x double> %94, %78
%96 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %81, <2 x double> %92, <2 x double> %95)
%97 = bitcast double* %49 to <2 x double>*
%98 = load <2 x double>, <2 x double>* %97
%99 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %86, <2 x double> %98, <2 x double> %96)
%100 = bitcast double* %90 to <2 x double>*
store <2 x double> %99, <2 x double>* %100
%101 = getelementptr inbounds double, double* %0, i64 5
%102 = load double, double* %101
%103 = getelementptr inbounds double, double* %3, i64 5
store double %102, double* %103
%104 = getelementptr inbounds double, double* %0, i64 6
%105 = load double, double* %104
%106 = getelementptr inbounds double, double* %3, i64 6
store double %105, double* %106
%107 = getelementptr inbounds double, double* %0, i64 7
%108 = load double, double* %107
%109 = getelementptr inbounds double, double* %3, i64 7
store double %108, double* %109
ret void
}