void kernel(double *g1, double *g2) {
  int x, y, z;
  int halo_width = 1;
  for (z = halo_width; z < 32 - halo_width; ++z) {
    for (y = halo_width; y < 32 - halo_width; ++y) {
      for (x = halo_width; x < 32 - halo_width; ++x) {
        double v = g1[((x) + (y)*32 + (z)*32 * 32)] +
                   g1[((x + 1) + (y)*32 + (z)*32 * 32)] +
                   g1[((x - 1) + (y)*32 + (z)*32 * 32)] +
                   g1[((x) + (y + 1) * 32 + (z)*32 * 32)] +
                   g1[((x) + (y - 1) * 32 + (z)*32 * 32)] +
                   g1[((x) + (y)*32 + (z - 1) * 32 * 32)] +
                   g1[((x) + (y)*32 + (z + 1) * 32 * 32)];
        g2[((x) + (y)*32 + (z)*32 * 32)] = v;
      }
    }
  }
  return;
}


x86 assembly (clang, -O0):
<kernel>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    $0x20,%ecx
sub    -0x20(%rbp),%ecx
cmp    %ecx,%eax
jge    12be <kernel+0x1be>
mov    -0x20(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
mov    $0x20,%ecx
sub    -0x20(%rbp),%ecx
cmp    %ecx,%eax
jge    12ab <kernel+0x1ab>
mov    -0x20(%rbp),%eax
mov    %eax,-0x14(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x20,%ecx
sub    -0x20(%rbp),%ecx
cmp    %ecx,%eax
jge    1298 <kernel+0x198>
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
mov    -0x18(%rbp),%edx
shl    $0x5,%edx
add    %edx,%ecx
mov    -0x1c(%rbp),%edx
shl    $0x5,%edx
shl    $0x5,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
mov    -0x18(%rbp),%edx
shl    $0x5,%edx
add    %edx,%ecx
mov    -0x1c(%rbp),%edx
shl    $0x5,%edx
shl    $0x5,%edx
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x18(%rbp),%edx
shl    $0x5,%edx
add    %edx,%ecx
mov    -0x1c(%rbp),%edx
shl    $0x5,%edx
shl    $0x5,%edx
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
mov    -0x18(%rbp),%edx
add    $0x1,%edx
shl    $0x5,%edx
add    %edx,%ecx
mov    -0x1c(%rbp),%edx
shl    $0x5,%edx
shl    $0x5,%edx
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
mov    -0x18(%rbp),%edx
sub    $0x1,%edx
shl    $0x5,%edx
add    %edx,%ecx
mov    -0x1c(%rbp),%edx
shl    $0x5,%edx
shl    $0x5,%edx
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
mov    -0x18(%rbp),%edx
shl    $0x5,%edx
add    %edx,%ecx
mov    -0x1c(%rbp),%edx
sub    $0x1,%edx
shl    $0x5,%edx
shl    $0x5,%edx
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
mov    -0x18(%rbp),%edx
shl    $0x5,%edx
add    %edx,%ecx
mov    -0x1c(%rbp),%edx
add    $0x1,%edx
shl    $0x5,%edx
shl    $0x5,%edx
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
mov    -0x18(%rbp),%edx
shl    $0x5,%edx
add    %edx,%ecx
mov    -0x1c(%rbp),%edx
shl    $0x5,%edx
shl    $0x5,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    114b <kernel+0x4b>
jmp    129d <kernel+0x19d>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1132 <kernel+0x32>
jmp    12b0 <kernel+0x1b0>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1119 <kernel+0x19>
pop    %rbp
ret

x86 assembly (clang, -O2):
<kernel>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
lea    0x4108(%rdi),%rbx
lea    0x2108(%rsi),%rdx
mov    $0x1,%eax
xor    %r11d,%r11d
jmp    115f <kernel+0x5f>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    -0x18(%rsp),%rax
add    $0x1,%rax
add    $0x1,%r11
mov    -0x8(%rsp),%rbx
add    $0x2000,%rbx
mov    -0x10(%rsp),%rdx
add    $0x2000,%rdx
cmp    $0x1f,%rax
je     12c7 <kernel+0x1c7>
mov    %rax,-0x18(%rsp)
mov    %r11,%r13
shl    $0xa,%r13
lea    0x421(%r13),%r12
lea    0x43f(%r13),%r14
lea    0x21(%r13),%r15
add    $0x83f,%r13
mov    %rdx,-0x10(%rsp)
mov    %rbx,-0x8(%rsp)
xor    %r10d,%r10d
mov    $0x1,%r9d
jmp    11c0 <kernel+0xc0>
nopl   0x0(%rax)
add    $0x1,%r9
add    $0x1,%r10
add    $0x100,%rbx
add    $0x100,%rdx
cmp    $0x1f,%r9
je     1130 <kernel+0x30>
mov    %r10,%rbp
shl    $0x5,%rbp
mov    %r13,%r8
add    %rbp,%r8
lea    (%rdi,%r8,8),%rax
lea    (%r12,%rbp,1),%rcx
lea    (%rsi,%rcx,8),%rcx
cmp    %rax,%rcx
jae    1250 <kernel+0x150>
lea    (%r15,%rbp,1),%rax
lea    (%rdi,%rax,8),%rax
add    %r14,%rbp
lea    (%rsi,%rbp,8),%rcx
cmp    %rcx,%rax
jae    1250 <kernel+0x150>
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  -0x2000(%rbx,%rbp,8),%xmm0
ff ff
addsd  -0x1ff8(%rbx,%rbp,8),%xmm0
ff ff
addsd  -0x2008(%rbx,%rbp,8),%xmm0
ff ff
addsd  -0x1f00(%rbx,%rbp,8),%xmm0
ff ff
addsd  -0x2100(%rbx,%rbp,8),%xmm0
ff ff
addsd  -0x4000(%rbx,%rbp,8),%xmm0
ff ff
addsd  (%rbx,%rbp,8),%xmm0
movsd  %xmm0,(%rdx,%rbp,8)
add    $0x1,%rbp
cmp    $0x1e,%rbp
jne    1200 <kernel+0x100>
jmp    11a0 <kernel+0xa0>
nop
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movupd -0x2000(%rbx,%rbp,8),%xmm0
ff ff
movupd -0x1ff8(%rbx,%rbp,8),%xmm1
ff ff
addpd  %xmm0,%xmm1
movupd -0x2008(%rbx,%rbp,8),%xmm0
ff ff
addpd  %xmm1,%xmm0
movupd -0x1f00(%rbx,%rbp,8),%xmm1
ff ff
addpd  %xmm0,%xmm1
movupd -0x2100(%rbx,%rbp,8),%xmm0
ff ff
addpd  %xmm1,%xmm0
movupd -0x4000(%rbx,%rbp,8),%xmm1
ff ff
addpd  %xmm0,%xmm1
movupd (%rbx,%rbp,8),%xmm0
addpd  %xmm1,%xmm0
movupd %xmm0,(%rdx,%rbp,8)
add    $0x2,%rbp
cmp    $0x1e,%rbp
jne    1260 <kernel+0x160>
jmp    11a0 <kernel+0xa0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @kernel(double* nocapture %0, double* nocapture writeonly %1) {
br label %3
3: ; preds = %2, %112
%4 = phi i64 [ 0, %2 ], [ %115, %112 ]
%5 = phi i64 [ 1, %2 ], [ %113, %112 ]
%6 = shl nuw nsw i64 %4, 10
%7 = add i64 %6, 1057
%8 = add i64 %6, 1087
%9 = or i64 %6, 33
%10 = add i64 %6, 2111
%11 = shl nsw i64 %5, 10
br label %12
12: ; preds = %3, %108
%13 = phi i64 [ 0, %3 ], [ %111, %108 ]
%14 = phi i64 [ 1, %3 ], [ %109, %108 ]
%15 = shl nuw nsw i64 %13, 5
%16 = shl nsw i64 %14, 5
%17 = add nuw nsw i64 %16, %11
%18 = add nsw i64 %17, -1
%19 = add nuw nsw i64 %17, 32
%20 = add nsw i64 %17, -32
%21 = add i64 %10, %15
%22 = getelementptr double, double* %0, i64 %21
%23 = add i64 %9, %15
%24 = getelementptr double, double* %0, i64 %23
%25 = add i64 %8, %15
%26 = getelementptr double, double* %1, i64 %25
%27 = add i64 %7, %15
%28 = getelementptr double, double* %1, i64 %27
%29 = icmp ult double* %28, %22
%30 = icmp ult double* %24, %26
%31 = and i1 %29, %30
br i1 %31, label %75, label %32
32: ; preds = %12, %32
%33 = phi i64 [ %73, %32 ], [ 0, %12 ]
%34 = or i64 %33, 1
%35 = add nuw nsw i64 %34, %16
%36 = add nuw nsw i64 %35, %11
%37 = getelementptr inbounds double, double* %0, i64 %36
%38 = bitcast double* %37 to <2 x double>*
%39 = load <2 x double>, <2 x double>* %38
%40 = add i64 %33, 2
%41 = add nuw nsw i64 %17, %40
%42 = getelementptr inbounds double, double* %0, i64 %41
%43 = bitcast double* %42 to <2 x double>*
%44 = load <2 x double>, <2 x double>* %43
%45 = fadd <2 x double> %39, %44
%46 = add nuw nsw i64 %18, %34
%47 = getelementptr inbounds double, double* %0, i64 %46
%48 = bitcast double* %47 to <2 x double>*
%49 = load <2 x double>, <2 x double>* %48
%50 = fadd <2 x double> %45, %49
%51 = add nuw nsw i64 %19, %34
%52 = getelementptr inbounds double, double* %0, i64 %51
%53 = bitcast double* %52 to <2 x double>*
%54 = load <2 x double>, <2 x double>* %53
%55 = fadd <2 x double> %50, %54
%56 = add nuw nsw i64 %20, %34
%57 = getelementptr inbounds double, double* %0, i64 %56
%58 = bitcast double* %57 to <2 x double>*
%59 = load <2 x double>, <2 x double>* %58
%60 = fadd <2 x double> %55, %59
%61 = add nsw i64 %36, -1024
%62 = getelementptr inbounds double, double* %0, i64 %61
%63 = bitcast double* %62 to <2 x double>*
%64 = load <2 x double>, <2 x double>* %63
%65 = fadd <2 x double> %60, %64
%66 = add nuw nsw i64 %36, 1024
%67 = getelementptr inbounds double, double* %0, i64 %66
%68 = bitcast double* %67 to <2 x double>*
%69 = load <2 x double>, <2 x double>* %68
%70 = fadd <2 x double> %65, %69
%71 = getelementptr inbounds double, double* %1, i64 %36
%72 = bitcast double* %71 to <2 x double>*
store <2 x double> %70, <2 x double>* %72
%73 = add nuw i64 %33, 2
%74 = icmp eq i64 %73, 30
br i1 %74, label %108, label %32
75: ; preds = %12, %75
%76 = phi i64 [ %81, %75 ], [ 1, %12 ]
%77 = add nuw nsw i64 %76, %16
%78 = add nuw nsw i64 %77, %11
%79 = getelementptr inbounds double, double* %0, i64 %78
%80 = load double, double* %79
%81 = add nuw nsw i64 %76, 1
%82 = add nuw nsw i64 %17, %81
%83 = getelementptr inbounds double, double* %0, i64 %82
%84 = load double, double* %83
%85 = fadd double %80, %84
%86 = add nuw nsw i64 %18, %76
%87 = getelementptr inbounds double, double* %0, i64 %86
%88 = load double, double* %87
%89 = fadd double %85, %88
%90 = add nuw nsw i64 %19, %76
%91 = getelementptr inbounds double, double* %0, i64 %90
%92 = load double, double* %91
%93 = fadd double %89, %92
%94 = add nuw nsw i64 %20, %76
%95 = getelementptr inbounds double, double* %0, i64 %94
%96 = load double, double* %95
%97 = fadd double %93, %96
%98 = add nsw i64 %78, -1024
%99 = getelementptr inbounds double, double* %0, i64 %98
%100 = load double, double* %99
%101 = fadd double %97, %100
%102 = add nuw nsw i64 %78, 1024
%103 = getelementptr inbounds double, double* %0, i64 %102
%104 = load double, double* %103
%105 = fadd double %101, %104
%106 = getelementptr inbounds double, double* %1, i64 %78
store double %105, double* %106
%107 = icmp eq i64 %81, 31
br i1 %107, label %108, label %75
108: ; preds = %32, %75
%109 = add nuw nsw i64 %14, 1
%110 = icmp eq i64 %109, 31
%111 = add i64 %13, 1
br i1 %110, label %112, label %12
112: ; preds = %108
%113 = add nuw nsw i64 %5, 1
%114 = icmp eq i64 %113, 31
%115 = add i64 %4, 1
br i1 %114, label %116, label %3
116: ; preds = %112
ret void
}