void cftbsub(int n, double *a, double *w) {
  void cft1st(int n, double *a, double *w);
  void cftmdl(int n, int l, double *a, double *w);
  int j, j1, j2, j3, j4, j5, j6, j7, l;
  double wn4r, x0r, x0i, x1r, x1i, x2r, x2i, x3r, x3i, y0r, y0i, y1r, y1i, y2r,
      y2i, y3r, y3i, y4r, y4i, y5r, y5i, y6r, y6i, y7r, y7i;
  l = 2;
  if (n > 16) {
    cft1st(n, a, w);
    l = 16;
    while ((l << 3) < n) {
      cftmdl(n, l, a, w);
      l <<= 3;
    }
  }
  if ((l << 2) < n) {
    wn4r = w[2];
    for (j = 0; j < l; j += 2) {
      j1 = j + l;
      j2 = j1 + l;
      j3 = j2 + l;
      j4 = j3 + l;
      j5 = j4 + l;
      j6 = j5 + l;
      j7 = j6 + l;
      x0r = a[j] + a[j1];
      x0i = -a[j + 1] - a[j1 + 1];
      x1r = a[j] - a[j1];
      x1i = -a[j + 1] + a[j1 + 1];
      x2r = a[j2] + a[j3];
      x2i = a[j2 + 1] + a[j3 + 1];
      x3r = a[j2] - a[j3];
      x3i = a[j2 + 1] - a[j3 + 1];
      y0r = x0r + x2r;
      y0i = x0i - x2i;
      y2r = x0r - x2r;
      y2i = x0i + x2i;
      y1r = x1r - x3i;
      y1i = x1i - x3r;
      y3r = x1r + x3i;
      y3i = x1i + x3r;
      x0r = a[j4] + a[j5];
      x0i = a[j4 + 1] + a[j5 + 1];
      x1r = a[j4] - a[j5];
      x1i = a[j4 + 1] - a[j5 + 1];
      x2r = a[j6] + a[j7];
      x2i = a[j6 + 1] + a[j7 + 1];
      x3r = a[j6] - a[j7];
      x3i = a[j6 + 1] - a[j7 + 1];
      y4r = x0r + x2r;
      y4i = x0i + x2i;
      y6r = x0r - x2r;
      y6i = x0i - x2i;
      x0r = x1r - x3i;
      x0i = x1i + x3r;
      x2r = x1r + x3i;
      x2i = x1i - x3r;
      y5r = wn4r * (x0r - x0i);
      y5i = wn4r * (x0r + x0i);
      y7r = wn4r * (x2r - x2i);
      y7i = wn4r * (x2r + x2i);
      a[j1] = y1r + y5r;
      a[j1 + 1] = y1i - y5i;
      a[j5] = y1r - y5r;
      a[j5 + 1] = y1i + y5i;
      a[j3] = y3r - y7i;
      a[j3 + 1] = y3i - y7r;
      a[j7] = y3r + y7i;
      a[j7 + 1] = y3i + y7r;
      a[j] = y0r + y4r;
      a[j + 1] = y0i - y4i;
      a[j4] = y0r - y4r;
      a[j4 + 1] = y0i + y4i;
      a[j2] = y2r - y6i;
      a[j2 + 1] = y2i - y6r;
      a[j6] = y2r + y6i;
      a[j6 + 1] = y2i + y6r;
    }
  } else if ((l << 2) == n) {
    for (j = 0; j < l; j += 2) {
      j1 = j + l;
      j2 = j1 + l;
      j3 = j2 + l;
      x0r = a[j] + a[j1];
      x0i = -a[j + 1] - a[j1 + 1];
      x1r = a[j] - a[j1];
      x1i = -a[j + 1] + a[j1 + 1];
      x2r = a[j2] + a[j3];
      x2i = a[j2 + 1] + a[j3 + 1];
      x3r = a[j2] - a[j3];
      x3i = a[j2 + 1] - a[j3 + 1];
      a[j] = x0r + x2r;
      a[j + 1] = x0i - x2i;
      a[j2] = x0r - x2r;
      a[j2 + 1] = x0i + x2i;
      a[j1] = x1r - x3i;
      a[j1 + 1] = x1i - x3r;
      a[j3] = x1r + x3i;
      a[j3 + 1] = x1i + x3r;
    }
  } else {
    for (j = 0; j < l; j += 2) {
      j1 = j + l;
      x0r = a[j] - a[j1];
      x0i = -a[j + 1] + a[j1 + 1];
      a[j] += a[j1];
      a[j + 1] = -a[j + 1] - a[j1 + 1];
      a[j1] = x0r;
      a[j1 + 1] = x0i;
    }
  }
}


x86 assembly (clang, -O0):
<cftbsub>:
push   %rbp
mov    %rsp,%rbp
sub    $0x110,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x2,-0x3c(%rbp)
cmpl   $0x10,-0x4(%rbp)
jle    1193 <cftbsub+0x73>
mov    -0x4(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x18(%rbp),%rdx
call   1030 <cft1st@plt>
movl   $0x10,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
shl    $0x3,%eax
cmp    -0x4(%rbp),%eax
jge    118e <cftbsub+0x6e>
mov    -0x4(%rbp),%edi
mov    -0x3c(%rbp),%esi
mov    -0x10(%rbp),%rdx
mov    -0x18(%rbp),%rcx
call   1040 <cftmdl@plt>
mov    -0x3c(%rbp),%eax
shl    $0x3,%eax
mov    %eax,-0x3c(%rbp)
jmp    115e <cftbsub+0x3e>
jmp    1193 <cftbsub+0x73>
mov    -0x3c(%rbp),%eax
shl    $0x2,%eax
cmp    -0x4(%rbp),%eax
jge    180d <cftbsub+0x6ed>
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    1808 <cftbsub+0x6e8>
mov    -0x1c(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x50(%rbp),%xmm0
addsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x58(%rbp),%xmm0
subsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x50(%rbp),%xmm0
subsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0x58(%rbp),%xmm0
addsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0x60(%rbp),%xmm0
subsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
subsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0x60(%rbp),%xmm0
addsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x50(%rbp),%xmm0
addsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
movsd  -0x58(%rbp),%xmm0
addsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
movsd  -0x50(%rbp),%xmm0
subsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0xf0(%rbp)
ff
movsd  -0x58(%rbp),%xmm0
subsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0xf8(%rbp)
ff
movsd  -0x60(%rbp),%xmm0
subsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x50(%rbp)
movsd  -0x68(%rbp),%xmm0
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x60(%rbp),%xmm0
addsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x70(%rbp)
movsd  -0x68(%rbp),%xmm0
subsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm1
subsd  -0x58(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
movsd  -0x48(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm1
addsd  -0x58(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xe8(%rbp)
ff
movsd  -0x48(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm1
subsd  -0x78(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x100(%rbp)
ff
movsd  -0x48(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm1
addsd  -0x78(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x108(%rbp)
ff
movsd  -0xa0(%rbp),%xmm0
ff
addsd  -0xe0(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xa8(%rbp),%xmm0
ff
subsd  -0xe8(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xa0(%rbp),%xmm0
ff
subsd  -0xe0(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xa8(%rbp),%xmm0
ff
addsd  -0xe8(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xc0(%rbp),%xmm0
ff
subsd  -0x108(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xc8(%rbp),%xmm0
ff
subsd  -0x100(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xc0(%rbp),%xmm0
ff
addsd  -0x108(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xc8(%rbp),%xmm0
ff
addsd  -0x100(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x90(%rbp),%xmm0
ff
addsd  -0xd0(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x98(%rbp),%xmm0
ff
subsd  -0xd8(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x90(%rbp),%xmm0
ff
subsd  -0xd0(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x98(%rbp),%xmm0
ff
addsd  -0xd8(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xb0(%rbp),%xmm0
ff
subsd  -0xf8(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xb8(%rbp),%xmm0
ff
subsd  -0xf0(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xb0(%rbp),%xmm0
ff
addsd  -0xf8(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xb8(%rbp),%xmm0
ff
addsd  -0xf0(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x1c(%rbp)
jmp    11b7 <cftbsub+0x97>
jmp    1ba8 <cftbsub+0xa88>
mov    -0x3c(%rbp),%eax
shl    $0x2,%eax
cmp    -0x4(%rbp),%eax
jne    1a80 <cftbsub+0x960>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    1a7b <cftbsub+0x95b>
mov    -0x1c(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x50(%rbp),%xmm0
addsd  -0x70(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x58(%rbp),%xmm0
subsd  -0x78(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
subsd  -0x70(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x58(%rbp),%xmm0
addsd  -0x78(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x60(%rbp),%xmm0
subsd  -0x88(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x68(%rbp),%xmm0
subsd  -0x80(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x60(%rbp),%xmm0
addsd  -0x88(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x68(%rbp),%xmm0
addsd  -0x80(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x1c(%rbp)
jmp    1823 <cftbsub+0x703>
jmp    1ba3 <cftbsub+0xa83>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    1b9e <cftbsub+0xa7e>
mov    -0x1c(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x58(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x1c(%rbp)
jmp    1a87 <cftbsub+0x967>
jmp    1ba3 <cftbsub+0xa83>
jmp    1ba8 <cftbsub+0xa88>
add    $0x110,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cftbsub>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
mov    %rdx,%r14
mov    %rsi,%rbx
mov    %edi,%ebp
mov    $0x2,%r15d
cmp    $0x11,%edi
jl     118f <cftbsub+0x6f>
mov    %ebp,%edi
mov    %rbx,%rsi
mov    %r14,%rdx
call   1030 <cft1st@plt>
mov    $0x10,%r15d
cmp    $0x81,%ebp
jb     118f <cftbsub+0x6f>
mov    $0x10,%esi
mov    $0x80,%eax
cs nopw 0x0(%rax,%rax,1)
mov    %eax,%r15d
mov    %ebp,%edi
mov    %rbx,%rdx
mov    %r14,%rcx
call   1040 <cftmdl@plt>
lea    0x0(,%r15,8),%eax
mov    %r15d,%esi
cmp    %ebp,%eax
jl     1170 <cftbsub+0x50>
lea    0x0(,%r15,4),%eax
cmp    %ebp,%eax
jge    1514 <cftbsub+0x3f4>
test   %r15d,%r15d
jle    16fd <cftbsub+0x5dd>
movsd  0x10(%r14),%xmm0
movsd  %xmm0,0x50(%rsp)
mov    %r15d,%r12d
add    %r15d,%r15d
or     $0x1,%r15
mov    %r15,0x30(%rsp)
mov    %r12,%r8
shl    $0x20,%r8
lea    (%rbx,%r12,8),%rax
mov    %rax,0x40(%rsp)
mov    %r12,%rax
shl    $0x4,%rax
add    %rbx,%rax
mov    %rax,0x8(%rsp)
lea    (%r8,%r8,2),%r10
lea    (%r12,%r12,2),%rax
lea    (%rbx,%rax,8),%rcx
mov    %rcx,0x38(%rsp)
mov    %r12,%r11
shl    $0x22,%r11
imul   $0x38,%r12,%rcx
add    %rbx,%rcx
add    $0x8,%rcx
mov    %rcx,0x20(%rsp)
shl    $0x4,%rax
add    %rbx,%rax
add    $0x8,%rax
mov    %rax,0x18(%rsp)
lea    (%r12,%r12,4),%rax
lea    (%rbx,%rax,8),%rax
add    $0x8,%rax
mov    %rax,0x10(%rsp)
mov    %r12,0x48(%rsp)
shl    $0x5,%r12
add    %rbx,%r12
mov    %r12,0x28(%rsp)
xor    %esi,%esi
mov    0x40(%rsp),%r14
mov    0x38(%rsp),%r15
mov    0x28(%rsp),%r13
mov    0x10(%rsp),%r9
mov    0x18(%rsp),%rdi
mov    0x20(%rsp),%rcx
nopl   0x0(%rax)
movsd  (%r14,%rsi,8),%xmm3
movsd  (%rbx,%rsi,8),%xmm0
movsd  0x8(%rbx,%rsi,8),%xmm4
movapd %xmm0,%xmm7
addsd  %xmm3,%xmm7
movapd %xmm4,%xmm2
xorpd  0xd7b(%rip),%xmm2
mov    %r8,%rdx
sar    $0x1d,%rdx
or     $0x8,%rdx
movsd  (%rbx,%rdx,1),%xmm1
subsd  %xmm1,%xmm2
subsd  %xmm3,%xmm0
subsd  %xmm4,%xmm1
mov    0x8(%rsp),%rax
movsd  (%rax,%rsi,8),%xmm10
movsd  (%r15,%rsi,8),%xmm3
movapd %xmm10,%xmm5
addsd  %xmm3,%xmm5
mov    0x30(%rsp),%rax
add    %esi,%eax
movslq %eax,%r12
movsd  (%rbx,%r12,8),%xmm13
mov    %r10,%rbp
sar    $0x1d,%rbp
or     $0x8,%rbp
movsd  (%rbx,%rbp,1),%xmm6
movapd %xmm13,%xmm11
addsd  %xmm6,%xmm11
subsd  %xmm3,%xmm10
subsd  %xmm6,%xmm13
movapd %xmm7,%xmm14
addsd  %xmm5,%xmm14
movapd %xmm2,%xmm3
subsd  %xmm11,%xmm3
movapd %xmm3,0x60(%rsp)
subsd  %xmm5,%xmm7
movsd  %xmm7,0x58(%rsp)
addsd  %xmm2,%xmm11
movapd %xmm0,%xmm7
subsd  %xmm13,%xmm7
movapd %xmm1,%xmm12
subsd  %xmm10,%xmm12
addsd  %xmm0,%xmm13
addsd  %xmm1,%xmm10
movsd  0x0(%r13,%rsi,8),%xmm6
movsd  -0x8(%r9,%rsi,8),%xmm3
movsd  (%r9,%rsi,8),%xmm8
movapd %xmm6,%xmm2
addsd  %xmm3,%xmm2
mov    %r11,%rax
sar    $0x1d,%rax
or     $0x8,%rax
movsd  (%rbx,%rax,1),%xmm0
movapd %xmm0,%xmm5
addsd  %xmm8,%xmm5
subsd  %xmm3,%xmm6
subsd  %xmm8,%xmm0
movsd  -0x8(%rdi,%rsi,8),%xmm1
movsd  (%rdi,%rsi,8),%xmm15
movsd  -0x8(%rcx,%rsi,8),%xmm3
movapd %xmm1,%xmm9
addsd  %xmm3,%xmm9
subsd  %xmm3,%xmm1
movsd  (%rcx,%rsi,8),%xmm3
movapd %xmm15,%xmm8
addsd  %xmm3,%xmm8
subsd  %xmm3,%xmm15
movapd %xmm2,%xmm3
addsd  %xmm9,%xmm3
subsd  %xmm9,%xmm2
movapd %xmm6,%xmm9
subsd  %xmm15,%xmm9
addsd  %xmm6,%xmm15
movapd %xmm0,%xmm6
addsd  %xmm1,%xmm6
subsd  %xmm1,%xmm0
movapd %xmm9,%xmm1
subsd  %xmm6,%xmm1
addsd  %xmm9,%xmm6
movapd %xmm15,%xmm9
subsd  %xmm0,%xmm9
addsd  %xmm15,%xmm0
movsd  0x50(%rsp),%xmm15
mulsd  %xmm15,%xmm1
movapd %xmm7,%xmm4
addsd  %xmm1,%xmm4
movsd  %xmm4,(%r14,%rsi,8)
mulsd  %xmm15,%xmm6
movapd %xmm12,%xmm4
subsd  %xmm6,%xmm4
movsd  %xmm4,(%rbx,%rdx,1)
subsd  %xmm1,%xmm7
movsd  %xmm7,-0x8(%r9,%rsi,8)
addsd  %xmm12,%xmm6
mulsd  %xmm15,%xmm9
mulsd  %xmm15,%xmm0
movsd  %xmm6,(%r9,%rsi,8)
movapd %xmm13,%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%r15,%rsi,8)
movapd %xmm10,%xmm1
subsd  %xmm9,%xmm1
movsd  %xmm1,(%rbx,%rbp,1)
addsd  %xmm13,%xmm0
movsd  %xmm0,-0x8(%rcx,%rsi,8)
addsd  %xmm10,%xmm9
movapd %xmm5,%xmm0
addsd  %xmm8,%xmm0
movsd  %xmm9,(%rcx,%rsi,8)
movapd %xmm14,%xmm1
addsd  %xmm3,%xmm1
movsd  %xmm1,(%rbx,%rsi,8)
movapd 0x60(%rsp),%xmm4
movapd %xmm4,%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rbx,%rsi,8)
subsd  %xmm3,%xmm14
movsd  %xmm14,0x0(%r13,%rsi,8)
addsd  %xmm4,%xmm0
movsd  %xmm0,(%rbx,%rax,1)
subsd  %xmm8,%xmm5
movsd  0x58(%rsp),%xmm1
movapd %xmm1,%xmm0
subsd  %xmm5,%xmm0
mov    0x8(%rsp),%rax
movsd  %xmm0,(%rax,%rsi,8)
movapd %xmm11,%xmm0
subsd  %xmm2,%xmm0
movsd  %xmm0,(%rbx,%r12,8)
addsd  %xmm1,%xmm5
movsd  %xmm5,-0x8(%rdi,%rsi,8)
addsd  %xmm11,%xmm2
movsd  %xmm2,(%rdi,%rsi,8)
add    $0x2,%rsi
movabs $0x200000000,%rax
add    %rax,%r8
add    %rax,%r10
add    %rax,%r11
cmp    0x48(%rsp),%rsi
jb     1260 <cftbsub+0x140>
jmp    16fd <cftbsub+0x5dd>
jne    1669 <cftbsub+0x549>
test   %r15d,%r15d
jle    16fd <cftbsub+0x5dd>
movabs $0x200000000,%r8
mov    %r15d,%r9d
add    %r15d,%r15d
or     $0x1,%r15
mov    %r9,%rsi
shl    $0x20,%rsi
lea    (%rbx,%r9,8),%r10
mov    %r9,%r14
shl    $0x4,%r14
add    %rbx,%r14
lea    (%rsi,%rsi,2),%rax
lea    (%r9,%r9,2),%rcx
lea    (%rbx,%rcx,8),%r11
xor    %ecx,%ecx
movapd 0xa9d(%rip),%xmm8
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%r10,%rcx,8),%xmm5
movsd  (%rbx,%rcx,8),%xmm10
movsd  0x8(%rbx,%rcx,8),%xmm6
movapd %xmm10,%xmm4
addsd  %xmm5,%xmm4
movapd %xmm6,%xmm3
xorpd  %xmm8,%xmm3
mov    %rsi,%rdx
sar    $0x1d,%rdx
or     $0x8,%rdx
movsd  (%rbx,%rdx,1),%xmm2
subsd  %xmm2,%xmm3
subsd  %xmm5,%xmm10
subsd  %xmm6,%xmm2
movsd  (%r14,%rcx,8),%xmm5
movsd  (%r11,%rcx,8),%xmm9
movapd %xmm5,%xmm7
addsd  %xmm9,%xmm7
lea    (%r15,%rcx,1),%ebp
movslq %ebp,%rdi
movsd  (%rbx,%rdi,8),%xmm1
mov    %rax,%rbp
sar    $0x1d,%rbp
or     $0x8,%rbp
movsd  (%rbx,%rbp,1),%xmm6
movapd %xmm1,%xmm0
addsd  %xmm6,%xmm0
subsd  %xmm9,%xmm5
subsd  %xmm6,%xmm1
movapd %xmm4,%xmm6
addsd  %xmm7,%xmm6
movsd  %xmm6,(%rbx,%rcx,8)
movapd %xmm3,%xmm6
subsd  %xmm0,%xmm6
movsd  %xmm6,0x8(%rbx,%rcx,8)
subsd  %xmm7,%xmm4
movsd  %xmm4,(%r14,%rcx,8)
addsd  %xmm3,%xmm0
movsd  %xmm0,(%rbx,%rdi,8)
movapd %xmm10,%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%r10,%rcx,8)
movapd %xmm2,%xmm0
subsd  %xmm5,%xmm0
movsd  %xmm0,(%rbx,%rdx,1)
addsd  %xmm10,%xmm1
movsd  %xmm1,(%r11,%rcx,8)
addsd  %xmm2,%xmm5
movsd  %xmm5,(%rbx,%rbp,1)
add    $0x2,%rcx
add    %r8,%rsi
add    %r8,%rax
cmp    %r9,%rcx
jb     1570 <cftbsub+0x450>
jmp    16fd <cftbsub+0x5dd>
test   %r15d,%r15d
jle    16fd <cftbsub+0x5dd>
movabs $0x200000000,%rax
mov    %r15d,%ebp
mov    %rbp,%rcx
shl    $0x20,%rcx
lea    (%rbx,%rbp,8),%rsi
xor    %edi,%edi
movapd 0x96c(%rip),%xmm0
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rsi,%rdi,8),%xmm1
movsd  (%rbx,%rdi,8),%xmm2
movsd  0x8(%rbx,%rdi,8),%xmm3
movapd %xmm2,%xmm4
subsd  %xmm1,%xmm4
mov    %rcx,%rdx
sar    $0x1d,%rdx
or     $0x8,%rdx
movsd  (%rbx,%rdx,1),%xmm5
movapd %xmm5,%xmm6
subsd  %xmm3,%xmm6
addsd  %xmm1,%xmm2
movsd  %xmm2,(%rbx,%rdi,8)
xorpd  %xmm0,%xmm3
subsd  %xmm5,%xmm3
movsd  %xmm3,0x8(%rbx,%rdi,8)
movsd  %xmm4,(%rsi,%rdi,8)
movsd  %xmm6,(%rbx,%rdx,1)
add    $0x2,%rdi
add    %rax,%rcx
cmp    %rbp,%rdi
jb     16a0 <cftbsub+0x580>
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @cftbsub(i32 %0, double* %1, double* %2) {
%4 = alloca i32
%5 = alloca double*
%6 = alloca double*
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca double
%17 = alloca double
%18 = alloca double
%19 = alloca double
%20 = alloca double
%21 = alloca double
%22 = alloca double
%23 = alloca double
%24 = alloca double
%25 = alloca double
%26 = alloca double
%27 = alloca double
%28 = alloca double
%29 = alloca double
%30 = alloca double
%31 = alloca double
%32 = alloca double
%33 = alloca double
%34 = alloca double
%35 = alloca double
%36 = alloca double
%37 = alloca double
%38 = alloca double
%39 = alloca double
%40 = alloca double
store i32 %0, i32* %4
store double* %1, double** %5
store double* %2, double** %6
store i32 2, i32* %15
%41 = load i32, i32* %4
%42 = icmp sgt i32 %41, 16
br i1 %42, label %43, label %60
43: ; preds = %3
%44 = load i32, i32* %4
%45 = load double*, double** %5
%46 = load double*, double** %6
store i32 16, i32* %15
br label %47
47: ; preds = %52, %43
%48 = load i32, i32* %15
%49 = shl i32 %48, 3
%50 = load i32, i32* %4
%51 = icmp slt i32 %49, %50
br i1 %51, label %52, label %59
52: ; preds = %47
%53 = load i32, i32* %4
%54 = load i32, i32* %15
%55 = load double*, double** %5
%56 = load double*, double** %6
%57 = load i32, i32* %15
%58 = shl i32 %57, 3
store i32 %58, i32* %15
br label %47
59: ; preds = %47
br label %60
60: ; preds = %59, %3
%61 = load i32, i32* %15
%62 = shl i32 %61, 2
%63 = load i32, i32* %4
%64 = icmp slt i32 %62, %63
br i1 %64, label %65, label %481
65: ; preds = %60
%66 = load double*, double** %6
%67 = getelementptr inbounds double, double* %66, i64 2
%68 = load double, double* %67
store double %68, double* %16
store i32 0, i32* %7
br label %69
69: ; preds = %477, %65
%70 = load i32, i32* %7
%71 = load i32, i32* %15
%72 = icmp slt i32 %70, %71
br i1 %72, label %73, label %480
73: ; preds = %69
%74 = load i32, i32* %7
%75 = load i32, i32* %15
%76 = add nsw i32 %74, %75
store i32 %76, i32* %8
%77 = load i32, i32* %8
%78 = load i32, i32* %15
%79 = add nsw i32 %77, %78
store i32 %79, i32* %9
%80 = load i32, i32* %9
%81 = load i32, i32* %15
%82 = add nsw i32 %80, %81
store i32 %82, i32* %10
%83 = load i32, i32* %10
%84 = load i32, i32* %15
%85 = add nsw i32 %83, %84
store i32 %85, i32* %11
%86 = load i32, i32* %11
%87 = load i32, i32* %15
%88 = add nsw i32 %86, %87
store i32 %88, i32* %12
%89 = load i32, i32* %12
%90 = load i32, i32* %15
%91 = add nsw i32 %89, %90
store i32 %91, i32* %13
%92 = load i32, i32* %13
%93 = load i32, i32* %15
%94 = add nsw i32 %92, %93
store i32 %94, i32* %14
%95 = load double*, double** %5
%96 = load i32, i32* %7
%97 = sext i32 %96 to i64
%98 = getelementptr inbounds double, double* %95, i64 %97
%99 = load double, double* %98
%100 = load double*, double** %5
%101 = load i32, i32* %8
%102 = sext i32 %101 to i64
%103 = getelementptr inbounds double, double* %100, i64 %102
%104 = load double, double* %103
%105 = fadd double %99, %104
store double %105, double* %17
%106 = load double*, double** %5
%107 = load i32, i32* %7
%108 = add nsw i32 %107, 1
%109 = sext i32 %108 to i64
%110 = getelementptr inbounds double, double* %106, i64 %109
%111 = load double, double* %110
%112 = fneg double %111
%113 = load double*, double** %5
%114 = load i32, i32* %8
%115 = add nsw i32 %114, 1
%116 = sext i32 %115 to i64
%117 = getelementptr inbounds double, double* %113, i64 %116
%118 = load double, double* %117
%119 = fsub double %112, %118
store double %119, double* %18
%120 = load double*, double** %5
%121 = load i32, i32* %7
%122 = sext i32 %121 to i64
%123 = getelementptr inbounds double, double* %120, i64 %122
%124 = load double, double* %123
%125 = load double*, double** %5
%126 = load i32, i32* %8
%127 = sext i32 %126 to i64
%128 = getelementptr inbounds double, double* %125, i64 %127
%129 = load double, double* %128
%130 = fsub double %124, %129
store double %130, double* %19
%131 = load double*, double** %5
%132 = load i32, i32* %7
%133 = add nsw i32 %132, 1
%134 = sext i32 %133 to i64
%135 = getelementptr inbounds double, double* %131, i64 %134
%136 = load double, double* %135
%137 = fneg double %136
%138 = load double*, double** %5
%139 = load i32, i32* %8
%140 = add nsw i32 %139, 1
%141 = sext i32 %140 to i64
%142 = getelementptr inbounds double, double* %138, i64 %141
%143 = load double, double* %142
%144 = fadd double %137, %143
store double %144, double* %20
%145 = load double*, double** %5
%146 = load i32, i32* %9
%147 = sext i32 %146 to i64
%148 = getelementptr inbounds double, double* %145, i64 %147
%149 = load double, double* %148
%150 = load double*, double** %5
%151 = load i32, i32* %10
%152 = sext i32 %151 to i64
%153 = getelementptr inbounds double, double* %150, i64 %152
%154 = load double, double* %153
%155 = fadd double %149, %154
store double %155, double* %21
%156 = load double*, double** %5
%157 = load i32, i32* %9
%158 = add nsw i32 %157, 1
%159 = sext i32 %158 to i64
%160 = getelementptr inbounds double, double* %156, i64 %159
%161 = load double, double* %160
%162 = load double*, double** %5
%163 = load i32, i32* %10
%164 = add nsw i32 %163, 1
%165 = sext i32 %164 to i64
%166 = getelementptr inbounds double, double* %162, i64 %165
%167 = load double, double* %166
%168 = fadd double %161, %167
store double %168, double* %22
%169 = load double*, double** %5
%170 = load i32, i32* %9
%171 = sext i32 %170 to i64
%172 = getelementptr inbounds double, double* %169, i64 %171
%173 = load double, double* %172
%174 = load double*, double** %5
%175 = load i32, i32* %10
%176 = sext i32 %175 to i64
%177 = getelementptr inbounds double, double* %174, i64 %176
%178 = load double, double* %177
%179 = fsub double %173, %178
store double %179, double* %23
%180 = load double*, double** %5
%181 = load i32, i32* %9
%182 = add nsw i32 %181, 1
%183 = sext i32 %182 to i64
%184 = getelementptr inbounds double, double* %180, i64 %183
%185 = load double, double* %184
%186 = load double*, double** %5
%187 = load i32, i32* %10
%188 = add nsw i32 %187, 1
%189 = sext i32 %188 to i64
%190 = getelementptr inbounds double, double* %186, i64 %189
%191 = load double, double* %190
%192 = fsub double %185, %191
store double %192, double* %24
%193 = load double, double* %17
%194 = load double, double* %21
%195 = fadd double %193, %194
store double %195, double* %25
%196 = load double, double* %18
%197 = load double, double* %22
%198 = fsub double %196, %197
store double %198, double* %26
%199 = load double, double* %17
%200 = load double, double* %21
%201 = fsub double %199, %200
store double %201, double* %29
%202 = load double, double* %18
%203 = load double, double* %22
%204 = fadd double %202, %203
store double %204, double* %30
%205 = load double, double* %19
%206 = load double, double* %24
%207 = fsub double %205, %206
store double %207, double* %27
%208 = load double, double* %20
%209 = load double, double* %23
%210 = fsub double %208, %209
store double %210, double* %28
%211 = load double, double* %19
%212 = load double, double* %24
%213 = fadd double %211, %212
store double %213, double* %31
%214 = load double, double* %20
%215 = load double, double* %23
%216 = fadd double %214, %215
store double %216, double* %32
%217 = load double*, double** %5
%218 = load i32, i32* %11
%219 = sext i32 %218 to i64
%220 = getelementptr inbounds double, double* %217, i64 %219
%221 = load double, double* %220
%222 = load double*, double** %5
%223 = load i32, i32* %12
%224 = sext i32 %223 to i64
%225 = getelementptr inbounds double, double* %222, i64 %224
%226 = load double, double* %225
%227 = fadd double %221, %226
store double %227, double* %17
%228 = load double*, double** %5
%229 = load i32, i32* %11
%230 = add nsw i32 %229, 1
%231 = sext i32 %230 to i64
%232 = getelementptr inbounds double, double* %228, i64 %231
%233 = load double, double* %232
%234 = load double*, double** %5
%235 = load i32, i32* %12
%236 = add nsw i32 %235, 1
%237 = sext i32 %236 to i64
%238 = getelementptr inbounds double, double* %234, i64 %237
%239 = load double, double* %238
%240 = fadd double %233, %239
store double %240, double* %18
%241 = load double*, double** %5
%242 = load i32, i32* %11
%243 = sext i32 %242 to i64
%244 = getelementptr inbounds double, double* %241, i64 %243
%245 = load double, double* %244
%246 = load double*, double** %5
%247 = load i32, i32* %12
%248 = sext i32 %247 to i64
%249 = getelementptr inbounds double, double* %246, i64 %248
%250 = load double, double* %249
%251 = fsub double %245, %250
store double %251, double* %19
%252 = load double*, double** %5
%253 = load i32, i32* %11
%254 = add nsw i32 %253, 1
%255 = sext i32 %254 to i64
%256 = getelementptr inbounds double, double* %252, i64 %255
%257 = load double, double* %256
%258 = load double*, double** %5
%259 = load i32, i32* %12
%260 = add nsw i32 %259, 1
%261 = sext i32 %260 to i64
%262 = getelementptr inbounds double, double* %258, i64 %261
%263 = load double, double* %262
%264 = fsub double %257, %263
store double %264, double* %20
%265 = load double*, double** %5
%266 = load i32, i32* %13
%267 = sext i32 %266 to i64
%268 = getelementptr inbounds double, double* %265, i64 %267
%269 = load double, double* %268
%270 = load double*, double** %5
%271 = load i32, i32* %14
%272 = sext i32 %271 to i64
%273 = getelementptr inbounds double, double* %270, i64 %272
%274 = load double, double* %273
%275 = fadd double %269, %274
store double %275, double* %21
%276 = load double*, double** %5
%277 = load i32, i32* %13
%278 = add nsw i32 %277, 1
%279 = sext i32 %278 to i64
%280 = getelementptr inbounds double, double* %276, i64 %279
%281 = load double, double* %280
%282 = load double*, double** %5
%283 = load i32, i32* %14
%284 = add nsw i32 %283, 1
%285 = sext i32 %284 to i64
%286 = getelementptr inbounds double, double* %282, i64 %285
%287 = load double, double* %286
%288 = fadd double %281, %287
store double %288, double* %22
%289 = load double*, double** %5
%290 = load i32, i32* %13
%291 = sext i32 %290 to i64
%292 = getelementptr inbounds double, double* %289, i64 %291
%293 = load double, double* %292
%294 = load double*, double** %5
%295 = load i32, i32* %14
%296 = sext i32 %295 to i64
%297 = getelementptr inbounds double, double* %294, i64 %296
%298 = load double, double* %297
%299 = fsub double %293, %298
store double %299, double* %23
%300 = load double*, double** %5
%301 = load i32, i32* %13
%302 = add nsw i32 %301, 1
%303 = sext i32 %302 to i64
%304 = getelementptr inbounds double, double* %300, i64 %303
%305 = load double, double* %304
%306 = load double*, double** %5
%307 = load i32, i32* %14
%308 = add nsw i32 %307, 1
%309 = sext i32 %308 to i64
%310 = getelementptr inbounds double, double* %306, i64 %309
%311 = load double, double* %310
%312 = fsub double %305, %311
store double %312, double* %24
%313 = load double, double* %17
%314 = load double, double* %21
%315 = fadd double %313, %314
store double %315, double* %33
%316 = load double, double* %18
%317 = load double, double* %22
%318 = fadd double %316, %317
store double %318, double* %34
%319 = load double, double* %17
%320 = load double, double* %21
%321 = fsub double %319, %320
store double %321, double* %37
%322 = load double, double* %18
%323 = load double, double* %22
%324 = fsub double %322, %323
store double %324, double* %38
%325 = load double, double* %19
%326 = load double, double* %24
%327 = fsub double %325, %326
store double %327, double* %17
%328 = load double, double* %20
%329 = load double, double* %23
%330 = fadd double %328, %329
store double %330, double* %18
%331 = load double, double* %19
%332 = load double, double* %24
%333 = fadd double %331, %332
store double %333, double* %21
%334 = load double, double* %20
%335 = load double, double* %23
%336 = fsub double %334, %335
store double %336, double* %22
%337 = load double, double* %16
%338 = load double, double* %17
%339 = load double, double* %18
%340 = fsub double %338, %339
%341 = fmul double %337, %340
store double %341, double* %35
%342 = load double, double* %16
%343 = load double, double* %17
%344 = load double, double* %18
%345 = fadd double %343, %344
%346 = fmul double %342, %345
store double %346, double* %36
%347 = load double, double* %16
%348 = load double, double* %21
%349 = load double, double* %22
%350 = fsub double %348, %349
%351 = fmul double %347, %350
store double %351, double* %39
%352 = load double, double* %16
%353 = load double, double* %21
%354 = load double, double* %22
%355 = fadd double %353, %354
%356 = fmul double %352, %355
store double %356, double* %40
%357 = load double, double* %27
%358 = load double, double* %35
%359 = fadd double %357, %358
%360 = load double*, double** %5
%361 = load i32, i32* %8
%362 = sext i32 %361 to i64
%363 = getelementptr inbounds double, double* %360, i64 %362
store double %359, double* %363
%364 = load double, double* %28
%365 = load double, double* %36
%366 = fsub double %364, %365
%367 = load double*, double** %5
%368 = load i32, i32* %8
%369 = add nsw i32 %368, 1
%370 = sext i32 %369 to i64
%371 = getelementptr inbounds double, double* %367, i64 %370
store double %366, double* %371
%372 = load double, double* %27
%373 = load double, double* %35
%374 = fsub double %372, %373
%375 = load double*, double** %5
%376 = load i32, i32* %12
%377 = sext i32 %376 to i64
%378 = getelementptr inbounds double, double* %375, i64 %377
store double %374, double* %378
%379 = load double, double* %28
%380 = load double, double* %36
%381 = fadd double %379, %380
%382 = load double*, double** %5
%383 = load i32, i32* %12
%384 = add nsw i32 %383, 1
%385 = sext i32 %384 to i64
%386 = getelementptr inbounds double, double* %382, i64 %385
store double %381, double* %386
%387 = load double, double* %31
%388 = load double, double* %40
%389 = fsub double %387, %388
%390 = load double*, double** %5
%391 = load i32, i32* %10
%392 = sext i32 %391 to i64
%393 = getelementptr inbounds double, double* %390, i64 %392
store double %389, double* %393
%394 = load double, double* %32
%395 = load double, double* %39
%396 = fsub double %394, %395
%397 = load double*, double** %5
%398 = load i32, i32* %10
%399 = add nsw i32 %398, 1
%400 = sext i32 %399 to i64
%401 = getelementptr inbounds double, double* %397, i64 %400
store double %396, double* %401
%402 = load double, double* %31
%403 = load double, double* %40
%404 = fadd double %402, %403
%405 = load double*, double** %5
%406 = load i32, i32* %14
%407 = sext i32 %406 to i64
%408 = getelementptr inbounds double, double* %405, i64 %407
store double %404, double* %408
%409 = load double, double* %32
%410 = load double, double* %39
%411 = fadd double %409, %410
%412 = load double*, double** %5
%413 = load i32, i32* %14
%414 = add nsw i32 %413, 1
%415 = sext i32 %414 to i64
%416 = getelementptr inbounds double, double* %412, i64 %415
store double %411, double* %416
%417 = load double, double* %25
%418 = load double, double* %33
%419 = fadd double %417, %418
%420 = load double*, double** %5
%421 = load i32, i32* %7
%422 = sext i32 %421 to i64
%423 = getelementptr inbounds double, double* %420, i64 %422
store double %419, double* %423
%424 = load double, double* %26
%425 = load double, double* %34
%426 = fsub double %424, %425
%427 = load double*, double** %5
%428 = load i32, i32* %7
%429 = add nsw i32 %428, 1
%430 = sext i32 %429 to i64
%431 = getelementptr inbounds double, double* %427, i64 %430
store double %426, double* %431
%432 = load double, double* %25
%433 = load double, double* %33
%434 = fsub double %432, %433
%435 = load double*, double** %5
%436 = load i32, i32* %11
%437 = sext i32 %436 to i64
%438 = getelementptr inbounds double, double* %435, i64 %437
store double %434, double* %438
%439 = load double, double* %26
%440 = load double, double* %34
%441 = fadd double %439, %440
%442 = load double*, double** %5
%443 = load i32, i32* %11
%444 = add nsw i32 %443, 1
%445 = sext i32 %444 to i64
%446 = getelementptr inbounds double, double* %442, i64 %445
store double %441, double* %446
%447 = load double, double* %29
%448 = load double, double* %38
%449 = fsub double %447, %448
%450 = load double*, double** %5
%451 = load i32, i32* %9
%452 = sext i32 %451 to i64
%453 = getelementptr inbounds double, double* %450, i64 %452
store double %449, double* %453
%454 = load double, double* %30
%455 = load double, double* %37
%456 = fsub double %454, %455
%457 = load double*, double** %5
%458 = load i32, i32* %9
%459 = add nsw i32 %458, 1
%460 = sext i32 %459 to i64
%461 = getelementptr inbounds double, double* %457, i64 %460
store double %456, double* %461
%462 = load double, double* %29
%463 = load double, double* %38
%464 = fadd double %462, %463
%465 = load double*, double** %5
%466 = load i32, i32* %13
%467 = sext i32 %466 to i64
%468 = getelementptr inbounds double, double* %465, i64 %467
store double %464, double* %468
%469 = load double, double* %30
%470 = load double, double* %37
%471 = fadd double %469, %470
%472 = load double*, double** %5
%473 = load i32, i32* %13
%474 = add nsw i32 %473, 1
%475 = sext i32 %474 to i64
%476 = getelementptr inbounds double, double* %472, i64 %475
store double %471, double* %476
br label %477
477: ; preds = %73
%478 = load i32, i32* %7
%479 = add nsw i32 %478, 2
store i32 %479, i32* %7
br label %69
480: ; preds = %69
br label %743
481: ; preds = %60
%482 = load i32, i32* %15
%483 = shl i32 %482, 2
%484 = load i32, i32* %4
%485 = icmp eq i32 %483, %484
br i1 %485, label %486, label %663
486: ; preds = %481
store i32 0, i32* %7
br label %487
487: ; preds = %659, %486
%488 = load i32, i32* %7
%489 = load i32, i32* %15
%490 = icmp slt i32 %488, %489
br i1 %490, label %491, label %662
491: ; preds = %487
%492 = load i32, i32* %7
%493 = load i32, i32* %15
%494 = add nsw i32 %492, %493
store i32 %494, i32* %8
%495 = load i32, i32* %8
%496 = load i32, i32* %15
%497 = add nsw i32 %495, %496
store i32 %497, i32* %9
%498 = load i32, i32* %9
%499 = load i32, i32* %15
%500 = add nsw i32 %498, %499
store i32 %500, i32* %10
%501 = load double*, double** %5
%502 = load i32, i32* %7
%503 = sext i32 %502 to i64
%504 = getelementptr inbounds double, double* %501, i64 %503
%505 = load double, double* %504
%506 = load double*, double** %5
%507 = load i32, i32* %8
%508 = sext i32 %507 to i64
%509 = getelementptr inbounds double, double* %506, i64 %508
%510 = load double, double* %509
%511 = fadd double %505, %510
store double %511, double* %17
%512 = load double*, double** %5
%513 = load i32, i32* %7
%514 = add nsw i32 %513, 1
%515 = sext i32 %514 to i64
%516 = getelementptr inbounds double, double* %512, i64 %515
%517 = load double, double* %516
%518 = fneg double %517
%519 = load double*, double** %5
%520 = load i32, i32* %8
%521 = add nsw i32 %520, 1
%522 = sext i32 %521 to i64
%523 = getelementptr inbounds double, double* %519, i64 %522
%524 = load double, double* %523
%525 = fsub double %518, %524
store double %525, double* %18
%526 = load double*, double** %5
%527 = load i32, i32* %7
%528 = sext i32 %527 to i64
%529 = getelementptr inbounds double, double* %526, i64 %528
%530 = load double, double* %529
%531 = load double*, double** %5
%532 = load i32, i32* %8
%533 = sext i32 %532 to i64
%534 = getelementptr inbounds double, double* %531, i64 %533
%535 = load double, double* %534
%536 = fsub double %530, %535
store double %536, double* %19
%537 = load double*, double** %5
%538 = load i32, i32* %7
%539 = add nsw i32 %538, 1
%540 = sext i32 %539 to i64
%541 = getelementptr inbounds double, double* %537, i64 %540
%542 = load double, double* %541
%543 = fneg double %542
%544 = load double*, double** %5
%545 = load i32, i32* %8
%546 = add nsw i32 %545, 1
%547 = sext i32 %546 to i64
%548 = getelementptr inbounds double, double* %544, i64 %547
%549 = load double, double* %548
%550 = fadd double %543, %549
store double %550, double* %20
%551 = load double*, double** %5
%552 = load i32, i32* %9
%553 = sext i32 %552 to i64
%554 = getelementptr inbounds double, double* %551, i64 %553
%555 = load double, double* %554
%556 = load double*, double** %5
%557 = load i32, i32* %10
%558 = sext i32 %557 to i64
%559 = getelementptr inbounds double, double* %556, i64 %558
%560 = load double, double* %559
%561 = fadd double %555, %560
store double %561, double* %21
%562 = load double*, double** %5
%563 = load i32, i32* %9
%564 = add nsw i32 %563, 1
%565 = sext i32 %564 to i64
%566 = getelementptr inbounds double, double* %562, i64 %565
%567 = load double, double* %566
%568 = load double*, double** %5
%569 = load i32, i32* %10
%570 = add nsw i32 %569, 1
%571 = sext i32 %570 to i64
%572 = getelementptr inbounds double, double* %568, i64 %571
%573 = load double, double* %572
%574 = fadd double %567, %573
store double %574, double* %22
%575 = load double*, double** %5
%576 = load i32, i32* %9
%577 = sext i32 %576 to i64
%578 = getelementptr inbounds double, double* %575, i64 %577
%579 = load double, double* %578
%580 = load double*, double** %5
%581 = load i32, i32* %10
%582 = sext i32 %581 to i64
%583 = getelementptr inbounds double, double* %580, i64 %582
%584 = load double, double* %583
%585 = fsub double %579, %584
store double %585, double* %23
%586 = load double*, double** %5
%587 = load i32, i32* %9
%588 = add nsw i32 %587, 1
%589 = sext i32 %588 to i64
%590 = getelementptr inbounds double, double* %586, i64 %589
%591 = load double, double* %590
%592 = load double*, double** %5
%593 = load i32, i32* %10
%594 = add nsw i32 %593, 1
%595 = sext i32 %594 to i64
%596 = getelementptr inbounds double, double* %592, i64 %595
%597 = load double, double* %596
%598 = fsub double %591, %597
store double %598, double* %24
%599 = load double, double* %17
%600 = load double, double* %21
%601 = fadd double %599, %600
%602 = load double*, double** %5
%603 = load i32, i32* %7
%604 = sext i32 %603 to i64
%605 = getelementptr inbounds double, double* %602, i64 %604
store double %601, double* %605
%606 = load double, double* %18
%607 = load double, double* %22
%608 = fsub double %606, %607
%609 = load double*, double** %5
%610 = load i32, i32* %7
%611 = add nsw i32 %610, 1
%612 = sext i32 %611 to i64
%613 = getelementptr inbounds double, double* %609, i64 %612
store double %608, double* %613
%614 = load double, double* %17
%615 = load double, double* %21
%616 = fsub double %614, %615
%617 = load double*, double** %5
%618 = load i32, i32* %9
%619 = sext i32 %618 to i64
%620 = getelementptr inbounds double, double* %617, i64 %619
store double %616, double* %620
%621 = load double, double* %18
%622 = load double, double* %22
%623 = fadd double %621, %622
%624 = load double*, double** %5
%625 = load i32, i32* %9
%626 = add nsw i32 %625, 1
%627 = sext i32 %626 to i64
%628 = getelementptr inbounds double, double* %624, i64 %627
store double %623, double* %628
%629 = load double, double* %19
%630 = load double, double* %24
%631 = fsub double %629, %630
%632 = load double*, double** %5
%633 = load i32, i32* %8
%634 = sext i32 %633 to i64
%635 = getelementptr inbounds double, double* %632, i64 %634
store double %631, double* %635
%636 = load double, double* %20
%637 = load double, double* %23
%638 = fsub double %636, %637
%639 = load double*, double** %5
%640 = load i32, i32* %8
%641 = add nsw i32 %640, 1
%642 = sext i32 %641 to i64
%643 = getelementptr inbounds double, double* %639, i64 %642
store double %638, double* %643
%644 = load double, double* %19
%645 = load double, double* %24
%646 = fadd double %644, %645
%647 = load double*, double** %5
%648 = load i32, i32* %10
%649 = sext i32 %648 to i64
%650 = getelementptr inbounds double, double* %647, i64 %649
store double %646, double* %650
%651 = load double, double* %20
%652 = load double, double* %23
%653 = fadd double %651, %652
%654 = load double*, double** %5
%655 = load i32, i32* %10
%656 = add nsw i32 %655, 1
%657 = sext i32 %656 to i64
%658 = getelementptr inbounds double, double* %654, i64 %657
store double %653, double* %658
br label %659
659: ; preds = %491
%660 = load i32, i32* %7
%661 = add nsw i32 %660, 2
store i32 %661, i32* %7
br label %487
662: ; preds = %487
br label %742
663: ; preds = %481
store i32 0, i32* %7
br label %664
664: ; preds = %738, %663
%665 = load i32, i32* %7
%666 = load i32, i32* %15
%667 = icmp slt i32 %665, %666
br i1 %667, label %668, label %741
668: ; preds = %664
%669 = load i32, i32* %7
%670 = load i32, i32* %15
%671 = add nsw i32 %669, %670
store i32 %671, i32* %8
%672 = load double*, double** %5
%673 = load i32, i32* %7
%674 = sext i32 %673 to i64
%675 = getelementptr inbounds double, double* %672, i64 %674
%676 = load double, double* %675
%677 = load double*, double** %5
%678 = load i32, i32* %8
%679 = sext i32 %678 to i64
%680 = getelementptr inbounds double, double* %677, i64 %679
%681 = load double, double* %680
%682 = fsub double %676, %681
store double %682, double* %17
%683 = load double*, double** %5
%684 = load i32, i32* %7
%685 = add nsw i32 %684, 1
%686 = sext i32 %685 to i64
%687 = getelementptr inbounds double, double* %683, i64 %686
%688 = load double, double* %687
%689 = fneg double %688
%690 = load double*, double** %5
%691 = load i32, i32* %8
%692 = add nsw i32 %691, 1
%693 = sext i32 %692 to i64
%694 = getelementptr inbounds double, double* %690, i64 %693
%695 = load double, double* %694
%696 = fadd double %689, %695
store double %696, double* %18
%697 = load double*, double** %5
%698 = load i32, i32* %8
%699 = sext i32 %698 to i64
%700 = getelementptr inbounds double, double* %697, i64 %699
%701 = load double, double* %700
%702 = load double*, double** %5
%703 = load i32, i32* %7
%704 = sext i32 %703 to i64
%705 = getelementptr inbounds double, double* %702, i64 %704
%706 = load double, double* %705
%707 = fadd double %706, %701
store double %707, double* %705
%708 = load double*, double** %5
%709 = load i32, i32* %7
%710 = add nsw i32 %709, 1
%711 = sext i32 %710 to i64
%712 = getelementptr inbounds double, double* %708, i64 %711
%713 = load double, double* %712
%714 = fneg double %713
%715 = load double*, double** %5
%716 = load i32, i32* %8
%717 = add nsw i32 %716, 1
%718 = sext i32 %717 to i64
%719 = getelementptr inbounds double, double* %715, i64 %718
%720 = load double, double* %719
%721 = fsub double %714, %720
%722 = load double*, double** %5
%723 = load i32, i32* %7
%724 = add nsw i32 %723, 1
%725 = sext i32 %724 to i64
%726 = getelementptr inbounds double, double* %722, i64 %725
store double %721, double* %726
%727 = load double, double* %17
%728 = load double*, double** %5
%729 = load i32, i32* %8
%730 = sext i32 %729 to i64
%731 = getelementptr inbounds double, double* %728, i64 %730
store double %727, double* %731
%732 = load double, double* %18
%733 = load double*, double** %5
%734 = load i32, i32* %8
%735 = add nsw i32 %734, 1
%736 = sext i32 %735 to i64
%737 = getelementptr inbounds double, double* %733, i64 %736
store double %732, double* %737
br label %738
738: ; preds = %668
%739 = load i32, i32* %7
%740 = add nsw i32 %739, 2
store i32 %740, i32* %7
br label %664
741: ; preds = %664
br label %742
742: ; preds = %741, %662
br label %743
743: ; preds = %742, %480
ret void
}

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