void bansol(double *e, double *y, double *c, int m, int n) {
  int i, k, nm1, mi, m2p1;
  double d;
  m2p1 = 2 * m + 1;
  nm1 = n - 1;
  if (m == 0) {
    for (i = 1; i <= n; i++) {
      c[i - 1] = y[i - 1] / e[(i - 1) * m2p1 + m];
    }
  } else {
    if (m == 1) {
      c[0] = y[0];
      for (i = 2; i <= n; i++) {
        c[i - 1] = y[i - 1] - e[(i - 1) * m2p1 - 1 + m] * c[i - 2];
      }
      c[n - 1] = c[n - 1] / e[(n - 1) * m2p1 + m];
      for (i = nm1; i >= 1; i--) {
        c[i - 1] = (c[i - 1] - e[(i - 1) * m2p1 + 1 + m] * c[i]) /
                   e[(i - 1) * m2p1 + m];
      }
    } else {
      c[0] = y[0];
      for (i = 2; i <= n; i++) {
        mi = (m < (i - 1)) ? m : i - 1;
        d = y[i - 1];
        for (k = 1; k <= mi; k++) {
          d -= e[(i - 1) * m2p1 - k + m] * c[i - k - 1];
        }
        c[i - 1] = d;
      }
      c[n - 1] /= e[(n - 1) * m2p1 + m];
      for (i = nm1; i >= 1; i--) {
        mi = (m < (n - i)) ? m : n - i;
        d = c[i - 1];
        for (k = 1; k <= mi; k++) {
          d -= e[(i - 1) * m2p1 + k + m] * c[i + k - 1];
        }
        c[i - 1] = d / e[(i - 1) * m2p1 + m];
      }
    }
  }
}


x86 assembly (clang, -O0):
<bansol>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    -0x1c(%rbp),%eax
shl    $0x1,%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x20(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x0,-0x1c(%rbp)
jne    1199 <bansol+0x99>
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x20(%rbp),%eax
jg     1194 <bansol+0x94>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x34(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    113d <bansol+0x3d>
jmp    1578 <bansol+0x478>
cmpl   $0x1,-0x1c(%rbp)
jne    132a <bansol+0x22a>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x2,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x20(%rbp),%eax
jg     1245 <bansol+0x145>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x34(%rbp),%ecx
sub    $0x1,%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x2,%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    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11ba <bansol+0xba>
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x34(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x1,-0x24(%rbp)
jl     1325 <bansol+0x225>
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x34(%rbp),%ecx
add    $0x1,%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%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    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x34(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    1288 <bansol+0x188>
jmp    1573 <bansol+0x473>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x2,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x20(%rbp),%eax
jg     142f <bansol+0x32f>
mov    -0x1c(%rbp),%eax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1369 <bansol+0x269>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x44(%rbp)
jmp    1372 <bansol+0x272>
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x30(%rbp),%eax
jg     140a <bansol+0x30a>
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x34(%rbp),%ecx
sub    -0x28(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    -0x28(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1396 <bansol+0x296>
movsd  -0x40(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1341 <bansol+0x241>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x34(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x1,-0x24(%rbp)
jl     156e <bansol+0x46e>
mov    -0x1c(%rbp),%eax
mov    -0x20(%rbp),%ecx
sub    -0x24(%rbp),%ecx
cmp    %ecx,%eax
jge    148f <bansol+0x38f>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x48(%rbp)
jmp    1498 <bansol+0x398>
mov    -0x20(%rbp),%eax
sub    -0x24(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x48(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x30(%rbp),%eax
jg     1530 <bansol+0x430>
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x34(%rbp),%ecx
add    -0x28(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x28(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    14bc <bansol+0x3bc>
movsd  -0x40(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x34(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    1469 <bansol+0x369>
jmp    1573 <bansol+0x473>
jmp    1578 <bansol+0x478>
pop    %rbp
ret

x86 assembly (clang, -O2):
<bansol>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
lea    (%rcx,%rcx,1),%r10d
add    $0x1,%r10d
lea    -0x1(%r8),%r15d
cmp    $0x1,%ecx
mov    %r10d,-0x28(%rsp)
je     1148 <bansol+0x48>
test   %ecx,%ecx
jne    126f <bansol+0x16f>
test   %r8d,%r8d
jle    15d6 <bansol+0x4d6>
movslq %r10d,%r9
mov    %r8d,%r11d
cmp    $0x1,%r8d
jne    1575 <bansol+0x475>
xor    %ebp,%ebp
jmp    15bc <bansol+0x4bc>
movsd  (%rsi),%xmm0
movsd  %xmm0,(%rdx)
cmp    $0x2,%r8d
jl     11f2 <bansol+0xf2>
movslq %r10d,%r14
lea    0x1(%r8),%r9d
mov    $0x1,%ebp
cmp    $0x3,%r9d
je     11cf <bansol+0xcf>
lea    -0x2(%r9),%r11
lea    -0x1(%r14),%r10
mov    %r10,%rax
shl    $0x4,%rax
add    $0x10,%rax
and    $0xfffffffffffffffe,%r11
neg    %r11
mov    $0x1,%ebp
mov    %rdi,%rcx
xchg   %ax,%ax
mulsd  0x8(%rcx,%r10,8),%xmm0
movsd  (%rsi,%rbp,8),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rbp,8)
movsd  0x8(%rsi,%rbp,8),%xmm0
mulsd  (%rcx,%rax,1),%xmm1
add    %rax,%rcx
subsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%rdx,%rbp,8)
lea    (%r11,%rbp,1),%rbx
add    $0x2,%rbx
add    $0x2,%rbp
cmp    $0x1,%rbx
jne    1190 <bansol+0x90>
test   $0x1,%r9b
mov    -0x28(%rsp),%r10d
je     11f2 <bansol+0xf2>
movsd  (%rsi,%rbp,8),%xmm1
imul   %rbp,%r14
mulsd  (%rdi,%r14,8),%xmm0
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rbp,8)
movslq %r15d,%rax
movsd  (%rdx,%rax,8),%xmm0
imul   %r10d,%r15d
movslq %r15d,%rcx
divsd  0x8(%rdi,%rcx,8),%xmm0
movsd  %xmm0,(%rdx,%rax,8)
cmp    $0x2,%r8d
jl     15d6 <bansol+0x4d6>
mov    %r8d,%eax
movslq %r10d,%rcx
lea    -0x2(%rax),%rsi
imul   %rcx,%rsi
lea    (%rdi,%rsi,8),%rsi
add    $0x10,%rsi
shl    $0x3,%rcx
neg    %rcx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mulsd  (%rsi),%xmm0
movsd  -0x10(%rdx,%rax,8),%xmm1
subsd  %xmm0,%xmm1
divsd  -0x8(%rsi),%xmm1
movsd  %xmm1,-0x10(%rdx,%rax,8)
add    $0xffffffffffffffff,%rax
add    %rcx,%rsi
movapd %xmm1,%xmm0
cmp    $0x1,%rax
ja     1240 <bansol+0x140>
jmp    15d6 <bansol+0x4d6>
mov    %r15d,-0x20(%rsp)
movsd  (%rsi),%xmm0
movsd  %xmm0,(%rdx)
mov    %r8,-0x30(%rsp)
cmp    $0x2,%r8d
jge    1453 <bansol+0x353>
mov    -0x20(%rsp),%esi
movslq %esi,%rax
imul   %r10d,%esi
add    %ecx,%esi
movslq %esi,%rsi
movsd  (%rdx,%rax,8),%xmm0
divsd  (%rdi,%rsi,8),%xmm0
movsd  %xmm0,(%rdx,%rax,8)
mov    -0x30(%rsp),%r9
cmp    $0x2,%r9d
jl     15d6 <bansol+0x4d6>
lea    (%rcx,%rcx,1),%ebp
movabs $0xffffffff00000000,%r11
ff ff ff
mov    %r9d,%r14d
mov    %r10d,%esi
lea    -0x1(%r14),%r10
movslq %esi,%rax
mov    %rax,-0x8(%rsp)
movslq %ecx,%r12
mov    %r14,-0x28(%rsp)
shl    $0x20,%r14
lea    -0x2(%r9),%eax
imul   %esi,%eax
add    %ecx,%eax
mov    $0xfffffffe,%esi
sub    %ebp,%esi
add    $0x1,%rsi
mov    %rsi,-0x20(%rsp)
mov    $0x1,%ebp
xor    %esi,%esi
movabs $0x200000000,%r13
mov    %r12,-0x10(%rsp)
jmp    134f <bansol+0x24f>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %r10,%rbx
imul   -0x8(%rsp),%rbx
mov    -0x38(%rsp),%esi
add    %r12,%rbx
divsd  (%rdi,%rbx,8),%xmm0
movsd  %xmm0,(%rdx,%r10,8)
add    $0x1,%ebp
add    %r11,%r14
add    -0x20(%rsp),%rax
cmp    %r9d,%ebp
je     15d6 <bansol+0x4d6>
mov    %r10,%r15
cmp    %ebp,%ecx
mov    %ebp,%r8d
cmovl  %ecx,%r8d
add    $0x1,%esi
cmp    %esi,%ecx
mov    %esi,-0x38(%rsp)
mov    %esi,%ebx
cmovl  %ecx,%ebx
mov    -0x28(%rsp),%rsi
sub    %r10,%rsi
cmp    %r12,%rsi
cmovg  %ecx,%esi
add    $0xffffffffffffffff,%r10
movsd  -0x8(%rdx,%r15,8),%xmm0
test   %esi,%esi
jle    1320 <bansol+0x220>
add    $0x1,%ebx
lea    -0x1(%rbx),%rsi
mov    $0x1,%r9d
cmp    $0x2,%ebx
mov    %rsi,-0x18(%rsp)
je     140e <bansol+0x30e>
add    $0x1,%r8d
add    $0xffffffffffffffff,%r8
and    $0xfffffffffffffffe,%r8
neg    %r8
mov    $0x1,%r9d
mov    %r14,%rbx
cs nopw 0x0(%rax,%rax,1)
nop
lea    (%rax,%r9,1),%rsi
movslq %esi,%rsi
movsd  (%rdi,%rsi,8),%xmm1
lea    (%rbx,%r11,1),%r12
sar    $0x1d,%r12
mulsd  (%rdx,%r12,1),%xmm1
subsd  %xmm1,%xmm0
add    $0x1,%esi
movslq %esi,%rsi
movsd  (%rdi,%rsi,8),%xmm1
mov    %rbx,%rsi
sar    $0x1d,%rsi
mulsd  (%rdx,%rsi,1),%xmm1
subsd  %xmm1,%xmm0
add    %r13,%rbx
lea    (%r8,%r9,1),%rsi
add    $0x2,%rsi
add    $0x2,%r9
cmp    $0x1,%rsi
jne    13c0 <bansol+0x2c0>
mov    %r10,%rbx
imul   -0x8(%rsp),%rbx
testb  $0x1,-0x18(%rsp)
je     1444 <bansol+0x344>
lea    (%rbx,%rcx,1),%r8d
mov    $0xffffffff,%esi
add    %rsi,%r15
add    %r9d,%r8d
movslq %r8d,%rsi
movsd  (%rdi,%rsi,8),%xmm1
add    %r9d,%r15d
movslq %r15d,%rsi
mulsd  (%rdx,%rsi,8),%xmm1
subsd  %xmm1,%xmm0
mov    -0x30(%rsp),%r9
mov    -0x10(%rsp),%r12
jmp    1329 <bansol+0x229>
movslq %ecx,%rax
mov    %rax,-0x38(%rsp)
mov    -0x30(%rsp),%rax
add    $0x1,%eax
mov    %rax,-0x18(%rsp)
lea    (%rcx,%rcx,2),%r12d
mov    $0x2,%r14d
xor    %r15d,%r15d
mov    $0x1,%r13d
mov    %rdx,%rbx
jmp    14a1 <bansol+0x3a1>
movsd  %xmm0,-0x8(%rdx,%r14,8)
add    $0x1,%r14
add    %r10d,%r12d
add    $0x8,%rbx
add    $0x1,%r13d
cmp    -0x18(%rsp),%r14
je     128b <bansol+0x18b>
cmp    %r13d,%ecx
mov    %r13d,%r9d
cmovl  %ecx,%r9d
add    $0x1,%r15d
cmp    %r15d,%ecx
mov    %r15d,%eax
cmovl  %ecx,%eax
lea    -0x1(%r14),%r8
cmp    -0x38(%rsp),%r8
mov    %r8d,%ebp
cmovg  %ecx,%ebp
movsd  -0x8(%rsi,%r14,8),%xmm0
test   %ebp,%ebp
jle    1480 <bansol+0x380>
add    $0x1,%eax
lea    -0x1(%rax),%r11
mov    $0x1,%ebp
cmp    $0x2,%eax
je     1545 <bansol+0x445>
add    $0x1,%r9d
add    $0xffffffffffffffff,%r9
and    $0xfffffffffffffffe,%r9
mov    $0xffffffffffffffff,%rbp
mov    %r12d,%r10d
nopl   0x0(%rax)
movslq %r10d,%r10
movsd  (%rdi,%r10,8),%xmm1
mulsd  0x8(%rbx,%rbp,8),%xmm1
subsd  %xmm1,%xmm0
lea    -0x1(%r10),%eax
cltq
movsd  (%rdi,%rax,8),%xmm1
mulsd  (%rbx,%rbp,8),%xmm1
subsd  %xmm1,%xmm0
add    $0xfffffffe,%r10d
lea    (%r9,%rbp,1),%rax
add    $0xfffffffffffffffe,%rax
add    $0xfffffffffffffffe,%rbp
cmp    $0xffffffffffffffff,%rax
jne    1500 <bansol+0x400>
neg    %rbp
mov    -0x28(%rsp),%r10d
test   $0x1,%r11b
je     1480 <bansol+0x380>
imul   %r10d,%r8d
add    %ecx,%r8d
sub    %ebp,%r8d
movslq %r8d,%rax
movsd  (%rdi,%rax,8),%xmm1
not    %rbp
add    %r14,%rbp
mulsd  (%rdx,%rbp,8),%xmm1
subsd  %xmm1,%xmm0
jmp    1480 <bansol+0x380>
mov    %r11d,%r10d
and    $0xfffffffe,%r10d
lea    -0x1(%r9),%rax
mov    %rax,%rbx
shl    $0x4,%rbx
add    $0x10,%rbx
xor    %ebp,%ebp
mov    %rdi,%rcx
movsd  (%rsi,%rbp,8),%xmm0
divsd  (%rcx),%xmm0
movsd  %xmm0,(%rdx,%rbp,8)
movsd  0x8(%rsi,%rbp,8),%xmm0
divsd  0x8(%rcx,%rax,8),%xmm0
movsd  %xmm0,0x8(%rdx,%rbp,8)
add    $0x2,%rbp
add    %rbx,%rcx
cmp    %rbp,%r10
jne    1590 <bansol+0x490>
test   $0x1,%r11b
je     15d6 <bansol+0x4d6>
movsd  (%rsi,%rbp,8),%xmm0
imul   %rbp,%r9
divsd  (%rdi,%r9,8),%xmm0
movsd  %xmm0,(%rdx,%rbp,8)
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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