void deconvolve3_rows(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 (width == 1) {
    q = a + 2.0 * b;
    for (i = 0; i < height; i++)
      data[i * rowstride] /= q;
    return;
  }
  if (width == 2) {
    q = a * (a + 2.0 * b);
    for (i = 0; i < height; i++) {
      row = data + i * rowstride;
      buffer[0] = (a + b) / q * row[0] - b / q * row[1];
      row[1] = (a + b) / q * row[1] - b / q * row[0];
      row[0] = buffer[0];
    }
    return;
  }
  buffer[0] = a + b;
  for (j = 1; j < width - 1; j++) {
    q = b / buffer[j - 1];
    buffer[j] = a - q * b;
    data[j] -= q * data[j - 1];
  }
  q = b / buffer[j - 1];
  buffer[j] = a + b * (1.0 - q);
  data[j] -= q * data[j - 1];
  data[j] /= buffer[j];
  do {
    j--;
    data[j] = (data[j] - b * data[j + 1]) / buffer[j];
  } while (j > 0);
  for (i = 1; i < height; i++) {
    row = data + i * rowstride;
    for (j = 1; j < width - 1; j++)
      row[j] -= b * row[j - 1] / buffer[j - 1];
    row[j] -= b * row[j - 1] / buffer[j - 1];
    row[j] /= buffer[j];
    do {
      j--;
      row[j] = (row[j] - b * row[j + 1]) / buffer[j];
    } while (j > 0);
  }
}


x86 assembly (clang, -O0):
<deconvolve3_rows>:
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_rows+0x33>
cmpl   $0x0,-0x4(%rbp)
jne    1138 <deconvolve3_rows+0x38>
jmp    162c <deconvolve3_rows+0x52c>
cmpl   $0x1,-0x4(%rbp)
jne    11a0 <deconvolve3_rows+0xa0>
movsd  -0x28(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm0
addsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    119b <deconvolve3_rows+0x9b>
movsd  -0x40(%rbp),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x44(%rbp),%ecx
imul   -0xc(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1160 <deconvolve3_rows+0x60>
jmp    162c <deconvolve3_rows+0x52c>
cmpl   $0x2,-0x4(%rbp)
jne    12bc <deconvolve3_rows+0x1bc>
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,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    12b7 <deconvolve3_rows+0x1b7>
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)
movsd  -0x28(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
divsd  -0x40(%rbp),%xmm0
mov    -0x38(%rbp),%rax
movsd  (%rax),%xmm2
movsd  -0x30(%rbp),%xmm1
divsd  -0x40(%rbp),%xmm1
mov    -0x38(%rbp),%rax
mulsd  0x8(%rax),%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    -0x38(%rbp),%rax
movsd  0x8(%rax),%xmm2
movsd  -0x30(%rbp),%xmm1
divsd  -0x40(%rbp),%xmm1
mov    -0x38(%rbp),%rax
mulsd  (%rax),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x38(%rbp),%rax
movsd  %xmm0,0x8(%rax)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x38(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    11d1 <deconvolve3_rows+0xd1>
jmp    162c <deconvolve3_rows+0x52c>
movsd  -0x28(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x48(%rbp)
mov    -0x48(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1393 <deconvolve3_rows+0x293>
movsd  -0x30(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x48(%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 -0x48(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x40(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x18(%rbp),%rax
movslq -0x48(%rbp),%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    12d5 <deconvolve3_rows+0x1d5>
movsd  -0x30(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x48(%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  0xc3f(%rip),%xmm2
subsd  -0x40(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x40(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x18(%rbp),%rax
movslq -0x48(%rbp),%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    -0x20(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%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    $0xffffffff,%eax
mov    %eax,-0x48(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x30(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    $0x1,%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 -0x48(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
cmpl   $0x0,-0x48(%rbp)
jg     1446 <deconvolve3_rows+0x346>
movl   $0x1,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    162c <deconvolve3_rows+0x52c>
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   $0x1,-0x48(%rbp)
mov    -0x48(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1547 <deconvolve3_rows+0x447>
movsd  -0x30(%rbp),%xmm1
mov    -0x38(%rbp),%rax
mov    -0x48(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
mov    -0x48(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm1
mov    -0x38(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    14e9 <deconvolve3_rows+0x3e9>
movsd  -0x30(%rbp),%xmm1
mov    -0x38(%rbp),%rax
mov    -0x48(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
mov    -0x48(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm1
mov    -0x38(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
movslq -0x48(%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    $0xffffffff,%eax
mov    %eax,-0x48(%rbp)
mov    -0x38(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x30(%rbp),%xmm0
mov    -0x38(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    $0x1,%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 -0x48(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
cmpl   $0x0,-0x48(%rbp)
jg     15a9 <deconvolve3_rows+0x4a9>
jmp    161e <deconvolve3_rows+0x51e>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    14bd <deconvolve3_rows+0x3bd>
pop    %rbp
ret

x86 assembly (clang, -O2):
<deconvolve3_rows>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %edi,%edi
je     145e <deconvolve3_rows+0x35e>
test   %esi,%esi
je     145e <deconvolve3_rows+0x35e>
cmp    $0x2,%edi
je     114e <deconvolve3_rows+0x4e>
cmp    $0x1,%edi
jne    11db <deconvolve3_rows+0xdb>
test   %esi,%esi
jle    145e <deconvolve3_rows+0x35e>
addsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
movslq %edx,%rax
mov    %esi,%r8d
cmp    $0x1,%esi
jne    1406 <deconvolve3_rows+0x306>
xor    %esi,%esi
jmp    1446 <deconvolve3_rows+0x346>
test   %esi,%esi
jle    145e <deconvolve3_rows+0x35e>
movapd %xmm1,%xmm2
addsd  %xmm1,%xmm2
addsd  %xmm0,%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm1,%xmm0
divsd  %xmm2,%xmm0
xorpd  0xe8a(%rip),%xmm1
divsd  %xmm2,%xmm1
movslq %edx,%rax
mov    %esi,%edx
add    $0x8,%rcx
shl    $0x3,%rax
nopw   0x0(%rax,%rax,1)
movsd  (%rcx),%xmm2
mulsd  %xmm1,%xmm2
movsd  -0x8(%rcx),%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm2,%xmm3
movsd  %xmm3,(%r8)
movsd  -0x8(%rcx),%xmm2
mulsd  %xmm1,%xmm2
movsd  (%rcx),%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm2,%xmm3
movsd  %xmm3,(%rcx)
movsd  (%r8),%xmm2
movsd  %xmm2,-0x8(%rcx)
add    %rax,%rcx
add    $0xffffffffffffffff,%rdx
jne    1190 <deconvolve3_rows+0x90>
jmp    145e <deconvolve3_rows+0x35e>
movapd %xmm0,%xmm3
addsd  %xmm1,%xmm3
movsd  %xmm3,(%r8)
lea    -0x1(%rdi),%r10d
movapd 0xe0c(%rip),%xmm2
xorpd  %xmm1,%xmm2
mov    $0x1,%eax
cmp    $0x3,%edi
jl     1251 <deconvolve3_rows+0x151>
mov    %r10d,%eax
add    $0xfffffffffffffffe,%rax
xor    %ebx,%ebx
nopl   0x0(%rax,%rax,1)
movapd %xmm2,%xmm4
divsd  %xmm3,%xmm4
movapd %xmm4,%xmm3
mulsd  %xmm1,%xmm3
addsd  %xmm0,%xmm3
movsd  %xmm3,0x8(%r8,%rbx,8)
mulsd  (%rcx,%rbx,8),%xmm4
addsd  0x8(%rcx,%rbx,8),%xmm4
movsd  %xmm4,0x8(%rcx,%rbx,8)
cmp    %rbx,%rax
je     124e <deconvolve3_rows+0x14e>
movsd  0x8(%r8,%rbx,8),%xmm3
add    $0x1,%rbx
jmp    1210 <deconvolve3_rows+0x110>
mov    %r10d,%eax
movslq %eax,%rbp
movapd %xmm1,%xmm3
divsd  -0x8(%r8,%rbp,8),%xmm3
movsd  0xda9(%rip),%xmm4
subsd  %xmm3,%xmm4
mulsd  %xmm1,%xmm4
addsd  %xmm0,%xmm4
mov    %ebp,%eax
movsd  %xmm4,(%r8,%rax,8)
movsd  (%rcx,%rax,8),%xmm0
mulsd  -0x8(%rcx,%rbp,8),%xmm3
subsd  %xmm3,%xmm0
movsd  %xmm0,(%rcx,%rax,8)
divsd  (%r8,%rax,8),%xmm0
movsd  %xmm0,(%rcx,%rax,8)
add    $0x1,%rax
xchg   %ax,%ax
mulsd  %xmm2,%xmm0
addsd  -0x10(%rcx,%rax,8),%xmm0
divsd  -0x10(%r8,%rax,8),%xmm0
movsd  %xmm0,-0x10(%rcx,%rax,8)
add    $0xffffffffffffffff,%rax
cmp    $0x1,%rax
ja     12a0 <deconvolve3_rows+0x1a0>
cmp    $0x2,%esi
jl     145e <deconvolve3_rows+0x35e>
movslq %edx,%r14
mov    %esi,%r11d
mov    %r10d,%r9d
add    $0xffffffffffffffff,%r9
lea    (%rcx,%r14,8),%rax
add    $0x8,%rax
lea    0x0(,%r14,8),%r15
mov    %r9,%r12
and    $0xfffffffffffffffe,%r12
neg    %r12
lea    (%rcx,%r14,8),%rdx
mov    $0x1,%r13d
nopl   0x0(%rax,%rax,1)
mov    %r13,%rsi
imul   %r14,%rsi
lea    (%rcx,%rsi,8),%rbx
mov    $0x1,%ebp
cmp    $0x3,%edi
jl     139a <deconvolve3_rows+0x29a>
movsd  (%rbx),%xmm0
mov    $0x1,%esi
cmp    $0x2,%r10d
je     1375 <deconvolve3_rows+0x275>
mov    $0x1,%esi
nopl   (%rax)
mulsd  %xmm1,%xmm0
divsd  -0x8(%r8,%rsi,8),%xmm0
movsd  -0x8(%rax,%rsi,8),%xmm3
subsd  %xmm0,%xmm3
movsd  (%rax,%rsi,8),%xmm0
movsd  %xmm3,-0x8(%rax,%rsi,8)
mulsd  %xmm1,%xmm3
divsd  (%r8,%rsi,8),%xmm3
subsd  %xmm3,%xmm0
movsd  %xmm0,(%rax,%rsi,8)
lea    (%r12,%rsi,1),%rbp
add    $0x2,%rbp
add    $0x2,%rsi
cmp    $0x1,%rbp
jne    1330 <deconvolve3_rows+0x230>
mov    %r10d,%ebp
test   $0x1,%r9b
je     139a <deconvolve3_rows+0x29a>
mulsd  %xmm1,%xmm0
divsd  -0x8(%r8,%rsi,8),%xmm0
movsd  (%rbx,%rsi,8),%xmm3
subsd  %xmm0,%xmm3
movsd  %xmm3,(%rbx,%rsi,8)
mov    %r10d,%ebp
movslq %ebp,%rsi
movsd  -0x8(%rbx,%rsi,8),%xmm3
mulsd  %xmm1,%xmm3
divsd  -0x8(%r8,%rsi,8),%xmm3
mov    %esi,%esi
movsd  (%rbx,%rsi,8),%xmm0
subsd  %xmm3,%xmm0
movsd  %xmm0,(%rbx,%rsi,8)
divsd  (%r8,%rsi,8),%xmm0
movsd  %xmm0,(%rbx,%rsi,8)
add    $0x1,%rsi
nopl   (%rax)
mulsd  %xmm2,%xmm0
addsd  -0x10(%rdx,%rsi,8),%xmm0
divsd  -0x10(%r8,%rsi,8),%xmm0
movsd  %xmm0,-0x10(%rdx,%rsi,8)
add    $0xffffffffffffffff,%rsi
cmp    $0x1,%rsi
ja     13d0 <deconvolve3_rows+0x2d0>
add    $0x1,%r13
add    %r15,%rax
add    %r15,%rdx
cmp    %r11,%r13
jne    1300 <deconvolve3_rows+0x200>
jmp    145e <deconvolve3_rows+0x35e>
mov    %r8d,%edi
and    $0xfffffffe,%edi
mov    %rax,%rbx
shl    $0x4,%rbx
xor    %esi,%esi
mov    %rcx,%rdx
nopl   0x0(%rax,%rax,1)
movsd  (%rdx),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rdx)
movsd  (%rdx,%rax,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rdx,%rax,8)
add    $0x2,%rsi
add    %rbx,%rdx
cmp    %rsi,%rdi
jne    1420 <deconvolve3_rows+0x320>
test   $0x1,%r8b
je     145e <deconvolve3_rows+0x35e>
imul   %rax,%rsi
movsd  (%rcx,%rsi,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rcx,%rsi,8)
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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