void convxy(double out[], int xdim, int ydim, double filtx[], double filty[],
            int fxdim, int fydim, int xoff, int yoff, double buff[]) {
  int x, y, k;
  for (y = 0; y < ydim; y++) {
    for (x = 0; x < xdim; x++) {
      buff[x] = out[x + y * xdim];
      if (!((buff[x]) * (buff[x]) >= 0.))
        buff[x] = 0.0;
    }
    for (x = 0; x < xdim; x++) {
      double sum1 = 0.0;
      int fstart, fend;
      fstart = ((x - xoff >= xdim) ? x - xdim - xoff + 1 : 0);
      fend = ((x - (xoff + fxdim) < 0) ? x - xoff + 1 : fxdim);
      for (k = fstart; k < fend; k++)
        sum1 += buff[x - xoff - k] * filtx[k];
      out[x + y * xdim] = sum1;
    }
  }
  for (x = 0; x < xdim; x++) {
    for (y = 0; y < ydim; y++)
      buff[y] = out[x + y * xdim];
    for (y = 0; y < ydim; y++) {
      double sum1 = 0.0;
      int fstart, fend;
      fstart = ((y - yoff >= ydim) ? y - ydim - yoff + 1 : 0);
      fend = ((y - (yoff + fydim) < 0) ? y - yoff + 1 : fydim);
      for (k = fstart; k < fend; k++)
        sum1 += buff[y - yoff - k] * filty[k];
      out[y * xdim + x] = sum1;
    }
  }
}


x86 assembly (clang, -O0):
<convxy>:
push   %rbp
mov    %rsp,%rbp
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9d,-0x24(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    12cd <convxy+0x1cd>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11bc <convxy+0xbc>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x28(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x28(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x28(%rbp),%rax
movslq -0x28(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jae    11a9 <convxy+0xa9>
mov    0x28(%rbp),%rax
movslq -0x28(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    11ae <convxy+0xae>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1141 <convxy+0x41>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12ba <convxy+0x1ba>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x28(%rbp),%eax
sub    0x18(%rbp),%eax
cmp    -0xc(%rbp),%eax
jl     11fa <convxy+0xfa>
mov    -0x28(%rbp),%eax
sub    -0xc(%rbp),%eax
sub    0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    1204 <convxy+0x104>
xor    %eax,%eax
mov    %eax,-0x54(%rbp)
jmp    1204 <convxy+0x104>
mov    -0x54(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x28(%rbp),%eax
mov    0x18(%rbp),%ecx
add    -0x24(%rbp),%ecx
sub    %ecx,%eax
cmp    $0x0,%eax
jge    122f <convxy+0x12f>
mov    -0x28(%rbp),%eax
sub    0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    1235 <convxy+0x135>
mov    -0x24(%rbp),%eax
mov    %eax,-0x58(%rbp)
mov    -0x58(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x40(%rbp),%eax
jge    128f <convxy+0x18f>
mov    0x28(%rbp),%rax
mov    -0x28(%rbp),%ecx
sub    0x18(%rbp),%ecx
sub    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1241 <convxy+0x141>
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    11c3 <convxy+0xc3>
jmp    12bf <convxy+0x1bf>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    112e <convxy+0x2e>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1435 <convxy+0x335>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1326 <convxy+0x226>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    12e7 <convxy+0x1e7>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1422 <convxy+0x322>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x2c(%rbp),%eax
sub    0x20(%rbp),%eax
cmp    -0x10(%rbp),%eax
jl     1364 <convxy+0x264>
mov    -0x2c(%rbp),%eax
sub    -0x10(%rbp),%eax
sub    0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    136e <convxy+0x26e>
xor    %eax,%eax
mov    %eax,-0x5c(%rbp)
jmp    136e <convxy+0x26e>
mov    -0x5c(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x2c(%rbp),%eax
mov    0x20(%rbp),%ecx
add    0x10(%rbp),%ecx
sub    %ecx,%eax
cmp    $0x0,%eax
jge    1399 <convxy+0x299>
mov    -0x2c(%rbp),%eax
sub    0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x60(%rbp)
jmp    139f <convxy+0x29f>
mov    0x10(%rbp),%eax
mov    %eax,-0x60(%rbp)
mov    -0x60(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x4c(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x50(%rbp),%eax
jge    13f9 <convxy+0x2f9>
mov    0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    0x20(%rbp),%ecx
sub    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    13ab <convxy+0x2ab>
movsd  -0x48(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    132d <convxy+0x22d>
jmp    1427 <convxy+0x327>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    12d4 <convxy+0x1d4>
pop    %rbp
ret

x86 assembly (clang, -O2):
<convxy>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x20,%rsp
mov    %r9,-0x60(%rsp)
mov    0x70(%rsp),%rax
movabs $0x200000000,%r10
movabs $0x100000000,%r11
mov    %edx,%r9d
mov    %rsi,-0x80(%rsp)
mov    %esi,%r12d
mov    %rdx,-0x40(%rsp)
test   %edx,%edx
mov    %rdi,-0x18(%rsp)
mov    %r12,-0x20(%rsp)
jle    151d <convxy+0x41d>
mov    0x60(%rsp),%ebx
mov    -0x80(%rsp),%rbp
mov    $0x1,%edx
sub    %ebp,%edx
sub    %ebx,%edx
mov    %rdx,-0x28(%rsp)
mov    -0x60(%rsp),%rsi
add    %ebx,%esi
movslq %esi,%rdx
mov    %rdx,-0x30(%rsp)
movslq %ebx,%rdx
mov    %rdx,-0x38(%rsp)
movslq %ebp,%rdx
lea    (%rax,%r12,8),%rsi
mov    %rsi,0x10(%rsp)
mov    %r12d,%esi
and    $0xfffffffc,%esi
mov    %rsi,0x8(%rsp)
add    $0xfffffffffffffffc,%rsi
mov    %rsi,(%rsp)
mov    %rsi,%rbp
shr    $0x2,%rbp
add    $0x1,%rbp
mov    %rbp,%rsi
mov    %rbp,-0x8(%rsp)
and    $0xfffffffffffffffe,%rbp
mov    %rbp,-0x10(%rsp)
mov    %r12,%rsi
neg    %rsi
mov    %rsi,-0x50(%rsp)
lea    0x30(%rdi),%rsi
mov    %rdx,-0x58(%rsp)
lea    0x0(,%rdx,8),%rdx
mov    %rdx,-0x48(%rsp)
lea    0x8(%rdi),%rbx
xorpd  %xmm0,%xmm0
xor    %edx,%edx
mov    %rdx,-0x70(%rsp)
mov    %r9,0x18(%rsp)
jmp    1224 <convxy+0x124>
xchg   %ax,%ax
mov    -0x70(%rsp),%rbp
add    $0x1,%rbp
mov    -0x68(%rsp),%rsi
mov    -0x48(%rsp),%rdx
add    %rdx,%rsi
mov    -0x78(%rsp),%rbx
add    %rdx,%rbx
mov    0x18(%rsp),%r9
mov    %rbp,%rdx
mov    %rbp,-0x70(%rsp)
cmp    %r9,%rbp
je     151d <convxy+0x41d>
mov    %rbx,-0x78(%rsp)
mov    %rsi,-0x68(%rsp)
cmpl   $0x0,-0x80(%rsp)
jle    11f0 <convxy+0xf0>
mov    -0x70(%rsp),%r14
imul   -0x58(%rsp),%r14
cmpl   $0x4,-0x80(%rsp)
jb     1267 <convxy+0x167>
lea    (%r14,%r12,1),%rsi
lea    (%rdi,%rsi,8),%rsi
cmp    %rax,%rsi
jbe    1300 <convxy+0x200>
lea    (%rdi,%r14,8),%rsi
cmp    0x10(%rsp),%rsi
jae    1300 <convxy+0x200>
xor    %ebp,%ebp
mov    -0x78(%rsp),%r9
mov    %rbp,%rdx
test   $0x1,%r12b
je     12a1 <convxy+0x1a1>
lea    (%r14,%rbp,1),%rsi
movsd  (%rdi,%rsi,8),%xmm1
movapd %xmm1,%xmm2
mulsd  %xmm1,%xmm2
xorpd  %xmm3,%xmm3
cmpnlesd %xmm2,%xmm3
andnpd %xmm1,%xmm3
movlpd %xmm3,(%rax,%rbp,8)
mov    %rbp,%rdx
or     $0x1,%rdx
not    %rbp
cmp    -0x50(%rsp),%rbp
je     1410 <convxy+0x310>
nop
movsd  -0x8(%r9,%rdx,8),%xmm1
movapd %xmm1,%xmm2
mulsd  %xmm1,%xmm2
xorpd  %xmm3,%xmm3
cmpnlesd %xmm2,%xmm3
andnpd %xmm1,%xmm3
movlpd %xmm3,(%rax,%rdx,8)
movsd  (%r9,%rdx,8),%xmm1
movapd %xmm1,%xmm2
mulsd  %xmm1,%xmm2
xorpd  %xmm3,%xmm3
cmpnlesd %xmm2,%xmm3
andnpd %xmm1,%xmm3
movlpd %xmm3,0x8(%rax,%rdx,8)
add    $0x2,%rdx
cmp    %rdx,%r12
jne    12b0 <convxy+0x1b0>
jmp    1410 <convxy+0x310>
cmpq   $0x0,(%rsp)
mov    -0x68(%rsp),%rsi
mov    -0x78(%rsp),%r9
je     150b <convxy+0x40b>
mov    -0x10(%rsp),%rdx
xor    %ebp,%ebp
nopl   0x0(%rax)
movupd -0x30(%rsi,%rbp,8),%xmm1
movupd -0x20(%rsi,%rbp,8),%xmm2
movapd %xmm1,%xmm3
mulpd  %xmm1,%xmm3
movapd %xmm2,%xmm4
mulpd  %xmm2,%xmm4
xorpd  %xmm5,%xmm5
cmplepd %xmm3,%xmm5
andpd  %xmm1,%xmm5
xorpd  %xmm1,%xmm1
cmplepd %xmm4,%xmm1
andpd  %xmm2,%xmm1
movupd %xmm5,(%rax,%rbp,8)
movupd %xmm1,0x10(%rax,%rbp,8)
movupd -0x10(%rsi,%rbp,8),%xmm1
movupd (%rsi,%rbp,8),%xmm2
movapd %xmm1,%xmm3
mulpd  %xmm1,%xmm3
movapd %xmm2,%xmm4
mulpd  %xmm2,%xmm4
xorpd  %xmm5,%xmm5
cmplepd %xmm3,%xmm5
andpd  %xmm1,%xmm5
xorpd  %xmm1,%xmm1
cmplepd %xmm4,%xmm1
andpd  %xmm2,%xmm1
movupd %xmm5,0x20(%rax,%rbp,8)
movupd %xmm1,0x30(%rax,%rbp,8)
add    $0x8,%rbp
add    $0xfffffffffffffffe,%rdx
jne    1320 <convxy+0x220>
testb  $0x1,-0x8(%rsp)
je     13fb <convxy+0x2fb>
lea    (%r14,%rbp,1),%rsi
movupd (%rdi,%rsi,8),%xmm1
movupd 0x10(%rdi,%rsi,8),%xmm2
movapd %xmm1,%xmm3
mulpd  %xmm1,%xmm3
movapd %xmm2,%xmm4
mulpd  %xmm2,%xmm4
xorpd  %xmm5,%xmm5
cmplepd %xmm3,%xmm5
andpd  %xmm1,%xmm5
xorpd  %xmm1,%xmm1
cmplepd %xmm4,%xmm1
andpd  %xmm2,%xmm1
movupd %xmm5,(%rax,%rbp,8)
movupd %xmm1,0x10(%rax,%rbp,8)
mov    0x8(%rsp),%rdx
mov    %rdx,%rbp
cmp    %r12,%rdx
jne    126e <convxy+0x16e>
nopl   0x0(%rax)
cmpl   $0x0,-0x80(%rsp)
jle    11f0 <convxy+0xf0>
xor    %r9d,%r9d
jmp    1436 <convxy+0x336>
lea    (%r9,%r14,1),%rdx
movsd  %xmm1,(%rdi,%rdx,8)
add    $0x1,%r9
cmp    %r12,%r9
je     11f0 <convxy+0xf0>
mov    %r9,%r15
sub    -0x38(%rsp),%r15
mov    -0x28(%rsp),%rdx
lea    (%rdx,%r9,1),%esi
cmp    -0x58(%rsp),%r15
mov    $0x0,%edx
cmovl  %edx,%esi
lea    0x1(%r15),%edx
cmp    -0x30(%rsp),%r9
cmovge -0x60(%rsp),%edx
xorpd  %xmm1,%xmm1
cmp    %edx,%esi
jge    1420 <convxy+0x320>
movslq %esi,%rbp
mov    %edx,%r12d
sub    %esi,%r12d
xorpd  %xmm1,%xmm1
mov    %rbp,%r13
test   $0x1,%r12b
je     1498 <convxy+0x398>
mov    %r15,%rsi
sub    %rbp,%rsi
movsd  (%rax,%rsi,8),%xmm1
mulsd  (%rcx,%rbp,8),%xmm1
addsd  %xmm0,%xmm1
lea    0x1(%rbp),%r13
movslq %edx,%rbx
not    %rbp
add    %rbx,%rbp
mov    -0x18(%rsp),%rdi
mov    -0x20(%rsp),%r12
je     1420 <convxy+0x320>
mov    %r13,%rbp
shl    $0x20,%rbp
nopl   0x0(%rax,%rax,1)
mov    %rbp,%rsi
sar    $0x20,%rsi
mov    %r15,%rdx
sub    %rsi,%rdx
movsd  (%rax,%rdx,8),%xmm2
mulsd  (%rcx,%r13,8),%xmm2
addsd  %xmm1,%xmm2
lea    (%r11,%rbp,1),%rdx
sar    $0x20,%rdx
mov    %r15,%rsi
sub    %rdx,%rsi
movsd  (%rax,%rsi,8),%xmm1
mulsd  0x8(%rcx,%r13,8),%xmm1
addsd  %xmm2,%xmm1
add    $0x2,%r13
add    %r10,%rbp
cmp    %r13,%rbx
jne    14c0 <convxy+0x3c0>
jmp    1420 <convxy+0x320>
xor    %ebp,%ebp
testb  $0x1,-0x8(%rsp)
jne    13b7 <convxy+0x2b7>
jmp    13fb <convxy+0x2fb>
cmpl   $0x0,-0x80(%rsp)
jle    178b <convxy+0x68b>
mov    0x68(%rsp),%ecx
mov    0x58(%rsp),%ebx
mov    -0x40(%rsp),%rsi
mov    $0x1,%ebp
sub    %esi,%ebp
sub    %ecx,%ebp
mov    %rbp,-0x58(%rsp)
lea    (%rcx,%rbx,1),%edx
movslq %edx,%rdx
mov    %rdx,-0x60(%rsp)
movslq %ecx,%rcx
mov    %rcx,-0x28(%rsp)
movslq %esi,%rcx
mov    %rcx,-0x30(%rsp)
lea    -0x1(%r9),%rcx
mov    %rcx,-0x38(%rsp)
mov    %r9d,%ecx
and    $0x3,%ecx
mov    %rcx,-0x70(%rsp)
mov    %r9d,%ecx
and    $0xfffffffc,%ecx
mov    %rcx,-0x80(%rsp)
mov    %r12,%rcx
shl    $0x4,%rcx
mov    %rcx,-0x68(%rsp)
mov    %r12,%rcx
shl    $0x5,%rcx
mov    %rcx,-0x48(%rsp)
lea    0x0(,%r12,8),%rcx
mov    %rcx,-0x78(%rsp)
lea    (%rcx,%rcx,2),%rcx
mov    %rcx,-0x50(%rsp)
xorpd  %xmm0,%xmm0
mov    %rdi,%rbx
xor    %r13d,%r13d
jmp    15d1 <convxy+0x4d1>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r13
add    $0x8,%rbx
cmp    %r12,%r13
je     178b <convxy+0x68b>
cmpl   $0x0,-0x40(%rsp)
jle    15c0 <convxy+0x4c0>
cmpq   $0x3,-0x38(%rsp)
jae    15f0 <convxy+0x4f0>
mov    %rbx,%r15
xor    %esi,%esi
jmp    1647 <convxy+0x547>
nopw   0x0(%rax,%rax,1)
mov    %rbx,%r15
mov    %rbx,%rcx
xor    %esi,%esi
mov    -0x80(%rsp),%r14
mov    -0x68(%rsp),%rbx
mov    -0x48(%rsp),%rdx
mov    -0x50(%rsp),%rbp
nopl   0x0(%rax)
movsd  (%rcx),%xmm1
movsd  %xmm1,(%rax,%rsi,8)
movsd  (%rcx,%r12,8),%xmm1
movsd  %xmm1,0x8(%rax,%rsi,8)
movsd  (%rcx,%rbx,1),%xmm1
movsd  %xmm1,0x10(%rax,%rsi,8)
movsd  (%rcx,%rbp,1),%xmm1
movsd  %xmm1,0x18(%rax,%rsi,8)
add    $0x4,%rsi
add    %rdx,%rcx
cmp    %rsi,%r14
jne    1610 <convxy+0x510>
mov    -0x70(%rsp),%rdx
test   %rdx,%rdx
mov    -0x78(%rsp),%rbx
je     1685 <convxy+0x585>
lea    (%rax,%rsi,8),%rcx
imul   %r12,%rsi
add    %r13,%rsi
lea    (%rdi,%rsi,8),%rsi
xor    %ebp,%ebp
nopw   0x0(%rax,%rax,1)
movsd  (%rsi),%xmm1
movsd  %xmm1,(%rcx,%rbp,8)
add    $0x1,%rbp
add    %rbx,%rsi
cmp    %rbp,%rdx
jne    1670 <convxy+0x570>
cmpl   $0x0,-0x40(%rsp)
mov    $0x0,%edx
mov    %r15,%rbx
jle    15c0 <convxy+0x4c0>
xor    %r14d,%r14d
jmp    16c6 <convxy+0x5c6>
nopl   (%rax)
mov    %r14,%rcx
mov    -0x20(%rsp),%r12
imul   %r12,%rcx
add    %r13,%rcx
mov    -0x18(%rsp),%rdi
movsd  %xmm1,(%rdi,%rcx,8)
add    $0x1,%r14
cmp    %r9,%r14
je     15c0 <convxy+0x4c0>
mov    %r14,%r12
sub    -0x28(%rsp),%r12
mov    -0x58(%rsp),%rcx
lea    (%rcx,%r14,1),%esi
cmp    -0x30(%rsp),%r12
cmovl  %edx,%esi
lea    0x1(%r12),%edi
cmp    -0x60(%rsp),%r14
mov    0x58(%rsp),%ecx
cmovge %ecx,%edi
xorpd  %xmm1,%xmm1
cmp    %edi,%esi
jge    16a0 <convxy+0x5a0>
movslq %esi,%rbp
mov    %edi,%ecx
sub    %esi,%ecx
xorpd  %xmm1,%xmm1
mov    %rbp,%rsi
test   $0x1,%cl
je     1724 <convxy+0x624>
mov    %r12,%rcx
sub    %rbp,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  (%r8,%rbp,8),%xmm1
addsd  %xmm0,%xmm1
lea    0x1(%rbp),%rsi
movslq %edi,%rdi
not    %rbp
add    %rdi,%rbp
je     16a0 <convxy+0x5a0>
mov    %rsi,%r15
shl    $0x20,%r15
nopw   0x0(%rax,%rax,1)
mov    %r15,%rcx
sar    $0x20,%rcx
mov    %r12,%rbp
sub    %rcx,%rbp
movsd  (%rax,%rbp,8),%xmm2
mulsd  (%r8,%rsi,8),%xmm2
addsd  %xmm1,%xmm2
lea    (%r15,%r11,1),%rcx
sar    $0x20,%rcx
mov    %r12,%rbp
sub    %rcx,%rbp
movsd  (%rax,%rbp,8),%xmm1
mulsd  0x8(%r8,%rsi,8),%xmm1
addsd  %xmm2,%xmm1
add    $0x2,%rsi
add    %r10,%r15
cmp    %rsi,%rdi
jne    1740 <convxy+0x640>
jmp    16a0 <convxy+0x5a0>
add    $0x20,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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