void seismic_exec(float **VEL, float **PPF, float **APF, float **NPF,
                  float *seismicPulseVector, int spPosX, int spPosY, int xDim,
                  int yDim, int timeSteps) {
  int i, j;
  int t;
  for (t = 0; t < timeSteps; t++) {
    APF[spPosY][spPosX] += seismicPulseVector[t];
    for (i = 2; i < (yDim - 2); i++) {
      for (j = 2; j < (xDim - 2); j++) {
        NPF[i][j] = 2.0f * APF[i][j] - PPF[i][j] +
                    VEL[i][j] * (16.0f * (APF[i][j - 1] + APF[i][j + 1] +
                                          APF[i - 1][j] + APF[i + 1][j]) -
                                 1.0f * (APF[i][j - 2] + APF[i][j + 2] +
                                         APF[i - 2][j] + APF[i + 2][j]) -
                                 60.0f * APF[i][j]);
      }
    }
    PPF = APF;
    APF = NPF;
    NPF = PPF;
  }
}


x86 assembly (clang, -O0):
<seismic_exec>:
push   %rbp
mov    %rsp,%rbp
mov    0x28(%rbp),%eax
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    0x28(%rbp),%eax
jge    1368 <seismic_exec+0x268>
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq 0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
movl   $0x2,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    0x20(%rbp),%ecx
sub    $0x2,%ecx
cmp    %ecx,%eax
jge    1342 <seismic_exec+0x242>
movl   $0x2,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    0x18(%rbp),%ecx
sub    $0x2,%ecx
cmp    %ecx,%eax
jge    132f <seismic_exec+0x22f>
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
addss  %xmm1,%xmm1
addss  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm3
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm3
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm3
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm3
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x34(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm4
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x34(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm4
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm4
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm4
movss  0xd4a(%rip),%xmm2
mulss  %xmm4,%xmm2
movd   %xmm2,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm2
movss  0xd2d(%rip),%xmm4
mulss  %xmm4,%xmm3
addss  %xmm2,%xmm3
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  0xd04(%rip),%xmm4
mulss  %xmm4,%xmm2
addss  %xmm3,%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1181 <seismic_exec+0x81>
jmp    1334 <seismic_exec+0x234>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1169 <seismic_exec+0x69>
mov    -0x18(%rbp),%rax
mov    %rax,-0x10(%rbp)
mov    -0x20(%rbp),%rax
mov    %rax,-0x18(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x20(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    112f <seismic_exec+0x2f>
pop    %rbp
ret

x86 assembly (clang, -O2):
<seismic_exec>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x30,%rsp
mov    %r8,(%rsp)
mov    %rcx,-0x50(%rsp)
mov    %rsi,-0x38(%rsp)
mov    %rdi,0x28(%rsp)
mov    0x80(%rsp),%eax
test   %eax,%eax
jle    151a <seismic_exec+0x41a>
mov    0x78(%rsp),%ecx
mov    0x70(%rsp),%esi
movslq 0x68(%rsp),%rdi
mov    %rdi,-0x8(%rsp)
movslq %r9d,%rdi
mov    %rdi,-0x10(%rsp)
lea    -0x2(%rcx),%edi
lea    -0x2(%rsi),%ecx
mov    %eax,%eax
mov    %rax,-0x18(%rsp)
lea    0x2(%rcx),%rax
mov    %rax,0x20(%rsp)
mov    %rcx,-0x28(%rsp)
lea    -0x2(%rcx),%r15
mov    %r15,%rax
and    $0xfffffffffffffffc,%rax
mov    %rax,0x18(%rsp)
add    $0x2,%rax
mov    %rax,0x10(%rsp)
xor    %ecx,%ecx
movss  0xe98(%rip),%xmm0
movss  0xe94(%rip),%xmm1
movaps 0xe68(%rip),%xmm8
movaps 0xe70(%rip),%xmm9
mov    %rdi,-0x40(%rsp)
mov    %r15,-0x48(%rsp)
jmp    11d8 <seismic_exec+0xd8>
nopl   0x0(%rax)
mov    0x8(%rsp),%rcx
add    $0x1,%rcx
mov    -0x30(%rsp),%rax
mov    %rax,-0x38(%rsp)
mov    -0x50(%rsp),%rdx
mov    %rax,-0x50(%rsp)
cmp    -0x18(%rsp),%rcx
je     151a <seismic_exec+0x41a>
mov    (%rsp),%rax
mov    %rcx,0x8(%rsp)
movss  (%rax,%rcx,4),%xmm4
mov    -0x8(%rsp),%rax
mov    %rdx,-0x30(%rsp)
mov    (%rdx,%rax,8),%rax
mov    -0x10(%rsp),%rcx
addss  (%rax,%rcx,4),%xmm4
movss  %xmm4,(%rax,%rcx,4)
mov    0x78(%rsp),%eax
cmp    $0x5,%eax
jl     11b0 <seismic_exec+0xb0>
mov    $0x2,%eax
jmp    122f <seismic_exec+0x12f>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%rax
mov    %rax,%rcx
mov    %rcx,%rax
cmp    %rdi,%rcx
je     11b0 <seismic_exec+0xb0>
cmp    $0x4,%esi
jle    1220 <seismic_exec+0x120>
lea    0x1(%rax),%rcx
mov    %rcx,-0x20(%rsp)
lea    0x2(%rax),%ecx
mov    -0x30(%rsp),%rbp
mov    0x0(%rbp,%rax,8),%r12
mov    -0x38(%rsp),%rdx
mov    (%rdx,%rax,8),%rsi
mov    0x28(%rsp),%rdx
mov    (%rdx,%rax,8),%r9
mov    -0x10(%rbp,%rax,8),%r14
mov    -0x8(%rbp,%rax,8),%r10
mov    0x8(%rbp,%rax,8),%rdi
mov    0x0(%rbp,%rcx,8),%r11
mov    -0x50(%rsp),%rcx
mov    (%rcx,%rax,8),%r13
mov    $0x2,%r8d
cmp    $0x4,%r15
jae    1340 <seismic_exec+0x240>
lea    (%r12,%r8,4),%rbx
lea    0x0(,%r8,4),%r12
add    %r13,%r12
lea    (%r11,%r8,4),%r11
lea    (%r14,%r8,4),%r14
lea    (%rdi,%r8,4),%rdi
lea    (%r10,%r8,4),%rdx
lea    (%r9,%r8,4),%rbp
lea    (%rsi,%r8,4),%rcx
mov    -0x28(%rsp),%rax
sub    %r8,%rax
xor    %esi,%esi
nopw   0x0(%rax,%rax,1)
movss  (%rbx,%rsi,4),%xmm2
movaps %xmm2,%xmm3
addss  %xmm2,%xmm3
subss  (%rcx,%rsi,4),%xmm3
movss  -0x8(%rbx,%rsi,4),%xmm4
movss  -0x4(%rbx,%rsi,4),%xmm5
addss  0x4(%rbx,%rsi,4),%xmm5
addss  (%rdx,%rsi,4),%xmm5
addss  (%rdi,%rsi,4),%xmm5
addss  0x8(%rbx,%rsi,4),%xmm4
addss  (%r14,%rsi,4),%xmm4
addss  (%r11,%rsi,4),%xmm4
mulss  %xmm0,%xmm5
subss  %xmm4,%xmm5
mulss  %xmm1,%xmm2
addss  %xmm5,%xmm2
mulss  0x0(%rbp,%rsi,4),%xmm2
addss  %xmm3,%xmm2
movss  %xmm2,(%r12,%rsi,4)
add    $0x1,%rsi
cmp    %rsi,%rax
jne    12c0 <seismic_exec+0x1c0>
mov    0x70(%rsp),%esi
mov    -0x40(%rsp),%rdi
mov    -0x20(%rsp),%rcx
jmp    1227 <seismic_exec+0x127>
nopl   0x0(%rax,%rax,1)
lea    0x8(%r13),%rcx
mov    -0x28(%rsp),%r15
lea    0x0(,%r15,4),%rax
add    %r13,%rax
mov    0x20(%rsp),%rdx
lea    (%r12,%rdx,4),%rdx
lea    0x8(%rsi),%rbx
mov    %rbx,-0x70(%rsp)
lea    (%rsi,%r15,4),%rbx
lea    0x8(%r9),%rbp
mov    %rbp,-0x60(%rsp)
lea    (%r9,%r15,4),%rbp
mov    %rbp,-0x78(%rsp)
lea    0x8(%r10),%rbp
mov    %rbp,-0x58(%rsp)
lea    (%r10,%r15,4),%rbp
mov    %rbp,-0x68(%rsp)
cmp    %rdx,%rcx
setb   -0x79(%rsp)
cmp    %rax,%r12
setb   -0x7a(%rsp)
cmp    %rbx,%rcx
setb   -0x7b(%rsp)
cmp    %rax,-0x70(%rsp)
lea    0x8(%rdi),%rdx
setb   -0x70(%rsp)
cmp    -0x78(%rsp),%rcx
setb   -0x78(%rsp)
cmp    %rax,-0x60(%rsp)
lea    (%rdi,%r15,4),%rbp
setb   -0x60(%rsp)
cmp    -0x68(%rsp),%rcx
setb   -0x68(%rsp)
cmp    %rax,-0x58(%rsp)
lea    0x8(%r14),%rbx
setb   -0x58(%rsp)
cmp    %rbp,%rcx
setb   -0x7c(%rsp)
cmp    %rax,%rdx
lea    (%r14,%r15,4),%rbp
setb   -0x7d(%rsp)
cmp    %rbp,%rcx
setb   %bpl
cmp    %rax,%rbx
lea    (%r11,%r15,4),%rbx
setb   %dl
cmp    %rbx,%rcx
lea    0x8(%r11),%rbx
setb   %cl
cmp    %rax,%rbx
setb   %al
mov    -0x7a(%rsp),%bl
test   %bl,-0x79(%rsp)
jne    1510 <seismic_exec+0x410>
mov    %edx,%ebx
mov    -0x7b(%rsp),%dl
and    -0x70(%rsp),%dl
mov    -0x48(%rsp),%r15
jne    1289 <seismic_exec+0x189>
mov    -0x60(%rsp),%dl
and    %dl,-0x78(%rsp)
jne    1289 <seismic_exec+0x189>
mov    -0x58(%rsp),%dl
and    %dl,-0x68(%rsp)
jne    1289 <seismic_exec+0x189>
mov    -0x7d(%rsp),%dl
and    %dl,-0x7c(%rsp)
jne    1289 <seismic_exec+0x189>
and    %bl,%bpl
jne    1289 <seismic_exec+0x189>
and    %al,%cl
jne    1289 <seismic_exec+0x189>
xor    %eax,%eax
mov    0x18(%rsp),%rcx
xchg   %ax,%ax
movups 0x8(%rsi,%rax,4),%xmm5
movups (%r12,%rax,4),%xmm6
movups 0x4(%r12,%rax,4),%xmm7
movups 0x8(%r12,%rax,4),%xmm4
movups 0xc(%r12,%rax,4),%xmm2
addps  %xmm7,%xmm2
movaps %xmm4,%xmm7
addps  %xmm4,%xmm7
subps  %xmm5,%xmm7
movups 0x8(%r9,%rax,4),%xmm5
movups 0x8(%r10,%rax,4),%xmm3
addps  %xmm2,%xmm3
movups 0x8(%rdi,%rax,4),%xmm2
addps  %xmm3,%xmm2
movups 0x10(%r12,%rax,4),%xmm3
addps  %xmm6,%xmm3
movups 0x8(%r14,%rax,4),%xmm6
addps  %xmm3,%xmm6
movups 0x8(%r11,%rax,4),%xmm3
addps  %xmm6,%xmm3
mulps  %xmm8,%xmm2
subps  %xmm3,%xmm2
mulps  %xmm9,%xmm4
addps  %xmm2,%xmm4
mulps  %xmm5,%xmm4
addps  %xmm7,%xmm4
movups %xmm4,0x8(%r13,%rax,4)
add    $0x4,%rax
cmp    %rax,%rcx
jne    1480 <seismic_exec+0x380>
mov    0x10(%rsp),%r8
cmp    %rcx,%r15
jne    1289 <seismic_exec+0x189>
jmp    1328 <seismic_exec+0x228>
mov    -0x48(%rsp),%r15
jmp    1289 <seismic_exec+0x189>
add    $0x30,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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