void horn_dt(image1, image2, dt, nx, ny) float *dt, *image1, *image2;
int nx, ny;
{
  int index, i, j;
  for (i = 0; i < ny - 1; i++) {
    for (j = 0; j < nx - 1; j++) {
      index = (i * nx) + j;
      dt[index] = ((image2[index] + image2[index + 1] + image2[index + nx] +
                    image2[index + nx + 1]) /
                   4.0) -
                  ((image1[index] + image1[index + 1] + image1[index + nx] +
                    image1[index + nx + 1]) /
                   4.0);
    }
  }
  for (j = 0; j < nx; j++) {
    dt[(ny - 1) * nx + j] = 0.0;
  }
  for (i = 0; i < ny; i++) {
    dt[(i * nx) + (nx - 1)] = 0.0;
  }
}


x86 assembly (clang, -O0):
<horn_dt>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1236 <horn_dt+0x136>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1223 <horn_dt+0x123>
mov    -0x28(%rbp),%eax
imul   -0x1c(%rbp),%eax
add    -0x2c(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe5a(%rip),%xmm1
divsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
movsd  0xe04(%rip),%xmm2
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1136 <horn_dt+0x36>
jmp    1228 <horn_dt+0x128>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    111e <horn_dt+0x1e>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1273 <horn_dt+0x173>
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    123d <horn_dt+0x13d>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    12b2 <horn_dt+0x1b2>
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
mov    -0x1c(%rbp),%edx
sub    $0x1,%edx
add    %edx,%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    127a <horn_dt+0x17a>
pop    %rbp
ret

x86 assembly (clang, -O2):
<horn_dt>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x88,%rsp
mov    %ecx,%ebx
mov    %rdx,0x18(%rsp)
mov    %rsi,0x28(%rsp)
mov    %rdi,0x20(%rsp)
lea    -0x1(%r8),%eax
mov    %eax,0xc(%rsp)
movslq %ecx,%r13
mov    %r8,0x58(%rsp)
cmp    $0x2,%r8d
jl     145b <horn_dt+0x34b>
lea    -0x1(%r13),%eax
mov    0xc(%rsp),%ecx
mov    %rcx,0x80(%rsp)
lea    (%rax,%r13,1),%rcx
add    $0x1,%rcx
mov    %rcx,0x70(%rsp)
lea    0x1(%rax),%rcx
mov    %rcx,0x68(%rsp)
mov    %rax,0x30(%rsp)
and    $0xfffffffc,%eax
mov    %rax,0x60(%rsp)
lea    0x0(,%r13,4),%r11
mov    0x20(%rsp),%rcx
lea    0x4(%rcx),%r15
lea    0x1(%r13),%rax
mov    %rax,0x50(%rsp)
lea    0x4(%rcx,%r13,4),%r8
mov    0x28(%rsp),%rax
lea    (%rax,%r13,4),%r14
add    $0x4,%r14
lea    0x4(%rax),%r9
xor    %r10d,%r10d
movsd  0xe56(%rip),%xmm0
movapd 0xe3e(%rip),%xmm1
mov    %rcx,0x48(%rsp)
mov    %rax,0x40(%rsp)
mov    0x18(%rsp),%r12
mov    %r13,0x10(%rsp)
mov    %ebx,0x3c(%rsp)
jmp    1214 <horn_dt+0x104>
nopl   0x0(%rax)
add    $0x1,%r10
add    %r11,%r12
add    %r11,%r15
add    %r11,%r8
add    %r11,%r14
add    %r11,%r9
add    %r13,0x50(%rsp)
add    %r11,0x40(%rsp)
add    %r11,0x48(%rsp)
cmp    0x80(%rsp),%r10
mov    0x3c(%rsp),%ebx
je     145b <horn_dt+0x34b>
cmp    $0x2,%ebx
jl     11e0 <horn_dt+0xd0>
mov    0x30(%rsp),%rcx
cmp    $0x4,%ecx
jae    12d0 <horn_dt+0x1c0>
xor    %edx,%edx
mov    0x30(%rsp),%rbx
sub    %rdx,%rbx
mov    0x50(%rsp),%rax
lea    (%rdx,%rax,1),%rbp
mov    0x20(%rsp),%rax
lea    (%rax,%rbp,4),%rcx
lea    (%r12,%rdx,4),%rax
mov    0x40(%rsp),%rsi
lea    (%rsi,%rdx,4),%rdi
mov    0x48(%rsp),%rsi
lea    (%rsi,%rdx,4),%rsi
mov    0x28(%rsp),%rdx
lea    (%rdx,%rbp,4),%rdx
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  (%rdi,%rbp,4),%xmm2
addss  0x4(%rdi,%rbp,4),%xmm2
addss  -0x4(%rdx,%rbp,4),%xmm2
addss  (%rdx,%rbp,4),%xmm2
cvtss2sd %xmm2,%xmm2
movss  (%rsi,%rbp,4),%xmm3
addss  0x4(%rsi,%rbp,4),%xmm3
addss  -0x4(%rcx,%rbp,4),%xmm3
mulsd  %xmm0,%xmm2
addss  (%rcx,%rbp,4),%xmm3
cvtss2sd %xmm3,%xmm3
mulsd  %xmm0,%xmm3
subsd  %xmm3,%xmm2
cvtsd2ss %xmm2,%xmm2
movss  %xmm2,(%rax,%rbp,4)
add    $0x1,%rbp
cmp    %rbp,%rbx
jne    1270 <horn_dt+0x160>
jmp    11e0 <horn_dt+0xd0>
nopw   0x0(%rax,%rax,1)
mov    %r10,%rax
imul   %r13,%rax
mov    %r13,%rdx
lea    (%rax,%rcx,1),%r13
lea    (%rax,%rdx,1),%rcx
mov    %rcx,0x78(%rsp)
mov    0x70(%rsp),%rdx
lea    (%rdx,%rax,1),%rsi
mov    0x68(%rsp),%rdx
lea    (%rdx,%rax,1),%rdi
mov    0x18(%rsp),%rbx
lea    (%rbx,%rax,4),%rbp
lea    (%rbx,%r13,4),%rbx
mov    0x28(%rsp),%rdx
lea    (%rdx,%rcx,4),%r13
lea    (%rdx,%rsi,4),%rcx
cmp    %rcx,%rbp
lea    (%rdx,%rdi,4),%rcx
setb   0xb(%rsp)
cmp    %rbx,%r13
setb   0xa(%rsp)
cmp    %rcx,%rbp
lea    (%rdx,%rax,4),%rcx
setb   %r13b
cmp    %rbx,%rcx
mov    0x20(%rsp),%rdx
lea    (%rdx,%rsi,4),%rcx
setb   0x9(%rsp)
cmp    %rcx,%rbp
mov    0x78(%rsp),%rcx
lea    (%rdx,%rcx,4),%rcx
setb   %sil
cmp    %rbx,%rcx
lea    (%rdx,%rdi,4),%rcx
setb   %dil
cmp    %rcx,%rbp
lea    (%rdx,%rax,4),%rcx
setb   %al
cmp    %rbx,%rcx
setb   %bl
mov    0xa(%rsp),%cl
test   %cl,0xb(%rsp)
jne    144f <horn_dt+0x33f>
and    0x9(%rsp),%r13b
jne    144f <horn_dt+0x33f>
and    %dil,%sil
mov    0x10(%rsp),%r13
jne    1227 <horn_dt+0x117>
mov    $0x0,%edx
and    %bl,%al
jne    1229 <horn_dt+0x119>
xor    %r13d,%r13d
mov    0x60(%rsp),%rax
nopw   0x0(%rax,%rax,1)
movups -0x4(%r9,%r13,4),%xmm2
movups (%r9,%r13,4),%xmm3
addps  %xmm2,%xmm3
movups -0x4(%r14,%r13,4),%xmm2
addps  %xmm3,%xmm2
movups (%r14,%r13,4),%xmm3
addps  %xmm2,%xmm3
cvtps2pd %xmm3,%xmm2
movhlps %xmm3,%xmm3
cvtps2pd %xmm3,%xmm3
mulpd  %xmm1,%xmm2
movups -0x4(%r15,%r13,4),%xmm4
movups (%r15,%r13,4),%xmm5
addps  %xmm4,%xmm5
movups -0x4(%r8,%r13,4),%xmm4
addps  %xmm5,%xmm4
movups (%r8,%r13,4),%xmm5
addps  %xmm4,%xmm5
cvtps2pd %xmm5,%xmm4
mulpd  %xmm1,%xmm3
movhlps %xmm5,%xmm5
cvtps2pd %xmm5,%xmm5
mulpd  %xmm1,%xmm4
subpd  %xmm4,%xmm2
mulpd  %xmm1,%xmm5
subpd  %xmm5,%xmm3
cvtpd2ps %xmm2,%xmm2
cvtpd2ps %xmm3,%xmm3
unpcklpd %xmm3,%xmm2
movupd %xmm2,(%r12,%r13,4)
add    $0x4,%r13
cmp    %r13,%rax
jne    13b0 <horn_dt+0x2a0>
mov    %rax,%rdx
cmp    0x30(%rsp),%rax
mov    0x10(%rsp),%r13
je     11e0 <horn_dt+0xd0>
jmp    1229 <horn_dt+0x119>
xor    %edx,%edx
mov    0x10(%rsp),%r13
jmp    1229 <horn_dt+0x119>
test   %ebx,%ebx
mov    0x18(%rsp),%r14
jle    147e <horn_dt+0x36e>
mov    0xc(%rsp),%eax
imul   %ebx,%eax
cltq
lea    (%r14,%rax,4),%rdi
mov    %ebx,%edx
shl    $0x2,%rdx
xor    %esi,%esi
call   1030 <memset@plt>
mov    0x58(%rsp),%rax
test   %eax,%eax
jle    152f <horn_dt+0x41f>
add    $0xffffffff,%ebx
movslq %ebx,%r8
mov    %eax,%r9d
lea    -0x1(%r9),%rcx
mov    %r9d,%eax
and    $0x3,%eax
cmp    $0x3,%rcx
jae    14a8 <horn_dt+0x398>
xor    %esi,%esi
jmp    1509 <horn_dt+0x3f9>
and    $0xfffffffc,%r9d
lea    (%r14,%r13,8),%rdi
lea    0x0(,%r8,4),%rbp
mov    %r13,%rbx
shl    $0x4,%rbx
lea    (%r14,%r13,4),%rcx
lea    0x0(,%r13,2),%rdx
add    %r13,%rdx
lea    (%r14,%rdx,4),%rdx
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movl   $0x0,(%r14,%rbp,1)
movl   $0x0,(%rcx,%rbp,1)
movl   $0x0,(%rdi,%rbp,1)
movl   $0x0,(%rdx,%rbp,1)
add    $0x4,%rsi
add    %rbx,%rbp
cmp    %rsi,%r9
jne    14e0 <horn_dt+0x3d0>
test   %rax,%rax
je     152f <horn_dt+0x41f>
imul   %r13,%rsi
add    %r8,%rsi
lea    (%r14,%rsi,4),%rcx
shl    $0x2,%r13
nopl   (%rax)
movl   $0x0,(%rcx)
add    %r13,%rcx
add    $0xffffffffffffffff,%rax
jne    1520 <horn_dt+0x410>
add    $0x88,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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