void horn_dx(image1, image2, dx, nx, ny) float *dx, *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;
      dx[index] = ((image1[index + 1] + image2[index + 1] +
                    image1[index + nx + 1] + image2[index + nx + 1]) /
                   4.0) -
                  ((image1[index] + image2[index] + image1[index + nx] +
                    image2[index + nx]) /
                   4.0);
    }
  }
  for (j = 0; j < nx; j++) {
    dx[(ny - 1) * nx + j] = 0.0;
  }
  for (i = 0; i < ny; i++) {
    dx[(i * nx) + (nx - 1)] = 0.0;
  }
}


x86 assembly (clang, -O0):
<horn_dx>:
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_dx+0x136>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1223 <horn_dx+0x123>
mov    -0x28(%rbp),%eax
imul   -0x1c(%rbp),%eax
add    -0x2c(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%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    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x1c(%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
add    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe52(%rip),%xmm1
divsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%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    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x1c(%rbp),%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_dx+0x36>
jmp    1228 <horn_dx+0x128>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    111e <horn_dx+0x1e>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1273 <horn_dx+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_dx+0x13d>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    12b2 <horn_dx+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_dx+0x17a>
pop    %rbp
ret

x86 assembly (clang, -O2):
<horn_dx>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x88,%rsp
mov    %ecx,%ebx
mov    %rdx,0x10(%rsp)
mov    %rsi,0x28(%rsp)
mov    %rdi,0x20(%rsp)
lea    -0x1(%r8),%eax
mov    %eax,0x4(%rsp)
movslq %ecx,%r13
mov    %r8,0x58(%rsp)
cmp    $0x2,%r8d
jl     146b <horn_dx+0x35b>
lea    -0x1(%r13),%eax
mov    0x4(%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),%r8
mov    0x28(%rsp),%rax
lea    0x4(%rax),%r15
mov    0x20(%rsp),%rcx
lea    0x4(%rcx),%r11
lea    0x1(%r13),%rdx
mov    %rdx,0x50(%rsp)
lea    (%rax,%r13,4),%r14
add    $0x4,%r14
lea    (%rcx,%r13,4),%r9
add    $0x4,%r9
xor    %r10d,%r10d
movsd  0xe53(%rip),%xmm0
movapd 0xe3b(%rip),%xmm1
mov    %rcx,0x48(%rsp)
mov    %rax,0x40(%rsp)
mov    0x10(%rsp),%r12
mov    %r13,0x8(%rsp)
mov    %ebx,0x3c(%rsp)
jmp    1214 <horn_dx+0x104>
nop
add    $0x1,%r10
add    %r8,%r12
add    %r8,%r15
add    %r8,%r11
add    %r8,%r14
add    %r8,%r9
add    %r8,0x40(%rsp)
add    %r8,0x48(%rsp)
add    %r13,0x50(%rsp)
cmp    0x80(%rsp),%r10
mov    0x3c(%rsp),%ebx
je     146b <horn_dx+0x35b>
cmp    $0x2,%ebx
jl     11e0 <horn_dx+0xd0>
mov    0x30(%rsp),%rax
cmp    $0x4,%eax
jae    12d0 <horn_dx+0x1c0>
xor    %eax,%eax
mov    0x30(%rsp),%rbx
sub    %rax,%rbx
lea    (%r12,%rax,4),%rcx
mov    0x40(%rsp),%rdx
lea    (%rdx,%rax,4),%rdi
mov    0x48(%rsp),%rdx
lea    (%rdx,%rax,4),%rsi
add    0x50(%rsp),%rax
mov    0x20(%rsp),%rdx
lea    (%rdx,%rax,4),%rdx
mov    0x28(%rsp),%rbp
lea    0x0(%rbp,%rax,4),%rbp
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movss  0x4(%rsi,%rax,4),%xmm2
addss  0x4(%rdi,%rax,4),%xmm2
addss  (%rdx,%rax,4),%xmm2
movss  (%rsi,%rax,4),%xmm3
addss  0x0(%rbp,%rax,4),%xmm2
cvtss2sd %xmm2,%xmm2
addss  (%rdi,%rax,4),%xmm3
addss  -0x4(%rdx,%rax,4),%xmm3
mulsd  %xmm0,%xmm2
addss  -0x4(%rbp,%rax,4),%xmm3
cvtss2sd %xmm3,%xmm3
mulsd  %xmm0,%xmm3
subsd  %xmm3,%xmm2
cvtsd2ss %xmm2,%xmm2
movss  %xmm2,(%rcx,%rax,4)
add    $0x1,%rax
cmp    %rax,%rbx
jne    1270 <horn_dx+0x160>
jmp    11e0 <horn_dx+0xd0>
nopl   0x0(%rax,%rax,1)
mov    %r10,%rcx
imul   %r13,%rcx
mov    %r13,%rdx
lea    (%rcx,%rax,1),%r13
lea    (%rcx,%rdx,1),%rax
mov    %rax,0x18(%rsp)
mov    0x70(%rsp),%rsi
add    %rcx,%rsi
mov    0x68(%rsp),%rdi
lea    (%rdi,%rcx,1),%rdx
mov    %rdx,0x78(%rsp)
mov    0x10(%rsp),%rbx
lea    (%rbx,%rcx,4),%rbp
lea    (%rbx,%r13,4),%rbx
mov    0x20(%rsp),%rdi
lea    (%rdi,%rax,4),%r13
lea    (%rdi,%rsi,4),%rax
cmp    %rax,%rbp
lea    (%rdi,%rdx,4),%rax
setb   0x3(%rsp)
cmp    %rbx,%r13
setb   0x2(%rsp)
cmp    %rax,%rbp
lea    (%rdi,%rcx,4),%rax
setb   %dil
cmp    %rbx,%rax
mov    0x28(%rsp),%r13
lea    0x0(%r13,%rsi,4),%rax
setb   %sil
cmp    %rax,%rbp
mov    0x18(%rsp),%rax
lea    0x0(%r13,%rax,4),%rax
setb   %dl
cmp    %rbx,%rax
mov    0x78(%rsp),%rax
lea    0x0(%r13,%rax,4),%rax
setb   0x18(%rsp)
cmp    %rax,%rbp
lea    0x0(%r13,%rcx,4),%rax
setb   %cl
cmp    %rbx,%rax
setb   %bl
mov    0x2(%rsp),%al
test   %al,0x3(%rsp)
jne    145f <horn_dx+0x34f>
and    %sil,%dil
jne    145f <horn_dx+0x34f>
and    0x18(%rsp),%dl
mov    0x8(%rsp),%r13
jne    1227 <horn_dx+0x117>
mov    $0x0,%eax
and    %bl,%cl
jne    1229 <horn_dx+0x119>
xor    %r13d,%r13d
mov    0x60(%rsp),%rcx
cs nopw 0x0(%rax,%rax,1)
nop
movups -0x4(%r11,%r13,4),%xmm2
movups (%r11,%r13,4),%xmm3
movups -0x4(%r15,%r13,4),%xmm4
addps  %xmm2,%xmm4
movups (%r15,%r13,4),%xmm2
addps  %xmm3,%xmm2
movups -0x4(%r9,%r13,4),%xmm3
addps  %xmm4,%xmm3
movups (%r9,%r13,4),%xmm4
addps  %xmm2,%xmm4
movups (%r14,%r13,4),%xmm2
addps  %xmm4,%xmm2
cvtps2pd %xmm2,%xmm4
movups -0x4(%r14,%r13,4),%xmm5
movhlps %xmm2,%xmm2
cvtps2pd %xmm2,%xmm2
addps  %xmm3,%xmm5
mulpd  %xmm1,%xmm4
cvtps2pd %xmm5,%xmm3
mulpd  %xmm1,%xmm2
movhlps %xmm5,%xmm5
cvtps2pd %xmm5,%xmm5
mulpd  %xmm1,%xmm3
subpd  %xmm3,%xmm4
mulpd  %xmm1,%xmm5
subpd  %xmm5,%xmm2
cvtpd2ps %xmm4,%xmm3
cvtpd2ps %xmm2,%xmm2
unpcklpd %xmm2,%xmm3
movupd %xmm3,(%r12,%r13,4)
add    $0x4,%r13
cmp    %r13,%rcx
jne    13c0 <horn_dx+0x2b0>
mov    %rcx,%rax
cmp    0x30(%rsp),%rcx
mov    0x8(%rsp),%r13
je     11e0 <horn_dx+0xd0>
jmp    1229 <horn_dx+0x119>
xor    %eax,%eax
mov    0x8(%rsp),%r13
jmp    1229 <horn_dx+0x119>
test   %ebx,%ebx
mov    0x10(%rsp),%r14
jle    148e <horn_dx+0x37e>
mov    0x4(%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    153f <horn_dx+0x42f>
add    $0xffffffff,%ebx
movslq %ebx,%r8
mov    %eax,%r9d
lea    -0x1(%r9),%rcx
mov    %r9d,%eax
and    $0x3,%eax
cmp    $0x3,%rcx
jae    14b8 <horn_dx+0x3a8>
xor    %esi,%esi
jmp    1519 <horn_dx+0x409>
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    14f0 <horn_dx+0x3e0>
test   %rax,%rax
je     153f <horn_dx+0x42f>
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    1530 <horn_dx+0x420>
add    $0x88,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @horn_dx(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** %6
%32 = load i32, i32* %11
%33 = add nsw i32 %32, 1
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds float, float* %31, i64 %34
%36 = load float, float* %35
%37 = load float*, float** %7
%38 = load i32, i32* %11
%39 = add nsw i32 %38, 1
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds float, float* %37, i64 %40
%42 = load float, float* %41
%43 = fadd float %36, %42
%44 = load float*, float** %6
%45 = load i32, i32* %11
%46 = load i32, i32* %9
%47 = add nsw i32 %45, %46
%48 = add nsw i32 %47, 1
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds float, float* %44, i64 %49
%51 = load float, float* %50
%52 = fadd float %43, %51
%53 = load float*, float** %7
%54 = load i32, i32* %11
%55 = load i32, i32* %9
%56 = add nsw i32 %54, %55
%57 = add nsw i32 %56, 1
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds float, float* %53, i64 %58
%60 = load float, float* %59
%61 = fadd float %52, %60
%62 = fpext float %61 to double
%63 = fdiv double %62, 4.000000e+00
%64 = load float*, float** %6
%65 = load i32, i32* %11
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds float, float* %64, i64 %66
%68 = load float, float* %67
%69 = load float*, float** %7
%70 = load i32, i32* %11
%71 = sext i32 %70 to i64
%72 = getelementptr inbounds float, float* %69, i64 %71
%73 = load float, float* %72
%74 = fadd float %68, %73
%75 = load float*, float** %6
%76 = load i32, i32* %11
%77 = load i32, i32* %9
%78 = add nsw i32 %76, %77
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds float, float* %75, i64 %79
%81 = load float, float* %80
%82 = fadd float %74, %81
%83 = load float*, float** %7
%84 = load i32, i32* %11
%85 = load i32, i32* %9
%86 = add nsw i32 %84, %85
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds float, float* %83, i64 %87
%89 = load float, float* %88
%90 = fadd float %82, %89
%91 = fpext float %90 to double
%92 = fdiv double %91, 4.000000e+00
%93 = fsub double %63, %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_dx(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* %0, i64 %26
%28 = add i64 %15, %22
%29 = getelementptr float, float* %0, i64 %28
%30 = getelementptr float, float* %0, i64 %22
%31 = add i64 %16, %22
%32 = getelementptr float, float* %0, i64 %31
%33 = getelementptr float, float* %1, i64 %26
%34 = getelementptr float, float* %1, i64 %28
%35 = getelementptr float, float* %1, i64 %22
%36 = getelementptr float, float* %1, 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 = add nsw i64 %57, 1
%59 = getelementptr inbounds float, float* %0, i64 %58
%60 = bitcast float* %59 to <4 x float>*
%61 = load <4 x float>, <4 x float>* %60
%62 = getelementptr inbounds float, float* %1, i64 %58
%63 = bitcast float* %62 to <4 x float>*
%64 = load <4 x float>, <4 x float>* %63
%65 = fadd <4 x float> %61, %64
%66 = add nsw i64 %57, %11
%67 = add nsw i64 %66, 1
%68 = getelementptr inbounds float, float* %0, i64 %67
%69 = bitcast float* %68 to <4 x float>*
%70 = load <4 x float>, <4 x float>* %69
%71 = fadd <4 x float> %65, %70
%72 = getelementptr inbounds float, float* %1, i64 %67
%73 = bitcast float* %72 to <4 x float>*
%74 = load <4 x float>, <4 x float>* %73
%75 = fadd <4 x float> %71, %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* %1, i64 %57
%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* %1, i64 %66
%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 = add nsw i64 %115, 1
%117 = getelementptr inbounds float, float* %0, i64 %116
%118 = load float, float* %117
%119 = getelementptr inbounds float, float* %1, i64 %116
%120 = load float, float* %119
%121 = fadd float %118, %120
%122 = add nsw i64 %115, %11
%123 = add nsw i64 %122, 1
%124 = getelementptr inbounds float, float* %0, i64 %123
%125 = load float, float* %124
%126 = fadd float %121, %125
%127 = getelementptr inbounds float, float* %1, i64 %123
%128 = load float, float* %127
%129 = fadd float %126, %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* %1, i64 %115
%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* %1, i64 %122
%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
}