void ReynoldsAveragedNavierStokes_kEpsilon_3D_Update(
    const int nPoints, const double nu, const double c_mu, const double *k,
    const double *grad_k, const double *epsilon, double *mom_u_diff_ten,
    double *mom_v_diff_ten, double *mom_w_diff_ten, double *mom_uv_diff_ten,
    double *mom_uw_diff_ten, double *mom_vu_diff_ten, double *mom_vw_diff_ten,
    double *mom_wu_diff_ten, double *mom_wv_diff_ten, double *mom_u_source,
    double *mom_v_source, double *mom_w_source) {
  int i;
  double nu_t, eval, kval;
  const double twoThirds = 2.0 / 3.0;
  const double div_zero = 1.0e-6;
  for (i = 0; i < nPoints; i++) {
    eval = epsilon[i] >= 0.0 ? epsilon[i] : 0.0;
    kval = k[i] >= 0.0 ? k[i] : 0.0;
    nu_t = c_mu * kval * kval / (eval + div_zero);
    mom_u_diff_ten[i * 9 + 0] += 2.0 * nu_t;
    mom_u_diff_ten[i * 9 + 4] += nu_t;
    mom_u_diff_ten[i * 9 + 8] += nu_t;
    mom_uv_diff_ten[i * 9 + 3] += nu_t;
    mom_uw_diff_ten[i * 9 + 6] += nu_t;
    mom_v_diff_ten[i * 9 + 0] += nu_t;
    mom_v_diff_ten[i * 9 + 4] += 2.0 * nu_t;
    mom_v_diff_ten[i * 9 + 8] += nu_t;
    mom_vu_diff_ten[i * 9 + 1] += nu_t;
    mom_vw_diff_ten[i * 9 + 7] += nu_t;
    mom_w_diff_ten[i * 9 + 0] += nu_t;
    mom_w_diff_ten[i * 9 + 4] += nu_t;
    mom_w_diff_ten[i * 9 + 8] += 2.0 * nu_t;
    mom_wu_diff_ten[i * 9 + 2] += nu_t;
    mom_wv_diff_ten[i * 9 + 5] += nu_t;
    mom_u_source[i] += twoThirds * grad_k[i * 3 + 0];
    mom_v_source[i] += twoThirds * grad_k[i * 3 + 1];
    mom_w_source[i] += twoThirds * grad_k[i * 3 + 2];
  }
}


x86 assembly (clang, -O0):
<ReynoldsAveragedNavierStokes_kEpsilon_3D_Update>:
push   %rbp
mov    %rsp,%rbp
mov    0x58(%rbp),%rax
mov    0x50(%rbp),%rax
mov    0x48(%rbp),%rax
mov    0x40(%rbp),%rax
mov    0x38(%rbp),%rax
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
mov    %rsi,-0x20(%rbp)
mov    %rdx,-0x28(%rbp)
mov    %rcx,-0x30(%rbp)
mov    %r8,-0x38(%rbp)
mov    %r9,-0x40(%rbp)
movsd  0xeb3(%rip),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  0xe9e(%rip),%xmm0
movsd  %xmm0,-0x70(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    14a7 <ReynoldsAveragedNavierStokes_kEpsilon_3D_Update+0x3a7>
mov    -0x30(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     11ab <ReynoldsAveragedNavierStokes_kEpsilon_3D_Update+0xab>
mov    -0x30(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    11b8 <ReynoldsAveragedNavierStokes_kEpsilon_3D_Update+0xb8>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    11b8 <ReynoldsAveragedNavierStokes_kEpsilon_3D_Update+0xb8>
movsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     11f3 <ReynoldsAveragedNavierStokes_kEpsilon_3D_Update+0xf3>
mov    -0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
jmp    1200 <ReynoldsAveragedNavierStokes_kEpsilon_3D_Update+0x100>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x80(%rbp)
jmp    1200 <ReynoldsAveragedNavierStokes_kEpsilon_3D_Update+0x100>
movsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x18(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
movsd  0xddf(%rip),%xmm1
addsd  -0x58(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
mov    -0x38(%rbp),%rax
imul   $0x9,-0x44(%rbp),%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
addsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    -0x38(%rbp),%rax
imul   $0x9,-0x44(%rbp),%ecx
add    $0x4,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    -0x38(%rbp),%rax
imul   $0x9,-0x44(%rbp),%ecx
add    $0x8,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    0x18(%rbp),%rax
imul   $0x9,-0x44(%rbp),%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    0x20(%rbp),%rax
imul   $0x9,-0x44(%rbp),%ecx
add    $0x6,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    -0x40(%rbp),%rax
imul   $0x9,-0x44(%rbp),%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    -0x40(%rbp),%rax
imul   $0x9,-0x44(%rbp),%ecx
add    $0x4,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
addsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    -0x40(%rbp),%rax
imul   $0x9,-0x44(%rbp),%ecx
add    $0x8,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    0x28(%rbp),%rax
imul   $0x9,-0x44(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    0x30(%rbp),%rax
imul   $0x9,-0x44(%rbp),%ecx
add    $0x7,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    0x10(%rbp),%rax
imul   $0x9,-0x44(%rbp),%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    0x10(%rbp),%rax
imul   $0x9,-0x44(%rbp),%ecx
add    $0x4,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    0x10(%rbp),%rax
imul   $0x9,-0x44(%rbp),%ecx
add    $0x8,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
addsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    0x38(%rbp),%rax
imul   $0x9,-0x44(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    0x40(%rbp),%rax
imul   $0x9,-0x44(%rbp),%ecx
add    $0x5,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
imul   $0x3,-0x44(%rbp),%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x48(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  0xbe6(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
imul   $0x3,-0x44(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x50(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  0xbb1(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
imul   $0x3,-0x44(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x58(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  0xb7c(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    116e <ReynoldsAveragedNavierStokes_kEpsilon_3D_Update+0x6e>
pop    %rbp
ret

x86 assembly (clang, -O2):
<ReynoldsAveragedNavierStokes_kEpsilon_3D_Update>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %edi,%edi
jle    1313 <ReynoldsAveragedNavierStokes_kEpsilon_3D_Update+0x213>
mov    0x50(%rsp),%r15
mov    0x48(%rsp),%rbx
mov    0x40(%rsp),%rbp
mov    0x38(%rsp),%r10
mov    %edi,%eax
shl    $0x3,%rax
lea    (%rax,%rax,8),%rdi
xor    %r12d,%r12d
movsd  0xec5(%rip),%xmm0
movsd  0xec5(%rip),%xmm2
xor    %r13d,%r13d
xor    %r11d,%r11d
nopl   0x0(%rax)
movsd  (%rcx,%r12,1),%xmm3
xorpd  %xmm4,%xmm4
cmplesd %xmm3,%xmm4
andpd  %xmm3,%xmm4
movsd  (%rsi,%r12,1),%xmm3
xorpd  %xmm5,%xmm5
cmplesd %xmm3,%xmm5
andpd  %xmm3,%xmm5
movapd %xmm5,%xmm3
mulsd  %xmm1,%xmm3
mulsd  %xmm5,%xmm3
addsd  %xmm0,%xmm4
divsd  %xmm4,%xmm3
movapd %xmm3,%xmm4
addsd  %xmm3,%xmm4
movsd  (%r8,%r13,1),%xmm5
addsd  %xmm4,%xmm5
movsd  %xmm5,(%r8,%r13,1)
movsd  0x20(%r8,%r13,1),%xmm5
addsd  %xmm3,%xmm5
movsd  %xmm5,0x20(%r8,%r13,1)
movsd  0x40(%r8,%r13,1),%xmm5
addsd  %xmm3,%xmm5
movsd  %xmm5,0x40(%r8,%r13,1)
movsd  0x18(%rbp,%r13,1),%xmm5
addsd  %xmm3,%xmm5
movsd  %xmm5,0x18(%rbp,%r13,1)
movsd  0x30(%rbx,%r13,1),%xmm5
addsd  %xmm3,%xmm5
movsd  %xmm5,0x30(%rbx,%r13,1)
movsd  (%r9,%r13,1),%xmm5
addsd  %xmm3,%xmm5
movsd  %xmm5,(%r9,%r13,1)
movsd  0x20(%r9,%r13,1),%xmm5
addsd  %xmm4,%xmm5
movsd  %xmm5,0x20(%r9,%r13,1)
movsd  0x40(%r9,%r13,1),%xmm5
addsd  %xmm3,%xmm5
movsd  %xmm5,0x40(%r9,%r13,1)
movsd  0x8(%r15,%r13,1),%xmm5
addsd  %xmm3,%xmm5
movsd  %xmm5,0x8(%r15,%r13,1)
mov    0x58(%rsp),%rax
movsd  0x38(%rax,%r13,1),%xmm5
addsd  %xmm3,%xmm5
movsd  %xmm5,0x38(%rax,%r13,1)
movsd  (%r10,%r13,1),%xmm5
addsd  %xmm3,%xmm5
movsd  %xmm5,(%r10,%r13,1)
movsd  0x20(%r10,%r13,1),%xmm5
addsd  %xmm3,%xmm5
movsd  %xmm5,0x20(%r10,%r13,1)
addsd  0x40(%r10,%r13,1),%xmm4
movsd  %xmm4,0x40(%r10,%r13,1)
mov    0x60(%rsp),%rax
movsd  0x10(%rax,%r13,1),%xmm4
addsd  %xmm3,%xmm4
movsd  %xmm4,0x10(%rax,%r13,1)
mov    0x68(%rsp),%rax
addsd  0x28(%rax,%r13,1),%xmm3
movsd  %xmm3,0x28(%rax,%r13,1)
mov    %r11d,%r14d
movsd  (%rdx,%r14,8),%xmm3
mulsd  %xmm2,%xmm3
mov    0x70(%rsp),%rax
addsd  (%rax,%r12,1),%xmm3
movsd  %xmm3,(%rax,%r12,1)
lea    0x1(%r11),%eax
movsd  (%rdx,%rax,8),%xmm3
mulsd  %xmm2,%xmm3
mov    0x78(%rsp),%rax
addsd  (%rax,%r12,1),%xmm3
movsd  %xmm3,(%rax,%r12,1)
lea    0x2(%r11),%eax
movsd  (%rdx,%rax,8),%xmm3
mulsd  %xmm2,%xmm3
mov    0x80(%rsp),%rax
addsd  (%rax,%r12,1),%xmm3
movsd  %xmm3,(%rax,%r12,1)
add    $0x3,%r11d
add    $0x48,%r13
add    $0x8,%r12
cmp    %r13,%rdi
jne    1150 <ReynoldsAveragedNavierStokes_kEpsilon_3D_Update+0x50>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @ReynoldsAveragedNavierStokes_kEpsilon_3D_Update(i32 %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, double* %13, double* %14, double* %15, double* %16, double* %17) {
%19 = alloca i32
%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 double*
%35 = alloca double*
%36 = alloca double*
%37 = alloca i32
%38 = alloca double
%39 = alloca double
%40 = alloca double
%41 = alloca double
%42 = alloca double
store i32 %0, i32* %19
store double %1, double* %20
store double %2, double* %21
store double* %3, double** %22
store double* %4, double** %23
store double* %5, double** %24
store double* %6, double** %25
store double* %7, double** %26
store double* %8, double** %27
store double* %9, double** %28
store double* %10, double** %29
store double* %11, double** %30
store double* %12, double** %31
store double* %13, double** %32
store double* %14, double** %33
store double* %15, double** %34
store double* %16, double** %35
store double* %17, double** %36
store double 0x3FE5555555555555, double* %41
store double 0x3EB0C6F7A0B5ED8D, double* %42
store i32 0, i32* %37
br label %43
43: ; preds = %260, %18
%44 = load i32, i32* %37
%45 = load i32, i32* %19
%46 = icmp slt i32 %44, %45
br i1 %46, label %47, label %263
47: ; preds = %43
%48 = load double*, double** %24
%49 = load i32, i32* %37
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds double, double* %48, i64 %50
%52 = load double, double* %51
%53 = fcmp oge double %52, 0.000000e+00
br i1 %53, label %54, label %60
54: ; preds = %47
%55 = load double*, double** %24
%56 = load i32, i32* %37
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds double, double* %55, i64 %57
%59 = load double, double* %58
br label %61
60: ; preds = %47
br label %61
61: ; preds = %60, %54
%62 = phi double [ %59, %54 ], [ 0.000000e+00, %60 ]
store double %62, double* %39
%63 = load double*, double** %22
%64 = load i32, i32* %37
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds double, double* %63, i64 %65
%67 = load double, double* %66
%68 = fcmp oge double %67, 0.000000e+00
br i1 %68, label %69, label %75
69: ; preds = %61
%70 = load double*, double** %22
%71 = load i32, i32* %37
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds double, double* %70, i64 %72
%74 = load double, double* %73
br label %76
75: ; preds = %61
br label %76
76: ; preds = %75, %69
%77 = phi double [ %74, %69 ], [ 0.000000e+00, %75 ]
store double %77, double* %40
%78 = load double, double* %21
%79 = load double, double* %40
%80 = fmul double %78, %79
%81 = load double, double* %40
%82 = fmul double %80, %81
%83 = load double, double* %39
%84 = fadd double %83, 0x3EB0C6F7A0B5ED8D
%85 = fdiv double %82, %84
store double %85, double* %38
%86 = load double, double* %38
%87 = load double*, double** %25
%88 = load i32, i32* %37
%89 = mul nsw i32 %88, 9
%90 = add nsw i32 %89, 0
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds double, double* %87, i64 %91
%93 = load double, double* %92
%94 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %86, double %93)
store double %94, double* %92
%95 = load double, double* %38
%96 = load double*, double** %25
%97 = load i32, i32* %37
%98 = mul nsw i32 %97, 9
%99 = add nsw i32 %98, 4
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds double, double* %96, i64 %100
%102 = load double, double* %101
%103 = fadd double %102, %95
store double %103, double* %101
%104 = load double, double* %38
%105 = load double*, double** %25
%106 = load i32, i32* %37
%107 = mul nsw i32 %106, 9
%108 = add nsw i32 %107, 8
%109 = sext i32 %108 to i64
%110 = getelementptr inbounds double, double* %105, i64 %109
%111 = load double, double* %110
%112 = fadd double %111, %104
store double %112, double* %110
%113 = load double, double* %38
%114 = load double*, double** %28
%115 = load i32, i32* %37
%116 = mul nsw i32 %115, 9
%117 = add nsw i32 %116, 3
%118 = sext i32 %117 to i64
%119 = getelementptr inbounds double, double* %114, i64 %118
%120 = load double, double* %119
%121 = fadd double %120, %113
store double %121, double* %119
%122 = load double, double* %38
%123 = load double*, double** %29
%124 = load i32, i32* %37
%125 = mul nsw i32 %124, 9
%126 = add nsw i32 %125, 6
%127 = sext i32 %126 to i64
%128 = getelementptr inbounds double, double* %123, i64 %127
%129 = load double, double* %128
%130 = fadd double %129, %122
store double %130, double* %128
%131 = load double, double* %38
%132 = load double*, double** %26
%133 = load i32, i32* %37
%134 = mul nsw i32 %133, 9
%135 = add nsw i32 %134, 0
%136 = sext i32 %135 to i64
%137 = getelementptr inbounds double, double* %132, i64 %136
%138 = load double, double* %137
%139 = fadd double %138, %131
store double %139, double* %137
%140 = load double, double* %38
%141 = load double*, double** %26
%142 = load i32, i32* %37
%143 = mul nsw i32 %142, 9
%144 = add nsw i32 %143, 4
%145 = sext i32 %144 to i64
%146 = getelementptr inbounds double, double* %141, i64 %145
%147 = load double, double* %146
%148 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %140, double %147)
store double %148, double* %146
%149 = load double, double* %38
%150 = load double*, double** %26
%151 = load i32, i32* %37
%152 = mul nsw i32 %151, 9
%153 = add nsw i32 %152, 8
%154 = sext i32 %153 to i64
%155 = getelementptr inbounds double, double* %150, i64 %154
%156 = load double, double* %155
%157 = fadd double %156, %149
store double %157, double* %155
%158 = load double, double* %38
%159 = load double*, double** %30
%160 = load i32, i32* %37
%161 = mul nsw i32 %160, 9
%162 = add nsw i32 %161, 1
%163 = sext i32 %162 to i64
%164 = getelementptr inbounds double, double* %159, i64 %163
%165 = load double, double* %164
%166 = fadd double %165, %158
store double %166, double* %164
%167 = load double, double* %38
%168 = load double*, double** %31
%169 = load i32, i32* %37
%170 = mul nsw i32 %169, 9
%171 = add nsw i32 %170, 7
%172 = sext i32 %171 to i64
%173 = getelementptr inbounds double, double* %168, i64 %172
%174 = load double, double* %173
%175 = fadd double %174, %167
store double %175, double* %173
%176 = load double, double* %38
%177 = load double*, double** %27
%178 = load i32, i32* %37
%179 = mul nsw i32 %178, 9
%180 = add nsw i32 %179, 0
%181 = sext i32 %180 to i64
%182 = getelementptr inbounds double, double* %177, i64 %181
%183 = load double, double* %182
%184 = fadd double %183, %176
store double %184, double* %182
%185 = load double, double* %38
%186 = load double*, double** %27
%187 = load i32, i32* %37
%188 = mul nsw i32 %187, 9
%189 = add nsw i32 %188, 4
%190 = sext i32 %189 to i64
%191 = getelementptr inbounds double, double* %186, i64 %190
%192 = load double, double* %191
%193 = fadd double %192, %185
store double %193, double* %191
%194 = load double, double* %38
%195 = load double*, double** %27
%196 = load i32, i32* %37
%197 = mul nsw i32 %196, 9
%198 = add nsw i32 %197, 8
%199 = sext i32 %198 to i64
%200 = getelementptr inbounds double, double* %195, i64 %199
%201 = load double, double* %200
%202 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %194, double %201)
store double %202, double* %200
%203 = load double, double* %38
%204 = load double*, double** %32
%205 = load i32, i32* %37
%206 = mul nsw i32 %205, 9
%207 = add nsw i32 %206, 2
%208 = sext i32 %207 to i64
%209 = getelementptr inbounds double, double* %204, i64 %208
%210 = load double, double* %209
%211 = fadd double %210, %203
store double %211, double* %209
%212 = load double, double* %38
%213 = load double*, double** %33
%214 = load i32, i32* %37
%215 = mul nsw i32 %214, 9
%216 = add nsw i32 %215, 5
%217 = sext i32 %216 to i64
%218 = getelementptr inbounds double, double* %213, i64 %217
%219 = load double, double* %218
%220 = fadd double %219, %212
store double %220, double* %218
%221 = load double*, double** %23
%222 = load i32, i32* %37
%223 = mul nsw i32 %222, 3
%224 = add nsw i32 %223, 0
%225 = sext i32 %224 to i64
%226 = getelementptr inbounds double, double* %221, i64 %225
%227 = load double, double* %226
%228 = load double*, double** %34
%229 = load i32, i32* %37
%230 = sext i32 %229 to i64
%231 = getelementptr inbounds double, double* %228, i64 %230
%232 = load double, double* %231
%233 = call double @llvm.fmuladd.f64(double 0x3FE5555555555555, double %227, double %232)
store double %233, double* %231
%234 = load double*, double** %23
%235 = load i32, i32* %37
%236 = mul nsw i32 %235, 3
%237 = add nsw i32 %236, 1
%238 = sext i32 %237 to i64
%239 = getelementptr inbounds double, double* %234, i64 %238
%240 = load double, double* %239
%241 = load double*, double** %35
%242 = load i32, i32* %37
%243 = sext i32 %242 to i64
%244 = getelementptr inbounds double, double* %241, i64 %243
%245 = load double, double* %244
%246 = call double @llvm.fmuladd.f64(double 0x3FE5555555555555, double %240, double %245)
store double %246, double* %244
%247 = load double*, double** %23
%248 = load i32, i32* %37
%249 = mul nsw i32 %248, 3
%250 = add nsw i32 %249, 2
%251 = sext i32 %250 to i64
%252 = getelementptr inbounds double, double* %247, i64 %251
%253 = load double, double* %252
%254 = load double*, double** %36
%255 = load i32, i32* %37
%256 = sext i32 %255 to i64
%257 = getelementptr inbounds double, double* %254, i64 %256
%258 = load double, double* %257
%259 = call double @llvm.fmuladd.f64(double 0x3FE5555555555555, double %253, double %258)
store double %259, double* %257
br label %260
260: ; preds = %76
%261 = load i32, i32* %37
%262 = add nsw i32 %261, 1
store i32 %262, i32* %37
br label %43
263: ; preds = %43
ret void
}

LLVM IR (-O2):
define void @ReynoldsAveragedNavierStokes_kEpsilon_3D_Update(i32 %0, double %1, double %2, double* nocapture %3, double* nocapture %4, double* nocapture %5, double* nocapture %6, double* nocapture %7, double* nocapture %8, double* nocapture %9, double* nocapture %10, double* nocapture %11, double* nocapture %12, double* nocapture %13, double* nocapture %14, double* nocapture %15, double* nocapture %16, double* nocapture %17) {
%19 = icmp sgt i32 %0, 0
br i1 %19, label %20, label %114
20: ; preds = %18
%21 = zext i32 %0 to i64
br label %22
22: ; preds = %20, %22
%23 = phi i64 [ 0, %20 ], [ %112, %22 ]
%24 = getelementptr inbounds double, double* %5, i64 %23
%25 = load double, double* %24
%26 = fcmp oge double %25, 0.000000e+00
%27 = select i1 %26, double %25, double 0.000000e+00
%28 = getelementptr inbounds double, double* %3, i64 %23
%29 = load double, double* %28
%30 = fcmp oge double %29, 0.000000e+00
%31 = select i1 %30, double %29, double 0.000000e+00
%32 = fmul double %31, %2
%33 = fmul double %31, %32
%34 = fadd double %27, 0x3EB0C6F7A0B5ED8D
%35 = fdiv double %33, %34
%36 = mul nuw nsw i64 %23, 9
%37 = getelementptr inbounds double, double* %6, i64 %36
%38 = load double, double* %37
%39 = tail call double @llvm.fmuladd.f64(double %35, double 2.000000e+00, double %38)
store double %39, double* %37
%40 = add nuw nsw i64 %36, 4
%41 = getelementptr inbounds double, double* %6, i64 %40
%42 = load double, double* %41
%43 = fadd double %42, %35
store double %43, double* %41
%44 = add nuw nsw i64 %36, 8
%45 = getelementptr inbounds double, double* %6, i64 %44
%46 = load double, double* %45
%47 = fadd double %46, %35
store double %47, double* %45
%48 = add nuw nsw i64 %36, 3
%49 = getelementptr inbounds double, double* %9, i64 %48
%50 = load double, double* %49
%51 = fadd double %35, %50
store double %51, double* %49
%52 = add nuw nsw i64 %36, 6
%53 = getelementptr inbounds double, double* %10, i64 %52
%54 = load double, double* %53
%55 = fadd double %35, %54
store double %55, double* %53
%56 = getelementptr inbounds double, double* %7, i64 %36
%57 = load double, double* %56
%58 = fadd double %35, %57
store double %58, double* %56
%59 = getelementptr inbounds double, double* %7, i64 %40
%60 = load double, double* %59
%61 = tail call double @llvm.fmuladd.f64(double %35, double 2.000000e+00, double %60)
store double %61, double* %59
%62 = getelementptr inbounds double, double* %7, i64 %44
%63 = load double, double* %62
%64 = fadd double %35, %63
store double %64, double* %62
%65 = add nuw nsw i64 %36, 1
%66 = getelementptr inbounds double, double* %11, i64 %65
%67 = load double, double* %66
%68 = fadd double %35, %67
store double %68, double* %66
%69 = add nuw nsw i64 %36, 7
%70 = getelementptr inbounds double, double* %12, i64 %69
%71 = load double, double* %70
%72 = fadd double %35, %71
store double %72, double* %70
%73 = getelementptr inbounds double, double* %8, i64 %36
%74 = load double, double* %73
%75 = fadd double %35, %74
store double %75, double* %73
%76 = getelementptr inbounds double, double* %8, i64 %40
%77 = load double, double* %76
%78 = fadd double %35, %77
store double %78, double* %76
%79 = getelementptr inbounds double, double* %8, i64 %44
%80 = load double, double* %79
%81 = tail call double @llvm.fmuladd.f64(double %35, double 2.000000e+00, double %80)
store double %81, double* %79
%82 = add nuw nsw i64 %36, 2
%83 = getelementptr inbounds double, double* %13, i64 %82
%84 = load double, double* %83
%85 = fadd double %35, %84
store double %85, double* %83
%86 = add nuw nsw i64 %36, 5
%87 = getelementptr inbounds double, double* %14, i64 %86
%88 = load double, double* %87
%89 = fadd double %35, %88
store double %89, double* %87
%90 = trunc i64 %23 to i32
%91 = mul nsw i32 %90, 3
%92 = zext i32 %91 to i64
%93 = getelementptr inbounds double, double* %4, i64 %92
%94 = load double, double* %93
%95 = getelementptr inbounds double, double* %15, i64 %23
%96 = load double, double* %95
%97 = tail call double @llvm.fmuladd.f64(double %94, double 0x3FE5555555555555, double %96)
store double %97, double* %95
%98 = add nuw nsw i32 %91, 1
%99 = zext i32 %98 to i64
%100 = getelementptr inbounds double, double* %4, i64 %99
%101 = load double, double* %100
%102 = getelementptr inbounds double, double* %16, i64 %23
%103 = load double, double* %102
%104 = tail call double @llvm.fmuladd.f64(double %101, double 0x3FE5555555555555, double %103)
store double %104, double* %102
%105 = add nuw nsw i32 %91, 2
%106 = zext i32 %105 to i64
%107 = getelementptr inbounds double, double* %4, i64 %106
%108 = load double, double* %107
%109 = getelementptr inbounds double, double* %17, i64 %23
%110 = load double, double* %109
%111 = tail call double @llvm.fmuladd.f64(double %108, double 0x3FE5555555555555, double %110)
store double %111, double* %109
%112 = add nuw nsw i64 %23, 1
%113 = icmp eq i64 %112, %21
br i1 %113, label %114, label %22
114: ; preds = %22, %18
ret void
}