int wave13pt(int nx, int ny, int ns, const double c0, const double c1,
             const double c2, double *w0, double *w1, double *w2) {
  int i, j, k;
  unsigned int szarray = (unsigned int)nx * ny * ns;
#pragma acc data present(w0 [0:szarray], w1 [0:szarray], w2 [0:szarray])
  {
#pragma acc kernels
    {
#pragma acc loop independent gang(64)
      for (k = 2; k < ns - 2; k++) {
#pragma acc loop independent gang(124) vector(2)
        for (j = 2; j < ny - 2; j++) {
#pragma acc loop independent vector(128)
          for (i = 2; i < nx - 2; i++) {
            (w2[(i) + nx * (j) + nx * ny * (k)]) =
                c0 * (w1[(i) + nx * (j) + nx * ny * (k)]) -
                (w0[(i) + nx * (j) + nx * ny * (k)]) +
                c1 * ((w1[(i + 1) + nx * (j) + nx * ny * (k)]) +
                      (w1[(i - 1) + nx * (j) + nx * ny * (k)]) +
                      (w1[(i) + nx * (j + 1) + nx * ny * (k)]) +
                      (w1[(i) + nx * (j - 1) + nx * ny * (k)]) +
                      (w1[(i) + nx * (j) + nx * ny * (k + 1)]) +
                      (w1[(i) + nx * (j) + nx * ny * (k - 1)])) +
                c2 * ((w1[(i + 2) + nx * (j) + nx * ny * (k)]) +
                      (w1[(i - 2) + nx * (j) + nx * ny * (k)]) +
                      (w1[(i) + nx * (j + 2) + nx * ny * (k)]) +
                      (w1[(i) + nx * (j - 2) + nx * ny * (k)]) +
                      (w1[(i) + nx * (j) + nx * ny * (k + 2)]) +
                      (w1[(i) + nx * (j) + nx * ny * (k - 2)]));
          }
        }
      }
    }
  }
  return 0;
}


x86 assembly (clang, -O0):
<wave13pt>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
movsd  %xmm2,-0x28(%rbp)
mov    %rcx,-0x30(%rbp)
mov    %r8,-0x38(%rbp)
mov    %r9,-0x40(%rbp)
mov    -0x4(%rbp),%eax
imul   -0x8(%rbp),%eax
imul   -0xc(%rbp),%eax
mov    %eax,-0x50(%rbp)
movl   $0x2,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x2,%ecx
cmp    %ecx,%eax
jge    144f <wave13pt+0x34f>
movl   $0x2,-0x48(%rbp)
mov    -0x48(%rbp),%eax
mov    -0x8(%rbp),%ecx
sub    $0x2,%ecx
cmp    %ecx,%eax
jge    143c <wave13pt+0x33c>
movl   $0x2,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x2,%ecx
cmp    %ecx,%eax
jge    1429 <wave13pt+0x329>
movsd  -0x18(%rbp),%xmm0
mov    -0x38(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    -0x4(%rbp),%edx
imul   -0x48(%rbp),%edx
add    %edx,%ecx
mov    -0x4(%rbp),%edx
imul   -0x8(%rbp),%edx
imul   -0x4c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x30(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    -0x4(%rbp),%edx
imul   -0x48(%rbp),%edx
add    %edx,%ecx
mov    -0x4(%rbp),%edx
imul   -0x8(%rbp),%edx
imul   -0x4c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0x20(%rbp),%xmm1
mov    -0x38(%rbp),%rax
mov    -0x44(%rbp),%ecx
add    $0x1,%ecx
mov    -0x4(%rbp),%edx
imul   -0x48(%rbp),%edx
add    %edx,%ecx
mov    -0x4(%rbp),%edx
imul   -0x8(%rbp),%edx
imul   -0x4c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x38(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x4(%rbp),%edx
imul   -0x48(%rbp),%edx
add    %edx,%ecx
mov    -0x4(%rbp),%edx
imul   -0x8(%rbp),%edx
imul   -0x4c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm2
mov    -0x38(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    -0x4(%rbp),%edx
mov    -0x48(%rbp),%esi
add    $0x1,%esi
imul   %esi,%edx
add    %edx,%ecx
mov    -0x4(%rbp),%edx
imul   -0x8(%rbp),%edx
imul   -0x4c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm2
mov    -0x38(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    -0x4(%rbp),%edx
mov    -0x48(%rbp),%esi
sub    $0x1,%esi
imul   %esi,%edx
add    %edx,%ecx
mov    -0x4(%rbp),%edx
imul   -0x8(%rbp),%edx
imul   -0x4c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm2
mov    -0x38(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    -0x4(%rbp),%edx
imul   -0x48(%rbp),%edx
add    %edx,%ecx
mov    -0x4(%rbp),%edx
imul   -0x8(%rbp),%edx
mov    -0x4c(%rbp),%esi
add    $0x1,%esi
imul   %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm2
mov    -0x38(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    -0x4(%rbp),%edx
imul   -0x48(%rbp),%edx
add    %edx,%ecx
mov    -0x4(%rbp),%edx
imul   -0x8(%rbp),%edx
mov    -0x4c(%rbp),%esi
sub    $0x1,%esi
imul   %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x28(%rbp),%xmm0
mov    -0x38(%rbp),%rax
mov    -0x44(%rbp),%ecx
add    $0x2,%ecx
mov    -0x4(%rbp),%edx
imul   -0x48(%rbp),%edx
add    %edx,%ecx
mov    -0x4(%rbp),%edx
imul   -0x8(%rbp),%edx
imul   -0x4c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x38(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x2,%ecx
mov    -0x4(%rbp),%edx
imul   -0x48(%rbp),%edx
add    %edx,%ecx
mov    -0x4(%rbp),%edx
imul   -0x8(%rbp),%edx
imul   -0x4c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm2
mov    -0x38(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    -0x4(%rbp),%edx
mov    -0x48(%rbp),%esi
add    $0x2,%esi
imul   %esi,%edx
add    %edx,%ecx
mov    -0x4(%rbp),%edx
imul   -0x8(%rbp),%edx
imul   -0x4c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm2
mov    -0x38(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    -0x4(%rbp),%edx
mov    -0x48(%rbp),%esi
sub    $0x2,%esi
imul   %esi,%edx
add    %edx,%ecx
mov    -0x4(%rbp),%edx
imul   -0x8(%rbp),%edx
imul   -0x4c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm2
mov    -0x38(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    -0x4(%rbp),%edx
imul   -0x48(%rbp),%edx
add    %edx,%ecx
mov    -0x4(%rbp),%edx
imul   -0x8(%rbp),%edx
mov    -0x4c(%rbp),%esi
add    $0x2,%esi
imul   %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm2
mov    -0x38(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    -0x4(%rbp),%edx
imul   -0x48(%rbp),%edx
add    %edx,%ecx
mov    -0x4(%rbp),%edx
imul   -0x8(%rbp),%edx
mov    -0x4c(%rbp),%esi
sub    $0x2,%esi
imul   %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x40(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    -0x4(%rbp),%edx
imul   -0x48(%rbp),%edx
add    %edx,%ecx
mov    -0x4(%rbp),%edx
imul   -0x8(%rbp),%edx
imul   -0x4c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    116d <wave13pt+0x6d>
jmp    142e <wave13pt+0x32e>
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1155 <wave13pt+0x55>
jmp    1441 <wave13pt+0x341>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    113d <wave13pt+0x3d>
xor    %eax,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<wave13pt>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %r9,-0x70(%rsp)
mov    %rcx,-0x78(%rsp)
cmp    $0x5,%edx
jl     142d <wave13pt+0x32d>
add    $0xfffffffe,%edx
mov    %esi,%eax
imul   %edi,%eax
lea    -0x2(%rsi),%ecx
lea    -0x2(%rdi),%ebx
movslq %edi,%rbp
mov    %rbp,-0x40(%rsp)
cltq
mov    %rax,-0x60(%rsp)
lea    0x10(%r8),%rax
mov    %rax,-0x68(%rsp)
addq   $0x10,-0x78(%rsp)
addq   $0x10,-0x70(%rsp)
add    $0xfffffffffffffffe,%rbx
mov    %rbx,-0x38(%rsp)
mov    $0x2,%ebp
movabs $0x100000000,%rbx
mov    %rdx,-0x50(%rsp)
mov    %rsi,-0x58(%rsp)
mov    %rdi,-0x28(%rsp)
mov    %rcx,-0x30(%rsp)
jmp    11a3 <wave13pt+0xa3>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %rbp,%rax
add    $0x1,%rax
mov    %rax,%rbp
cmp    %rdx,%rax
je     142d <wave13pt+0x32d>
cmp    $0x4,%esi
jle    1190 <wave13pt+0x90>
mov    %rbp,%rax
mov    %rbp,%r14
mov    -0x60(%rsp),%rdx
imul   %rdx,%r14
lea    0x1(%rbp),%r9
mov    %r9,-0x48(%rsp)
imul   %rdx,%r9
lea    -0x1(%rbp),%r10
imul   %rdx,%r10
lea    0x2(%rbp),%r11
imul   %rdx,%r11
add    $0xfffffffffffffffe,%rbp
imul   %rdx,%rbp
mov    -0x68(%rsp),%rdx
lea    (%rdx,%rbp,8),%rbp
mov    %rbp,0x18(%rsp)
lea    (%rdx,%r11,8),%rbp
mov    %rbp,0x10(%rsp)
lea    (%rdx,%r10,8),%rbp
mov    %rbp,0x8(%rsp)
lea    (%rdx,%r9,8),%rdx
mov    %rdx,(%rsp)
imul   %esi,%eax
lea    -0x2(%rax),%edx
mov    %rdx,-0x8(%rsp)
lea    0x2(%rax),%edx
mov    %rdx,-0x10(%rsp)
lea    -0x1(%rax),%edx
mov    %rdx,-0x18(%rsp)
mov    %rax,0x28(%rsp)
add    $0x1,%eax
mov    %eax,-0x7c(%rsp)
mov    %r14,0x20(%rsp)
lea    0x3(%r14),%rax
mov    %rax,-0x20(%rsp)
mov    $0x2,%r9d
jmp    1260 <wave13pt+0x160>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x1,%r9
mov    %r9,%rax
cmp    %rcx,%rax
je     1412 <wave13pt+0x312>
cmp    $0x4,%edi
jle    1250 <wave13pt+0x150>
mov    %r9,%rsi
imul   -0x40(%rsp),%rsi
lea    0x1(%r9),%rax
mov    %rax,0x30(%rsp)
mov    0x18(%rsp),%rax
lea    (%rax,%rsi,8),%rax
mov    %rax,0x50(%rsp)
mov    0x10(%rsp),%rax
lea    (%rax,%rsi,8),%rax
mov    %rax,0x48(%rsp)
mov    0x8(%rsp),%rax
lea    (%rax,%rsi,8),%rax
mov    %rax,0x40(%rsp)
mov    (%rsp),%rax
lea    (%rax,%rsi,8),%rax
mov    %rax,0x38(%rsp)
mov    0x20(%rsp),%rax
add    %rsi,%rax
mov    -0x68(%rsp),%rcx
lea    (%rcx,%rax,8),%r14
mov    -0x78(%rsp),%rdx
lea    (%rdx,%rax,8),%rbp
mov    -0x70(%rsp),%rdx
lea    (%rdx,%rax,8),%r10
mov    -0x8(%rsp),%rax
lea    (%rax,%r9,1),%r12d
imul   %edi,%r12d
add    $0x2,%r12d
mov    -0x10(%rsp),%rax
lea    (%rax,%r9,1),%r13d
imul   %edi,%r13d
add    $0x2,%r13d
mov    0x28(%rsp),%rax
lea    (%rax,%r9,1),%ecx
imul   %edi,%ecx
mov    -0x18(%rsp),%rax
lea    (%rax,%r9,1),%r15d
imul   %edi,%r15d
add    $0x2,%r15d
add    -0x7c(%rsp),%r9d
imul   %edi,%r9d
add    $0x2,%r9d
add    -0x20(%rsp),%esi
shl    $0x20,%rsi
xor    %eax,%eax
mov    -0x38(%rsp),%rdi
nopl   0x0(%rax)
mov    %rsi,%rdx
sar    $0x1d,%rdx
movsd  (%r8,%rdx,1),%xmm3
lea    (%rcx,%rax,1),%r11
lea    0x1(%r11),%edx
movslq %edx,%rdx
addsd  (%r8,%rdx,8),%xmm3
lea    (%r9,%rax,1),%edx
movslq %edx,%rdx
addsd  (%r8,%rdx,8),%xmm3
lea    (%r15,%rax,1),%edx
movslq %edx,%rdx
addsd  (%r8,%rdx,8),%xmm3
lea    0x4(%r11),%edx
movslq %edx,%rdx
movsd  (%r8,%rdx,8),%xmm4
movslq %r11d,%rdx
addsd  (%r8,%rdx,8),%xmm4
lea    0x0(%r13,%rax,1),%edx
movslq %edx,%rdx
addsd  (%r8,%rdx,8),%xmm4
lea    (%r12,%rax,1),%edx
movslq %edx,%rdx
addsd  (%r8,%rdx,8),%xmm4
movsd  (%r14,%rax,8),%xmm5
mulsd  %xmm0,%xmm5
subsd  0x0(%rbp,%rax,8),%xmm5
mov    0x38(%rsp),%rdx
addsd  (%rdx,%rax,8),%xmm3
mov    0x40(%rsp),%rdx
addsd  (%rdx,%rax,8),%xmm3
mulsd  %xmm1,%xmm3
addsd  %xmm5,%xmm3
mov    0x48(%rsp),%rdx
addsd  (%rdx,%rax,8),%xmm4
mov    0x50(%rsp),%rdx
addsd  (%rdx,%rax,8),%xmm4
mulsd  %xmm2,%xmm4
addsd  %xmm3,%xmm4
movsd  %xmm4,(%r10,%rax,8)
add    $0x1,%rax
add    %rbx,%rsi
cmp    %rax,%rdi
jne    1330 <wave13pt+0x230>
mov    0x30(%rsp),%rax
mov    %rax,%r9
mov    -0x28(%rsp),%rdi
mov    -0x30(%rsp),%rcx
cmp    %rcx,%rax
jne    1260 <wave13pt+0x160>
mov    -0x48(%rsp),%rax
mov    %rax,%rbp
mov    -0x50(%rsp),%rdx
mov    -0x58(%rsp),%rsi
cmp    %rdx,%rax
jne    11a3 <wave13pt+0xa3>
xor    %eax,%eax
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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