void dfct(int n, double *a, double *t, int *ip, double *w) {
  void makewt(int nw, int *ip, double *w);
  void makect(int nc, int *ip, double *c);
  void cftfsub(int n, double *a, int *ip, int nw, double *w);
  void rftfsub(int n, double *a, int nc, double *c);
  void dctsub(int n, double *a, int nc, double *c);
  int j, k, l, m, mh, nw, nc;
  double xr, xi, yr, yi;
  nw = ip[0];
  if (n > (nw << 3)) {
    nw = n >> 3;
    makewt(nw, ip, w);
  }
  nc = ip[1];
  if (n > (nc << 1)) {
    nc = n >> 1;
    makect(nc, ip, w + nw);
  }
  m = n >> 1;
  yi = a[m];
  xi = a[0] + a[n];
  a[0] -= a[n];
  t[0] = xi - yi;
  t[m] = xi + yi;
  if (n > 2) {
    mh = m >> 1;
    for (j = 1; j < mh; j++) {
      k = m - j;
      xr = a[j] - a[n - j];
      xi = a[j] + a[n - j];
      yr = a[k] - a[n - k];
      yi = a[k] + a[n - k];
      a[j] = xr;
      a[k] = yr;
      t[j] = xi - yi;
      t[k] = xi + yi;
    }
    t[mh] = a[mh] + a[n - mh];
    a[mh] -= a[n - mh];
    dctsub(m, a, nc, w + nw);
    if (m > 4) {
      cftfsub(m, a, ip, nw, w);
      rftfsub(m, a, nc, w + nw);
    } else if (m == 4) {
      cftfsub(m, a, ip, nw, w);
    }
    a[n - 1] = a[0] - a[1];
    a[1] = a[0] + a[1];
    for (j = m - 2; j >= 2; j -= 2) {
      a[2 * j + 1] = a[j] + a[j + 1];
      a[2 * j - 1] = a[j] - a[j + 1];
    }
    l = 2;
    m = mh;
    while (m >= 2) {
      dctsub(m, t, nc, w + nw);
      if (m > 4) {
        cftfsub(m, t, ip, nw, w);
        rftfsub(m, t, nc, w + nw);
      } else if (m == 4) {
        cftfsub(m, t, ip, nw, w);
      }
      a[n - l] = t[0] - t[1];
      a[l] = t[0] + t[1];
      k = 0;
      for (j = 2; j < m; j += 2) {
        k += l << 2;
        a[k - l] = t[j] - t[j + 1];
        a[k + l] = t[j] + t[j + 1];
      }
      l <<= 1;
      mh = m >> 1;
      for (j = 0; j < mh; j++) {
        k = m - j;
        t[j] = t[m + k] - t[m + j];
        t[k] = t[m + k] + t[m + j];
      }
      t[mh] = t[m + mh];
      m = mh;
    }
    a[l] = t[0];
    a[n] = t[2] - t[1];
    a[0] = t[2] + t[1];
  } else {
    a[1] = a[0];
    a[2] = t[0];
    a[0] = t[1];
  }
}


x86 assembly (clang, -O0):
<dfct>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x20(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x40(%rbp)
mov    -0x4(%rbp),%eax
mov    -0x40(%rbp),%ecx
shl    $0x3,%ecx
cmp    %ecx,%eax
jle    119e <dfct+0x4e>
mov    -0x4(%rbp),%eax
sar    $0x3,%eax
mov    %eax,-0x40(%rbp)
mov    -0x40(%rbp),%edi
mov    -0x20(%rbp),%rsi
mov    -0x28(%rbp),%rdx
call   1040 <makewt@plt>
mov    -0x20(%rbp),%rax
mov    0x4(%rax),%eax
mov    %eax,-0x44(%rbp)
mov    -0x4(%rbp),%eax
mov    -0x44(%rbp),%ecx
shl    $0x1,%ecx
cmp    %ecx,%eax
jle    11dd <dfct+0x8d>
mov    -0x4(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%edi
mov    -0x20(%rbp),%rsi
mov    -0x28(%rbp),%rdx
movslq -0x40(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdx
call   1060 <makect@plt>
mov    -0x4(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x38(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x4(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
movsd  -0x58(%rbp),%xmm0
subsd  -0x68(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x58(%rbp),%xmm0
addsd  -0x68(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
cmpl   $0x2,-0x4(%rbp)
jle    17aa <dfct+0x65a>
mov    -0x38(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x3c(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    1377 <dfct+0x227>
mov    -0x38(%rbp),%eax
sub    -0x2c(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    -0x30(%rbp),%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    -0x30(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x50(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x60(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x58(%rbp),%xmm0
subsd  -0x68(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x58(%rbp),%xmm0
addsd  -0x68(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1272 <dfct+0x122>
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x44(%rbp),%edx
mov    -0x28(%rbp),%rcx
movslq -0x40(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rcx
call   1030 <dctsub@plt>
cmpl   $0x4,-0x38(%rbp)
jle    142d <dfct+0x2dd>
mov    -0x38(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x20(%rbp),%rdx
mov    -0x40(%rbp),%ecx
mov    -0x28(%rbp),%r8
call   1070 <cftfsub@plt>
mov    -0x38(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x44(%rbp),%edx
mov    -0x28(%rbp),%rcx
movslq -0x40(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rcx
call   1050 <rftfsub@plt>
jmp    1453 <dfct+0x303>
cmpl   $0x4,-0x38(%rbp)
jne    144e <dfct+0x2fe>
mov    -0x38(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x20(%rbp),%rdx
mov    -0x40(%rbp),%ecx
mov    -0x28(%rbp),%r8
call   1070 <cftfsub@plt>
jmp    1453 <dfct+0x303>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
subsd  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
addsd  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x8(%rax)
mov    -0x38(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x2,-0x2c(%rbp)
jl     1519 <dfct+0x3c9>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
shl    $0x1,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x2c(%rbp)
jmp    1499 <dfct+0x349>
movl   $0x2,-0x34(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x38(%rbp)
cmpl   $0x2,-0x38(%rbp)
jl     1757 <dfct+0x607>
mov    -0x38(%rbp),%edi
mov    -0x18(%rbp),%rsi
mov    -0x44(%rbp),%edx
mov    -0x28(%rbp),%rcx
movslq -0x40(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rcx
call   1030 <dctsub@plt>
cmpl   $0x4,-0x38(%rbp)
jle    1592 <dfct+0x442>
mov    -0x38(%rbp),%edi
mov    -0x18(%rbp),%rsi
mov    -0x20(%rbp),%rdx
mov    -0x40(%rbp),%ecx
mov    -0x28(%rbp),%r8
call   1070 <cftfsub@plt>
mov    -0x38(%rbp),%edi
mov    -0x18(%rbp),%rsi
mov    -0x44(%rbp),%edx
mov    -0x28(%rbp),%rcx
movslq -0x40(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rcx
call   1050 <rftfsub@plt>
jmp    15b8 <dfct+0x468>
cmpl   $0x4,-0x38(%rbp)
jne    15b3 <dfct+0x463>
mov    -0x38(%rbp),%edi
mov    -0x18(%rbp),%rsi
mov    -0x20(%rbp),%rdx
mov    -0x40(%rbp),%ecx
mov    -0x28(%rbp),%r8
call   1070 <cftfsub@plt>
jmp    15b8 <dfct+0x468>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
subsd  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
addsd  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x30(%rbp)
movl   $0x2,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x38(%rbp),%eax
jge    168f <dfct+0x53f>
mov    -0x34(%rbp),%eax
shl    $0x2,%eax
add    -0x30(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x2c(%rbp)
jmp    1607 <dfct+0x4b7>
mov    -0x34(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x38(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x3c(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    172d <dfct+0x5dd>
mov    -0x38(%rbp),%eax
sub    -0x2c(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    16a8 <dfct+0x558>
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x38(%rbp)
jmp    1526 <dfct+0x3d6>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x18(%rbp),%rax
subsd  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x4(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x18(%rbp),%rax
addsd  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,(%rax)
jmp    17dd <dfct+0x68d>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x8(%rax)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x10(%rax)
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,(%rax)
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dfct>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %rcx,%r15
mov    %rdx,%rbx
mov    %rsi,%r14
mov    %edi,%ebp
mov    (%rcx),%r10d
lea    0x0(,%r10,8),%eax
cmp    %edi,%eax
mov    %r8,0x10(%rsp)
jge    119a <dfct+0x4a>
mov    %ebp,%r12d
sar    $0x3,%r12d
mov    %r12d,%edi
mov    %r15,%rsi
mov    %r8,%rdx
call   1040 <makewt@plt>
mov    0x10(%rsp),%r8
mov    %r12,%r10
mov    %r15,0x18(%rsp)
mov    0x4(%r15),%ecx
lea    (%rcx,%rcx,1),%eax
mov    %ebp,%r15d
sar    %r15d
cmp    %ebp,%eax
jge    11d2 <dfct+0x82>
movslq %r10d,%rax
lea    (%r8,%rax,8),%rdx
mov    %r15d,%edi
mov    0x18(%rsp),%rsi
mov    %r10,%r13
call   1060 <makect@plt>
mov    0x10(%rsp),%r8
mov    %r13,%r10
mov    %r15d,%ecx
movslq %r15d,%rdi
movsd  (%r14,%rdi,8),%xmm0
movsd  (%r14),%xmm1
movslq %ebp,%r13
movsd  (%r14,%r13,8),%xmm2
movapd %xmm1,%xmm3
addsd  %xmm2,%xmm3
subsd  %xmm2,%xmm1
movsd  %xmm1,(%r14)
movapd %xmm3,%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx)
addsd  %xmm0,%xmm3
movsd  %xmm3,(%rbx,%rdi,8)
cmp    $0x3,%ebp
jl     1368 <dfct+0x218>
mov    %r15,0x8(%rsp)
mov    %rcx,0x30(%rsp)
mov    %ebp,%r12d
shr    $0x2,%r12d
mov    %rbp,0x38(%rsp)
cmp    $0x8,%ebp
jb     12d4 <dfct+0x184>
cmp    $0x3,%r12d
mov    $0x2,%ebp
cmovae %r12d,%ebp
lea    (%r14,%r13,8),%rax
mov    $0x1,%ecx
mov    $0x1,%edx
sub    %rdi,%rdx
lea    (%r14,%rdi,8),%rsi
add    $0xfffffffffffffff8,%rsi
lea    (%rbx,%rdi,8),%r11
add    $0xfffffffffffffff8,%r11
lea    -0x8(,%rbp,8),%r9
ff
xor    %ebp,%ebp
nop
movsd  (%r14,%rcx,8),%xmm0
movsd  -0x8(%rax,%rbp,1),%xmm1
movapd %xmm0,%xmm2
subsd  %xmm1,%xmm2
addsd  %xmm1,%xmm0
movsd  (%rsi,%rbp,1),%xmm1
movsd  (%rax,%rdx,8),%xmm3
movapd %xmm1,%xmm4
subsd  %xmm3,%xmm4
addsd  %xmm3,%xmm1
movsd  %xmm2,(%r14,%rcx,8)
movsd  %xmm4,(%rsi,%rbp,1)
movapd %xmm0,%xmm2
subsd  %xmm1,%xmm2
movsd  %xmm2,(%rbx,%rcx,8)
addsd  %xmm0,%xmm1
movsd  %xmm1,(%r11,%rbp,1)
add    $0x1,%rcx
add    $0xfffffffffffffff8,%rbp
add    $0x1,%rdx
mov    %r9,%rdi
add    %rbp,%rdi
jne    1270 <dfct+0x120>
mov    %r13,0x40(%rsp)
mov    %r12d,%eax
movsd  (%r14,%rax,8),%xmm0
mov    0x38(%rsp),%r15
mov    %r15d,%ecx
sub    %r12d,%ecx
movslq %ecx,%rcx
addsd  (%r14,%rcx,8),%xmm0
movsd  %xmm0,(%rbx,%rax,8)
movsd  (%r14,%rax,8),%xmm0
subsd  (%r14,%rcx,8),%xmm0
movsd  %xmm0,(%r14,%rax,8)
movslq %r10d,%rax
lea    (%r8,%rax,8),%r13
mov    0x8(%rsp),%rbp
mov    %ebp,%edi
mov    %r14,%rsi
mov    0x30(%rsp),%rdx
mov    %r13,%rcx
mov    %r10,0x20(%rsp)
call   1030 <dctsub@plt>
cmp    $0xa,%r15d
mov    %r13,0x50(%rsp)
jl     1387 <dfct+0x237>
mov    %ebp,%edi
mov    %r14,%rsi
mov    0x18(%rsp),%rdx
mov    0x20(%rsp),%rcx
mov    0x10(%rsp),%r8
call   1070 <cftfsub@plt>
mov    %ebp,%edi
mov    %r14,%rsi
mov    0x30(%rsp),%rdx
mov    %r13,%rcx
call   1050 <rftfsub@plt>
jmp    13a8 <dfct+0x258>
movsd  (%r14),%xmm0
movsd  %xmm0,0x8(%r14)
movsd  (%rbx),%xmm0
movsd  %xmm0,0x10(%r14)
movsd  0x8(%rbx),%xmm0
jmp    1819 <dfct+0x6c9>
cmp    $0x4,%ebp
jne    13a8 <dfct+0x258>
mov    $0x4,%edi
mov    %r14,%rsi
mov    0x18(%rsp),%rdx
mov    0x20(%rsp),%rcx
mov    0x10(%rsp),%r8
call   1070 <cftfsub@plt>
movsd  (%r14),%xmm0
subsd  0x8(%r14),%xmm0
lea    -0x1(%r15),%eax
movsd  %xmm0,(%r14,%rax,8)
movsd  (%r14),%xmm0
addsd  0x8(%r14),%xmm0
movsd  %xmm0,0x8(%r14)
cmp    $0x8,%r15d
jl     1428 <dfct+0x2d8>
mov    0x8(%rsp),%rcx
mov    %ecx,%eax
add    $0x2,%eax
add    %ecx,%ecx
add    $0xfffffffb,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
lea    -0x4(%rax),%edx
movsd  (%r14,%rdx,8),%xmm0
lea    -0x3(%rax),%edx
movsd  (%r14,%rdx,8),%xmm1
addsd  %xmm0,%xmm1
lea    0x2(%rcx),%esi
movsd  %xmm1,(%r14,%rsi,8)
subsd  (%r14,%rdx,8),%xmm0
mov    %ecx,%edx
movsd  %xmm0,(%r14,%rdx,8)
add    $0xfffffffe,%eax
add    $0xfffffffc,%ecx
cmp    $0x5,%eax
jg     13f0 <dfct+0x2a0>
mov    $0x2,%eax
cmp    $0x8,%r15d
jb     17f0 <dfct+0x6a0>
lea    0x8(%rbx),%rax
mov    %rax,0x48(%rsp)
movl   $0x2,0x8(%rsp)
jmp    1470 <dfct+0x320>
nopw   0x0(%rax,%rax,1)
lea    (%r12,%r13,1),%eax
movsd  (%rbx,%rax,8),%xmm0
mov    %r12d,%eax
movsd  %xmm0,(%rbx,%rax,8)
cmp    $0x3,%r13d
mov    0x38(%rsp),%r15
jbe    17eb <dfct+0x69b>
mov    %r12d,%r13d
mov    %r12d,%edi
mov    %rbx,%rsi
mov    0x30(%rsp),%rbp
mov    %ebp,%edx
mov    0x50(%rsp),%r12
mov    %r12,%rcx
call   1030 <dctsub@plt>
cmp    $0x5,%r13d
jb     14d0 <dfct+0x380>
mov    %r13d,%edi
mov    %rbx,%rsi
mov    0x18(%rsp),%rdx
mov    0x20(%rsp),%rcx
mov    0x10(%rsp),%r8
call   1070 <cftfsub@plt>
mov    %r13d,%edi
mov    %rbx,%rsi
mov    %ebp,%edx
mov    %r12,%rcx
call   1050 <rftfsub@plt>
mov    0x8(%rsp),%r12d
jmp    14f7 <dfct+0x3a7>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
cmp    $0x4,%r13d
mov    0x8(%rsp),%r12d
jne    14f7 <dfct+0x3a7>
mov    $0x4,%edi
mov    %rbx,%rsi
mov    0x18(%rsp),%rdx
mov    0x20(%rsp),%rcx
mov    0x10(%rsp),%r8
call   1070 <cftfsub@plt>
movsd  (%rbx),%xmm0
subsd  0x8(%rbx),%xmm0
mov    %r15d,%eax
sub    %r12d,%eax
cltq
movsd  %xmm0,(%r14,%rax,8)
movsd  (%rbx),%xmm0
addsd  0x8(%rbx),%xmm0
movslq %r12d,%rax
movsd  %xmm0,(%r14,%rax,8)
mov    %r13d,%ecx
cmp    $0x3,%r13d
jb     157c <dfct+0x42c>
lea    0x0(,%rax,4),%edx
movslq %edx,%rbp
lea    (%r14,%rax,8),%rdx
shl    $0x3,%rax
mov    %r14,%rsi
sub    %rax,%rsi
mov    $0x2,%edi
mov    %rbp,%rax
nopl   0x0(%rax)
movsd  (%rbx,%rdi,8),%xmm0
subsd  0x8(%rbx,%rdi,8),%xmm0
movsd  %xmm0,(%rsi,%rax,8)
movsd  (%rbx,%rdi,8),%xmm0
addsd  0x8(%rbx,%rdi,8),%xmm0
movsd  %xmm0,(%rdx,%rax,8)
add    $0x2,%rdi
add    %rbp,%rax
cmp    %rcx,%rdi
jb     1550 <dfct+0x400>
add    %r12d,%r12d
mov    %r12d,0x8(%rsp)
mov    %r13d,%r12d
shr    %r12d
cmp    $0x1,%r13d
jbe    17da <dfct+0x68a>
cmp    $0x2,%r12d
mov    $0x1,%esi
cmovae %r12d,%esi
cmp    $0x2,%esi
jae    1670 <dfct+0x520>
xor    %eax,%eax
mov    %rax,%rdx
test   $0x1,%sil
je     15e7 <dfct+0x497>
mov    %rcx,%rdx
sub    %rax,%rdx
lea    (%rdx,%rcx,1),%rdi
movsd  (%rbx,%rdi,8),%xmm0
lea    (%rax,%rcx,1),%rbp
movsd  (%rbx,%rbp,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rbx,%rax,8)
addsd  (%rbx,%rdi,8),%xmm1
movsd  %xmm1,(%rbx,%rdx,8)
mov    %rax,%rdx
or     $0x1,%rdx
not    %rax
add    %rsi,%rax
je     1450 <dfct+0x300>
lea    (%rdx,%rcx,1),%rax
add    $0x1,%rax
mov    %rcx,%rdi
sub    %rdx,%rdi
shl    $0x4,%rcx
lea    0x0(,%rdx,8),%rbp
sub    %rbp,%rcx
movsd  (%rbx,%rcx,1),%xmm0
movsd  -0x8(%rbx,%rax,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rbx,%rdx,8)
addsd  (%rbx,%rcx,1),%xmm1
movsd  %xmm1,(%rbx,%rdi,8)
movsd  -0x8(%rbx,%rcx,1),%xmm0
movsd  (%rbx,%rax,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%rbx,%rdx,8)
addsd  -0x8(%rbx,%rcx,1),%xmm1
movsd  %xmm1,-0x8(%rbx,%rdi,8)
add    $0x2,%rdx
add    $0x2,%rax
add    $0xfffffffffffffffe,%rdi
add    $0xfffffffffffffff0,%rcx
cmp    %rdx,%rsi
jne    1610 <dfct+0x4c0>
jmp    1450 <dfct+0x300>
nopl   0x0(%rax)
lea    -0x1(%rsi),%rax
mov    $0x8,%edx
mul    %rdx
lea    (%rbx,%rcx,8),%r8
seto   %dl
mov    %r8,%rdi
sub    %rax,%rdi
cmp    %r8,%rdi
ja     15aa <dfct+0x45a>
test   %dl,%dl
jne    15aa <dfct+0x45a>
mov    %rcx,%rbp
shl    $0x4,%rbp
lea    (%rbx,%rbp,1),%r10
mov    %r10,%rdi
sub    %rax,%rdi
cmp    %r10,%rdi
ja     15aa <dfct+0x45a>
mov    $0x0,%eax
test   %dl,%dl
jne    15ac <dfct+0x45c>
mov    0x48(%rsp),%rdx
cmp    %rbx,%rdx
seta   0x4(%rsp)
lea    (%rbx,%rsi,8),%rdi
lea    (%rcx,%rsi,1),%rax
lea    (%rbx,%rax,8),%r11
mov    %r11,0x28(%rsp)
mov    %rcx,%rax
sub    %rsi,%rax
lea    (%rdx,%rax,8),%r9
lea    (%rdx,%rcx,8),%r15
lea    (%rcx,%rcx,1),%rax
sub    %rsi,%rax
add    %rdx,%rbp
cmp    %rbx,%r11
seta   0x7(%rsp)
lea    (%rdx,%rax,8),%rdx
cmp    %rdi,%r8
setb   0x6(%rsp)
cmp    %rbx,%r15
seta   %al
cmp    %rdi,%r9
setb   0x5(%rsp)
cmp    %rbx,%rbp
seta   %r11b
cmp    %rdi,%rdx
setb   %dil
cmp    0x28(%rsp),%r9
setb   0x28(%rsp)
cmp    %rbp,%r9
setb   %r9b
cmp    %r15,%rdx
setb   %bpl
mov    0x6(%rsp),%dl
test   %dl,0x7(%rsp)
jne    15aa <dfct+0x45a>
and    0x5(%rsp),%al
jne    15aa <dfct+0x45a>
and    %dil,%r11b
jne    15aa <dfct+0x45a>
mov    0x28(%rsp),%al
and    %al,0x4(%rsp)
jne    15aa <dfct+0x45a>
mov    $0x0,%eax
and    %bpl,%r9b
jne    15ac <dfct+0x45c>
mov    %esi,%eax
and    $0x1ffffffe,%eax
mov    $0xfffffffffffffff8,%rbp
xor    %edx,%edx
nopl   0x0(%rax)
movupd (%r10,%rbp,1),%xmm0
shufpd $0x1,%xmm0,%xmm0
movupd (%r8,%rdx,8),%xmm1
subpd  %xmm1,%xmm0
movupd %xmm0,(%rbx,%rdx,8)
movupd (%r10,%rbp,1),%xmm0
shufpd $0x1,%xmm1,%xmm1
addpd  %xmm0,%xmm1
movupd %xmm1,(%r8,%rbp,1)
add    $0x2,%rdx
add    $0xfffffffffffffff0,%rbp
cmp    %rdx,%rax
jne    1790 <dfct+0x640>
cmp    %rsi,%rax
je     1450 <dfct+0x300>
jmp    15ac <dfct+0x45c>
add    %r12d,%r13d
movsd  (%rbx,%r13,8),%xmm0
mov    %r12d,%eax
movsd  %xmm0,(%rbx,%rax,8)
movslq 0x8(%rsp),%rax
movsd  (%rbx),%xmm0
movsd  %xmm0,(%r14,%rax,8)
movsd  0x10(%rbx),%xmm0
subsd  0x8(%rbx),%xmm0
mov    0x40(%rsp),%rax
movsd  %xmm0,(%r14,%rax,8)
movsd  0x10(%rbx),%xmm0
addsd  0x8(%rbx),%xmm0
movsd  %xmm0,(%r14)
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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