void jacobi_(int nx, int ny, double c0, double c1, double c2, double *w0,
             double *w1) {
  int i, j;
  int i00, im10, ip10, i0m1, i0p1;
  int im1m1, im1p1, ip1m1, ip1p1;
#pragma acc data present(w0 [0:nx * ny], w1 [0:nx * ny])
  {
#pragma acc kernels
    {
#pragma acc loop independent gang(2046)
      for (j = 1; j < ny - 1; j++) {
#pragma acc loop independent gang(16) vector(128)
        for (i = 1; i < nx - 1; i++) {
          i00 = ((i) + (nx)*j);
          im10 = ((i - 1) + (nx)*j);
          ip10 = ((i + 1) + (nx)*j);
          i0m1 = ((i) + (nx)*j - 1);
          i0p1 = ((i) + (nx)*j + 1);
          im1m1 = ((i - 1) + (nx)*j - 1);
          im1p1 = ((i - 1) + (nx)*j + 1);
          ip1m1 = ((i + 1) + (nx)*j - 1);
          ip1p1 = ((i + 1) + (nx)*j + 1);
          w1[i00] = c0 * w0[i00] +
                    c1 * (w0[im10] + w0[i0m1] + w0[ip10] + w0[i0p1]) +
                    c2 * (w0[im1m1] + w0[im1p1] + w0[ip1m1] + w0[ip1p1]);
        }
      }
    }
  }
}


x86 assembly (clang, -O0):
<jacobi_>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
movsd  %xmm2,-0x20(%rbp)
mov    %rdx,-0x28(%rbp)
mov    %rcx,-0x30(%rbp)
movl   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%eax
mov    -0x8(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    12c2 <jacobi_+0x1c2>
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    12af <jacobi_+0x1af>
mov    -0x34(%rbp),%eax
mov    -0x4(%rbp),%ecx
imul   -0x38(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x34(%rbp),%eax
sub    $0x1,%eax
mov    -0x4(%rbp),%ecx
imul   -0x38(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x40(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    -0x4(%rbp),%ecx
imul   -0x38(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x44(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x4(%rbp),%ecx
imul   -0x38(%rbp),%ecx
add    %ecx,%eax
sub    $0x1,%eax
mov    %eax,-0x48(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x4(%rbp),%ecx
imul   -0x38(%rbp),%ecx
add    %ecx,%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x34(%rbp),%eax
sub    $0x1,%eax
mov    -0x4(%rbp),%ecx
imul   -0x38(%rbp),%ecx
add    %ecx,%eax
sub    $0x1,%eax
mov    %eax,-0x50(%rbp)
mov    -0x34(%rbp),%eax
sub    $0x1,%eax
mov    -0x4(%rbp),%ecx
imul   -0x38(%rbp),%ecx
add    %ecx,%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    -0x4(%rbp),%ecx
imul   -0x38(%rbp),%ecx
add    %ecx,%eax
sub    $0x1,%eax
mov    %eax,-0x58(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    -0x4(%rbp),%ecx
imul   -0x38(%rbp),%ecx
add    %ecx,%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
movsd  -0x10(%rbp),%xmm1
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x18(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm3
mov    -0x28(%rbp),%rax
movslq -0x48(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm3
mov    -0x28(%rbp),%rax
movslq -0x44(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm3
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm3
mulsd  %xmm3,%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x20(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x28(%rbp),%rax
movslq -0x54(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm2
mov    -0x28(%rbp),%rax
movslq -0x58(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm2
mov    -0x28(%rbp),%rax
movslq -0x5c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1140 <jacobi_+0x40>
jmp    12b4 <jacobi_+0x1b4>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1128 <jacobi_+0x28>
pop    %rbp
ret

x86 assembly (clang, -O2):
<jacobi_>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %rcx,-0x78(%rsp)
mov    %rsi,-0x68(%rsp)
cmp    $0x3,%esi
jl     1479 <jacobi_+0x379>
mov    -0x68(%rsp),%rax
add    $0xffffffff,%eax
mov    %rax,-0x68(%rsp)
lea    -0x1(%rdi),%r13d
movslq %edi,%rcx
lea    -0x2(%r13),%rax
mov    %rax,0x8(%rsp)
lea    0x1(%rcx),%rax
mov    %rax,(%rsp)
lea    (%rcx,%r13,1),%rax
mov    %rax,-0x8(%rsp)
lea    0x2(%rcx,%r13,1),%rax
mov    %rax,-0x10(%rsp)
mov    %ecx,%r15d
lea    -0x1(%r13),%rax
lea    (%rdx,%r13,8),%rsi
add    $0xfffffffffffffff8,%rsi
mov    %rsi,-0x18(%rsp)
mov    %rax,%rsi
and    $0xfffffffffffffffe,%rsi
mov    %rsi,-0x20(%rsp)
mov    %rax,-0x48(%rsp)
or     $0x1,%rax
mov    %rax,-0x28(%rsp)
movapd %xmm1,%xmm8
unpcklpd %xmm1,%xmm8
movapd %xmm0,%xmm9
unpcklpd %xmm0,%xmm9
movapd %xmm2,%xmm5
unpcklpd %xmm2,%xmm5
lea    0x18(%rdx,%rcx,8),%r10
lea    0x0(,%rcx,8),%rax
mov    %rax,0x10(%rsp)
mov    -0x78(%rsp),%rax
lea    0x8(%rax,%rcx,8),%r14
lea    -0x2(%rcx),%eax
mov    %rax,-0x50(%rsp)
mov    $0x1,%r12d
xor    %eax,%eax
mov    %rax,-0x60(%rsp)
mov    %rcx,-0x38(%rsp)
mov    %rcx,-0x58(%rsp)
mov    %r15,-0x40(%rsp)
mov    %r13,-0x70(%rsp)
mov    %rdi,-0x30(%rsp)
jmp    123b <jacobi_+0x13b>
nop
add    $0x1,%r12
addq   $0x1,-0x60(%rsp)
mov    0x10(%rsp),%rcx
add    %rcx,%r10
mov    -0x40(%rsp),%rax
add    %rax,%r13
add    %rax,%r15
add    %rcx,%r14
mov    -0x30(%rsp),%rdi
mov    -0x50(%rsp),%rax
add    %edi,%eax
mov    %rax,-0x50(%rsp)
mov    -0x58(%rsp),%rax
add    -0x38(%rsp),%rax
mov    %rax,-0x58(%rsp)
cmp    -0x68(%rsp),%r12
je     1479 <jacobi_+0x379>
cmp    $0x3,%edi
jl     11f0 <jacobi_+0xf0>
mov    $0x1,%r9d
cmpq   $0x2,-0x48(%rsp)
jae    12f0 <jacobi_+0x1f0>
mov    -0x50(%rsp),%rax
lea    (%rax,%r9,1),%esi
lea    (%r9,%r13,1),%r11d
mov    -0x70(%rsp),%r8
sub    %r9,%r8
add    -0x58(%rsp),%r9
mov    -0x78(%rsp),%rax
lea    (%rax,%r9,8),%rax
lea    (%rdx,%r9,8),%rbp
xor    %edi,%edi
nopl   0x0(%rax,%rax,1)
lea    (%r11,%rdi,1),%ebx
lea    (%rsi,%rdi,1),%ecx
movslq %ebx,%rbx
movsd  (%rdx,%rbx,8),%xmm3
addsd  -0x8(%rbp,%rdi,8),%xmm3
movsd  0x0(%rbp,%rdi,8),%xmm4
movsd  0x8(%rbp,%rdi,8),%xmm6
addsd  %xmm6,%xmm3
addsd  %xmm6,%xmm3
mulsd  %xmm1,%xmm3
movapd %xmm0,%xmm6
mulsd  %xmm4,%xmm6
addsd  %xmm3,%xmm6
movslq %ecx,%rcx
movsd  (%rdx,%rcx,8),%xmm3
addsd  %xmm4,%xmm3
addsd  %xmm4,%xmm3
addsd  0x10(%rbp,%rdi,8),%xmm3
mulsd  %xmm2,%xmm3
addsd  %xmm6,%xmm3
movsd  %xmm3,(%rax,%rdi,8)
add    $0x1,%rdi
cmp    %rdi,%r8
jne    1280 <jacobi_+0x180>
jmp    11f0 <jacobi_+0xf0>
nopw   0x0(%rax,%rax,1)
mov    -0x60(%rsp),%rax
mov    -0x30(%rsp),%rcx
imul   %ecx,%eax
lea    (%rax,%rcx,1),%esi
mov    0x8(%rsp),%rcx
lea    (%rsi,%rcx,1),%edi
cmp    %esi,%edi
jl     1252 <jacobi_+0x152>
mov    %rcx,%rsi
shr    $0x20,%rsi
jne    1252 <jacobi_+0x152>
add    -0x70(%rsp),%eax
lea    (%rax,%rcx,1),%edi
cmp    %eax,%edi
jl     1252 <jacobi_+0x152>
test   %rsi,%rsi
jne    1252 <jacobi_+0x152>
mov    -0x60(%rsp),%rbp
mov    %rbp,%rax
mov    -0x38(%rsp),%rbx
imul   %rbx,%rax
mov    (%rsp),%rcx
lea    (%rcx,%rax,1),%r11
mov    -0x8(%rsp),%rcx
lea    (%rcx,%rax,1),%rdi
mov    -0x10(%rsp),%rcx
lea    (%rcx,%rax,1),%r8
add    %rbx,%rax
imul   -0x40(%rsp),%ebp
mov    -0x70(%rsp),%rcx
lea    (%rcx,%rbp,1),%esi
add    %ebp,%ebx
mov    -0x78(%rsp),%rcx
lea    (%rcx,%r11,8),%r11
lea    (%rcx,%rdi,8),%rdi
lea    (%rdx,%rax,8),%rax
lea    (%rdx,%r8,8),%r8
movslq %esi,%rsi
cmp    %r8,%r11
setb   -0x79(%rsp)
cmp    %rdi,%rax
lea    (%rdx,%rsi,8),%r8
mov    -0x18(%rsp),%rcx
lea    (%rcx,%rsi,8),%rax
movslq %ebx,%rsi
setb   %bpl
cmp    %rax,%r11
setb   %al
cmp    %rdi,%r8
lea    (%rcx,%rsi,8),%r8
setb   %cl
cmp    %r8,%r11
lea    (%rdx,%rsi,8),%rsi
setb   %r8b
cmp    %rdi,%rsi
setb   %dil
test   %bpl,-0x79(%rsp)
jne    1252 <jacobi_+0x152>
and    %cl,%al
jne    1252 <jacobi_+0x152>
and    %dil,%r8b
jne    1252 <jacobi_+0x152>
xor    %eax,%eax
mov    -0x20(%rsp),%rcx
xchg   %ax,%ax
lea    (%r15,%rax,1),%esi
lea    (%rax,%r13,1),%edi
movslq %esi,%rsi
movupd (%rdx,%rsi,8),%xmm6
movupd -0x18(%r10,%rax,8),%xmm7
addpd  %xmm6,%xmm7
movupd -0x10(%r10,%rax,8),%xmm6
movupd -0x8(%r10,%rax,8),%xmm3
movupd (%r10,%rax,8),%xmm4
addpd  %xmm3,%xmm7
addpd  %xmm3,%xmm7
mulpd  %xmm8,%xmm7
movapd %xmm9,%xmm3
mulpd  %xmm6,%xmm3
addpd  %xmm7,%xmm3
movslq %edi,%rsi
movupd (%rdx,%rsi,8),%xmm7
addpd  %xmm6,%xmm7
addpd  %xmm6,%xmm7
addpd  %xmm4,%xmm7
mulpd  %xmm5,%xmm7
addpd  %xmm3,%xmm7
movupd %xmm7,(%r14,%rax,8)
add    $0x2,%rax
cmp    %rax,%rcx
jne    13f0 <jacobi_+0x2f0>
mov    -0x28(%rsp),%r9
cmp    %rcx,-0x48(%rsp)
je     11f0 <jacobi_+0xf0>
jmp    1252 <jacobi_+0x152>
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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