void Periodic_YZ_flux(double *phi_fv, double *u_buff, int nX, int nY, int nZ,
                      int X_low, int X_high) {
  int iY;
  int iZ;
  double u_face;
  for (iY = 0; iY < nY; iY++)
    for (iZ = 0; iZ < nZ; iZ++) {
      u_face =
          0.5 *
          (*(u_buff + X_low + iY * nX + iZ * nX * nY + (0) * nX * nY * nZ) +
           *(u_buff + (X_low + 1) + iY * nX + iZ * nX * nY +
             (0) * nX * nY * nZ));
      if (u_face > 0)
        *(phi_fv + iZ * nX * nY + iY * nX + X_low + (0) * nX * nY * nZ) =
            *(phi_fv + iZ * nX * nY + iY * nX + X_high + (0) * nX * nY * nZ);
      else
        *(phi_fv + iZ * nX * nY + iY * nX + X_high + (0) * nX * nY * nZ) =
            *(phi_fv + iZ * nX * nY + iY * nX + X_low + (0) * nX * nY * nZ);
    }
}


x86 assembly (clang, -O0):
<Periodic_YZ_flux>:
push   %rbp
mov    %rsp,%rbp
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
mov    %r9d,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    134c <Periodic_YZ_flux+0x24c>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1339 <Periodic_YZ_flux+0x239>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    -0x24(%rbp),%ecx
imul   -0x14(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    -0x28(%rbp),%ecx
imul   -0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
imul   $0x0,-0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    -0x24(%rbp),%ecx
imul   -0x14(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    -0x28(%rbp),%ecx
imul   -0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
imul   $0x0,-0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm1
movsd  0xe1e(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1294 <Periodic_YZ_flux+0x194>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    -0x24(%rbp),%ecx
imul   -0x14(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq 0x10(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
imul   $0x0,-0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    -0x24(%rbp),%ecx
imul   -0x14(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
imul   $0x0,-0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1326 <Periodic_YZ_flux+0x226>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    -0x24(%rbp),%ecx
imul   -0x14(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
imul   $0x0,-0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    -0x24(%rbp),%ecx
imul   -0x14(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq 0x10(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
imul   $0x0,-0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    132b <Periodic_YZ_flux+0x22b>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1137 <Periodic_YZ_flux+0x37>
jmp    133e <Periodic_YZ_flux+0x23e>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1124 <Periodic_YZ_flux+0x24>
pop    %rbp
ret

x86 assembly (clang, -O2):
<Periodic_YZ_flux>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %ecx,%ecx
jle    11ae <Periodic_YZ_flux+0xae>
movslq %r9d,%r12
lea    (%rsi,%r12,8),%r9
movslq 0x38(%rsp),%r13
movslq %edx,%r10
mov    %ecx,%r14d
mov    %r8d,%r11d
mov    %r10d,%ebx
imul   %ecx,%ebx
xor    %r15d,%r15d
movsd  0xec8(%rip),%xmm0
xorpd  %xmm1,%xmm1
jmp    1149 <Periodic_YZ_flux+0x49>
xchg   %ax,%ax
add    $0x1,%r15
cmp    %r14,%r15
je     11ae <Periodic_YZ_flux+0xae>
test   %r8d,%r8d
jle    1140 <Periodic_YZ_flux+0x40>
mov    %r15,%rax
imul   %r10,%rax
lea    (%r9,%rax,8),%rdx
lea    (%rdi,%rax,8),%rcx
xor    %ebp,%ebp
mov    %r11,%rax
jmp    1184 <Periodic_YZ_flux+0x84>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rsi,%r12,8),%xmm2
movsd  %xmm2,(%rsi,%r13,8)
add    %ebx,%ebp
add    $0xffffffffffffffff,%rax
je     1140 <Periodic_YZ_flux+0x40>
movslq %ebp,%rsi
movsd  (%rdx,%rsi,8),%xmm2
addsd  0x8(%rdx,%rsi,8),%xmm2
mulsd  %xmm0,%xmm2
lea    (%rcx,%rsi,8),%rsi
ucomisd %xmm1,%xmm2
jbe    1170 <Periodic_YZ_flux+0x70>
movsd  (%rsi,%r13,8),%xmm2
movsd  %xmm2,(%rsi,%r12,8)
jmp    117c <Periodic_YZ_flux+0x7c>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @Periodic_YZ_flux(double* nocapture %0, double* nocapture %1, i32 %2, i32 %3, i32 %4, i32 %5, i32 %6) {
%8 = icmp sgt i32 %3, 0
br i1 %8, label %9, label %57
9: ; preds = %7
%10 = icmp sgt i32 %4, 0
%11 = sext i32 %5 to i64
%12 = getelementptr inbounds double, double* %1, i64 %11
%13 = mul i32 %3, %2
%14 = add nsw i32 %5, 1
%15 = sext i32 %14 to i64
%16 = getelementptr inbounds double, double* %1, i64 %15
%17 = sext i32 %6 to i64
%18 = sext i32 %2 to i64
%19 = zext i32 %3 to i64
%20 = zext i32 %4 to i64
br label %21
21: ; preds = %9, %54
%22 = phi i64 [ 0, %9 ], [ %55, %54 ]
br i1 %10, label %23, label %54
23: ; preds = %21
%24 = mul nsw i64 %22, %18
%25 = getelementptr inbounds double, double* %12, i64 %24
%26 = getelementptr inbounds double, double* %16, i64 %24
%27 = getelementptr inbounds double, double* %0, i64 %24
%28 = getelementptr inbounds double, double* %0, i64 %24
br label %29
29: ; preds = %23, %51
%30 = phi i64 [ 0, %23 ], [ %52, %51 ]
%31 = trunc i64 %30 to i32
%32 = mul i32 %13, %31
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds double, double* %25, i64 %33
%35 = load double, double* %34
%36 = getelementptr inbounds double, double* %26, i64 %33
%37 = load double, double* %36
%38 = fadd double %35, %37
%39 = fmul double %38, 5.000000e-01
%40 = fcmp ogt double %39, 0.000000e+00
br i1 %40, label %41, label %46
41: ; preds = %29
%42 = getelementptr inbounds double, double* %28, i64 %33
%43 = getelementptr inbounds double, double* %42, i64 %17
%44 = load double, double* %43
%45 = getelementptr inbounds double, double* %42, i64 %11
store double %44, double* %45
br label %51
46: ; preds = %29
%47 = getelementptr inbounds double, double* %27, i64 %33
%48 = getelementptr inbounds double, double* %47, i64 %11
%49 = load double, double* %48
%50 = getelementptr inbounds double, double* %47, i64 %17
store double %49, double* %50
br label %51
51: ; preds = %41, %46
%52 = add nuw nsw i64 %30, 1
%53 = icmp eq i64 %52, %20
br i1 %53, label %54, label %29
54: ; preds = %51, %21
%55 = add nuw nsw i64 %22, 1
%56 = icmp eq i64 %55, %19
br i1 %56, label %57, label %21
57: ; preds = %54, %7
ret void
}