void sift_compute_gradient(const float *im, float *im_x, float *im_y, int w,
                           int h) {
  const float *p_in;
  const float *p_in_p;
  const float *p_in_m;
  float *p_out;
  p_in = &im[1];
  p_out = &im_y[1];
  p_in_p = p_in + 1;
  p_in_m = p_in - 1;
  for (int i = 1; i < h * w - 1; i++) {
    *p_out = (*p_in_p - *p_in_m) * 0.5;
    p_out++;
    p_in_m++;
    p_in_p++;
  }
  for (int i = 0; i < h; i++) {
    im_y[i * w] = im[i * w + 1] - im[i * w];
    im_y[i * w + (w - 1)] = im[i * w + (w - 1)] - im[i * w + (w - 2)];
  }
  p_in = &im[w];
  p_out = &im_x[w];
  p_in_p = p_in + w;
  p_in_m = p_in - w;
  for (int i = w; i < h * w - w; i++) {
    *p_out = (*p_in_p - *p_in_m) * 0.5;
    p_out++;
    p_in_m++;
    p_in_p++;
  }
  for (int i = 0; i < w; i++) {
    im_x[i] = im[i + w] - im[i];
    im_x[(h - 1) * w + i] = im[(h - 1) * w + i] - im[(h - 1 - 1) * w + i];
  }
}


x86 assembly (clang, -O0):
<sift_compute_gradient>:
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)
mov    -0x8(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x28(%rbp)
mov    -0x18(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x40(%rbp)
mov    -0x28(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x30(%rbp)
mov    -0x28(%rbp),%rax
add    $0xfffffffffffffffc,%rax
mov    %rax,-0x38(%rbp)
movl   $0x1,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    11c1 <sift_compute_gradient+0xc1>
mov    -0x30(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x38(%rbp),%rax
subss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe81(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x40(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x40(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x40(%rbp)
mov    -0x38(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x38(%rbp)
mov    -0x30(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x30(%rbp)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    114e <sift_compute_gradient+0x4e>
movl   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    126f <sift_compute_gradient+0x16f>
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
mov    -0x1c(%rbp),%edx
sub    $0x1,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
mov    -0x1c(%rbp),%edx
sub    $0x2,%edx
add    %edx,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
mov    -0x1c(%rbp),%edx
sub    $0x1,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    11c8 <sift_compute_gradient+0xc8>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x28(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x1c(%rbp),%rdx
xor    %ecx,%ecx
sub    %rdx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
mov    -0x1c(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
cmp    %ecx,%eax
jge    1339 <sift_compute_gradient+0x239>
mov    -0x30(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x38(%rbp),%rax
subss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xd09(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x40(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x40(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x40(%rbp)
mov    -0x38(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x38(%rbp)
mov    -0x30(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x30(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    12c6 <sift_compute_gradient+0x1c6>
movl   $0x0,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    13d4 <sift_compute_gradient+0x2d4>
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x50(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x50(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
sub    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x50(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x50(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1340 <sift_compute_gradient+0x240>
pop    %rbp
ret

x86 assembly (clang, -O2):
<sift_compute_gradient>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %r8d,%r9d
imul   %ecx,%r9d
cmp    $0x3,%r9d
jl     132c <sift_compute_gradient+0x22c>
lea    0x8(%rdi),%r10
lea    -0x3(%r9),%r15d
mov    $0x1,%eax
cmp    $0x7,%r15d
jb     1154 <sift_compute_gradient+0x54>
lea    0x4(%rdx),%rbp
lea    (%rdi,%r15,4),%rbx
add    $0xc,%rbx
cmp    %rbx,%rbp
jae    124f <sift_compute_gradient+0x14f>
lea    (%rdx,%r15,4),%rbp
add    $0x8,%rbp
cmp    %rdi,%rbp
jbe    124f <sift_compute_gradient+0x14f>
mov    %rdx,%r11
mov    %rdi,%r12
mov    %eax,%ebp
not    %ebp
add    %r9d,%ebp
mov    %r9d,%r14d
sub    %eax,%r14d
add    $0xfffffffe,%r14d
test   $0x3,%bpl
je     11bc <sift_compute_gradient+0xbc>
mov    %eax,%ebx
not    %bl
add    %r9b,%bl
movzbl %bl,%r15d
and    $0x3,%r15d
shl    $0x2,%r15
xor    %ebp,%ebp
movss  0xe82(%rip),%xmm0
xchg   %ax,%ax
movss  (%r10,%rbp,1),%xmm1
subss  (%r12,%rbp,1),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,0x4(%r11,%rbp,1)
add    $0x1,%eax
add    $0x4,%rbp
cmp    %ebp,%r15d
jne    1190 <sift_compute_gradient+0x90>
add    %rbp,%r11
add    %rbp,%r12
add    %rbp,%r10
cmp    $0x3,%r14d
jb     132c <sift_compute_gradient+0x22c>
not    %eax
add    %r9d,%eax
xor    %ebp,%ebp
movss  0xe3b(%rip),%xmm0
cs nopw 0x0(%rax,%rax,1)
nop
movss  (%r10,%rbp,4),%xmm1
subss  (%r12,%rbp,4),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,0x4(%r11,%rbp,4)
movss  0x4(%r10,%rbp,4),%xmm1
subss  0x4(%r12,%rbp,4),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,0x8(%r11,%rbp,4)
movss  0x8(%r10,%rbp,4),%xmm1
subss  0x8(%r12,%rbp,4),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,0xc(%r11,%rbp,4)
movss  0xc(%r10,%rbp,4),%xmm1
subss  0xc(%r12,%rbp,4),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,0x10(%r11,%rbp,4)
add    $0x4,%rbp
cmp    %ebp,%eax
jne    11e0 <sift_compute_gradient+0xe0>
jmp    132c <sift_compute_gradient+0x22c>
add    $0x1,%r15
mov    %r15,%r14
and    $0xfffffffffffffff8,%r14
lea    -0x8(%r14),%rax
mov    %rax,%r11
shr    $0x3,%r11
add    $0x1,%r11
test   %rax,%rax
je     171d <sift_compute_gradient+0x61d>
mov    %r11,%rbp
and    $0xfffffffffffffffe,%rbp
xor    %eax,%eax
movaps 0xd7e(%rip),%xmm0
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movups (%rdi,%rax,4),%xmm1
movups 0x8(%rdi,%rax,4),%xmm2
subps  %xmm1,%xmm2
movups 0x10(%rdi,%rax,4),%xmm1
movups 0x18(%rdi,%rax,4),%xmm3
subps  %xmm1,%xmm3
mulps  %xmm0,%xmm2
mulps  %xmm0,%xmm3
movups %xmm2,0x4(%rdx,%rax,4)
movups %xmm3,0x14(%rdx,%rax,4)
movups 0x20(%rdi,%rax,4),%xmm1
movups 0x28(%rdi,%rax,4),%xmm2
subps  %xmm1,%xmm2
movups 0x30(%rdi,%rax,4),%xmm1
movups 0x38(%rdi,%rax,4),%xmm3
subps  %xmm1,%xmm3
mulps  %xmm0,%xmm2
mulps  %xmm0,%xmm3
movups %xmm2,0x24(%rdx,%rax,4)
movups %xmm3,0x34(%rdx,%rax,4)
add    $0x10,%rax
add    $0xfffffffffffffffe,%rbp
jne    1290 <sift_compute_gradient+0x190>
test   $0x1,%r11b
je     1323 <sift_compute_gradient+0x223>
movups (%rdi,%rax,4),%xmm0
movups 0x8(%rdi,%rax,4),%xmm1
subps  %xmm0,%xmm1
movups 0x10(%rdi,%rax,4),%xmm0
movups 0x18(%rdi,%rax,4),%xmm2
subps  %xmm0,%xmm2
movaps 0xced(%rip),%xmm0
mulps  %xmm0,%xmm1
mulps  %xmm0,%xmm2
movups %xmm1,0x4(%rdx,%rax,4)
movups %xmm2,0x14(%rdx,%rax,4)
cmp    %r14,%r15
jne    16f3 <sift_compute_gradient+0x5f3>
movslq %ecx,%r10
test   %r8d,%r8d
jle    1443 <sift_compute_gradient+0x343>
mov    %r8d,%r15d
xor    %eax,%eax
cmp    $0x4,%r8d
jb     13fc <sift_compute_gradient+0x2fc>
cmp    $0x1,%ecx
jne    13fc <sift_compute_gradient+0x2fc>
lea    (%rdx,%r15,4),%r11
lea    -0x1(%r15),%rbp
imul   %r10,%rbp
lea    (%rdx,%rbp,4),%rbp
cmp    %rdx,%rbp
mov    %rbp,%rbx
cmova  %rdx,%rbx
lea    -0x4(%rdi),%r14
lea    0x4(%rdi,%r15,4),%rax
cmovb  %rdx,%rbp
add    $0x4,%rbp
cmp    %rdx,%rax
seta   %r12b
cmp    %r11,%r14
setb   %r13b
cmp    %rax,%rbx
setb   %bl
cmp    %rbp,%r14
setb   %r11b
xor    %eax,%eax
test   %r13b,%r12b
jne    13fc <sift_compute_gradient+0x2fc>
and    %r11b,%bl
jne    13fc <sift_compute_gradient+0x2fc>
mov    %r15d,%eax
and    $0xfffffffc,%eax
mov    %r10,%r11
shl    $0x4,%r11
lea    -0x2(%rcx),%ebp
lea    0x0(,%rcx,4),%r14d
xor    %r13d,%r13d
mov    %rax,%r12
movups (%rdi,%r13,1),%xmm0
movups 0x4(%rdi,%r13,1),%xmm1
subps  %xmm0,%xmm1
movups %xmm1,(%rdx,%r13,1)
lea    0x1(%rbp),%ebx
movslq %ebx,%rbx
movups (%rdi,%rbx,4),%xmm0
movslq %ebp,%rbp
movups (%rdi,%rbp,4),%xmm1
subps  %xmm1,%xmm0
movups %xmm0,(%rdx,%rbx,4)
add    %r11,%r13
add    %r14d,%ebp
add    $0xfffffffffffffffc,%r12
jne    13c0 <sift_compute_gradient+0x2c0>
cmp    %r15,%rax
je     1443 <sift_compute_gradient+0x343>
lea    0x1(%rax),%ebp
imul   %ecx,%ebp
add    $0xfffffffe,%ebp
sub    %rax,%r15
imul   %r10,%rax
nopl   0x0(%rax)
movss  0x4(%rdi,%rax,4),%xmm0
subss  (%rdi,%rax,4),%xmm0
movss  %xmm0,(%rdx,%rax,4)
lea    0x1(%rbp),%ebx
movslq %ebx,%rbx
movss  (%rdi,%rbx,4),%xmm0
movslq %ebp,%rbp
subss  (%rdi,%rbp,4),%xmm0
movss  %xmm0,(%rdx,%rbx,4)
add    %ecx,%ebp
add    %r10,%rax
add    $0xffffffffffffffff,%r15
jne    1410 <sift_compute_gradient+0x310>
sub    %ecx,%r9d
cmp    %ecx,%r9d
jle    168a <sift_compute_gradient+0x58a>
lea    (%rdi,%r10,4),%rax
lea    (%rax,%r10,4),%rdx
lea    (%rsi,%r10,4),%rax
lea    -0x2(%r8),%r15d
imul   %ecx,%r15d
add    $0xffffffff,%r15d
cmp    $0x7,%r15d
jae    1543 <sift_compute_gradient+0x443>
mov    %ecx,%r14d
mov    %rdi,%rbx
mov    %r9d,%r15d
sub    %r14d,%r15d
mov    %r14d,%r11d
not    %r11d
add    %r9d,%r11d
and    $0x3,%r15d
je     14c7 <sift_compute_gradient+0x3c7>
neg    %r15d
xor    %ebp,%ebp
movss  0xb77(%rip),%xmm0
nopl   0x0(%rax)
movss  (%rdx),%xmm1
subss  (%rbx),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,(%rax)
add    $0x4,%rax
add    $0x4,%rbx
add    $0x4,%rdx
add    $0xffffffff,%ebp
cmp    %ebp,%r15d
jne    14a0 <sift_compute_gradient+0x3a0>
sub    %ebp,%r14d
cmp    $0x3,%r11d
jb     168a <sift_compute_gradient+0x58a>
sub    %r14d,%r9d
xor    %ebp,%ebp
movss  0xb32(%rip),%xmm0
xchg   %ax,%ax
movss  (%rdx,%rbp,4),%xmm1
subss  (%rbx,%rbp,4),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,(%rax,%rbp,4)
movss  0x4(%rdx,%rbp,4),%xmm1
subss  0x4(%rbx,%rbp,4),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,0x4(%rax,%rbp,4)
movss  0x8(%rdx,%rbp,4),%xmm1
subss  0x8(%rbx,%rbp,4),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,0x8(%rax,%rbp,4)
movss  0xc(%rdx,%rbp,4),%xmm1
subss  0xc(%rbx,%rbp,4),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,0xc(%rax,%rbp,4)
add    $0x4,%rbp
cmp    %ebp,%r9d
jne    14e0 <sift_compute_gradient+0x3e0>
jmp    168a <sift_compute_gradient+0x58a>
lea    (%r10,%r15,1),%rbp
lea    (%rsi,%rbp,4),%rbp
add    $0x4,%rbp
lea    (%rdi,%r15,4),%r11
add    $0x4,%r11
lea    (%r15,%r10,2),%rbx
lea    (%rdi,%rbx,4),%rbx
add    $0x4,%rbx
cmp    %r11,%rax
setb   %r14b
cmp    %rdi,%rbp
seta   %r12b
cmp    %rbx,%rax
setb   %bl
cmp    %rbp,%rdx
setb   %r11b
test   %r12b,%r14b
jne    1471 <sift_compute_gradient+0x371>
and    %r11b,%bl
jne    1471 <sift_compute_gradient+0x371>
add    $0x1,%r15
mov    %r15,%rbp
and    $0xfffffffffffffff8,%rbp
mov    %rbp,-0x8(%rsp)
add    $0xfffffffffffffff8,%rbp
mov    %rbp,%r14
shr    $0x3,%r14
add    $0x1,%r14
test   %rbp,%rbp
je     172e <sift_compute_gradient+0x62e>
mov    %r14,%rbx
and    $0xfffffffffffffffe,%rbx
lea    (%rsi,%r10,4),%r12
add    $0x30,%r12
lea    0x30(%rdi),%r13
lea    (%rdi,%r10,8),%r11
add    $0x30,%r11
xor    %ebp,%ebp
movaps 0xa24(%rip),%xmm0
nopl   0x0(%rax)
movups -0x30(%r11,%rbp,4),%xmm1
movups -0x20(%r11,%rbp,4),%xmm2
movups -0x30(%r13,%rbp,4),%xmm3
subps  %xmm3,%xmm1
movups -0x20(%r13,%rbp,4),%xmm3
subps  %xmm3,%xmm2
mulps  %xmm0,%xmm1
mulps  %xmm0,%xmm2
movups %xmm1,-0x30(%r12,%rbp,4)
movups %xmm2,-0x20(%r12,%rbp,4)
movups -0x10(%r11,%rbp,4),%xmm1
movups (%r11,%rbp,4),%xmm2
movups -0x10(%r13,%rbp,4),%xmm3
subps  %xmm3,%xmm1
movups 0x0(%r13,%rbp,4),%xmm3
subps  %xmm3,%xmm2
mulps  %xmm0,%xmm1
mulps  %xmm0,%xmm2
movups %xmm1,-0x10(%r12,%rbp,4)
movups %xmm2,(%r12,%rbp,4)
add    $0x10,%rbp
add    $0xfffffffffffffffe,%rbx
jne    15e0 <sift_compute_gradient+0x4e0>
test   $0x1,%r14b
je     167c <sift_compute_gradient+0x57c>
movups (%rdx,%rbp,4),%xmm0
movups 0x10(%rdx,%rbp,4),%xmm1
movups (%rdi,%rbp,4),%xmm2
subps  %xmm2,%xmm0
movups 0x10(%rdi,%rbp,4),%xmm2
subps  %xmm2,%xmm1
movaps 0x993(%rip),%xmm2
mulps  %xmm2,%xmm0
mulps  %xmm2,%xmm1
movups %xmm0,(%rax,%rbp,4)
movups %xmm1,0x10(%rax,%rbp,4)
mov    -0x8(%rsp),%rbp
cmp    %rbp,%r15
jne    1708 <sift_compute_gradient+0x608>
test   %ecx,%ecx
jle    16e8 <sift_compute_gradient+0x5e8>
lea    -0x1(%r8),%eax
imul   %ecx,%eax
add    $0xfffffffe,%r8d
imul   %ecx,%r8d
cltq
movslq %r8d,%rdx
mov    %ecx,%r8d
lea    (%rdi,%rdx,4),%rcx
lea    (%rdi,%rax,4),%rdx
lea    (%rsi,%rax,4),%rbp
lea    (%rdi,%r10,4),%rbx
xor    %eax,%eax
nopw   0x0(%rax,%rax,1)
movss  (%rbx,%rax,4),%xmm0
subss  (%rdi,%rax,4),%xmm0
movss  %xmm0,(%rsi,%rax,4)
movss  (%rdx,%rax,4),%xmm0
subss  (%rcx,%rax,4),%xmm0
movss  %xmm0,0x0(%rbp,%rax,4)
add    $0x1,%rax
cmp    %rax,%r8
jne    16c0 <sift_compute_gradient+0x5c0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
lea    (%rdx,%r14,4),%r11
lea    0x1(%r14),%eax
lea    (%rdi,%r14,4),%r12
lea    (%r10,%r14,4),%r10
jmp    115a <sift_compute_gradient+0x5a>
lea    (%rcx,%rbp,1),%r14d
lea    (%rax,%rbp,4),%rax
lea    (%rdi,%rbp,4),%rbx
lea    (%rdx,%rbp,4),%rdx
jmp    1477 <sift_compute_gradient+0x377>
xor    %eax,%eax
test   $0x1,%r11b
jne    12f3 <sift_compute_gradient+0x1f3>
jmp    1323 <sift_compute_gradient+0x223>
xor    %ebp,%ebp
test   $0x1,%r14b
jne    164e <sift_compute_gradient+0x54e>
jmp    167c <sift_compute_gradient+0x57c>

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

LLVM IR (-O2):
define void @sift_compute_gradient(float* nocapture %0, float* nocapture writeonly %1, float* nocapture writeonly %2, i32 %3, i32 %4) {
%6 = mul i32 %4, %3
%7 = add i32 %6, -1
%8 = icmp sgt i32 %6, 2
br i1 %8, label %9, label %149
9: ; preds = %5
%10 = getelementptr inbounds float, float* %0, i64 2
%11 = add i32 %6, -3
%12 = zext i32 %11 to i64
%13 = add nuw nsw i64 %12, 1
%14 = icmp ult i32 %11, 7
br i1 %14, label %116, label %15
15: ; preds = %9
%16 = getelementptr float, float* %2, i64 1
%17 = add i32 %6, -3
%18 = zext i32 %17 to i64
%19 = add nuw nsw i64 %18, 2
%20 = getelementptr float, float* %2, i64 %19
%21 = add nuw nsw i64 %18, 3
%22 = getelementptr float, float* %0, i64 %21
%23 = icmp ult float* %16, %22
%24 = icmp ugt float* %20, %0
%25 = and i1 %23, %24
br i1 %25, label %116, label %26
26: ; preds = %15
%27 = and i64 %13, 8589934584
%28 = getelementptr float, float* %2, i64 %27
%29 = trunc i64 %27 to i32
%30 = or i32 %29, 1
%31 = getelementptr float, float* %0, i64 %27
%32 = getelementptr float, float* %10, i64 %27
%33 = add nsw i64 %27, -8
%34 = lshr exact i64 %33, 3
%35 = add nuw nsw i64 %34, 1
%36 = and i64 %35, 1
%37 = icmp eq i64 %33, 0
br i1 %37, label %89, label %38
38: ; preds = %26
%39 = and i64 %35, 4611686018427387902
%40 = getelementptr float, float* %2, i64 1
%41 = getelementptr float, float* %2, i64 1
br label %42
42: ; preds = %42, %38
%43 = phi i64 [ 0, %38 ], [ %86, %42 ]
%44 = phi i64 [ 0, %38 ], [ %87, %42 ]
%45 = getelementptr float, float* %0, i64 %43
%46 = getelementptr float, float* %10, i64 %43
%47 = getelementptr inbounds float, float* %40, i64 %43
%48 = bitcast float* %46 to <4 x float>*
%49 = load <4 x float>, <4 x float>* %48
%50 = getelementptr float, float* %46, i64 4
%51 = bitcast float* %50 to <4 x float>*
%52 = load <4 x float>, <4 x float>* %51
%53 = bitcast float* %45 to <4 x float>*
%54 = load <4 x float>, <4 x float>* %53
%55 = getelementptr float, float* %45, i64 4
%56 = bitcast float* %55 to <4 x float>*
%57 = load <4 x float>, <4 x float>* %56
%58 = fsub <4 x float> %49, %54
%59 = fsub <4 x float> %52, %57
%60 = fmul <4 x float> %58, <float 5.000000e-01, float 5.000000e-01, float 5.000000e-01, float 5.000000e-01>
%61 = fmul <4 x float> %59, <float 5.000000e-01, float 5.000000e-01, float 5.000000e-01, float 5.000000e-01>
%62 = bitcast float* %47 to <4 x float>*
store <4 x float> %60, <4 x float>* %62
%63 = getelementptr inbounds float, float* %47, i64 4
%64 = bitcast float* %63 to <4 x float>*
store <4 x float> %61, <4 x float>* %64
%65 = or i64 %43, 8
%66 = getelementptr float, float* %0, i64 %65
%67 = getelementptr float, float* %10, i64 %65
%68 = getelementptr inbounds float, float* %41, i64 %65
%69 = bitcast float* %67 to <4 x float>*
%70 = load <4 x float>, <4 x float>* %69
%71 = getelementptr float, float* %67, i64 4
%72 = bitcast float* %71 to <4 x float>*
%73 = load <4 x float>, <4 x float>* %72
%74 = bitcast float* %66 to <4 x float>*
%75 = load <4 x float>, <4 x float>* %74
%76 = getelementptr float, float* %66, i64 4
%77 = bitcast float* %76 to <4 x float>*
%78 = load <4 x float>, <4 x float>* %77
%79 = fsub <4 x float> %70, %75
%80 = fsub <4 x float> %73, %78
%81 = fmul <4 x float> %79, <float 5.000000e-01, float 5.000000e-01, float 5.000000e-01, float 5.000000e-01>
%82 = fmul <4 x float> %80, <float 5.000000e-01, float 5.000000e-01, float 5.000000e-01, float 5.000000e-01>
%83 = bitcast float* %68 to <4 x float>*
store <4 x float> %81, <4 x float>* %83
%84 = getelementptr inbounds float, float* %68, i64 4
%85 = bitcast float* %84 to <4 x float>*
store <4 x float> %82, <4 x float>* %85
%86 = add nuw i64 %43, 16
%87 = add i64 %44, 2
%88 = icmp eq i64 %87, %39
br i1 %88, label %89, label %42
89: ; preds = %42, %26
%90 = phi i64 [ 0, %26 ], [ %86, %42 ]
%91 = icmp eq i64 %36, 0
br i1 %91, label %114, label %92
92: ; preds = %89
%93 = getelementptr float, float* %2, i64 1
%94 = getelementptr float, float* %0, i64 %90
%95 = getelementptr float, float* %10, i64 %90
%96 = getelementptr inbounds float, float* %93, i64 %90
%97 = bitcast float* %95 to <4 x float>*
%98 = load <4 x float>, <4 x float>* %97
%99 = getelementptr float, float* %95, i64 4
%100 = bitcast float* %99 to <4 x float>*
%101 = load <4 x float>, <4 x float>* %100
%102 = bitcast float* %94 to <4 x float>*
%103 = load <4 x float>, <4 x float>* %102
%104 = getelementptr float, float* %94, i64 4
%105 = bitcast float* %104 to <4 x float>*
%106 = load <4 x float>, <4 x float>* %105
%107 = fsub <4 x float> %98, %103
%108 = fsub <4 x float> %101, %106
%109 = fmul <4 x float> %107, <float 5.000000e-01, float 5.000000e-01, float 5.000000e-01, float 5.000000e-01>
%110 = fmul <4 x float> %108, <float 5.000000e-01, float 5.000000e-01, float 5.000000e-01, float 5.000000e-01>
%111 = bitcast float* %96 to <4 x float>*
store <4 x float> %109, <4 x float>* %111
%112 = getelementptr inbounds float, float* %96, i64 4
%113 = bitcast float* %112 to <4 x float>*
store <4 x float> %110, <4 x float>* %113
br label %114
114: ; preds = %89, %92
%115 = icmp eq i64 %13, %27
br i1 %115, label %149, label %116
116: ; preds = %15, %9, %114
%117 = phi float* [ %2, %15 ], [ %2, %9 ], [ %28, %114 ]
%118 = phi i32 [ 1, %15 ], [ 1, %9 ], [ %30, %114 ]
%119 = phi float* [ %0, %15 ], [ %0, %9 ], [ %31, %114 ]
%120 = phi float* [ %10, %15 ], [ %10, %9 ], [ %32, %114 ]
%121 = xor i32 %118, -1
%122 = add i32 %6, %121
%123 = add i32 %6, -2
%124 = sub i32 %123, %118
%125 = and i32 %122, 3
%126 = icmp eq i32 %125, 0
br i1 %126, label %143, label %127
127: ; preds = %116, %127
%128 = phi float* [ %133, %127 ], [ %117, %116 ]
%129 = phi i32 [ %140, %127 ], [ %118, %116 ]
%130 = phi float* [ %138, %127 ], [ %119, %116 ]
%131 = phi float* [ %139, %127 ], [ %120, %116 ]
%132 = phi i32 [ %141, %127 ], [ 0, %116 ]
%133 = getelementptr inbounds float, float* %128, i64 1
%134 = load float, float* %131
%135 = load float, float* %130
%136 = fsub float %134, %135
%137 = fmul float %136, 5.000000e-01
store float %137, float* %133
%138 = getelementptr inbounds float, float* %130, i64 1
%139 = getelementptr inbounds float, float* %131, i64 1
%140 = add nuw nsw i32 %129, 1
%141 = add i32 %132, 1
%142 = icmp eq i32 %141, %125
br i1 %142, label %143, label %127
143: ; preds = %127, %116
%144 = phi float* [ %117, %116 ], [ %133, %127 ]
%145 = phi i32 [ %118, %116 ], [ %140, %127 ]
%146 = phi float* [ %119, %116 ], [ %138, %127 ]
%147 = phi float* [ %120, %116 ], [ %139, %127 ]
%148 = icmp ult i32 %124, 3
br i1 %148, label %149, label %213
149: ; preds = %143, %213, %114, %5
%150 = icmp sgt i32 %4, 0
%151 = sext i32 %3 to i64
br i1 %150, label %152, label %248
152: ; preds = %149
%153 = zext i32 %4 to i64
%154 = icmp ugt i32 %4, 3
%155 = icmp eq i32 %3, 1
%156 = and i1 %154, %155
br i1 %156, label %157, label %211
157: ; preds = %152
%158 = getelementptr float, float* %2, i64 %153
%159 = getelementptr float, float* %0, i64 -1
%160 = add nuw nsw i64 %153, 1
%161 = getelementptr float, float* %0, i64 %160
%162 = add nsw i64 %153, -1
%163 = mul i64 %162, %151
%164 = getelementptr float, float* %2, i64 %163
%165 = icmp ugt float* %164, %2
%166 = select i1 %165, float* %2, float* %164
%167 = icmp ult float* %164, %2
%168 = select i1 %167, float* %2, float* %164
%169 = getelementptr float, float* %168, i64 1
%170 = icmp ugt float* %161, %2
%171 = icmp ult float* %159, %158
%172 = and i1 %170, %171
%173 = icmp ult float* %166, %161
%174 = icmp ult float* %159, %169
%175 = and i1 %173, %174
%176 = or i1 %172, %175
br i1 %176, label %211, label %177
177: ; preds = %157
%178 = and i64 %153, 4294967292
br label %179
179: ; preds = %179, %177
%180 = phi i64 [ 0, %177 ], [ %207, %179 ]
%181 = mul nsw i64 %180, %151
%182 = or i64 %181, 1
%183 = getelementptr inbounds float, float* %0, i64 %182
%184 = bitcast float* %183 to <4 x float>*
%185 = load <4 x float>, <4 x float>* %184
%186 = getelementptr inbounds float, float* %0, i64 %181
%187 = bitcast float* %186 to <4 x float>*
%188 = load <4 x float>, <4 x float>* %187
%189 = fsub <4 x float> %185, %188
%190 = getelementptr inbounds float, float* %2, i64 %181
%191 = bitcast float* %190 to <4 x float>*
store <4 x float> %189, <4 x float>* %191
%192 = trunc i64 %181 to i32
%193 = add i32 %192, %3
%194 = add i32 %193, -1
%195 = sext i32 %194 to i64
%196 = getelementptr inbounds float, float* %0, i64 %195
%197 = bitcast float* %196 to <4 x float>*
%198 = load <4 x float>, <4 x float>* %197
%199 = add i32 %193, -2
%200 = sext i32 %199 to i64
%201 = getelementptr inbounds float, float* %0, i64 %200
%202 = bitcast float* %201 to <4 x float>*
%203 = load <4 x float>, <4 x float>* %202
%204 = fsub <4 x float> %198, %203
%205 = getelementptr inbounds float, float* %2, i64 %195
%206 = bitcast float* %205 to <4 x float>*
store <4 x float> %204, <4 x float>* %206
%207 = add nuw i64 %180, 4
%208 = icmp eq i64 %207, %178
br i1 %208, label %209, label %179
209: ; preds = %179
%210 = icmp eq i64 %178, %153
br i1 %210, label %248, label %211
211: ; preds = %157, %152, %209
%212 = phi i64 [ 0, %157 ], [ 0, %152 ], [ %178, %209 ]
br label %401
213: ; preds = %143, %213
%214 = phi float* [ %239, %213 ], [ %144, %143 ]
%215 = phi i32 [ %246, %213 ], [ %145, %143 ]
%216 = phi float* [ %244, %213 ], [ %146, %143 ]
%217 = phi float* [ %245, %213 ], [ %147, %143 ]
%218 = getelementptr inbounds float, float* %214, i64 1
%219 = load float, float* %217
%220 = load float, float* %216
%221 = fsub float %219, %220
%222 = fmul float %221, 5.000000e-01
store float %222, float* %218
%223 = getelementptr inbounds float, float* %216, i64 1
%224 = getelementptr inbounds float, float* %217, i64 1
%225 = getelementptr inbounds float, float* %214, i64 2
%226 = load float, float* %224
%227 = load float, float* %223
%228 = fsub float %226, %227
%229 = fmul float %228, 5.000000e-01
store float %229, float* %225
%230 = getelementptr inbounds float, float* %216, i64 2
%231 = getelementptr inbounds float, float* %217, i64 2
%232 = getelementptr inbounds float, float* %214, i64 3
%233 = load float, float* %231
%234 = load float, float* %230
%235 = fsub float %233, %234
%236 = fmul float %235, 5.000000e-01
store float %236, float* %232
%237 = getelementptr inbounds float, float* %216, i64 3
%238 = getelementptr inbounds float, float* %217, i64 3
%239 = getelementptr inbounds float, float* %214, i64 4
%240 = load float, float* %238
%241 = load float, float* %237
%242 = fsub float %240, %241
%243 = fmul float %242, 5.000000e-01
store float %243, float* %239
%244 = getelementptr inbounds float, float* %216, i64 4
%245 = getelementptr inbounds float, float* %217, i64 4
%246 = add nuw nsw i32 %215, 4
%247 = icmp eq i32 %246, %7
br i1 %247, label %149, label %213
248: ; preds = %401, %209, %149
%249 = sub i32 %6, %3
%250 = icmp sgt i32 %249, %3
br i1 %250, label %251, label %425
251: ; preds = %248
%252 = getelementptr inbounds float, float* %0, i64 %151
%253 = getelementptr float, float* %252, i64 %151
%254 = getelementptr float, float* %1, i64 %151
%255 = add i32 %4, -2
%256 = mul i32 %255, %3
%257 = add i32 %256, -1
%258 = zext i32 %257 to i64
%259 = add nuw nsw i64 %258, 1
%260 = icmp ult i32 %257, 7
br i1 %260, label %369, label %261
261: ; preds = %251
%262 = add i32 %4, -2
%263 = mul i32 %262, %3
%264 = add i32 %263, -1
%265 = zext i32 %264 to i64
%266 = add nsw i64 %151, %265
%267 = add nsw i64 %266, 1
%268 = getelementptr float, float* %1, i64 %267
%269 = add nuw nsw i64 %265, 1
%270 = getelementptr float, float* %0, i64 %269
%271 = shl nsw i64 %151, 1
%272 = add nsw i64 %271, %265
%273 = add nsw i64 %272, 1
%274 = getelementptr float, float* %0, i64 %273
%275 = icmp ult float* %254, %270
%276 = icmp ugt float* %268, %0
%277 = and i1 %275, %276
%278 = icmp ult float* %254, %274
%279 = icmp ult float* %253, %268
%280 = and i1 %278, %279
%281 = or i1 %277, %280
br i1 %281, label %369, label %282
282: ; preds = %261
%283 = and i64 %259, 8589934584
%284 = trunc i64 %283 to i32
%285 = add i32 %284, %3
%286 = getelementptr float, float* %254, i64 %283
%287 = getelementptr float, float* %0, i64 %283
%288 = getelementptr float, float* %253, i64 %283
%289 = add nsw i64 %283, -8
%290 = lshr exact i64 %289, 3
%291 = add nuw nsw i64 %290, 1
%292 = and i64 %291, 1
%293 = icmp eq i64 %289, 0
br i1 %293, label %343, label %294
294: ; preds = %282
%295 = and i64 %291, 4611686018427387902
br label %296
296: ; preds = %296, %294
%297 = phi i64 [ 0, %294 ], [ %340, %296 ]
%298 = phi i64 [ 0, %294 ], [ %341, %296 ]
%299 = getelementptr float, float* %254, i64 %297
%300 = getelementptr float, float* %0, i64 %297
%301 = getelementptr float, float* %253, i64 %297
%302 = bitcast float* %301 to <4 x float>*
%303 = load <4 x float>, <4 x float>* %302
%304 = getelementptr float, float* %301, i64 4
%305 = bitcast float* %304 to <4 x float>*
%306 = load <4 x float>, <4 x float>* %305
%307 = bitcast float* %300 to <4 x float>*
%308 = load <4 x float>, <4 x float>* %307
%309 = getelementptr float, float* %300, i64 4
%310 = bitcast float* %309 to <4 x float>*
%311 = load <4 x float>, <4 x float>* %310
%312 = fsub <4 x float> %303, %308
%313 = fsub <4 x float> %306, %311
%314 = fmul <4 x float> %312, <float 5.000000e-01, float 5.000000e-01, float 5.000000e-01, float 5.000000e-01>
%315 = fmul <4 x float> %313, <float 5.000000e-01, float 5.000000e-01, float 5.000000e-01, float 5.000000e-01>
%316 = bitcast float* %299 to <4 x float>*
store <4 x float> %314, <4 x float>* %316
%317 = getelementptr float, float* %299, i64 4
%318 = bitcast float* %317 to <4 x float>*
store <4 x float> %315, <4 x float>* %318
%319 = or i64 %297, 8
%320 = getelementptr float, float* %254, i64 %319
%321 = getelementptr float, float* %0, i64 %319
%322 = getelementptr float, float* %253, i64 %319
%323 = bitcast float* %322 to <4 x float>*
%324 = load <4 x float>, <4 x float>* %323
%325 = getelementptr float, float* %322, i64 4
%326 = bitcast float* %325 to <4 x float>*
%327 = load <4 x float>, <4 x float>* %326
%328 = bitcast float* %321 to <4 x float>*
%329 = load <4 x float>, <4 x float>* %328
%330 = getelementptr float, float* %321, i64 4
%331 = bitcast float* %330 to <4 x float>*
%332 = load <4 x float>, <4 x float>* %331
%333 = fsub <4 x float> %324, %329
%334 = fsub <4 x float> %327, %332
%335 = fmul <4 x float> %333, <float 5.000000e-01, float 5.000000e-01, float 5.000000e-01, float 5.000000e-01>
%336 = fmul <4 x float> %334, <float 5.000000e-01, float 5.000000e-01, float 5.000000e-01, float 5.000000e-01>
%337 = bitcast float* %320 to <4 x float>*
store <4 x float> %335, <4 x float>* %337
%338 = getelementptr float, float* %320, i64 4
%339 = bitcast float* %338 to <4 x float>*
store <4 x float> %336, <4 x float>* %339
%340 = add nuw i64 %297, 16
%341 = add i64 %298, 2
%342 = icmp eq i64 %341, %295
br i1 %342, label %343, label %296
343: ; preds = %296, %282
%344 = phi i64 [ 0, %282 ], [ %340, %296 ]
%345 = icmp eq i64 %292, 0
br i1 %345, label %367, label %346
346: ; preds = %343
%347 = getelementptr float, float* %254, i64 %344
%348 = getelementptr float, float* %0, i64 %344
%349 = getelementptr float, float* %253, i64 %344
%350 = bitcast float* %349 to <4 x float>*
%351 = load <4 x float>, <4 x float>* %350
%352 = getelementptr float, float* %349, i64 4
%353 = bitcast float* %352 to <4 x float>*
%354 = load <4 x float>, <4 x float>* %353
%355 = bitcast float* %348 to <4 x float>*
%356 = load <4 x float>, <4 x float>* %355
%357 = getelementptr float, float* %348, i64 4
%358 = bitcast float* %357 to <4 x float>*
%359 = load <4 x float>, <4 x float>* %358
%360 = fsub <4 x float> %351, %356
%361 = fsub <4 x float> %354, %359
%362 = fmul <4 x float> %360, <float 5.000000e-01, float 5.000000e-01, float 5.000000e-01, float 5.000000e-01>
%363 = fmul <4 x float> %361, <float 5.000000e-01, float 5.000000e-01, float 5.000000e-01, float 5.000000e-01>
%364 = bitcast float* %347 to <4 x float>*
store <4 x float> %362, <4 x float>* %364
%365 = getelementptr float, float* %347, i64 4
%366 = bitcast float* %365 to <4 x float>*
store <4 x float> %363, <4 x float>* %366
br label %367
367: ; preds = %343, %346
%368 = icmp eq i64 %259, %283
br i1 %368, label %425, label %369
369: ; preds = %261, %251, %367
%370 = phi i32 [ %3, %261 ], [ %3, %251 ], [ %285, %367 ]
%371 = phi float* [ %254, %261 ], [ %254, %251 ], [ %286, %367 ]
%372 = phi float* [ %0, %261 ], [ %0, %251 ], [ %287, %367 ]
%373 = phi float* [ %253, %261 ], [ %253, %251 ], [ %288, %367 ]
%374 = sub i32 %249, %370
%375 = xor i32 %370, -1
%376 = add i32 %249, %375
%377 = and i32 %374, 3
%378 = icmp eq i32 %377, 0
br i1 %378, label %395, label %379
379: ; preds = %369, %379
%380 = phi i32 [ %392, %379 ], [ %370, %369 ]
%381 = phi float* [ %389, %379 ], [ %371, %369 ]
%382 = phi float* [ %390, %379 ], [ %372, %369 ]
%383 = phi float* [ %391, %379 ], [ %373, %369 ]
%384 = phi i32 [ %393, %379 ], [ 0, %369 ]
%385 = load float, float* %383
%386 = load float, float* %382
%387 = fsub float %385, %386
%388 = fmul float %387, 5.000000e-01
store float %388, float* %381
%389 = getelementptr inbounds float, float* %381, i64 1
%390 = getelementptr inbounds float, float* %382, i64 1
%391 = getelementptr inbounds float, float* %383, i64 1
%392 = add nsw i32 %380, 1
%393 = add i32 %384, 1
%394 = icmp eq i32 %393, %377
br i1 %394, label %395, label %379
395: ; preds = %379, %369
%396 = phi i32 [ %370, %369 ], [ %392, %379 ]
%397 = phi float* [ %371, %369 ], [ %389, %379 ]
%398 = phi float* [ %372, %369 ], [ %390, %379 ]
%399 = phi float* [ %373, %369 ], [ %391, %379 ]
%400 = icmp ult i32 %376, 3
br i1 %400, label %425, label %435
401: ; preds = %211, %401
%402 = phi i64 [ %423, %401 ], [ %212, %211 ]
%403 = mul nsw i64 %402, %151
%404 = add nsw i64 %403, 1
%405 = getelementptr inbounds float, float* %0, i64 %404
%406 = load float, float* %405
%407 = getelementptr inbounds float, float* %0, i64 %403
%408 = load float, float* %407
%409 = fsub float %406, %408
%410 = getelementptr inbounds float, float* %2, i64 %403
store float %409, float* %410
%411 = trunc i64 %403 to i32
%412 = add i32 %411, %3
%413 = add i32 %412, -1
%414 = sext i32 %413 to i64
%415 = getelementptr inbounds float, float* %0, i64 %414
%416 = load float, float* %415
%417 = add i32 %412, -2
%418 = sext i32 %417 to i64
%419 = getelementptr inbounds float, float* %0, i64 %418
%420 = load float, float* %419
%421 = fsub float %416, %420
%422 = getelementptr inbounds float, float* %2, i64 %414
store float %421, float* %422
%423 = add nuw nsw i64 %402, 1
%424 = icmp eq i64 %423, %153
br i1 %424, label %248, label %401
425: ; preds = %395, %435, %367, %248
%426 = icmp sgt i32 %3, 0
br i1 %426, label %427, label %470
427: ; preds = %425
%428 = add nsw i32 %4, -1
%429 = mul nsw i32 %428, %3
%430 = add nsw i32 %4, -2
%431 = mul nsw i32 %430, %3
%432 = sext i32 %429 to i64
%433 = sext i32 %431 to i64
%434 = zext i32 %3 to i64
br label %471
435: ; preds = %395, %435
%436 = phi i32 [ %468, %435 ], [ %396, %395 ]
%437 = phi float* [ %465, %435 ], [ %397, %395 ]
%438 = phi float* [ %466, %435 ], [ %398, %395 ]
%439 = phi float* [ %467, %435 ], [ %399, %395 ]
%440 = load float, float* %439
%441 = load float, float* %438
%442 = fsub float %440, %441
%443 = fmul float %442, 5.000000e-01
store float %443, float* %437
%444 = getelementptr inbounds float, float* %437, i64 1
%445 = getelementptr inbounds float, float* %438, i64 1
%446 = getelementptr inbounds float, float* %439, i64 1
%447 = load float, float* %446
%448 = load float, float* %445
%449 = fsub float %447, %448
%450 = fmul float %449, 5.000000e-01
store float %450, float* %444
%451 = getelementptr inbounds float, float* %437, i64 2
%452 = getelementptr inbounds float, float* %438, i64 2
%453 = getelementptr inbounds float, float* %439, i64 2
%454 = load float, float* %453
%455 = load float, float* %452
%456 = fsub float %454, %455
%457 = fmul float %456, 5.000000e-01
store float %457, float* %451
%458 = getelementptr inbounds float, float* %437, i64 3
%459 = getelementptr inbounds float, float* %438, i64 3
%460 = getelementptr inbounds float, float* %439, i64 3
%461 = load float, float* %460
%462 = load float, float* %459
%463 = fsub float %461, %462
%464 = fmul float %463, 5.000000e-01
store float %464, float* %458
%465 = getelementptr inbounds float, float* %437, i64 4
%466 = getelementptr inbounds float, float* %438, i64 4
%467 = getelementptr inbounds float, float* %439, i64 4
%468 = add nsw i32 %436, 4
%469 = icmp eq i32 %468, %249
br i1 %469, label %425, label %435
470: ; preds = %471, %425
ret void
471: ; preds = %427, %471
%472 = phi i64 [ 0, %427 ], [ %488, %471 ]
%473 = add nsw i64 %472, %151
%474 = getelementptr inbounds float, float* %0, i64 %473
%475 = load float, float* %474
%476 = getelementptr inbounds float, float* %0, i64 %472
%477 = load float, float* %476
%478 = fsub float %475, %477
%479 = getelementptr inbounds float, float* %1, i64 %472
store float %478, float* %479
%480 = add nsw i64 %472, %432
%481 = getelementptr inbounds float, float* %0, i64 %480
%482 = load float, float* %481
%483 = add nsw i64 %472, %433
%484 = getelementptr inbounds float, float* %0, i64 %483
%485 = load float, float* %484
%486 = fsub float %482, %485
%487 = getelementptr inbounds float, float* %1, i64 %480
store float %486, float* %487
%488 = add nuw nsw i64 %472, 1
%489 = icmp eq i64 %488, %434
br i1 %489, label %470, label %471
}