void deconvolve3_columns(int width, int height, int rowstride, double *data,
                         double *buffer, double a, double b) {
  double *row;
  double q;
  int i, j;
  if (!height || !width)
    return;
  if (height == 1) {
    q = a + 2.0 * b;
    for (j = 0; j < width; j++)
      data[j] /= q;
    return;
  }
  if (height == 2) {
    q = a * (a + 2.0 * b);
    for (j = 0; j < width; j++) {
      buffer[0] = (a + b) / q * data[j] - b / q * data[rowstride + j];
      data[rowstride + j] = (a + b) / q * data[rowstride + j] - b / q * data[j];
      data[j] = buffer[0];
    }
    return;
  }
  buffer[0] = a + b;
  for (i = 1; i < height - 1; i++) {
    q = b / buffer[i - 1];
    buffer[i] = a - q * b;
    row = data + (i - 1) * rowstride;
    for (j = 0; j < width; j++)
      row[rowstride + j] -= q * row[j];
  }
  q = b / buffer[i - 1];
  buffer[i] = a + b * (1.0 - q);
  row = data + (i - 1) * rowstride;
  for (j = 0; j < width; j++)
    row[rowstride + j] -= q * row[j];
  row += rowstride;
  for (j = 0; j < width; j++)
    row[j] /= buffer[i];
  do {
    i--;
    row = data + i * rowstride;
    for (j = 0; j < width; j++)
      row[j] = (row[j] - b * row[rowstride + j]) / buffer[i];
  } while (i > 0);
}


x86 assembly (clang, -O0):
<deconvolve3_columns>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
movsd  %xmm0,-0x28(%rbp)
movsd  %xmm1,-0x30(%rbp)
cmpl   $0x0,-0x8(%rbp)
je     1133 <deconvolve3_columns+0x33>
cmpl   $0x0,-0x4(%rbp)
jne    1138 <deconvolve3_columns+0x38>
jmp    15ba <deconvolve3_columns+0x4ba>
cmpl   $0x1,-0x8(%rbp)
jne    119a <deconvolve3_columns+0x9a>
movsd  -0x28(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm0
addsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1195 <deconvolve3_columns+0x95>
movsd  -0x40(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1160 <deconvolve3_columns+0x60>
jmp    15ba <deconvolve3_columns+0x4ba>
cmpl   $0x2,-0x8(%rbp)
jne    12c7 <deconvolve3_columns+0x1c7>
movsd  -0x28(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm2
movsd  -0x30(%rbp),%xmm1
addsd  %xmm1,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    12c2 <deconvolve3_columns+0x1c2>
movsd  -0x28(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
divsd  -0x40(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x30(%rbp),%xmm1
divsd  -0x40(%rbp),%xmm1
mov    -0x18(%rbp),%rax
mov    -0xc(%rbp),%ecx
add    -0x48(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x28(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
divsd  -0x40(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0xc(%rbp),%ecx
add    -0x48(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x30(%rbp),%xmm1
divsd  -0x40(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x48(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0xc(%rbp),%ecx
add    -0x48(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    11cb <deconvolve3_columns+0xcb>
jmp    15ba <deconvolve3_columns+0x4ba>
movsd  -0x28(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    -0x8(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    13e0 <deconvolve3_columns+0x2e0>
movsd  -0x30(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x28(%rbp),%xmm1
movsd  -0x40(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
imul   -0xc(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    13cd <deconvolve3_columns+0x2cd>
movsd  -0x40(%rbp),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x38(%rbp),%rax
mov    -0xc(%rbp),%ecx
add    -0x48(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    136b <deconvolve3_columns+0x26b>
jmp    13d2 <deconvolve3_columns+0x2d2>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    12e0 <deconvolve3_columns+0x1e0>
movsd  -0x30(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x28(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm0
movsd  0xbf2(%rip),%xmm2
subsd  -0x40(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
imul   -0xc(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    14ad <deconvolve3_columns+0x3ad>
movsd  -0x40(%rbp),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x38(%rbp),%rax
mov    -0xc(%rbp),%ecx
add    -0x48(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    144b <deconvolve3_columns+0x34b>
mov    -0xc(%rbp),%ecx
mov    -0x38(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1506 <deconvolve3_columns+0x406>
mov    -0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x38(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    14c9 <deconvolve3_columns+0x3c9>
jmp    150b <deconvolve3_columns+0x40b>
mov    -0x44(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x44(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x44(%rbp),%ecx
imul   -0xc(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    15ab <deconvolve3_columns+0x4ab>
mov    -0x38(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x30(%rbp),%xmm0
mov    -0x38(%rbp),%rax
mov    -0xc(%rbp),%ecx
add    -0x48(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1534 <deconvolve3_columns+0x434>
jmp    15b0 <deconvolve3_columns+0x4b0>
cmpl   $0x0,-0x44(%rbp)
jg     150b <deconvolve3_columns+0x40b>
pop    %rbp
ret

x86 assembly (clang, -O2):
<deconvolve3_columns>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %edi,%edi
je     15e9 <deconvolve3_columns+0x4e9>
test   %esi,%esi
je     15e9 <deconvolve3_columns+0x4e9>
mov    %rcx,%r10
cmp    $0x2,%esi
je     114d <deconvolve3_columns+0x4d>
cmp    $0x1,%esi
jne    11e1 <deconvolve3_columns+0xe1>
test   %edi,%edi
jle    15e9 <deconvolve3_columns+0x4e9>
addsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
mov    %edi,%eax
cmp    $0x1,%edi
jne    154f <deconvolve3_columns+0x44f>
xor    %ecx,%ecx
jmp    15d0 <deconvolve3_columns+0x4d0>
test   %edi,%edi
jle    15e9 <deconvolve3_columns+0x4e9>
movapd %xmm1,%xmm2
addsd  %xmm1,%xmm2
addsd  %xmm0,%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm1,%xmm0
divsd  %xmm2,%xmm0
xorpd  0xe8b(%rip),%xmm1
divsd  %xmm2,%xmm1
movslq %edx,%rcx
mov    %edi,%eax
lea    (%r10,%rcx,8),%rcx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rcx,%rdx,8),%xmm2
mulsd  %xmm1,%xmm2
movsd  (%r10,%rdx,8),%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm2,%xmm3
movsd  %xmm3,(%r8)
movsd  (%r10,%rdx,8),%xmm2
mulsd  %xmm1,%xmm2
movsd  (%rcx,%rdx,8),%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm2,%xmm3
movsd  %xmm3,(%rcx,%rdx,8)
movsd  (%r8),%xmm2
movsd  %xmm2,(%r10,%rdx,8)
add    $0x1,%rdx
cmp    %rdx,%rax
jne    1190 <deconvolve3_columns+0x90>
jmp    15e9 <deconvolve3_columns+0x4e9>
movapd %xmm0,%xmm4
addsd  %xmm1,%xmm4
movsd  %xmm4,(%r8)
mov    $0x1,%eax
mov    %edi,%r11d
cmp    $0x3,%esi
mov    %r10,-0x70(%rsp)
mov    %edi,-0x74(%rsp)
jl     14f1 <deconvolve3_columns+0x3f1>
add    $0xffffffff,%esi
movapd 0xdec(%rip),%xmm8
xorpd  %xmm1,%xmm8
mov    %edx,-0x48(%rsp)
movslq %edx,%r9
mov    %esi,-0x50(%rsp)
mov    %esi,%r13d
lea    (%r9,%r11,1),%rax
mov    %rax,-0x60(%rsp)
mov    %r11d,%eax
and    $0xfffffffc,%eax
mov    %rax,-0x28(%rsp)
add    $0xfffffffffffffffc,%rax
mov    %rax,-0x30(%rsp)
mov    %rax,%rcx
shr    $0x2,%rcx
add    $0x1,%rcx
mov    %rcx,%rax
mov    %rcx,-0x38(%rsp)
and    $0xfffffffffffffffe,%rcx
mov    %rcx,-0x40(%rsp)
mov    %r11,%rax
neg    %rax
mov    %rax,-0x58(%rsp)
lea    0x30(%r10),%rbp
lea    0x0(,%r9,8),%r15
lea    (%r10,%r9,8),%r14
add    $0x8,%r14
mov    $0x1,%edx
mov    %r10,%rcx
xor    %r10d,%r10d
mov    %r9,-0x20(%rsp)
mov    %r13,-0x68(%rsp)
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movapd %xmm8,%xmm9
divsd  %xmm4,%xmm9
movapd %xmm9,%xmm4
mulsd  %xmm1,%xmm4
addsd  %xmm0,%xmm4
movsd  %xmm4,(%r8,%rdx,8)
test   %edi,%edi
jle    14c0 <deconvolve3_columns+0x3c0>
lea    -0x1(%rdx),%rax
imul   %r9,%rax
mov    %edi,%esi
mov    -0x70(%rsp),%rdi
lea    (%rdi,%rax,8),%r12
cmp    $0x4,%esi
jb     1316 <deconvolve3_columns+0x216>
mov    %r10,%rsi
imul   %r9,%rsi
lea    (%rsi,%r9,1),%rbx
lea    (%rsi,%r11,1),%rax
lea    (%rdi,%rbx,8),%rbx
lea    (%rdi,%rax,8),%rax
cmp    %rax,%rbx
jae    139d <deconvolve3_columns+0x29d>
mov    -0x60(%rsp),%rax
add    %rsi,%rax
lea    (%rdi,%rax,8),%rax
lea    (%rdi,%rsi,8),%rsi
cmp    %rax,%rsi
jae    139d <deconvolve3_columns+0x29d>
xor    %esi,%esi
mov    -0x74(%rsp),%edi
mov    %rsi,%rbx
test   $0x1,%r11b
je     1347 <deconvolve3_columns+0x247>
lea    (%rsi,%r9,1),%rax
movsd  (%r12,%rsi,8),%xmm2
mulsd  %xmm9,%xmm2
addsd  (%r12,%rax,8),%xmm2
movsd  %xmm2,(%r12,%rax,8)
mov    %rsi,%rbx
or     $0x1,%rbx
not    %rsi
cmp    -0x58(%rsp),%rsi
je     14c0 <deconvolve3_columns+0x3c0>
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rcx,%rbx,8),%xmm2
mulsd  %xmm9,%xmm2
addsd  -0x8(%r14,%rbx,8),%xmm2
movsd  %xmm2,-0x8(%r14,%rbx,8)
movsd  0x8(%rcx,%rbx,8),%xmm2
mulsd  %xmm9,%xmm2
addsd  (%r14,%rbx,8),%xmm2
movsd  %xmm2,(%r14,%rbx,8)
add    $0x2,%rbx
cmp    %rbx,%r11
jne    1360 <deconvolve3_columns+0x260>
jmp    14c0 <deconvolve3_columns+0x3c0>
movapd %xmm9,%xmm4
unpcklpd %xmm9,%xmm4
cmpq   $0x0,-0x30(%rsp)
je     144e <deconvolve3_columns+0x34e>
lea    (%r15,%rbp,1),%rbx
mov    -0x40(%rsp),%r9
xor    %r13d,%r13d
mov    -0x74(%rsp),%edi
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movupd -0x30(%rbp,%r13,8),%xmm5
movupd -0x20(%rbp,%r13,8),%xmm6
mulpd  %xmm4,%xmm5
movupd -0x30(%rbx,%r13,8),%xmm7
addpd  %xmm5,%xmm7
movupd -0x20(%rbx,%r13,8),%xmm5
movupd -0x10(%rbx,%r13,8),%xmm2
movupd (%rbx,%r13,8),%xmm3
mulpd  %xmm4,%xmm6
addpd  %xmm5,%xmm6
movupd %xmm7,-0x30(%rbx,%r13,8)
movupd %xmm6,-0x20(%rbx,%r13,8)
movupd -0x10(%rbp,%r13,8),%xmm5
movupd 0x0(%rbp,%r13,8),%xmm6
mulpd  %xmm4,%xmm5
addpd  %xmm2,%xmm5
mulpd  %xmm4,%xmm6
addpd  %xmm3,%xmm6
movupd %xmm5,-0x10(%rbx,%r13,8)
movupd %xmm6,(%rbx,%r13,8)
add    $0x8,%r13
add    $0xfffffffffffffffe,%r9
jne    13d0 <deconvolve3_columns+0x2d0>
jmp    1455 <deconvolve3_columns+0x355>
xor    %r13d,%r13d
mov    -0x74(%rsp),%edi
testb  $0x1,-0x38(%rsp)
mov    -0x20(%rsp),%r9
je     149b <deconvolve3_columns+0x39b>
movupd (%r12,%r13,8),%xmm2
movupd 0x10(%r12,%r13,8),%xmm3
add    %r9,%r13
movupd (%r12,%r13,8),%xmm5
movupd 0x10(%r12,%r13,8),%xmm6
mulpd  %xmm4,%xmm2
addpd  %xmm5,%xmm2
mulpd  %xmm4,%xmm3
addpd  %xmm6,%xmm3
movupd %xmm2,(%r12,%r13,8)
movupd %xmm3,0x10(%r12,%r13,8)
mov    -0x28(%rsp),%rax
mov    %rax,%rsi
cmp    %r11,%rax
mov    -0x68(%rsp),%r13
jne    131c <deconvolve3_columns+0x21c>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
lea    0x1(%rdx),%rax
cmp    %r13,%rax
je     14e4 <deconvolve3_columns+0x3e4>
movsd  (%r8,%rdx,8),%xmm4
add    $0x1,%r10
add    %r15,%rbp
add    %r15,%r14
add    %r15,%rcx
mov    %rax,%rdx
jmp    12a0 <deconvolve3_columns+0x1a0>
mov    -0x50(%rsp),%eax
mov    -0x70(%rsp),%r10
mov    -0x48(%rsp),%edx
movslq %eax,%rcx
movapd %xmm1,%xmm2
divsd  -0x8(%r8,%rcx,8),%xmm2
movsd  0xb09(%rip),%xmm3
subsd  %xmm2,%xmm3
mulsd  %xmm1,%xmm3
addsd  %xmm0,%xmm3
mov    %eax,%esi
movsd  %xmm3,(%r8,%rsi,8)
test   %edi,%edi
jle    1547 <deconvolve3_columns+0x447>
add    $0xffffffff,%eax
imul   %edx,%eax
movslq %eax,%r9
lea    (%r10,%r9,8),%r12
xorpd  0xacc(%rip),%xmm2
movslq %edx,%rbp
cmp    $0x4,%edi
jae    15f4 <deconvolve3_columns+0x4f4>
xor    %eax,%eax
jmp    174b <deconvolve3_columns+0x64b>
movslq %edx,%rbp
jmp    190b <deconvolve3_columns+0x80b>
mov    %eax,%ecx
and    $0xfffffffe,%ecx
movapd %xmm1,%xmm0
unpcklpd %xmm1,%xmm0
lea    -0x2(%rcx),%rsi
mov    %rsi,%rdx
shr    %rdx
add    $0x1,%rdx
test   %rsi,%rsi
je     16e5 <deconvolve3_columns+0x5e5>
mov    %rdx,%rdi
and    $0xfffffffffffffffe,%rdi
xor    %esi,%esi
nopl   0x0(%rax)
movupd (%r10,%rsi,8),%xmm2
movupd 0x10(%r10,%rsi,8),%xmm3
divpd  %xmm0,%xmm2
movupd %xmm2,(%r10,%rsi,8)
divpd  %xmm0,%xmm3
movupd %xmm3,0x10(%r10,%rsi,8)
add    $0x4,%rsi
add    $0xfffffffffffffffe,%rdi
jne    1580 <deconvolve3_columns+0x480>
test   $0x1,%dl
je     15c1 <deconvolve3_columns+0x4c1>
movupd (%r10,%rsi,8),%xmm2
divpd  %xmm0,%xmm2
movupd %xmm2,(%r10,%rsi,8)
cmp    %rax,%rcx
je     15e9 <deconvolve3_columns+0x4e9>
cs nopw 0x0(%rax,%rax,1)
movsd  (%r10,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%r10,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    15d0 <deconvolve3_columns+0x4d0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
lea    (%r9,%rbp,1),%rax
lea    (%r10,%rax,8),%rcx
lea    (%r9,%r11,1),%rdi
lea    (%r10,%rdi,8),%rdi
cmp    %rdi,%rcx
jae    1620 <deconvolve3_columns+0x520>
add    %r11,%rax
lea    (%r10,%rax,8),%rax
cmp    %rax,%r12
jae    1620 <deconvolve3_columns+0x520>
xor    %eax,%eax
mov    -0x74(%rsp),%edi
jmp    174b <deconvolve3_columns+0x64b>
mov    %r11d,%eax
and    $0xfffffffc,%eax
movapd %xmm2,%xmm0
unpcklpd %xmm2,%xmm0
lea    -0x4(%rax),%rcx
mov    %rcx,%r10
shr    $0x2,%r10
add    $0x1,%r10
test   %rcx,%rcx
je     16f5 <deconvolve3_columns+0x5f5>
mov    %r10,%rcx
and    $0xfffffffffffffffe,%rcx
mov    -0x70(%rsp),%rdi
lea    (%rdi,%r9,8),%r14
add    $0x30,%r14
mov    %rbp,%r15
lea    (%r12,%rbp,8),%rbx
add    $0x30,%rbx
xor    %ebp,%ebp
mov    -0x74(%rsp),%edi
nopl   0x0(%rax,%rax,1)
movupd -0x30(%r14,%rbp,8),%xmm3
movupd -0x20(%r14,%rbp,8),%xmm4
mulpd  %xmm0,%xmm3
movupd -0x30(%rbx,%rbp,8),%xmm5
addpd  %xmm3,%xmm5
movupd -0x20(%rbx,%rbp,8),%xmm3
movupd -0x10(%rbx,%rbp,8),%xmm6
movupd (%rbx,%rbp,8),%xmm7
mulpd  %xmm0,%xmm4
addpd  %xmm3,%xmm4
movupd %xmm5,-0x30(%rbx,%rbp,8)
movupd %xmm4,-0x20(%rbx,%rbp,8)
movupd -0x10(%r14,%rbp,8),%xmm3
movupd (%r14,%rbp,8),%xmm4
mulpd  %xmm0,%xmm3
addpd  %xmm6,%xmm3
mulpd  %xmm0,%xmm4
addpd  %xmm7,%xmm4
movupd %xmm3,-0x10(%rbx,%rbp,8)
movupd %xmm4,(%rbx,%rbp,8)
add    $0x8,%rbp
add    $0xfffffffffffffffe,%rcx
jne    1670 <deconvolve3_columns+0x570>
jmp    16fe <deconvolve3_columns+0x5fe>
xor    %esi,%esi
test   $0x1,%dl
jne    15b1 <deconvolve3_columns+0x4b1>
jmp    15c1 <deconvolve3_columns+0x4c1>
mov    %rbp,%r15
xor    %ebp,%ebp
mov    -0x74(%rsp),%edi
test   $0x1,%r10b
je     173e <deconvolve3_columns+0x63e>
movupd (%r12,%rbp,8),%xmm3
movupd 0x10(%r12,%rbp,8),%xmm4
add    %r15,%rbp
movupd (%r12,%rbp,8),%xmm5
movupd 0x10(%r12,%rbp,8),%xmm6
mulpd  %xmm0,%xmm3
addpd  %xmm5,%xmm3
mulpd  %xmm0,%xmm4
addpd  %xmm6,%xmm4
movupd %xmm3,(%r12,%rbp,8)
movupd %xmm4,0x10(%r12,%rbp,8)
cmp    %r11,%rax
mov    -0x70(%rsp),%r10
mov    %r15,%rbp
je     17c4 <deconvolve3_columns+0x6c4>
mov    %rax,%rcx
not    %rcx
test   $0x1,%r11b
je     1779 <deconvolve3_columns+0x679>
lea    (%rax,%rbp,1),%rdi
movsd  (%r12,%rax,8),%xmm0
mulsd  %xmm2,%xmm0
addsd  (%r12,%rdi,8),%xmm0
movsd  %xmm0,(%r12,%rdi,8)
mov    -0x74(%rsp),%edi
or     $0x1,%rax
add    %r11,%rcx
je     17c4 <deconvolve3_columns+0x6c4>
lea    (%r9,%rbp,1),%rcx
lea    (%r10,%rcx,8),%rcx
cs nopw 0x0(%rax,%rax,1)
movsd  (%r12,%rax,8),%xmm0
mulsd  %xmm2,%xmm0
addsd  (%rcx,%rax,8),%xmm0
movsd  %xmm0,(%rcx,%rax,8)
movsd  0x8(%r12,%rax,8),%xmm0
mulsd  %xmm2,%xmm0
addsd  0x8(%rcx,%rax,8),%xmm0
movsd  %xmm0,0x8(%rcx,%rax,8)
add    $0x2,%rax
cmp    %rax,%r11
jne    1790 <deconvolve3_columns+0x690>
test   %edi,%edi
jle    190b <deconvolve3_columns+0x80b>
lea    (%r8,%rsi,8),%r14
lea    (%r12,%rbp,8),%r10
cmp    $0x1,%edi
je     1800 <deconvolve3_columns+0x700>
lea    (%r9,%rbp,1),%rbx
lea    (%r8,%rsi,8),%rax
add    $0x8,%rax
cmp    %rax,%r10
jae    186e <deconvolve3_columns+0x76e>
lea    (%rbx,%r11,1),%rax
mov    -0x70(%rsp),%rcx
lea    (%rcx,%rax,8),%rax
cmp    %rax,%r14
jae    186e <deconvolve3_columns+0x76e>
xor    %ecx,%ecx
mov    %rcx,%rax
not    %rax
test   $0x1,%r11b
je     1823 <deconvolve3_columns+0x723>
movsd  (%r10,%rcx,8),%xmm0
divsd  (%r14),%xmm0
movsd  %xmm0,(%r10,%rcx,8)
or     $0x1,%rcx
mov    -0x70(%rsp),%r10
add    %r11,%rax
je     190b <deconvolve3_columns+0x80b>
add    %rbp,%r9
lea    (%r10,%r9,8),%rax
add    $0x8,%rax
nopl   0x0(%rax)
movsd  -0x8(%rax,%rcx,8),%xmm0
divsd  (%r14),%xmm0
movsd  (%rax,%rcx,8),%xmm2
movsd  %xmm0,-0x8(%rax,%rcx,8)
divsd  (%r14),%xmm2
movsd  %xmm2,(%rax,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%r11
jne    1840 <deconvolve3_columns+0x740>
jmp    190b <deconvolve3_columns+0x80b>
mov    %rbp,%r15
mov    %r11d,%ecx
and    $0xfffffffe,%ecx
lea    -0x2(%rcx),%rax
mov    %rax,%rdx
shr    %rdx
add    $0x1,%rdx
test   %rax,%rax
je     18da <deconvolve3_columns+0x7da>
mov    %rdx,%rbp
and    $0xfffffffffffffffe,%rbp
movsd  (%r14),%xmm0
unpcklpd %xmm0,%xmm0
mov    -0x70(%rsp),%rax
lea    (%rax,%rbx,8),%rax
add    $0x10,%rax
xor    %ebx,%ebx
nopl   0x0(%rax)
movupd -0x10(%rax,%rbx,8),%xmm2
movupd (%rax,%rbx,8),%xmm3
divpd  %xmm0,%xmm2
movupd %xmm2,-0x10(%rax,%rbx,8)
divpd  %xmm0,%xmm3
movupd %xmm3,(%rax,%rbx,8)
add    $0x4,%rbx
add    $0xfffffffffffffffe,%rbp
jne    18b0 <deconvolve3_columns+0x7b0>
jmp    18dc <deconvolve3_columns+0x7dc>
xor    %ebx,%ebx
test   $0x1,%dl
je     18fa <deconvolve3_columns+0x7fa>
movsd  (%r14),%xmm0
unpcklpd %xmm0,%xmm0
movupd (%r10,%rbx,8),%xmm2
divpd  %xmm0,%xmm2
movupd %xmm2,(%r10,%rbx,8)
cmp    %r11,%rcx
mov    %r15,%rbp
jne    1802 <deconvolve3_columns+0x702>
mov    -0x70(%rsp),%r10
xorpd  0x6ed(%rip),%xmm1
lea    -0x1(%rsi),%rax
mov    %rbp,%rcx
mov    %rax,-0x28(%rsp)
imul   %rax,%rcx
lea    (%rcx,%r11,1),%rax
mov    %rax,-0x38(%rsp)
mov    %rbp,%rax
imul   %rsi,%rax
lea    (%rax,%r11,1),%rdx
mov    %rdx,-0x48(%rsp)
mov    %r11d,%edx
and    $0xfffffffe,%edx
mov    %rdx,-0x50(%rsp)
add    $0xfffffffffffffffe,%rdx
mov    %rdx,-0x8(%rsp)
mov    %rdx,%rbx
shr    %rbx
add    $0x1,%rbx
pshufd $0x44,%xmm1,%xmm0
mov    %rbx,%rdx
mov    %rbx,-0x10(%rsp)
and    $0xfffffffffffffffe,%rbx
mov    %rbx,-0x18(%rsp)
mov    %r11,%rdx
neg    %rdx
mov    %rdx,-0x20(%rsp)
lea    (%r10,%rcx,8),%r12
add    $0x10,%r12
lea    0x0(,%rbp,8),%rdx
neg    %rdx
lea    (%r10,%rax,8),%r14
add    $0x10,%r14
mov    %rax,-0x40(%rsp)
lea    (%r10,%rax,8),%r13
add    $0x8,%r13
mov    %rcx,-0x30(%rsp)
lea    (%r10,%rcx,8),%r9
xor    %r15d,%r15d
mov    %rbp,-0x68(%rsp)
jmp    19dc <deconvolve3_columns+0x8dc>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%r15
add    %rdx,%r12
add    %rdx,%r14
add    %rdx,%r13
add    %rdx,%r9
cmpq   $0x1,-0x58(%rsp)
jle    15e9 <deconvolve3_columns+0x4e9>
mov    %rsi,%rax
mov    %rsi,-0x58(%rsp)
add    $0xffffffffffffffff,%rsi
test   %edi,%edi
jle    19c0 <deconvolve3_columns+0x8c0>
mov    %rsi,%rax
imul   %rbp,%rax
lea    (%r10,%rax,8),%rax
mov    %rax,-0x60(%rsp)
cmp    $0x1,%edi
jne    1a10 <deconvolve3_columns+0x910>
xor    %ebx,%ebx
jmp    1b70 <deconvolve3_columns+0xa70>
nopl   0x0(%rax,%rax,1)
mov    %rbp,%rcx
imul   %r15,%rcx
mov    -0x30(%rsp),%rdi
sub    %rcx,%rdi
mov    -0x38(%rsp),%rbx
sub    %rcx,%rbx
mov    -0x40(%rsp),%rbp
sub    %rcx,%rbp
mov    -0x48(%rsp),%rax
sub    %rcx,%rax
mov    -0x28(%rsp),%rcx
sub    %r15,%rcx
lea    (%r10,%rdi,8),%rdi
lea    (%r10,%rbx,8),%rbx
lea    (%r10,%rbp,8),%rbp
lea    (%r10,%rax,8),%rax
cmp    %rax,%rdi
setb   %al
cmp    %rbx,%rbp
mov    -0x58(%rsp),%rbp
lea    (%r8,%rbp,8),%rbp
setb   %r10b
cmp    %rbp,%rdi
lea    (%r8,%rcx,8),%rdi
setb   %cl
cmp    %rbx,%rdi
setb   %bpl
test   %r10b,%al
jne    1b51 <deconvolve3_columns+0xa51>
mov    $0x0,%ebx
and    %bpl,%cl
mov    -0x70(%rsp),%r10
mov    -0x74(%rsp),%edi
mov    -0x68(%rsp),%rbp
jne    1b70 <deconvolve3_columns+0xa70>
cmpq   $0x0,-0x8(%rsp)
je     1bfa <deconvolve3_columns+0xafa>
movsd  (%r8,%rsi,8),%xmm2
unpcklpd %xmm2,%xmm2
mov    -0x18(%rsp),%rbx
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
movupd -0x10(%r14,%rcx,8),%xmm3
mulpd  %xmm0,%xmm3
movupd -0x10(%r12,%rcx,8),%xmm4
addpd  %xmm3,%xmm4
movupd (%r12,%rcx,8),%xmm3
divpd  %xmm2,%xmm4
movupd %xmm4,-0x10(%r12,%rcx,8)
movupd (%r14,%rcx,8),%xmm4
mulpd  %xmm0,%xmm4
addpd  %xmm3,%xmm4
divpd  %xmm2,%xmm4
movupd %xmm4,(%r12,%rcx,8)
add    $0x4,%rcx
add    $0xfffffffffffffffe,%rbx
jne    1ac0 <deconvolve3_columns+0x9c0>
testb  $0x1,-0x10(%rsp)
je     1b3e <deconvolve3_columns+0xa3e>
mov    -0x60(%rsp),%rbx
movupd (%rbx,%rcx,8),%xmm2
lea    (%rcx,%rbp,1),%rax
movupd (%rbx,%rax,8),%xmm3
mulpd  %xmm0,%xmm3
addpd  %xmm2,%xmm3
movsd  (%r8,%rsi,8),%xmm2
unpcklpd %xmm2,%xmm2
divpd  %xmm2,%xmm3
movupd %xmm3,(%rbx,%rcx,8)
mov    -0x50(%rsp),%rax
mov    %rax,%rbx
cmp    %r11,%rax
je     19c0 <deconvolve3_columns+0x8c0>
jmp    1b70 <deconvolve3_columns+0xa70>
xor    %ebx,%ebx
mov    -0x70(%rsp),%r10
mov    -0x74(%rsp),%edi
mov    -0x68(%rsp),%rbp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %rbx,%rcx
test   $0x1,%r11b
je     1ba2 <deconvolve3_columns+0xaa2>
lea    (%rbx,%rbp,1),%rax
mov    -0x60(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  (%rcx,%rbx,8),%xmm2
divsd  (%r8,%rsi,8),%xmm2
movsd  %xmm2,(%rcx,%rbx,8)
mov    %rbx,%rcx
or     $0x1,%rcx
not    %rbx
cmp    -0x20(%rsp),%rbx
je     19c0 <deconvolve3_columns+0x8c0>
movsd  -0x8(%r13,%rcx,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  (%r9,%rcx,8),%xmm2
divsd  (%r8,%rsi,8),%xmm2
movsd  %xmm2,(%r9,%rcx,8)
movsd  0x0(%r13,%rcx,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  0x8(%r9,%rcx,8),%xmm2
divsd  (%r8,%rsi,8),%xmm2
movsd  %xmm2,0x8(%r9,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%r11
jne    1bb0 <deconvolve3_columns+0xab0>
jmp    19c0 <deconvolve3_columns+0x8c0>
xor    %ecx,%ecx
testb  $0x1,-0x10(%rsp)
jne    1b10 <deconvolve3_columns+0xa10>
jmp    1b3e <deconvolve3_columns+0xa3e>

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

LLVM IR (-O2):
define void @deconvolve3_columns(i32 %0, i32 %1, i32 %2, double* nocapture %3, double* nocapture %4, double %5, double %6) {
%8 = icmp ne i32 %1, 0
%9 = icmp ne i32 %0, 0
%10 = and i1 %9, %8
br i1 %10, label %11, label %621
11: ; preds = %7
switch i32 %1, label %93 [
i32 1, label %12
i32 2, label %66
]
12: ; preds = %11
%13 = tail call double @llvm.fmuladd.f64(double %6, double 2.000000e+00, double %5)
%14 = icmp sgt i32 %0, 0
br i1 %14, label %15, label %621
15: ; preds = %12
%16 = zext i32 %0 to i64
%17 = icmp eq i32 %0, 1
br i1 %17, label %57, label %18
18: ; preds = %15
%19 = and i64 %16, 4294967294
%20 = insertelement <2 x double> poison, double %13, i64 0
%21 = shufflevector <2 x double> %20, <2 x double> poison, <2 x i32> zeroinitializer
%22 = add nsw i64 %19, -2
%23 = lshr exact i64 %22, 1
%24 = add nuw i64 %23, 1
%25 = and i64 %24, 1
%26 = icmp eq i64 %22, 0
br i1 %26, label %46, label %27
27: ; preds = %18
%28 = and i64 %24, -2
br label %29
29: ; preds = %29, %27
%30 = phi i64 [ 0, %27 ], [ %43, %29 ]
%31 = phi i64 [ 0, %27 ], [ %44, %29 ]
%32 = getelementptr inbounds double, double* %3, i64 %30
%33 = bitcast double* %32 to <2 x double>*
%34 = load <2 x double>, <2 x double>* %33
%35 = fdiv <2 x double> %34, %21
%36 = bitcast double* %32 to <2 x double>*
store <2 x double> %35, <2 x double>* %36
%37 = or i64 %30, 2
%38 = getelementptr inbounds double, double* %3, i64 %37
%39 = bitcast double* %38 to <2 x double>*
%40 = load <2 x double>, <2 x double>* %39
%41 = fdiv <2 x double> %40, %21
%42 = bitcast double* %38 to <2 x double>*
store <2 x double> %41, <2 x double>* %42
%43 = add nuw i64 %30, 4
%44 = add nuw i64 %31, 2
%45 = icmp eq i64 %44, %28
br i1 %45, label %46, label %29
46: ; preds = %29, %18
%47 = phi i64 [ 0, %18 ], [ %43, %29 ]
%48 = icmp eq i64 %25, 0
br i1 %48, label %55, label %49
49: ; preds = %46
%50 = getelementptr inbounds double, double* %3, i64 %47
%51 = bitcast double* %50 to <2 x double>*
%52 = load <2 x double>, <2 x double>* %51
%53 = fdiv <2 x double> %52, %21
%54 = bitcast double* %50 to <2 x double>*
store <2 x double> %53, <2 x double>* %54
br label %55
55: ; preds = %46, %49
%56 = icmp eq i64 %19, %16
br i1 %56, label %621, label %57
57: ; preds = %15, %55
%58 = phi i64 [ 0, %15 ], [ %19, %55 ]
br label %59
59: ; preds = %57, %59
%60 = phi i64 [ %64, %59 ], [ %58, %57 ]
%61 = getelementptr inbounds double, double* %3, i64 %60
%62 = load double, double* %61
%63 = fdiv double %62, %13
store double %63, double* %61
%64 = add nuw nsw i64 %60, 1
%65 = icmp eq i64 %64, %16
br i1 %65, label %621, label %59
66: ; preds = %11
%67 = tail call double @llvm.fmuladd.f64(double %6, double 2.000000e+00, double %5)
%68 = fmul double %67, %5
%69 = icmp sgt i32 %0, 0
br i1 %69, label %70, label %621
70: ; preds = %66
%71 = fadd double %5, %6
%72 = fdiv double %71, %68
%73 = fneg double %6
%74 = fdiv double %73, %68
%75 = sext i32 %2 to i64
%76 = zext i32 %0 to i64
br label %77
77: ; preds = %70, %77
%78 = phi i64 [ 0, %70 ], [ %91, %77 ]
%79 = getelementptr inbounds double, double* %3, i64 %78
%80 = load double, double* %79
%81 = add nsw i64 %78, %75
%82 = getelementptr inbounds double, double* %3, i64 %81
%83 = load double, double* %82
%84 = fmul double %74, %83
%85 = tail call double @llvm.fmuladd.f64(double %72, double %80, double %84)
store double %85, double* %4
%86 = load double, double* %82
%87 = load double, double* %79
%88 = fmul double %74, %87
%89 = tail call double @llvm.fmuladd.f64(double %72, double %86, double %88)
store double %89, double* %82
%90 = load double, double* %4
store double %90, double* %79
%91 = add nuw nsw i64 %78, 1
%92 = icmp eq i64 %91, %76
br i1 %92, label %621, label %77
93: ; preds = %11
%94 = fadd double %5, %6
store double %94, double* %4
%95 = add i32 %1, -1
%96 = icmp sgt i32 %1, 2
br i1 %96, label %97, label %245
97: ; preds = %93
%98 = fneg double %6
%99 = icmp sgt i32 %0, 0
%100 = sext i32 %2 to i64
%101 = zext i32 %95 to i64
%102 = zext i32 %0 to i64
%103 = add nsw i64 %100, %102
%104 = and i64 %102, 4294967292
%105 = add nsw i64 %104, -4
%106 = lshr exact i64 %105, 2
%107 = add nuw nsw i64 %106, 1
%108 = icmp ult i32 %0, 4
%109 = and i64 %102, 4294967292
%110 = and i64 %107, 1
%111 = icmp eq i64 %105, 0
%112 = and i64 %107, 9223372036854775806
%113 = icmp eq i64 %110, 0
%114 = icmp eq i64 %109, %102
%115 = and i64 %102, 1
%116 = icmp eq i64 %115, 0
%117 = sub nsw i64 0, %102
br label %118
118: ; preds = %242, %97
%119 = phi i64 [ %244, %242 ], [ 0, %97 ]
%120 = phi double [ %243, %242 ], [ %94, %97 ]
%121 = phi i64 [ %240, %242 ], [ 1, %97 ]
%122 = mul i64 %119, %100
%123 = add i64 %122, %100
%124 = getelementptr double, double* %3, i64 %123
%125 = add i64 %103, %122
%126 = getelementptr double, double* %3, i64 %125
%127 = getelementptr double, double* %3, i64 %122
%128 = add i64 %122, %102
%129 = getelementptr double, double* %3, i64 %128
%130 = add nsw i64 %121, -1
%131 = fdiv double %98, %120
%132 = tail call double @llvm.fmuladd.f64(double %131, double %6, double %5)
%133 = getelementptr inbounds double, double* %4, i64 %121
store double %132, double* %133
%134 = mul nsw i64 %130, %100
%135 = getelementptr inbounds double, double* %3, i64 %134
br i1 %99, label %136, label %239
136: ; preds = %118
br i1 %108, label %208, label %137
137: ; preds = %136
%138 = icmp ult double* %124, %129
%139 = icmp ult double* %127, %126
%140 = and i1 %138, %139
br i1 %140, label %208, label %141
141: ; preds = %137
%142 = insertelement <2 x double> poison, double %131, i64 0
%143 = shufflevector <2 x double> %142, <2 x double> poison, <2 x i32> zeroinitializer
%144 = insertelement <2 x double> poison, double %131, i64 0
%145 = shufflevector <2 x double> %144, <2 x double> poison, <2 x i32> zeroinitializer
br i1 %111, label %187, label %146
146: ; preds = %141, %146
%147 = phi i64 [ %184, %146 ], [ 0, %141 ]
%148 = phi i64 [ %185, %146 ], [ 0, %141 ]
%149 = getelementptr inbounds double, double* %135, i64 %147
%150 = bitcast double* %149 to <2 x double>*
%151 = load <2 x double>, <2 x double>* %150
%152 = getelementptr inbounds double, double* %149, i64 2
%153 = bitcast double* %152 to <2 x double>*
%154 = load <2 x double>, <2 x double>* %153
%155 = add nsw i64 %147, %100
%156 = getelementptr inbounds double, double* %135, i64 %155
%157 = bitcast double* %156 to <2 x double>*
%158 = load <2 x double>, <2 x double>* %157
%159 = getelementptr inbounds double, double* %156, i64 2
%160 = bitcast double* %159 to <2 x double>*
%161 = load <2 x double>, <2 x double>* %160
%162 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %143, <2 x double> %151, <2 x double> %158)
%163 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %145, <2 x double> %154, <2 x double> %161)
%164 = bitcast double* %156 to <2 x double>*
store <2 x double> %162, <2 x double>* %164
%165 = bitcast double* %159 to <2 x double>*
store <2 x double> %163, <2 x double>* %165
%166 = or i64 %147, 4
%167 = getelementptr inbounds double, double* %135, i64 %166
%168 = bitcast double* %167 to <2 x double>*
%169 = load <2 x double>, <2 x double>* %168
%170 = getelementptr inbounds double, double* %167, i64 2
%171 = bitcast double* %170 to <2 x double>*
%172 = load <2 x double>, <2 x double>* %171
%173 = add nsw i64 %166, %100
%174 = getelementptr inbounds double, double* %135, i64 %173
%175 = bitcast double* %174 to <2 x double>*
%176 = load <2 x double>, <2 x double>* %175
%177 = getelementptr inbounds double, double* %174, i64 2
%178 = bitcast double* %177 to <2 x double>*
%179 = load <2 x double>, <2 x double>* %178
%180 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %143, <2 x double> %169, <2 x double> %176)
%181 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %145, <2 x double> %172, <2 x double> %179)
%182 = bitcast double* %174 to <2 x double>*
store <2 x double> %180, <2 x double>* %182
%183 = bitcast double* %177 to <2 x double>*
store <2 x double> %181, <2 x double>* %183
%184 = add nuw i64 %147, 8
%185 = add i64 %148, 2
%186 = icmp eq i64 %185, %112
br i1 %186, label %187, label %146
187: ; preds = %146, %141
%188 = phi i64 [ 0, %141 ], [ %184, %146 ]
br i1 %113, label %207, label %189
189: ; preds = %187
%190 = getelementptr inbounds double, double* %135, i64 %188
%191 = bitcast double* %190 to <2 x double>*
%192 = load <2 x double>, <2 x double>* %191
%193 = getelementptr inbounds double, double* %190, i64 2
%194 = bitcast double* %193 to <2 x double>*
%195 = load <2 x double>, <2 x double>* %194
%196 = add nsw i64 %188, %100
%197 = getelementptr inbounds double, double* %135, i64 %196
%198 = bitcast double* %197 to <2 x double>*
%199 = load <2 x double>, <2 x double>* %198
%200 = getelementptr inbounds double, double* %197, i64 2
%201 = bitcast double* %200 to <2 x double>*
%202 = load <2 x double>, <2 x double>* %201
%203 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %143, <2 x double> %192, <2 x double> %199)
%204 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %145, <2 x double> %195, <2 x double> %202)
%205 = bitcast double* %197 to <2 x double>*
store <2 x double> %203, <2 x double>* %205
%206 = bitcast double* %200 to <2 x double>*
store <2 x double> %204, <2 x double>* %206
br label %207
207: ; preds = %187, %189
br i1 %114, label %239, label %208
208: ; preds = %137, %136, %207
%209 = phi i64 [ 0, %137 ], [ 0, %136 ], [ %109, %207 ]
%210 = xor i64 %209, -1
br i1 %116, label %219, label %211
211: ; preds = %208
%212 = getelementptr inbounds double, double* %135, i64 %209
%213 = load double, double* %212
%214 = add nsw i64 %209, %100
%215 = getelementptr inbounds double, double* %135, i64 %214
%216 = load double, double* %215
%217 = tail call double @llvm.fmuladd.f64(double %131, double %213, double %216)
store double %217, double* %215
%218 = or i64 %209, 1
br label %219
219: ; preds = %211, %208
%220 = phi i64 [ %209, %208 ], [ %218, %211 ]
%221 = icmp eq i64 %210, %117
br i1 %221, label %239, label %222
222: ; preds = %219, %222
%223 = phi i64 [ %237, %222 ], [ %220, %219 ]
%224 = getelementptr inbounds double, double* %135, i64 %223
%225 = load double, double* %224
%226 = add nsw i64 %223, %100
%227 = getelementptr inbounds double, double* %135, i64 %226
%228 = load double, double* %227
%229 = tail call double @llvm.fmuladd.f64(double %131, double %225, double %228)
store double %229, double* %227
%230 = add nuw nsw i64 %223, 1
%231 = getelementptr inbounds double, double* %135, i64 %230
%232 = load double, double* %231
%233 = add nsw i64 %230, %100
%234 = getelementptr inbounds double, double* %135, i64 %233
%235 = load double, double* %234
%236 = tail call double @llvm.fmuladd.f64(double %131, double %232, double %235)
store double %236, double* %234
%237 = add nuw nsw i64 %223, 2
%238 = icmp eq i64 %237, %102
br i1 %238, label %239, label %222
239: ; preds = %219, %222, %207, %118
%240 = add nuw nsw i64 %121, 1
%241 = icmp eq i64 %240, %101
br i1 %241, label %245, label %242
242: ; preds = %239
%243 = load double, double* %133
%244 = add i64 %119, 1
br label %118
245: ; preds = %239, %93
%246 = phi i32 [ 1, %93 ], [ %95, %239 ]
%247 = add nsw i32 %246, -1
%248 = sext i32 %247 to i64
%249 = getelementptr inbounds double, double* %4, i64 %248
%250 = load double, double* %249
%251 = fdiv double %6, %250
%252 = fsub double 1.000000e+00, %251
%253 = tail call double @llvm.fmuladd.f64(double %6, double %252, double %5)
%254 = zext i32 %246 to i64
%255 = getelementptr inbounds double, double* %4, i64 %254
store double %253, double* %255
%256 = mul nsw i32 %247, %2
%257 = sext i32 %256 to i64
%258 = getelementptr double, double* %3, i64 %257
%259 = icmp sgt i32 %0, 0
br i1 %259, label %262, label %260
260: ; preds = %245
%261 = sext i32 %2 to i64
br label %463
262: ; preds = %245
%263 = fneg double %251
%264 = sext i32 %2 to i64
%265 = zext i32 %0 to i64
%266 = icmp ult i32 %0, 4
br i1 %266, label %354, label %267
267: ; preds = %262
%268 = add nsw i64 %264, %257
%269 = getelementptr double, double* %3, i64 %268
%270 = add nsw i64 %268, %265
%271 = getelementptr double, double* %3, i64 %270
%272 = add nsw i64 %257, %265
%273 = getelementptr double, double* %3, i64 %272
%274 = icmp ult double* %269, %273
%275 = icmp ult double* %258, %271
%276 = and i1 %274, %275
br i1 %276, label %354, label %277
277: ; preds = %267
%278 = and i64 %265, 4294967292
%279 = insertelement <2 x double> poison, double %263, i64 0
%280 = shufflevector <2 x double> %279, <2 x double> poison, <2 x i32> zeroinitializer
%281 = insertelement <2 x double> poison, double %263, i64 0
%282 = shufflevector <2 x double> %281, <2 x double> poison, <2 x i32> zeroinitializer
%283 = add nsw i64 %278, -4
%284 = lshr exact i64 %283, 2
%285 = add nuw nsw i64 %284, 1
%286 = and i64 %285, 1
%287 = icmp eq i64 %283, 0
br i1 %287, label %331, label %288
288: ; preds = %277
%289 = and i64 %285, 9223372036854775806
br label %290
290: ; preds = %290, %288
%291 = phi i64 [ 0, %288 ], [ %328, %290 ]
%292 = phi i64 [ 0, %288 ], [ %329, %290 ]
%293 = getelementptr inbounds double, double* %258, i64 %291
%294 = bitcast double* %293 to <2 x double>*
%295 = load <2 x double>, <2 x double>* %294
%296 = getelementptr inbounds double, double* %293, i64 2
%297 = bitcast double* %296 to <2 x double>*
%298 = load <2 x double>, <2 x double>* %297
%299 = add nsw i64 %291, %264
%300 = getelementptr inbounds double, double* %258, i64 %299
%301 = bitcast double* %300 to <2 x double>*
%302 = load <2 x double>, <2 x double>* %301
%303 = getelementptr inbounds double, double* %300, i64 2
%304 = bitcast double* %303 to <2 x double>*
%305 = load <2 x double>, <2 x double>* %304
%306 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %280, <2 x double> %295, <2 x double> %302)
%307 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %282, <2 x double> %298, <2 x double> %305)
%308 = bitcast double* %300 to <2 x double>*
store <2 x double> %306, <2 x double>* %308
%309 = bitcast double* %303 to <2 x double>*
store <2 x double> %307, <2 x double>* %309
%310 = or i64 %291, 4
%311 = getelementptr inbounds double, double* %258, i64 %310
%312 = bitcast double* %311 to <2 x double>*
%313 = load <2 x double>, <2 x double>* %312
%314 = getelementptr inbounds double, double* %311, i64 2
%315 = bitcast double* %314 to <2 x double>*
%316 = load <2 x double>, <2 x double>* %315
%317 = add nsw i64 %310, %264
%318 = getelementptr inbounds double, double* %258, i64 %317
%319 = bitcast double* %318 to <2 x double>*
%320 = load <2 x double>, <2 x double>* %319
%321 = getelementptr inbounds double, double* %318, i64 2
%322 = bitcast double* %321 to <2 x double>*
%323 = load <2 x double>, <2 x double>* %322
%324 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %280, <2 x double> %313, <2 x double> %320)
%325 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %282, <2 x double> %316, <2 x double> %323)
%326 = bitcast double* %318 to <2 x double>*
store <2 x double> %324, <2 x double>* %326
%327 = bitcast double* %321 to <2 x double>*
store <2 x double> %325, <2 x double>* %327
%328 = add nuw i64 %291, 8
%329 = add i64 %292, 2
%330 = icmp eq i64 %329, %289
br i1 %330, label %331, label %290
331: ; preds = %290, %277
%332 = phi i64 [ 0, %277 ], [ %328, %290 ]
%333 = icmp eq i64 %286, 0
br i1 %333, label %352, label %334
334: ; preds = %331
%335 = getelementptr inbounds double, double* %258, i64 %332
%336 = bitcast double* %335 to <2 x double>*
%337 = load <2 x double>, <2 x double>* %336
%338 = getelementptr inbounds double, double* %335, i64 2
%339 = bitcast double* %338 to <2 x double>*
%340 = load <2 x double>, <2 x double>* %339
%341 = add nsw i64 %332, %264
%342 = getelementptr inbounds double, double* %258, i64 %341
%343 = bitcast double* %342 to <2 x double>*
%344 = load <2 x double>, <2 x double>* %343
%345 = getelementptr inbounds double, double* %342, i64 2
%346 = bitcast double* %345 to <2 x double>*
%347 = load <2 x double>, <2 x double>* %346
%348 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %280, <2 x double> %337, <2 x double> %344)
%349 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %282, <2 x double> %340, <2 x double> %347)
%350 = bitcast double* %342 to <2 x double>*
store <2 x double> %348, <2 x double>* %350
%351 = bitcast double* %345 to <2 x double>*
store <2 x double> %349, <2 x double>* %351
br label %352
352: ; preds = %331, %334
%353 = icmp eq i64 %278, %265
br i1 %353, label %388, label %354
354: ; preds = %267, %262, %352
%355 = phi i64 [ 0, %267 ], [ 0, %262 ], [ %278, %352 ]
%356 = xor i64 %355, -1
%357 = and i64 %265, 1
%358 = icmp eq i64 %357, 0
br i1 %358, label %367, label %359
359: ; preds = %354
%360 = getelementptr inbounds double, double* %258, i64 %355
%361 = load double, double* %360
%362 = add nsw i64 %355, %264
%363 = getelementptr inbounds double, double* %258, i64 %362
%364 = load double, double* %363
%365 = tail call double @llvm.fmuladd.f64(double %263, double %361, double %364)
store double %365, double* %363
%366 = or i64 %355, 1
br label %367
367: ; preds = %359, %354
%368 = phi i64 [ %355, %354 ], [ %366, %359 ]
%369 = sub nsw i64 0, %265
%370 = icmp eq i64 %356, %369
br i1 %370, label %388, label %371
371: ; preds = %367, %371
%372 = phi i64 [ %386, %371 ], [ %368, %367 ]
%373 = getelementptr inbounds double, double* %258, i64 %372
%374 = load double, double* %373
%375 = add nsw i64 %372, %264
%376 = getelementptr inbounds double, double* %258, i64 %375
%377 = load double, double* %376
%378 = tail call double @llvm.fmuladd.f64(double %263, double %374, double %377)
store double %378, double* %376
%379 = add nuw nsw i64 %372, 1
%380 = getelementptr inbounds double, double* %258, i64 %379
%381 = load double, double* %380
%382 = add nsw i64 %379, %264
%383 = getelementptr inbounds double, double* %258, i64 %382
%384 = load double, double* %383
%385 = tail call double @llvm.fmuladd.f64(double %263, double %381, double %384)
store double %385, double* %383
%386 = add nuw nsw i64 %372, 2
%387 = icmp eq i64 %386, %265
br i1 %387, label %388, label %371
388: ; preds = %367, %371, %352
%389 = getelementptr double, double* %258, i64 %264
br i1 %259, label %390, label %463
390: ; preds = %388
%391 = zext i32 %0 to i64
%392 = icmp eq i32 %0, 1
br i1 %392, label %448, label %393
393: ; preds = %390
%394 = add nsw i64 %264, %257
%395 = add nsw i64 %394, %265
%396 = getelementptr double, double* %3, i64 %395
%397 = add nuw nsw i64 %254, 1
%398 = getelementptr double, double* %4, i64 %397
%399 = icmp ult double* %389, %398
%400 = icmp ult double* %255, %396
%401 = and i1 %399, %400
br i1 %401, label %448, label %402
402: ; preds = %393
%403 = and i64 %265, 4294967294
%404 = add nsw i64 %403, -2
%405 = lshr exact i64 %404, 1
%406 = add nuw i64 %405, 1
%407 = and i64 %406, 1
%408 = icmp eq i64 %404, 0
br i1 %408, label %434, label %409
409: ; preds = %402
%410 = and i64 %406, -2
%411 = load double, double* %255
%412 = insertelement <2 x double> poison, double %411, i64 0
%413 = shufflevector <2 x double> %412, <2 x double> poison, <2 x i32> zeroinitializer
%414 = load double, double* %255
%415 = insertelement <2 x double> poison, double %414, i64 0
%416 = shufflevector <2 x double> %415, <2 x double> poison, <2 x i32> zeroinitializer
br label %417
417: ; preds = %417, %409
%418 = phi i64 [ 0, %409 ], [ %431, %417 ]
%419 = phi i64 [ 0, %409 ], [ %432, %417 ]
%420 = getelementptr inbounds double, double* %389, i64 %418
%421 = bitcast double* %420 to <2 x double>*
%422 = load <2 x double>, <2 x double>* %421
%423 = fdiv <2 x double> %422, %413
%424 = bitcast double* %420 to <2 x double>*
store <2 x double> %423, <2 x double>* %424
%425 = or i64 %418, 2
%426 = getelementptr inbounds double, double* %389, i64 %425
%427 = bitcast double* %426 to <2 x double>*
%428 = load <2 x double>, <2 x double>* %427
%429 = fdiv <2 x double> %428, %416
%430 = bitcast double* %426 to <2 x double>*
store <2 x double> %429, <2 x double>* %430
%431 = add nuw i64 %418, 4
%432 = add nuw i64 %419, 2
%433 = icmp eq i64 %432, %410
br i1 %433, label %434, label %417
434: ; preds = %417, %402
%435 = phi i64 [ 0, %402 ], [ %431, %417 ]
%436 = icmp eq i64 %407, 0
br i1 %436, label %446, label %437
437: ; preds = %434
%438 = load double, double* %255
%439 = insertelement <2 x double> poison, double %438, i64 0
%440 = shufflevector <2 x double> %439, <2 x double> poison, <2 x i32> zeroinitializer
%441 = getelementptr inbounds double, double* %389, i64 %435
%442 = bitcast double* %441 to <2 x double>*
%443 = load <2 x double>, <2 x double>* %442
%444 = fdiv <2 x double> %443, %440
%445 = bitcast double* %441 to <2 x double>*
store <2 x double> %444, <2 x double>* %445
br label %446
446: ; preds = %434, %437
%447 = icmp eq i64 %403, %265
br i1 %447, label %463, label %448
448: ; preds = %393, %390, %446
%449 = phi i64 [ 0, %393 ], [ 0, %390 ], [ %403, %446 ]
%450 = xor i64 %449, -1
%451 = and i64 %265, 1
%452 = icmp eq i64 %451, 0
br i1 %452, label %459, label %453
453: ; preds = %448
%454 = load double, double* %255
%455 = getelementptr inbounds double, double* %389, i64 %449
%456 = load double, double* %455
%457 = fdiv double %456, %454
store double %457, double* %455
%458 = or i64 %449, 1
br label %459
459: ; preds = %453, %448
%460 = phi i64 [ %449, %448 ], [ %458, %453 ]
%461 = sub nsw i64 0, %265
%462 = icmp eq i64 %450, %461
br i1 %462, label %463, label %488
463: ; preds = %459, %488, %446, %260, %388
%464 = phi i64 [ %261, %260 ], [ %264, %388 ], [ %264, %446 ], [ %264, %488 ], [ %264, %459 ]
%465 = fneg double %6
%466 = zext i32 %0 to i64
%467 = add nsw i64 %254, -1
%468 = mul i64 %464, %467
%469 = add i64 %468, %466
%470 = mul nsw i64 %464, %254
%471 = add i64 %470, %466
%472 = and i64 %466, 4294967294
%473 = add nsw i64 %472, -2
%474 = lshr exact i64 %473, 1
%475 = add nuw i64 %474, 1
%476 = icmp eq i32 %0, 1
%477 = and i64 %466, 4294967294
%478 = insertelement <2 x double> poison, double %465, i64 0
%479 = shufflevector <2 x double> %478, <2 x double> poison, <2 x i32> zeroinitializer
%480 = and i64 %475, 1
%481 = icmp eq i64 %473, 0
%482 = and i64 %475, -2
%483 = icmp eq i64 %480, 0
%484 = icmp eq i64 %477, %466
%485 = and i64 %466, 1
%486 = icmp eq i64 %485, 0
%487 = sub nsw i64 0, %466
br label %501
488: ; preds = %459, %488
%489 = phi i64 [ %499, %488 ], [ %460, %459 ]
%490 = load double, double* %255
%491 = getelementptr inbounds double, double* %389, i64 %489
%492 = load double, double* %491
%493 = fdiv double %492, %490
store double %493, double* %491
%494 = add nuw nsw i64 %489, 1
%495 = load double, double* %255
%496 = getelementptr inbounds double, double* %389, i64 %494
%497 = load double, double* %496
%498 = fdiv double %497, %495
store double %498, double* %496
%499 = add nuw nsw i64 %489, 2
%500 = icmp eq i64 %499, %391
br i1 %500, label %463, label %488
501: ; preds = %463, %618
%502 = phi i64 [ 0, %463 ], [ %620, %618 ]
%503 = phi i64 [ %254, %463 ], [ %516, %618 ]
%504 = mul i64 %464, %502
%505 = sub i64 %468, %504
%506 = getelementptr double, double* %3, i64 %505
%507 = sub i64 %469, %504
%508 = getelementptr double, double* %3, i64 %507
%509 = sub i64 %470, %504
%510 = getelementptr double, double* %3, i64 %509
%511 = sub i64 %471, %504
%512 = getelementptr double, double* %3, i64 %511
%513 = sub i64 %467, %502
%514 = getelementptr double, double* %4, i64 %513
%515 = getelementptr double, double* %4, i64 %503
%516 = add nsw i64 %503, -1
%517 = mul nsw i64 %516, %464
%518 = getelementptr inbounds double, double* %3, i64 %517
br i1 %259, label %519, label %618
519: ; preds = %501
%520 = getelementptr inbounds double, double* %4, i64 %516
br i1 %476, label %581, label %521
521: ; preds = %519
%522 = icmp ult double* %506, %512
%523 = icmp ult double* %510, %508
%524 = and i1 %522, %523
%525 = icmp ult double* %506, %515
%526 = icmp ult double* %514, %508
%527 = and i1 %525, %526
%528 = or i1 %524, %527
br i1 %528, label %581, label %529
529: ; preds = %521
br i1 %481, label %564, label %530
530: ; preds = %529
%531 = load double, double* %520
%532 = insertelement <2 x double> poison, double %531, i64 0
%533 = shufflevector <2 x double> %532, <2 x double> poison, <2 x i32> zeroinitializer
%534 = load double, double* %520
%535 = insertelement <2 x double> poison, double %534, i64 0
%536 = shufflevector <2 x double> %535, <2 x double> poison, <2 x i32> zeroinitializer
br label %537
537: ; preds = %537, %530
%538 = phi i64 [ 0, %530 ], [ %561, %537 ]
%539 = phi i64 [ 0, %530 ], [ %562, %537 ]
%540 = getelementptr inbounds double, double* %518, i64 %538
%541 = bitcast double* %540 to <2 x double>*
%542 = load <2 x double>, <2 x double>* %541
%543 = add nsw i64 %538, %464
%544 = getelementptr inbounds double, double* %518, i64 %543
%545 = bitcast double* %544 to <2 x double>*
%546 = load <2 x double>, <2 x double>* %545
%547 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %479, <2 x double> %546, <2 x double> %542)
%548 = fdiv <2 x double> %547, %533
%549 = bitcast double* %540 to <2 x double>*
store <2 x double> %548, <2 x double>* %549
%550 = or i64 %538, 2
%551 = getelementptr inbounds double, double* %518, i64 %550
%552 = bitcast double* %551 to <2 x double>*
%553 = load <2 x double>, <2 x double>* %552
%554 = add nsw i64 %550, %464
%555 = getelementptr inbounds double, double* %518, i64 %554
%556 = bitcast double* %555 to <2 x double>*
%557 = load <2 x double>, <2 x double>* %556
%558 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %479, <2 x double> %557, <2 x double> %553)
%559 = fdiv <2 x double> %558, %536
%560 = bitcast double* %551 to <2 x double>*
store <2 x double> %559, <2 x double>* %560
%561 = add nuw i64 %538, 4
%562 = add nuw i64 %539, 2
%563 = icmp eq i64 %562, %482
br i1 %563, label %564, label %537
564: ; preds = %537, %529
%565 = phi i64 [ 0, %529 ], [ %561, %537 ]
br i1 %483, label %580, label %566
566: ; preds = %564
%567 = getelementptr inbounds double, double* %518, i64 %565
%568 = bitcast double* %567 to <2 x double>*
%569 = load <2 x double>, <2 x double>* %568
%570 = add nsw i64 %565, %464
%571 = getelementptr inbounds double, double* %518, i64 %570
%572 = bitcast double* %571 to <2 x double>*
%573 = load <2 x double>, <2 x double>* %572
%574 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %479, <2 x double> %573, <2 x double> %569)
%575 = load double, double* %520
%576 = insertelement <2 x double> poison, double %575, i64 0
%577 = shufflevector <2 x double> %576, <2 x double> poison, <2 x i32> zeroinitializer
%578 = fdiv <2 x double> %574, %577
%579 = bitcast double* %567 to <2 x double>*
store <2 x double> %578, <2 x double>* %579
br label %580
580: ; preds = %564, %566
br i1 %484, label %618, label %581
581: ; preds = %521, %519, %580
%582 = phi i64 [ 0, %521 ], [ 0, %519 ], [ %477, %580 ]
%583 = xor i64 %582, -1
br i1 %486, label %594, label %584
584: ; preds = %581
%585 = getelementptr inbounds double, double* %518, i64 %582
%586 = load double, double* %585
%587 = add nsw i64 %582, %464
%588 = getelementptr inbounds double, double* %518, i64 %587
%589 = load double, double* %588
%590 = tail call double @llvm.fmuladd.f64(double %465, double %589, double %586)
%591 = load double, double* %520
%592 = fdiv double %590, %591
store double %592, double* %585
%593 = or i64 %582, 1
br label %594
594: ; preds = %584, %581
%595 = phi i64 [ %582, %581 ], [ %593, %584 ]
%596 = icmp eq i64 %583, %487
br i1 %596, label %618, label %597
597: ; preds = %594, %597
%598 = phi i64 [ %616, %597 ], [ %595, %594 ]
%599 = getelementptr inbounds double, double* %518, i64 %598
%600 = load double, double* %599
%601 = add nsw i64 %598, %464
%602 = getelementptr inbounds double, double* %518, i64 %601
%603 = load double, double* %602
%604 = tail call double @llvm.fmuladd.f64(double %465, double %603, double %600)
%605 = load double, double* %520
%606 = fdiv double %604, %605
store double %606, double* %599
%607 = add nuw nsw i64 %598, 1
%608 = getelementptr inbounds double, double* %518, i64 %607
%609 = load double, double* %608
%610 = add nsw i64 %607, %464
%611 = getelementptr inbounds double, double* %518, i64 %610
%612 = load double, double* %611
%613 = tail call double @llvm.fmuladd.f64(double %465, double %612, double %609)
%614 = load double, double* %520
%615 = fdiv double %613, %614
store double %615, double* %608
%616 = add nuw nsw i64 %598, 2
%617 = icmp eq i64 %616, %466
br i1 %617, label %618, label %597
618: ; preds = %594, %597, %580, %501
%619 = icmp sgt i64 %503, 1
%620 = add i64 %502, 1
br i1 %619, label %501, label %621
621: ; preds = %77, %59, %618, %55, %66, %12, %7
ret void
}