void jacobi(double **grid, double **new, int gridSize, int iter) {
  for (int k = 0; k < iter; k++) {
    for (int i = 1; i < gridSize + 1; i++) {
      for (int j = 1; j < gridSize + 1; j++) {
        new[i][j] = (grid[i - 1][j] + grid[i + 1][j] + grid[i][j - 1] +
                     grid[i][j + 1]) *
                    0.25;
      }
    }
    for (int i = 1; i < gridSize + 1; i++) {
      for (int j = 1; j < gridSize + 1; j++) {
        grid[i][j] =
            (new[i - 1][j] + new[i + 1][j] + new[i][j - 1] + new[i][j + 1]) *
            0.25;
      }
    }
  }
}


x86 assembly (clang, -O0):
<jacobi>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    12ec <jacobi+0x1ec>
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    11ff <jacobi+0xff>
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    11ec <jacobi+0xec>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x24(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  0xe3b(%rip),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1144 <jacobi+0x44>
jmp    11f1 <jacobi+0xf1>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    112c <jacobi+0x2c>
movl   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    12d9 <jacobi+0x1d9>
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    12c6 <jacobi+0x1c6>
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  0xd61(%rip),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    121e <jacobi+0x11e>
jmp    12cb <jacobi+0x1cb>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1206 <jacobi+0x106>
jmp    12de <jacobi+0x1de>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1119 <jacobi+0x19>
pop    %rbp
ret

x86 assembly (clang, -O2):
<jacobi>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rdx,-0x40(%rsp)
mov    %ecx,-0x34(%rsp)
test   %ecx,%ecx
jle    1564 <jacobi+0x464>
mov    %rsi,%r8
mov    -0x40(%rsp),%rax
lea    0x1(%rax),%r13d
lea    0x1(%r13),%rax
mov    %rax,-0x28(%rsp)
lea    -0x1(%r13),%rcx
mov    %rcx,%r10
and    $0xfffffffffffffffc,%r10
lea    0x1(%r10),%rax
mov    %rax,-0x30(%rsp)
mov    %r13,%rax
neg    %rax
mov    %rax,-0x10(%rsp)
xor    %eax,%eax
movsd  0xeb7(%rip),%xmm0
movapd 0xe9f(%rip),%xmm1
mov    %rsi,-0x18(%rsp)
mov    %rdi,-0x20(%rsp)
mov    %rcx,-0x48(%rsp)
jmp    1192 <jacobi+0x92>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    -0x8(%rsp),%rax
add    $0x1,%eax
cmp    -0x34(%rsp),%eax
je     1564 <jacobi+0x464>
mov    %rax,-0x8(%rsp)
cmpl   $0x0,-0x40(%rsp)
jle    1180 <jacobi+0x80>
mov    $0x1,%r12d
jmp    11b9 <jacobi+0xb9>
cs nopw 0x0(%rax,%rax,1)
cmp    %r13,%r12
je     1380 <jacobi+0x280>
mov    -0x8(%rdi,%r12,8),%rbx
mov    (%rdi,%r12,8),%rax
mov    (%r8,%r12,8),%rbp
add    $0x1,%r12
mov    %r12d,%ecx
mov    (%rdi,%rcx,8),%r14
mov    $0x1,%edx
cmpq   $0x4,-0x48(%rsp)
jb     12e0 <jacobi+0x1e0>
lea    0x8(%rbp),%rcx
lea    0x0(,%r13,8),%rdx
add    %rbp,%rdx
lea    0x8(%rbx),%r8
lea    (%rbx,%r13,8),%r11
lea    0x8(%r14),%r15
lea    (%r14,%r13,8),%r9
mov    -0x28(%rsp),%rsi
lea    (%rax,%rsi,8),%rsi
cmp    %r11,%rcx
setb   %dil
cmp    %rdx,%r8
setb   %r8b
cmp    %r9,%rcx
setb   %r11b
cmp    %rdx,%r15
setb   %r15b
cmp    %rsi,%rcx
setb   %cl
cmp    %rdx,%rax
setb   %r9b
mov    $0x1,%edx
test   %r8b,%dil
mov    -0x20(%rsp),%rdi
mov    -0x18(%rsp),%r8
jne    12e0 <jacobi+0x1e0>
mov    $0x1,%edx
and    %r15b,%r11b
jne    12e0 <jacobi+0x1e0>
mov    $0x1,%edx
and    %r9b,%cl
jne    12e0 <jacobi+0x1e0>
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nop
movupd 0x8(%rbx,%rcx,8),%xmm2
movupd 0x18(%rbx,%rcx,8),%xmm3
movupd 0x8(%r14,%rcx,8),%xmm4
addpd  %xmm2,%xmm4
movupd 0x18(%r14,%rcx,8),%xmm2
addpd  %xmm3,%xmm2
movupd (%rax,%rcx,8),%xmm3
addpd  %xmm4,%xmm3
movupd 0x10(%rax,%rcx,8),%xmm4
movupd 0x20(%rax,%rcx,8),%xmm5
addpd  %xmm4,%xmm2
addpd  %xmm5,%xmm2
addpd  %xmm4,%xmm3
mulpd  %xmm1,%xmm3
mulpd  %xmm1,%xmm2
movupd %xmm3,0x8(%rbp,%rcx,8)
movupd %xmm2,0x18(%rbp,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%r10
jne    1270 <jacobi+0x170>
mov    -0x30(%rsp),%rdx
cmp    %r10,-0x48(%rsp)
je     11b0 <jacobi+0xb0>
mov    %rdx,%rcx
test   $0x1,%r13b
jne    130e <jacobi+0x20e>
movsd  (%rbx,%rdx,8),%xmm2
addsd  (%r14,%rdx,8),%xmm2
addsd  -0x8(%rax,%rdx,8),%xmm2
lea    0x1(%rdx),%rcx
addsd  0x8(%rax,%rdx,8),%xmm2
mulsd  %xmm0,%xmm2
movsd  %xmm2,0x0(%rbp,%rdx,8)
not    %rdx
cmp    -0x10(%rsp),%rdx
je     11b0 <jacobi+0xb0>
nopl   0x0(%rax)
movsd  (%rbx,%rcx,8),%xmm2
addsd  (%r14,%rcx,8),%xmm2
addsd  -0x8(%rax,%rcx,8),%xmm2
addsd  0x8(%rax,%rcx,8),%xmm2
mulsd  %xmm0,%xmm2
movsd  %xmm2,0x0(%rbp,%rcx,8)
movsd  0x8(%rbx,%rcx,8),%xmm2
addsd  0x8(%r14,%rcx,8),%xmm2
addsd  (%rax,%rcx,8),%xmm2
addsd  0x10(%rax,%rcx,8),%xmm2
mulsd  %xmm0,%xmm2
movsd  %xmm2,0x8(%rbp,%rcx,8)
lea    0x2(%rcx),%rdx
mov    %rdx,%rcx
cmp    %rdx,%r13
jne    1320 <jacobi+0x220>
jmp    11b0 <jacobi+0xb0>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
cmpl   $0x0,-0x40(%rsp)
jle    1180 <jacobi+0x80>
mov    $0x1,%r12d
jmp    13a9 <jacobi+0x2a9>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
cmp    %r13,%r12
je     1180 <jacobi+0x80>
mov    -0x8(%r8,%r12,8),%r15
mov    (%r8,%r12,8),%rax
mov    (%rdi,%r12,8),%rbp
add    $0x1,%r12
mov    %r12d,%ecx
mov    (%r8,%rcx,8),%rcx
mov    $0x1,%r14d
cmpq   $0x4,-0x48(%rsp)
jb     14d0 <jacobi+0x3d0>
lea    0x8(%rbp),%rsi
lea    0x0(,%r13,8),%rdi
add    %rbp,%rdi
lea    0x8(%r15),%r8
lea    (%r15,%r13,8),%rdx
lea    0x8(%rcx),%r9
lea    (%rcx,%r13,8),%r11
mov    -0x28(%rsp),%rbx
lea    (%rax,%rbx,8),%r14
cmp    %rdx,%rsi
mov    %r12,%rbx
setb   %r12b
cmp    %rdi,%r8
setb   %r8b
cmp    %r11,%rsi
setb   %dl
cmp    %rdi,%r9
setb   %r9b
cmp    %r14,%rsi
setb   %r11b
cmp    %rdi,%rax
setb   %sil
mov    $0x1,%r14d
test   %r8b,%r12b
mov    %rbx,%r12
mov    -0x20(%rsp),%rdi
mov    -0x18(%rsp),%r8
jne    14d0 <jacobi+0x3d0>
mov    $0x1,%r14d
and    %r9b,%dl
jne    14d0 <jacobi+0x3d0>
mov    $0x1,%r14d
and    %sil,%r11b
jne    14d0 <jacobi+0x3d0>
xor    %edx,%edx
nop
movupd 0x8(%r15,%rdx,8),%xmm2
movupd 0x18(%r15,%rdx,8),%xmm3
movupd 0x8(%rcx,%rdx,8),%xmm4
addpd  %xmm2,%xmm4
movupd 0x18(%rcx,%rdx,8),%xmm2
addpd  %xmm3,%xmm2
movupd (%rax,%rdx,8),%xmm3
addpd  %xmm4,%xmm3
movupd 0x10(%rax,%rdx,8),%xmm4
movupd 0x20(%rax,%rdx,8),%xmm5
addpd  %xmm4,%xmm2
addpd  %xmm5,%xmm2
addpd  %xmm4,%xmm3
mulpd  %xmm1,%xmm3
mulpd  %xmm1,%xmm2
movupd %xmm3,0x8(%rbp,%rdx,8)
movupd %xmm2,0x18(%rbp,%rdx,8)
add    $0x4,%rdx
cmp    %rdx,%r10
jne    1460 <jacobi+0x360>
mov    -0x30(%rsp),%r14
cmp    %r10,-0x48(%rsp)
je     13a0 <jacobi+0x2a0>
mov    %r14,%rdx
test   $0x1,%r13b
jne    1502 <jacobi+0x402>
movsd  (%r15,%r14,8),%xmm2
addsd  (%rcx,%r14,8),%xmm2
addsd  -0x8(%rax,%r14,8),%xmm2
lea    0x1(%r14),%rdx
addsd  0x8(%rax,%r14,8),%xmm2
mulsd  %xmm0,%xmm2
movsd  %xmm2,0x0(%rbp,%r14,8)
not    %r14
cmp    -0x10(%rsp),%r14
je     13a0 <jacobi+0x2a0>
movsd  (%r15,%rdx,8),%xmm2
addsd  (%rcx,%rdx,8),%xmm2
addsd  -0x8(%rax,%rdx,8),%xmm2
addsd  0x8(%rax,%rdx,8),%xmm2
mulsd  %xmm0,%xmm2
movsd  %xmm2,0x0(%rbp,%rdx,8)
movsd  0x8(%r15,%rdx,8),%xmm2
addsd  0x8(%rcx,%rdx,8),%xmm2
addsd  (%rax,%rdx,8),%xmm2
addsd  0x10(%rax,%rdx,8),%xmm2
mulsd  %xmm0,%xmm2
movsd  %xmm2,0x8(%rbp,%rdx,8)
lea    0x2(%rdx),%rbx
mov    %rbx,%rdx
cmp    %rbx,%r13
jne    1510 <jacobi+0x410>
jmp    13a0 <jacobi+0x2a0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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