void ReynoldsAveragedNavierStokes_kEpsilon_3D_Update_sd(
    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 * 3 + 0] += 2.0 * nu_t;
    mom_u_diff_ten[i * 3 + 1] += nu_t;
    mom_u_diff_ten[i * 3 + 2] += nu_t;
    mom_uv_diff_ten[i] += nu_t;
    mom_uw_diff_ten[i] += nu_t;
    mom_v_diff_ten[i * 3 + 0] += nu_t;
    mom_v_diff_ten[i * 3 + 1] += 2.0 * nu_t;
    mom_v_diff_ten[i * 3 + 2] += nu_t;
    mom_vu_diff_ten[i] += nu_t;
    mom_vw_diff_ten[i] += nu_t;
    mom_w_diff_ten[i * 3 + 0] += nu_t;
    mom_w_diff_ten[i * 3 + 1] += nu_t;
    mom_w_diff_ten[i * 3 + 2] += 2.0 * nu_t;
    mom_wu_diff_ten[i] += nu_t;
    mom_wv_diff_ten[i] += 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_sd>:
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    1483 <ReynoldsAveragedNavierStokes_kEpsilon_3D_Update_sd+0x383>
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_sd+0xab>
mov    -0x30(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    11b8 <ReynoldsAveragedNavierStokes_kEpsilon_3D_Update_sd+0xb8>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    11b8 <ReynoldsAveragedNavierStokes_kEpsilon_3D_Update_sd+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_sd+0xf3>
mov    -0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
jmp    1200 <ReynoldsAveragedNavierStokes_kEpsilon_3D_Update_sd+0x100>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x80(%rbp)
jmp    1200 <ReynoldsAveragedNavierStokes_kEpsilon_3D_Update_sd+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   $0x3,-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   $0x3,-0x44(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    -0x38(%rbp),%rax
imul   $0x3,-0x44(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    -0x40(%rbp),%rax
imul   $0x3,-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   $0x3,-0x44(%rbp),%ecx
add    $0x1,%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   $0x3,-0x44(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    0x28(%rbp),%rax
movslq -0x44(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    0x30(%rbp),%rax
movslq -0x44(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    0x10(%rbp),%rax
imul   $0x3,-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   $0x3,-0x44(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    0x10(%rbp),%rax
imul   $0x3,-0x44(%rbp),%ecx
add    $0x2,%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
movslq -0x44(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    0x40(%rbp),%rax
movslq -0x44(%rbp),%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  0xc0a(%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  0xbd5(%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  0xba0(%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_sd+0x6e>
pop    %rbp
ret

x86 assembly (clang, -O2):
<ReynoldsAveragedNavierStokes_kEpsilon_3D_Update_sd>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rcx,-0x8(%rsp)
test   %edi,%edi
jle    12dd <ReynoldsAveragedNavierStokes_kEpsilon_3D_Update_sd+0x1dd>
mov    0x60(%rsp),%r12
mov    0x58(%rsp),%r13
mov    0x50(%rsp),%r15
mov    0x48(%rsp),%rbx
mov    0x40(%rsp),%rbp
mov    0x38(%rsp),%r10
mov    %edi,%edi
shl    $0x3,%rdi
xor    %ecx,%ecx
movsd  0xebb(%rip),%xmm0
movsd  0xebb(%rip),%xmm2
xor    %r14d,%r14d
mov    -0x8(%rsp),%r11
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%r11,%rcx,1),%xmm3
xorpd  %xmm4,%xmm4
cmplesd %xmm3,%xmm4
andpd  %xmm3,%xmm4
movsd  (%rsi,%rcx,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,%r14,1),%xmm5
addsd  %xmm4,%xmm5
movsd  %xmm5,(%r8,%r14,1)
movupd 0x8(%r8,%r14,1),%xmm6
movapd %xmm3,%xmm5
unpcklpd %xmm3,%xmm5
addpd  %xmm5,%xmm6
movupd %xmm6,0x8(%r8,%r14,1)
movsd  0x0(%rbp,%rcx,1),%xmm6
addsd  %xmm3,%xmm6
movsd  %xmm6,0x0(%rbp,%rcx,1)
movsd  (%rbx,%rcx,1),%xmm6
addsd  %xmm3,%xmm6
movsd  %xmm6,(%rbx,%rcx,1)
movsd  (%r9,%r14,1),%xmm6
addsd  %xmm3,%xmm6
movsd  %xmm6,(%r9,%r14,1)
movsd  0x8(%r9,%r14,1),%xmm6
addsd  %xmm4,%xmm6
movsd  %xmm6,0x8(%r9,%r14,1)
movsd  0x10(%r9,%r14,1),%xmm6
addsd  %xmm3,%xmm6
movsd  %xmm6,0x10(%r9,%r14,1)
movsd  (%r15,%rcx,1),%xmm6
addsd  %xmm3,%xmm6
movsd  %xmm6,(%r15,%rcx,1)
movsd  0x0(%r13,%rcx,1),%xmm6
addsd  %xmm3,%xmm6
movsd  %xmm6,0x0(%r13,%rcx,1)
movupd (%r10,%r14,1),%xmm6
addpd  %xmm5,%xmm6
movupd %xmm6,(%r10,%r14,1)
addsd  0x10(%r10,%r14,1),%xmm4
movsd  %xmm4,0x10(%r10,%r14,1)
movsd  (%r12,%rcx,1),%xmm4
addsd  %xmm3,%xmm4
movsd  %xmm4,(%r12,%rcx,1)
mov    0x68(%rsp),%rax
addsd  (%rax,%rcx,1),%xmm3
movsd  %xmm3,(%rax,%rcx,1)
movsd  (%rdx,%r14,1),%xmm3
mulsd  %xmm2,%xmm3
mov    0x70(%rsp),%rax
addsd  (%rax,%rcx,1),%xmm3
movsd  %xmm3,(%rax,%rcx,1)
movsd  0x8(%rdx,%r14,1),%xmm3
mulsd  %xmm2,%xmm3
mov    0x78(%rsp),%rax
addsd  (%rax,%rcx,1),%xmm3
movsd  %xmm3,(%rax,%rcx,1)
movsd  0x10(%rdx,%r14,1),%xmm3
mulsd  %xmm2,%xmm3
mov    0x80(%rsp),%rax
addsd  (%rax,%rcx,1),%xmm3
movsd  %xmm3,(%rax,%rcx,1)
add    $0x18,%r14
add    $0x8,%rcx
cmp    %rcx,%rdi
jne    1160 <ReynoldsAveragedNavierStokes_kEpsilon_3D_Update_sd+0x60>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @ReynoldsAveragedNavierStokes_kEpsilon_3D_Update_sd(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 = %248, %18
%44 = load i32, i32* %37
%45 = load i32, i32* %19
%46 = icmp slt i32 %44, %45
br i1 %46, label %47, label %251
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, 3
%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, 3
%99 = add nsw i32 %98, 1
%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, 3
%108 = add nsw i32 %107, 2
%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 = sext i32 %115 to i64
%117 = getelementptr inbounds double, double* %114, i64 %116
%118 = load double, double* %117
%119 = fadd double %118, %113
store double %119, double* %117
%120 = load double, double* %38
%121 = load double*, double** %29
%122 = load i32, i32* %37
%123 = sext i32 %122 to i64
%124 = getelementptr inbounds double, double* %121, i64 %123
%125 = load double, double* %124
%126 = fadd double %125, %120
store double %126, double* %124
%127 = load double, double* %38
%128 = load double*, double** %26
%129 = load i32, i32* %37
%130 = mul nsw i32 %129, 3
%131 = add nsw i32 %130, 0
%132 = sext i32 %131 to i64
%133 = getelementptr inbounds double, double* %128, i64 %132
%134 = load double, double* %133
%135 = fadd double %134, %127
store double %135, double* %133
%136 = load double, double* %38
%137 = load double*, double** %26
%138 = load i32, i32* %37
%139 = mul nsw i32 %138, 3
%140 = add nsw i32 %139, 1
%141 = sext i32 %140 to i64
%142 = getelementptr inbounds double, double* %137, i64 %141
%143 = load double, double* %142
%144 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %136, double %143)
store double %144, double* %142
%145 = load double, double* %38
%146 = load double*, double** %26
%147 = load i32, i32* %37
%148 = mul nsw i32 %147, 3
%149 = add nsw i32 %148, 2
%150 = sext i32 %149 to i64
%151 = getelementptr inbounds double, double* %146, i64 %150
%152 = load double, double* %151
%153 = fadd double %152, %145
store double %153, double* %151
%154 = load double, double* %38
%155 = load double*, double** %30
%156 = load i32, i32* %37
%157 = sext i32 %156 to i64
%158 = getelementptr inbounds double, double* %155, i64 %157
%159 = load double, double* %158
%160 = fadd double %159, %154
store double %160, double* %158
%161 = load double, double* %38
%162 = load double*, double** %31
%163 = load i32, i32* %37
%164 = sext i32 %163 to i64
%165 = getelementptr inbounds double, double* %162, i64 %164
%166 = load double, double* %165
%167 = fadd double %166, %161
store double %167, double* %165
%168 = load double, double* %38
%169 = load double*, double** %27
%170 = load i32, i32* %37
%171 = mul nsw i32 %170, 3
%172 = add nsw i32 %171, 0
%173 = sext i32 %172 to i64
%174 = getelementptr inbounds double, double* %169, i64 %173
%175 = load double, double* %174
%176 = fadd double %175, %168
store double %176, double* %174
%177 = load double, double* %38
%178 = load double*, double** %27
%179 = load i32, i32* %37
%180 = mul nsw i32 %179, 3
%181 = add nsw i32 %180, 1
%182 = sext i32 %181 to i64
%183 = getelementptr inbounds double, double* %178, i64 %182
%184 = load double, double* %183
%185 = fadd double %184, %177
store double %185, double* %183
%186 = load double, double* %38
%187 = load double*, double** %27
%188 = load i32, i32* %37
%189 = mul nsw i32 %188, 3
%190 = add nsw i32 %189, 2
%191 = sext i32 %190 to i64
%192 = getelementptr inbounds double, double* %187, i64 %191
%193 = load double, double* %192
%194 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %186, double %193)
store double %194, double* %192
%195 = load double, double* %38
%196 = load double*, double** %32
%197 = load i32, i32* %37
%198 = sext i32 %197 to i64
%199 = getelementptr inbounds double, double* %196, i64 %198
%200 = load double, double* %199
%201 = fadd double %200, %195
store double %201, double* %199
%202 = load double, double* %38
%203 = load double*, double** %33
%204 = load i32, i32* %37
%205 = sext i32 %204 to i64
%206 = getelementptr inbounds double, double* %203, i64 %205
%207 = load double, double* %206
%208 = fadd double %207, %202
store double %208, double* %206
%209 = load double*, double** %23
%210 = load i32, i32* %37
%211 = mul nsw i32 %210, 3
%212 = add nsw i32 %211, 0
%213 = sext i32 %212 to i64
%214 = getelementptr inbounds double, double* %209, i64 %213
%215 = load double, double* %214
%216 = load double*, double** %34
%217 = load i32, i32* %37
%218 = sext i32 %217 to i64
%219 = getelementptr inbounds double, double* %216, i64 %218
%220 = load double, double* %219
%221 = call double @llvm.fmuladd.f64(double 0x3FE5555555555555, double %215, double %220)
store double %221, double* %219
%222 = load double*, double** %23
%223 = load i32, i32* %37
%224 = mul nsw i32 %223, 3
%225 = add nsw i32 %224, 1
%226 = sext i32 %225 to i64
%227 = getelementptr inbounds double, double* %222, i64 %226
%228 = load double, double* %227
%229 = load double*, double** %35
%230 = load i32, i32* %37
%231 = sext i32 %230 to i64
%232 = getelementptr inbounds double, double* %229, i64 %231
%233 = load double, double* %232
%234 = call double @llvm.fmuladd.f64(double 0x3FE5555555555555, double %228, double %233)
store double %234, double* %232
%235 = load double*, double** %23
%236 = load i32, i32* %37
%237 = mul nsw i32 %236, 3
%238 = add nsw i32 %237, 2
%239 = sext i32 %238 to i64
%240 = getelementptr inbounds double, double* %235, i64 %239
%241 = load double, double* %240
%242 = load double*, double** %36
%243 = load i32, i32* %37
%244 = sext i32 %243 to i64
%245 = getelementptr inbounds double, double* %242, i64 %244
%246 = load double, double* %245
%247 = call double @llvm.fmuladd.f64(double 0x3FE5555555555555, double %241, double %246)
store double %247, double* %245
br label %248
248: ; preds = %76
%249 = load i32, i32* %37
%250 = add nsw i32 %249, 1
store i32 %250, i32* %37
br label %43
251: ; preds = %43
ret void
}

LLVM IR (-O2):
define void @ReynoldsAveragedNavierStokes_kEpsilon_3D_Update_sd(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 %101
20: ; preds = %18
%21 = zext i32 %0 to i64
br label %22
22: ; preds = %20, %22
%23 = phi i64 [ 0, %20 ], [ %99, %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, 3
%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, 1
%41 = getelementptr inbounds double, double* %6, i64 %40
%42 = add nuw nsw i64 %36, 2
%43 = bitcast double* %41 to <2 x double>*
%44 = load <2 x double>, <2 x double>* %43
%45 = insertelement <2 x double> poison, double %35, i64 0
%46 = shufflevector <2 x double> %45, <2 x double> poison, <2 x i32> zeroinitializer
%47 = fadd <2 x double> %44, %46
%48 = bitcast double* %41 to <2 x double>*
store <2 x double> %47, <2 x double>* %48
%49 = getelementptr inbounds double, double* %9, i64 %23
%50 = load double, double* %49
%51 = fadd double %35, %50
store double %51, double* %49
%52 = getelementptr inbounds double, double* %10, i64 %23
%53 = load double, double* %52
%54 = fadd double %35, %53
store double %54, double* %52
%55 = getelementptr inbounds double, double* %7, i64 %36
%56 = load double, double* %55
%57 = fadd double %35, %56
store double %57, double* %55
%58 = getelementptr inbounds double, double* %7, i64 %40
%59 = load double, double* %58
%60 = tail call double @llvm.fmuladd.f64(double %35, double 2.000000e+00, double %59)
store double %60, double* %58
%61 = getelementptr inbounds double, double* %7, i64 %42
%62 = load double, double* %61
%63 = fadd double %35, %62
store double %63, double* %61
%64 = getelementptr inbounds double, double* %11, i64 %23
%65 = load double, double* %64
%66 = fadd double %35, %65
store double %66, double* %64
%67 = getelementptr inbounds double, double* %12, i64 %23
%68 = load double, double* %67
%69 = fadd double %35, %68
store double %69, double* %67
%70 = getelementptr inbounds double, double* %8, i64 %36
%71 = bitcast double* %70 to <2 x double>*
%72 = load <2 x double>, <2 x double>* %71
%73 = fadd <2 x double> %46, %72
%74 = bitcast double* %70 to <2 x double>*
store <2 x double> %73, <2 x double>* %74
%75 = getelementptr inbounds double, double* %8, i64 %42
%76 = load double, double* %75
%77 = tail call double @llvm.fmuladd.f64(double %35, double 2.000000e+00, double %76)
store double %77, double* %75
%78 = getelementptr inbounds double, double* %13, i64 %23
%79 = load double, double* %78
%80 = fadd double %35, %79
store double %80, double* %78
%81 = getelementptr inbounds double, double* %14, i64 %23
%82 = load double, double* %81
%83 = fadd double %35, %82
store double %83, double* %81
%84 = getelementptr inbounds double, double* %4, i64 %36
%85 = load double, double* %84
%86 = getelementptr inbounds double, double* %15, i64 %23
%87 = load double, double* %86
%88 = tail call double @llvm.fmuladd.f64(double %85, double 0x3FE5555555555555, double %87)
store double %88, double* %86
%89 = getelementptr inbounds double, double* %4, i64 %40
%90 = load double, double* %89
%91 = getelementptr inbounds double, double* %16, i64 %23
%92 = load double, double* %91
%93 = tail call double @llvm.fmuladd.f64(double %90, double 0x3FE5555555555555, double %92)
store double %93, double* %91
%94 = getelementptr inbounds double, double* %4, i64 %42
%95 = load double, double* %94
%96 = getelementptr inbounds double, double* %17, i64 %23
%97 = load double, double* %96
%98 = tail call double @llvm.fmuladd.f64(double %95, double 0x3FE5555555555555, double %97)
store double %98, double* %96
%99 = add nuw nsw i64 %23, 1
%100 = icmp eq i64 %99, %21
br i1 %100, label %101, label %22
101: ; preds = %22, %18
ret void
}