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 bitrv2(int n, int *ip, double *a);
  void cftfsub(int n, double *a, 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) {
      bitrv2(m, ip + 2, a);
      cftfsub(m, a, w);
      rftfsub(m, a, nc, w + nw);
    } else if (m == 4) {
      cftfsub(m, a, 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) {
        bitrv2(m, ip + 2, t);
        cftfsub(m, t, w);
        rftfsub(m, t, nc, w + nw);
      } else if (m == 4) {
        cftfsub(m, t, 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    11ae <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    11ed <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   1070 <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    17c6 <dfct+0x666>
mov    -0x38(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x3c(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    1387 <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    1282 <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    144a <dfct+0x2ea>
mov    -0x38(%rbp),%edi
mov    -0x20(%rbp),%rsi
add    $0x8,%rsi
mov    -0x10(%rbp),%rdx
call   1060 <bitrv2@plt>
mov    -0x38(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x28(%rbp),%rdx
call   1080 <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    1469 <dfct+0x309>
cmpl   $0x4,-0x38(%rbp)
jne    1464 <dfct+0x304>
mov    -0x38(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x28(%rbp),%rdx
call   1080 <cftfsub@plt>
jmp    1469 <dfct+0x309>
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     152f <dfct+0x3cf>
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    14af <dfct+0x34f>
movl   $0x2,-0x34(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x38(%rbp)
cmpl   $0x2,-0x38(%rbp)
jl     1773 <dfct+0x613>
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    15b5 <dfct+0x455>
mov    -0x38(%rbp),%edi
mov    -0x20(%rbp),%rsi
add    $0x8,%rsi
mov    -0x18(%rbp),%rdx
call   1060 <bitrv2@plt>
mov    -0x38(%rbp),%edi
mov    -0x18(%rbp),%rsi
mov    -0x28(%rbp),%rdx
call   1080 <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    15d4 <dfct+0x474>
cmpl   $0x4,-0x38(%rbp)
jne    15cf <dfct+0x46f>
mov    -0x38(%rbp),%edi
mov    -0x18(%rbp),%rsi
mov    -0x28(%rbp),%rdx
call   1080 <cftfsub@plt>
jmp    15d4 <dfct+0x474>
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    16ab <dfct+0x54b>
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    1623 <dfct+0x4c3>
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    1749 <dfct+0x5e9>
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    16c4 <dfct+0x564>
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    153c <dfct+0x3dc>
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    17f9 <dfct+0x699>
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    $0x48,%rsp
mov    %r8,%r12
mov    %rcx,%rbp
mov    %rdx,%rbx
mov    %rsi,%r14
mov    (%rcx),%r13d
lea    0x0(,%r13,8),%eax
cmp    %edi,%eax
mov    %rdi,0x8(%rsp)
jge    11a8 <dfct+0x48>
mov    %edi,%r13d
sar    $0x3,%r13d
mov    %r13d,%edi
mov    %rbp,%rsi
mov    %r12,%rdx
call   1040 <makewt@plt>
mov    0x8(%rsp),%rdi
mov    0x4(%rbp),%eax
mov    %rax,0x10(%rsp)
add    %eax,%eax
mov    %edi,%r15d
sar    %r15d
cmp    %edi,%eax
jge    11db <dfct+0x7b>
movslq %r13d,%rax
lea    (%r12,%rax,8),%rdx
mov    %r15d,%edi
mov    %rbp,%rsi
call   1070 <makect@plt>
mov    0x8(%rsp),%rdi
mov    %r15d,%eax
mov    %rax,0x10(%rsp)
movslq %r15d,%rax
movsd  (%r14,%rax,8),%xmm0
movsd  (%r14),%xmm1
movslq %edi,%rcx
movsd  (%r14,%rcx,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,%rax,8)
cmp    $0x3,%edi
jl     137a <dfct+0x21a>
mov    %rcx,0x28(%rsp)
mov    %rbp,0x20(%rsp)
mov    %r12,0x30(%rsp)
mov    %edi,%r12d
shr    $0x2,%r12d
cmp    $0x8,%edi
jb     12e6 <dfct+0x186>
cmp    $0x3,%r12d
mov    $0x2,%ebp
cmovae %r12d,%ebp
mov    0x28(%rsp),%rcx
lea    (%r14,%rcx,8),%r10
mov    $0x1,%ecx
mov    $0x1,%edx
sub    %rax,%rdx
lea    (%r14,%rax,8),%rsi
add    $0xfffffffffffffff8,%rsi
lea    (%rbx,%rax,8),%r8
add    $0xfffffffffffffff8,%r8
lea    -0x8(,%rbp,8),%r9
ff
xor    %eax,%eax
nopl   (%rax)
movsd  (%r14,%rcx,8),%xmm0
movsd  -0x8(%r10,%rax,1),%xmm1
movapd %xmm0,%xmm2
subsd  %xmm1,%xmm2
addsd  %xmm1,%xmm0
movsd  (%rsi,%rax,1),%xmm1
movsd  (%r10,%rdx,8),%xmm3
movapd %xmm1,%xmm4
subsd  %xmm3,%xmm4
addsd  %xmm3,%xmm1
movsd  %xmm2,(%r14,%rcx,8)
movsd  %xmm4,(%rsi,%rax,1)
movapd %xmm0,%xmm2
subsd  %xmm1,%xmm2
movsd  %xmm2,(%rbx,%rcx,8)
addsd  %xmm0,%xmm1
movsd  %xmm1,(%r8,%rax,1)
add    $0x1,%rcx
add    $0xfffffffffffffff8,%rax
add    $0x1,%rdx
mov    %r9,%rbp
add    %rax,%rbp
jne    1280 <dfct+0x120>
mov    %r12d,%eax
movsd  (%r14,%rax,8),%xmm0
mov    %edi,%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 %r13d,%rax
mov    0x30(%rsp),%r13
lea    0x0(%r13,%rax,8),%rbp
mov    %r15d,%edi
mov    %r14,%rsi
mov    0x10(%rsp),%rdx
mov    %rbp,%rcx
call   1030 <dctsub@plt>
mov    0x8(%rsp),%rdi
cmp    $0xa,%edi
mov    %rbp,0x40(%rsp)
jl     1399 <dfct+0x239>
mov    0x20(%rsp),%rax
lea    0x8(%rax),%rsi
mov    %r15d,%edi
mov    %r14,%rdx
call   1060 <bitrv2@plt>
mov    %r15d,%edi
mov    %r14,%rsi
mov    %r13,%rdx
call   1080 <cftfsub@plt>
mov    %r15d,%edi
mov    %r14,%rsi
mov    0x10(%rsp),%rdx
mov    %rbp,%rcx
call   1050 <rftfsub@plt>
jmp    13af <dfct+0x24f>
movsd  (%r14),%xmm0
movsd  %xmm0,0x8(%r14)
movsd  (%rbx),%xmm0
movsd  %xmm0,0x10(%r14)
movsd  0x8(%rbx),%xmm0
jmp    1819 <dfct+0x6b9>
cmp    $0x4,%r15d
jne    13b4 <dfct+0x254>
mov    $0x4,%edi
mov    %r14,%rsi
mov    %r13,%rdx
call   1080 <cftfsub@plt>
mov    0x8(%rsp),%rdi
movsd  (%r14),%xmm0
subsd  0x8(%r14),%xmm0
lea    -0x1(%rdi),%eax
movsd  %xmm0,(%r14,%rax,8)
movsd  (%r14),%xmm0
addsd  0x8(%r14),%xmm0
movsd  %xmm0,0x8(%r14)
cmp    $0x8,%edi
jl     1428 <dfct+0x2c8>
mov    %r15d,%eax
add    $0x2,%eax
lea    (%r15,%r15,1),%ecx
add    $0xfffffffb,%ecx
nopl   0x0(%rax,%rax,1)
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+0x290>
mov    $0x2,%eax
cmp    $0x8,%edi
jb     17f0 <dfct+0x690>
addq   $0x8,0x20(%rsp)
lea    0x8(%rbx),%rax
mov    %rax,0x38(%rsp)
movl   $0x2,0x4(%rsp)
jmp    1470 <dfct+0x310>
nop
lea    (%r12,%r15,1),%eax
movsd  (%rbx,%rax,8),%xmm0
mov    %r12d,%eax
movsd  %xmm0,(%rbx,%rax,8)
cmp    $0x3,%r15d
mov    0x30(%rsp),%r13
jbe    17eb <dfct+0x68b>
mov    %r12d,%r15d
mov    %r12d,%edi
mov    %rbx,%rsi
mov    0x10(%rsp),%rbp
mov    %ebp,%edx
mov    0x40(%rsp),%r12
mov    %r12,%rcx
call   1030 <dctsub@plt>
cmp    $0x5,%r15d
jb     14d0 <dfct+0x370>
mov    %r15d,%edi
mov    0x20(%rsp),%rsi
mov    %rbx,%rdx
call   1060 <bitrv2@plt>
mov    %r15d,%edi
mov    %rbx,%rsi
mov    %r13,%rdx
call   1080 <cftfsub@plt>
mov    %r15d,%edi
mov    %rbx,%rsi
mov    %ebp,%edx
mov    %r12,%rcx
call   1050 <rftfsub@plt>
mov    0x8(%rsp),%rax
mov    0x4(%rsp),%r12d
jmp    14f5 <dfct+0x395>
nopl   (%rax)
cmp    $0x4,%r15d
mov    0x8(%rsp),%rax
mov    0x4(%rsp),%r12d
jne    14f5 <dfct+0x395>
mov    $0x4,%edi
mov    %rbx,%rsi
mov    %r13,%rdx
call   1080 <cftfsub@plt>
mov    0x8(%rsp),%rax
movsd  (%rbx),%xmm0
subsd  0x8(%rbx),%xmm0
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    %r15d,%ecx
cmp    $0x3,%r15d
jb     157c <dfct+0x41c>
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
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
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+0x3f0>
add    %r12d,%r12d
mov    %r12d,0x4(%rsp)
mov    %r15d,%r12d
shr    %r12d
cmp    $0x1,%r15d
jbe    17da <dfct+0x67a>
cmp    $0x2,%r12d
mov    $0x1,%esi
cmovae %r12d,%esi
cmp    $0x2,%esi
jae    1670 <dfct+0x510>
xor    %eax,%eax
mov    %rax,%rdx
test   $0x1,%sil
je     15e7 <dfct+0x487>
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+0x2f0>
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+0x4b0>
jmp    1450 <dfct+0x2f0>
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+0x44a>
test   %dl,%dl
jne    15aa <dfct+0x44a>
mov    %rcx,%rbp
shl    $0x4,%rbp
lea    (%rbx,%rbp,1),%r10
mov    %r10,%rdi
sub    %rax,%rdi
cmp    %r10,%rdi
ja     15aa <dfct+0x44a>
mov    $0x0,%eax
test   %dl,%dl
jne    15ac <dfct+0x44c>
mov    0x38(%rsp),%rdx
cmp    %rbx,%rdx
seta   (%rsp)
lea    (%rbx,%rsi,8),%rdi
lea    (%rcx,%rsi,1),%rax
lea    (%rbx,%rax,8),%r11
mov    %r11,0x18(%rsp)
mov    %rcx,%rax
sub    %rsi,%rax
lea    (%rdx,%rax,8),%r9
lea    (%rdx,%rcx,8),%r13
lea    (%rcx,%rcx,1),%rax
sub    %rsi,%rax
add    %rdx,%rbp
cmp    %rbx,%r11
seta   0x3(%rsp)
lea    (%rdx,%rax,8),%rdx
cmp    %rdi,%r8
setb   0x2(%rsp)
cmp    %rbx,%r13
seta   %al
cmp    %rdi,%r9
setb   0x1(%rsp)
cmp    %rbx,%rbp
seta   %r11b
cmp    %rdi,%rdx
setb   %dil
cmp    0x18(%rsp),%r9
setb   0x18(%rsp)
cmp    %rbp,%r9
setb   %r9b
cmp    %r13,%rdx
setb   %bpl
mov    0x2(%rsp),%dl
test   %dl,0x3(%rsp)
jne    15aa <dfct+0x44a>
and    0x1(%rsp),%al
jne    15aa <dfct+0x44a>
and    %dil,%r11b
jne    15aa <dfct+0x44a>
mov    0x18(%rsp),%al
and    %al,(%rsp)
jne    15aa <dfct+0x44a>
mov    $0x0,%eax
and    %bpl,%r9b
jne    15ac <dfct+0x44c>
mov    %esi,%eax
and    $0x1ffffffe,%eax
mov    $0xfffffffffffffff8,%rbp
xor    %edx,%edx
nopw   0x0(%rax,%rax,1)
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+0x630>
cmp    %rsi,%rax
je     1450 <dfct+0x2f0>
jmp    15ac <dfct+0x44c>
add    %r12d,%r15d
movsd  (%rbx,%r15,8),%xmm0
mov    %r12d,%eax
movsd  %xmm0,(%rbx,%rax,8)
movslq 0x4(%rsp),%rax
movsd  (%rbx),%xmm0
movsd  %xmm0,(%r14,%rax,8)
movsd  0x10(%rbx),%xmm0
subsd  0x8(%rbx),%xmm0
mov    0x28(%rsp),%rax
movsd  %xmm0,(%r14,%rax,8)
movsd  0x10(%rbx),%xmm0
addsd  0x8(%rbx),%xmm0
movsd  %xmm0,(%r14)
add    $0x48,%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 %237
222: ; preds = %182
%223 = load i32, i32* %14
%224 = load i32*, i32** %9
%225 = getelementptr inbounds i32, i32* %224, i64 2
%226 = load double*, double** %7
%227 = load i32, i32* %14
%228 = load double*, double** %7
%229 = load double*, double** %10
%230 = load i32, i32* %14
%231 = load double*, double** %7
%232 = load i32, i32* %17
%233 = load double*, double** %10
%234 = load i32, i32* %16
%235 = sext i32 %234 to i64
%236 = getelementptr inbounds double, double* %233, i64 %235
br label %245
237: ; preds = %182
%238 = load i32, i32* %14
%239 = icmp eq i32 %238, 4
br i1 %239, label %240, label %244
240: ; preds = %237
%241 = load i32, i32* %14
%242 = load double*, double** %7
%243 = load double*, double** %10
br label %244
244: ; preds = %240, %237
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 %342
327: ; preds = %317
%328 = load i32, i32* %14
%329 = load i32*, i32** %9
%330 = getelementptr inbounds i32, i32* %329, i64 2
%331 = load double*, double** %8
%332 = load i32, i32* %14
%333 = load double*, double** %8
%334 = load double*, double** %10
%335 = load i32, i32* %14
%336 = load double*, double** %8
%337 = load i32, i32* %17
%338 = load double*, double** %10
%339 = load i32, i32* %16
%340 = sext i32 %339 to i64
%341 = getelementptr inbounds double, double* %338, i64 %340
br label %350
342: ; preds = %317
%343 = load i32, i32* %14
%344 = icmp eq i32 %343, 4
br i1 %344, label %345, label %349
345: ; preds = %342
%346 = load i32, i32* %14
%347 = load double*, double** %8
%348 = load double*, double** %10
br label %349
349: ; preds = %345, %342
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 %344
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 %83
81: ; preds = %65
%82 = getelementptr inbounds i32, i32* %3, i64 2
tail call void @bitrv2(i32 %17, i32* %82, double* %1)
tail call void @cftfsub(i32 %17, double* %1, double* %4)
tail call void @rftfsub(i32 %17, double* %1, i32 %22, double* %79)
br label %86
83: ; preds = %65
%84 = icmp eq i32 %17, 4
br i1 %84, label %85, label %86
85: ; preds = %83
tail call void @cftfsub(i32 4, double* %1, double* %4)
br label %86
86: ; preds = %83, %85, %81
%87 = load double, double* %1
%88 = getelementptr inbounds double, double* %1, i64 1
%89 = load double, double* %88
%90 = fsub double %87, %89
%91 = add nsw i32 %0, -1
%92 = zext i32 %91 to i64
%93 = getelementptr inbounds double, double* %1, i64 %92
store double %90, double* %93
%94 = load double, double* %1
%95 = load double, double* %88
%96 = fadd double %94, %95
store double %96, double* %88
%97 = icmp sgt i32 %0, 7
br i1 %97, label %109, label %98
98: ; preds = %109, %86
br i1 %38, label %99, label %332
99: ; preds = %98
%100 = getelementptr inbounds i32, i32* %3, i64 2
%101 = getelementptr inbounds double, double* %2, i64 1
%102 = getelementptr double, double* %2, i64 1
%103 = getelementptr double, double* %2, i64 1
%104 = getelementptr double, double* %2, i64 1
%105 = getelementptr double, double* %2, i64 1
%106 = icmp ugt double* %103, %2
%107 = getelementptr inbounds double, double* %2, i64 -1
%108 = getelementptr inbounds double, double* %2, i64 -1
br label %130
109: ; preds = %86, %109
%110 = phi i32 [ %111, %109 ], [ %17, %86 ]
%111 = add nsw i32 %110, -2
%112 = zext i32 %111 to i64
%113 = getelementptr inbounds double, double* %1, i64 %112
%114 = load double, double* %113
%115 = add nsw i32 %110, -1
%116 = zext i32 %115 to i64
%117 = getelementptr inbounds double, double* %1, i64 %116
%118 = load double, double* %117
%119 = fadd double %114, %118
%120 = shl nuw nsw i32 %111, 1
%121 = or i32 %120, 1
%122 = zext i32 %121 to i64
%123 = getelementptr inbounds double, double* %1, i64 %122
store double %119, double* %123
%124 = load double, double* %117
%125 = fsub double %114, %124
%126 = add nsw i32 %120, -1
%127 = zext i32 %126 to i64
%128 = getelementptr inbounds double, double* %1, i64 %127
store double %125, double* %128
%129 = icmp sgt i32 %110, 5
br i1 %129, label %109, label %98
130: ; preds = %99, %322
%131 = phi i32 [ %37, %99 ], [ %176, %322 ]
%132 = phi i32 [ 2, %99 ], [ %175, %322 ]
tail call void @dctsub(i32 %131, double* %2, i32 %22, double* %79)
%133 = icmp ugt i32 %131, 4
br i1 %133, label %134, label %135
134: ; preds = %130
tail call void @bitrv2(i32 %131, i32* %100, double* %2)
tail call void @cftfsub(i32 %131, double* %2, double* %4)
tail call void @rftfsub(i32 %131, double* %2, i32 %22, double* %79)
br label %138
135: ; preds = %130
%136 = icmp eq i32 %131, 4
br i1 %136, label %137, label %138
137: ; preds = %135
tail call void @cftfsub(i32 4, double* %2, double* %4)
br label %138
138: ; preds = %135, %137, %134
%139 = load double, double* %2
%140 = load double, double* %101
%141 = fsub double %139, %140
%142 = sub nsw i32 %0, %132
%143 = sext i32 %142 to i64
%144 = getelementptr inbounds double, double* %1, i64 %143
store double %141, double* %144
%145 = load double, double* %2
%146 = load double, double* %101
%147 = fadd double %145, %146
%148 = sext i32 %132 to i64
%149 = getelementptr inbounds double, double* %1, i64 %148
store double %147, double* %149
%150 = icmp ugt i32 %131, 2
br i1 %150, label %151, label %174
151: ; preds = %138
%152 = shl i32 %132, 2
%153 = zext i32 %131 to i64
%154 = sext i32 %152 to i64
br label %155
155: ; preds = %151, %155
%156 = phi i64 [ 0, %151 ], [ %158, %155 ]
%157 = phi i64 [ 2, %151 ], [ %172, %155 ]
%158 = add nsw i64 %156, %154
%159 = getelementptr inbounds double, double* %2, i64 %157
%160 = load double, double* %159
%161 = or i64 %157, 1
%162 = getelementptr inbounds double, double* %2, i64 %161
%163 = load double, double* %162
%164 = fsub double %160, %163
%165 = sub nsw i64 %158, %148
%166 = getelementptr inbounds double, double* %1, i64 %165
store double %164, double* %166
%167 = load double, double* %159
%168 = load double, double* %162
%169 = fadd double %167, %168
%170 = add nsw i64 %158, %148
%171 = getelementptr inbounds double, double* %1, i64 %170
store double %169, double* %171
%172 = add nuw nsw i64 %157, 2
%173 = icmp ult i64 %172, %153
br i1 %173, label %155, label %174
174: ; preds = %155, %138
%175 = shl i32 %132, 1
%176 = lshr i32 %131, 1
%177 = icmp ult i32 %131, 2
br i1 %177, label %178, label %185
178: ; preds = %174
%179 = add nuw nsw i32 %176, %131
%180 = zext i32 %179 to i64
%181 = getelementptr inbounds double, double* %2, i64 %180
%182 = load double, double* %181
%183 = zext i32 %176 to i64
%184 = getelementptr inbounds double, double* %2, i64 %183
store double %182, double* %184
br label %330
185: ; preds = %174
%186 = zext i32 %131 to i64
%187 = call i32 @llvm.umax.i32(i32 %176, i32 1)
%188 = zext i32 %187 to i64
%189 = icmp ult i32 %187, 2
br i1 %189, label %270, label %190
190: ; preds = %185
%191 = add nsw i64 %188, -1
%192 = getelementptr double, double* %2, i64 %186
%193 = call { i64, i1 } @llvm.umul.with.overflow.i64(i64 %191, i64 8)
%194 = extractvalue { i64, i1 } %193, 0
%195 = extractvalue { i64, i1 } %193, 1
%196 = bitcast double* %192 to i8*
%197 = sub i64 0, %194
%198 = getelementptr i8, i8* %196, i64 %197
%199 = icmp ugt i8* %198, %196
%200 = or i1 %199, %195
%201 = shl nuw nsw i64 %186, 1
%202 = getelementptr double, double* %2, i64 %201
%203 = call { i64, i1 } @llvm.umul.with.overflow.i64(i64 %191, i64 8)
%204 = extractvalue { i64, i1 } %203, 0
%205 = extractvalue { i64, i1 } %203, 1
%206 = bitcast double* %202 to i8*
%207 = sub i64 0, %204
%208 = getelementptr i8, i8* %206, i64 %207
%209 = icmp ugt i8* %208, %206
%210 = or i1 %209, %205
%211 = or i1 %200, %210
br i1 %211, label %270, label %212
212: ; preds = %190
%213 = getelementptr double, double* %2, i64 %188
%214 = getelementptr double, double* %2, i64 %186
%215 = add nuw nsw i64 %186, %188
%216 = getelementptr double, double* %2, i64 %215
%217 = sub nsw i64 %186, %188
%218 = getelementptr double, double* %102, i64 %217
%219 = getelementptr double, double* %103, i64 %186
%220 = shl nuw nsw i64 %186, 1
%221 = sub nsw i64 %220, %188
%222 = getelementptr double, double* %104, i64 %221
%223 = getelementptr double, double* %105, i64 %220
%224 = icmp ugt double* %216, %2
%225 = icmp ult double* %214, %213
%226 = and i1 %224, %225
%227 = icmp ugt double* %219, %2
%228 = icmp ult double* %218, %213
%229 = and i1 %227, %228
%230 = or i1 %226, %229
%231 = icmp ugt double* %223, %2
%232 = icmp ult double* %222, %213
%233 = and i1 %231, %232
%234 = or i1 %230, %233
%235 = icmp ult double* %218, %216
%236 = and i1 %106, %235
%237 = or i1 %234, %236
%238 = icmp ult double* %218, %223
%239 = icmp ult double* %222, %219
%240 = and i1 %238, %239
%241 = or i1 %237, %240
br i1 %241, label %270, label %242
242: ; preds = %212
%243 = and i64 %188, 536870910
br label %244
244: ; preds = %244, %242
%245 = phi i64 [ 0, %242 ], [ %266, %244 ]
%246 = sub nsw i64 %186, %245
%247 = add nsw i64 %246, %186
%248 = getelementptr inbounds double, double* %107, i64 %247
%249 = bitcast double* %248 to <2 x double>*
%250 = load <2 x double>, <2 x double>* %249
%251 = shufflevector <2 x double> %250, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%252 = add nuw nsw i64 %245, %186
%253 = getelementptr inbounds double, double* %2, i64 %252
%254 = bitcast double* %253 to <2 x double>*
%255 = load <2 x double>, <2 x double>* %254
%256 = fsub <2 x double> %251, %255
%257 = getelementptr inbounds double, double* %2, i64 %245
%258 = bitcast double* %257 to <2 x double>*
store <2 x double> %256, <2 x double>* %258
%259 = bitcast double* %248 to <2 x double>*
%260 = load <2 x double>, <2 x double>* %259
%261 = shufflevector <2 x double> %260, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%262 = fadd <2 x double> %255, %261
%263 = shufflevector <2 x double> %262, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%264 = getelementptr inbounds double, double* %108, i64 %246
%265 = bitcast double* %264 to <2 x double>*
store <2 x double> %263, <2 x double>* %265
%266 = add nuw i64 %245, 2
%267 = icmp eq i64 %266, %243
br i1 %267, label %268, label %244
268: ; preds = %244
%269 = icmp eq i64 %243, %188
br i1 %269, label %322, label %270
270: ; preds = %212, %190, %185, %268
%271 = phi i64 [ 0, %212 ], [ 0, %190 ], [ 0, %185 ], [ %243, %268 ]
%272 = xor i64 %271, -1
%273 = and i64 %188, 1
%274 = icmp eq i64 %273, 0
br i1 %274, label %289, label %275
275: ; preds = %270
%276 = sub nsw i64 %186, %271
%277 = add nsw i64 %276, %186
%278 = getelementptr inbounds double, double* %2, i64 %277
%279 = load double, double* %278
%280 = add nuw nsw i64 %271, %186
%281 = getelementptr inbounds double, double* %2, i64 %280
%282 = load double, double* %281
%283 = fsub double %279, %282
%284 = getelementptr inbounds double, double* %2, i64 %271
store double %283, double* %284
%285 = load double, double* %278
%286 = fadd double %282, %285
%287 = getelementptr inbounds double, double* %2, i64 %276
store double %286, double* %287
%288 = or i64 %271, 1
br label %289
289: ; preds = %275, %270
%290 = phi i64 [ %271, %270 ], [ %288, %275 ]
%291 = sub nsw i64 0, %188
%292 = icmp eq i64 %272, %291
br i1 %292, label %322, label %293
293: ; preds = %289, %293
%294 = phi i64 [ %320, %293 ], [ %290, %289 ]
%295 = sub nsw i64 %186, %294
%296 = add nsw i64 %295, %186
%297 = getelementptr inbounds double, double* %2, i64 %296
%298 = load double, double* %297
%299 = add nuw nsw i64 %294, %186
%300 = getelementptr inbounds double, double* %2, i64 %299
%301 = load double, double* %300
%302 = fsub double %298, %301
%303 = getelementptr inbounds double, double* %2, i64 %294
store double %302, double* %303
%304 = load double, double* %297
%305 = fadd double %301, %304
%306 = getelementptr inbounds double, double* %2, i64 %295
store double %305, double* %306
%307 = add nuw nsw i64 %294, 1
%308 = sub nsw i64 %186, %307
%309 = add nsw i64 %308, %186
%310 = getelementptr inbounds double, double* %2, i64 %309
%311 = load double, double* %310
%312 = add nuw nsw i64 %307, %186
%313 = getelementptr inbounds double, double* %2, i64 %312
%314 = load double, double* %313
%315 = fsub double %311, %314
%316 = getelementptr inbounds double, double* %2, i64 %307
store double %315, double* %316
%317 = load double, double* %310
%318 = fadd double %314, %317
%319 = getelementptr inbounds double, double* %2, i64 %308
store double %318, double* %319
%320 = add nuw nsw i64 %294, 2
%321 = icmp eq i64 %320, %188
br i1 %321, label %322, label %293
322: ; preds = %289, %293, %268
%323 = add nuw nsw i32 %176, %131
%324 = zext i32 %323 to i64
%325 = getelementptr inbounds double, double* %2, i64 %324
%326 = load double, double* %325
%327 = zext i32 %176 to i64
%328 = getelementptr inbounds double, double* %2, i64 %327
store double %326, double* %328
%329 = icmp ugt i32 %131, 3
br i1 %329, label %130, label %330
330: ; preds = %322, %178
%331 = sext i32 %175 to i64
br label %332
332: ; preds = %330, %98
%333 = phi i64 [ 2, %98 ], [ %331, %330 ]
%334 = load double, double* %2
%335 = getelementptr inbounds double, double* %1, i64 %333
store double %334, double* %335
%336 = getelementptr inbounds double, double* %2, i64 2
%337 = load double, double* %336
%338 = getelementptr inbounds double, double* %2, i64 1
%339 = load double, double* %338
%340 = fsub double %337, %339
store double %340, double* %28
%341 = load double, double* %336
%342 = load double, double* %338
%343 = fadd double %341, %342
br label %351
344: ; preds = %21
%345 = load double, double* %1
%346 = getelementptr inbounds double, double* %1, i64 1
store double %345, double* %346
%347 = load double, double* %2
%348 = getelementptr inbounds double, double* %1, i64 2
store double %347, double* %348
%349 = getelementptr inbounds double, double* %2, i64 1
%350 = load double, double* %349
br label %351
351: ; preds = %344, %332
%352 = phi double [ %350, %344 ], [ %343, %332 ]
store double %352, double* %1
ret void
}