void dfst(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 dstsub(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);
  }
  if (n > 2) {
    m = n >> 1;
    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[0] = a[mh] - a[n - mh];
    a[mh] += a[n - mh];
    a[0] = a[m];
    dstsub(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[1] - a[0];
    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) {
      dstsub(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[1] - t[0];
      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 = 1; j < mh; j++) {
        k = m - j;
        t[j] = t[m + k] + t[m + j];
        t[k] = t[m + k] - t[m + j];
      }
      t[0] = t[m + mh];
      m = mh;
    }
    a[l] = t[0];
  }
  a[0] = 0;
}


x86 assembly (clang, -O0):
<dfst>:
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 <dfst+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   1030 <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 <dfst+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>
cmpl   $0x2,-0x4(%rbp)
jle    172f <dfst+0x5df>
mov    -0x4(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x3c(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    1305 <dfst+0x1b5>
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
addsd  (%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
subsd  (%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
addsd  (%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
subsd  (%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
addsd  -0x68(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x58(%rbp),%xmm0
subsd  -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    1200 <dfst+0xb0>
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
subsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,(%rax)
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 <dstsub@plt>
cmpl   $0x4,-0x38(%rbp)
jle    13c7 <dfst+0x277>
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   1040 <rftfsub@plt>
jmp    13ed <dfst+0x29d>
cmpl   $0x4,-0x38(%rbp)
jne    13e8 <dfst+0x298>
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    13ed <dfst+0x29d>
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
subsd  (%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     14ca <dfst+0x37a>
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
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%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    1433 <dfst+0x2e3>
movl   $0x2,-0x34(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x38(%rbp)
cmpl   $0x2,-0x38(%rbp)
jl     171a <dfst+0x5ca>
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 <dstsub@plt>
cmpl   $0x4,-0x38(%rbp)
jle    1543 <dfst+0x3f3>
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   1040 <rftfsub@plt>
jmp    1569 <dfst+0x419>
cmpl   $0x4,-0x38(%rbp)
jne    1564 <dfst+0x414>
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    1569 <dfst+0x419>
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x18(%rbp),%rax
subsd  (%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    1657 <dfst+0x507>
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
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%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
subsd  (%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    15b8 <dfst+0x468>
mov    -0x34(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x38(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x3c(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    16f5 <dfst+0x5a5>
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
addsd  (%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
subsd  (%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    1670 <dfst+0x520>
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x38(%rbp)
jmp    14d7 <dfst+0x387>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dfst>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
mov    %r8,%r12
mov    %rcx,%r13
mov    %rdx,%rbx
mov    %rsi,%r14
mov    %edi,%r15d
mov    (%rcx),%ebp
lea    0x0(,%rbp,8),%eax
cmp    %edi,%eax
jge    118d <dfst+0x3d>
mov    %r15d,%ebp
sar    $0x3,%ebp
mov    %ebp,%edi
mov    %r13,%rsi
mov    %r12,%rdx
call   1030 <makewt@plt>
mov    0x4(%r13),%eax
mov    %rax,0x28(%rsp)
add    %eax,%eax
movslq %ebp,%rcx
cmp    %r15d,%eax
mov    %r12,0x8(%rsp)
jge    11c1 <dfst+0x71>
mov    %r15d,%edi
sar    %edi
lea    (%r12,%rcx,8),%rdx
mov    %rdi,0x28(%rsp)
mov    %r13,%rsi
mov    %rcx,%r12
call   1060 <makect@plt>
mov    %r12,%rcx
cmp    $0x3,%r15d
jl     184b <dfst+0x6fb>
mov    %rcx,0x38(%rsp)
mov    %rbp,0x20(%rsp)
mov    %r13,0x18(%rsp)
mov    %r15d,%edi
shr    %edi
mov    %r15d,%r12d
shr    $0x2,%r12d
mov    %edi,%r8d
mov    %r15,0x10(%rsp)
cmp    $0x8,%r15d
jb     12a6 <dfst+0x156>
mov    0x10(%rsp),%eax
cmp    $0x3,%r12d
mov    $0x2,%ebp
cmovae %r12d,%ebp
lea    (%r14,%rax,8),%r11
mov    $0x1,%edx
mov    $0x1,%esi
sub    %r8,%rsi
lea    (%r14,%r8,8),%rcx
add    $0xfffffffffffffff8,%rcx
lea    (%rbx,%r8,8),%r9
add    $0xfffffffffffffff8,%r9
lea    -0x8(,%rbp,8),%r10
ff
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%r14,%rdx,8),%xmm0
movsd  -0x8(%r11,%rax,1),%xmm1
movapd %xmm0,%xmm2
addsd  %xmm1,%xmm2
subsd  %xmm1,%xmm0
movsd  (%rcx,%rax,1),%xmm1
movsd  (%r11,%rsi,8),%xmm3
movapd %xmm1,%xmm4
addsd  %xmm3,%xmm4
subsd  %xmm3,%xmm1
movsd  %xmm2,(%r14,%rdx,8)
movsd  %xmm4,(%rcx,%rax,1)
movapd %xmm0,%xmm2
addsd  %xmm1,%xmm2
movsd  %xmm2,(%rbx,%rdx,8)
subsd  %xmm1,%xmm0
movsd  %xmm0,(%r9,%rax,1)
add    $0x1,%rdx
add    $0xfffffffffffffff8,%rax
add    $0x1,%rsi
mov    %r10,%rbp
add    %rax,%rbp
jne    1240 <dfst+0xf0>
mov    %r12d,%eax
movsd  (%r14,%rax,8),%xmm0
mov    0x10(%rsp),%r15
mov    %r15d,%ecx
sub    %r12d,%ecx
movslq %ecx,%rcx
subsd  (%r14,%rcx,8),%xmm0
movsd  %xmm0,(%rbx)
movsd  (%r14,%rcx,8),%xmm0
addsd  (%r14,%rax,8),%xmm0
movsd  %xmm0,(%r14,%rax,8)
movsd  (%r14,%r8,8),%xmm0
movsd  %xmm0,(%r14)
mov    0x38(%rsp),%rax
mov    0x8(%rsp),%rcx
lea    (%rcx,%rax,8),%r13
mov    %r14,%rsi
mov    0x28(%rsp),%rdx
mov    %r13,%rcx
mov    %rdi,%rbp
call   1050 <dstsub@plt>
cmp    $0xa,%r15d
mov    %r13,0x38(%rsp)
jl     133d <dfst+0x1ed>
mov    %ebp,%edi
mov    %r14,%rsi
mov    0x18(%rsp),%rdx
mov    0x20(%rsp),%rcx
mov    0x8(%rsp),%r8
call   1070 <cftfsub@plt>
mov    %ebp,%edi
mov    %r14,%rsi
mov    0x28(%rsp),%rdx
mov    %r13,%rcx
call   1040 <rftfsub@plt>
jmp    135e <dfst+0x20e>
cmp    $0x4,%ebp
mov    0x20(%rsp),%rcx
jne    135e <dfst+0x20e>
mov    $0x4,%edi
mov    %r14,%rsi
mov    0x18(%rsp),%rdx
mov    0x8(%rsp),%r8
call   1070 <cftfsub@plt>
mov    0x10(%rsp),%rcx
lea    -0x1(%rcx),%eax
movsd  0x8(%r14),%xmm0
subsd  (%r14),%xmm0
movsd  %xmm0,(%r14,%rax,8)
movsd  (%r14),%xmm0
addsd  0x8(%r14),%xmm0
movsd  %xmm0,0x8(%r14)
cmp    $0x8,%ecx
jb     1821 <dfst+0x6d1>
mov    %ebp,%eax
add    $0x2,%eax
lea    -0x5(%rbp,%rbp,1),%ecx
movapd 0xc5e(%rip),%xmm2
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
lea    -0x4(%rax),%edx
movsd  (%r14,%rdx,8),%xmm0
lea    -0x3(%rax),%edx
movapd %xmm0,%xmm1
subsd  (%r14,%rdx,8),%xmm1
lea    0x2(%rcx),%esi
movsd  %xmm1,(%r14,%rsi,8)
xorpd  %xmm2,%xmm0
subsd  (%r14,%rdx,8),%xmm0
mov    %ecx,%edx
movsd  %xmm0,(%r14,%rdx,8)
add    $0xfffffffe,%eax
add    $0xfffffffc,%ecx
cmp    $0x5,%eax
jg     13b0 <dfst+0x260>
cmpl   $0x7,0x10(%rsp)
jbe    1821 <dfst+0x6d1>
lea    -0x8(%rbx),%rax
mov    %rax,0x48(%rsp)
lea    0x8(%rbx),%rax
mov    %rax,0x58(%rsp)
lea    -0x10(%rbx),%rax
mov    %rax,0x50(%rsp)
mov    $0x2,%r15d
jmp    1437 <dfst+0x2e7>
nopw   0x0(%rax,%rax,1)
lea    (%r12,%r13,1),%eax
movsd  (%rbx,%rax,8),%xmm0
movsd  %xmm0,(%rbx)
cmp    $0x3,%r13d
jbe    1842 <dfst+0x6f2>
mov    %r12d,%r13d
mov    %r12d,%edi
mov    %rbx,%rsi
mov    0x28(%rsp),%rbp
mov    %ebp,%edx
mov    0x38(%rsp),%r12
mov    %r12,%rcx
call   1050 <dstsub@plt>
cmp    $0x5,%r13d
jb     1490 <dfst+0x340>
mov    %r13d,%edi
mov    %rbx,%rsi
mov    0x18(%rsp),%rdx
mov    0x20(%rsp),%rcx
mov    0x8(%rsp),%r8
call   1070 <cftfsub@plt>
mov    %r13d,%edi
mov    %rbx,%rsi
mov    %ebp,%edx
mov    %r12,%rcx
call   1040 <rftfsub@plt>
jmp    14ba <dfst+0x36a>
cs nopw 0x0(%rax,%rax,1)
cmp    $0x4,%r13d
movapd 0xb64(%rip),%xmm1
jne    14c2 <dfst+0x372>
mov    $0x4,%edi
mov    %rbx,%rsi
mov    0x18(%rsp),%rdx
mov    0x20(%rsp),%rcx
mov    0x8(%rsp),%r8
call   1070 <cftfsub@plt>
movapd 0xb3e(%rip),%xmm1
movsd  0x8(%rbx),%xmm0
subsd  (%rbx),%xmm0
mov    0x10(%rsp),%rax
sub    %r15d,%eax
cltq
movsd  %xmm0,(%r14,%rax,8)
movsd  (%rbx),%xmm0
addsd  0x8(%rbx),%xmm0
movslq %r15d,%rax
movsd  %xmm0,(%r14,%rax,8)
cmp    $0x2,%r13d
jbe    182c <dfst+0x6dc>
lea    0x0(,%r15,4),%edx
mov    %r13d,%ecx
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,%rax,1)
movsd  (%rbx,%rdi,8),%xmm0
xorpd  %xmm1,%xmm0
subsd  0x8(%rbx,%rdi,8),%xmm0
movsd  %xmm0,(%rsi,%rax,8)
movsd  (%rbx,%rdi,8),%xmm0
subsd  0x8(%rbx,%rdi,8),%xmm0
movsd  %xmm0,(%rdx,%rax,8)
add    $0x2,%rdi
add    %rbp,%rax
cmp    %rcx,%rdi
jb     1520 <dfst+0x3d0>
add    %r15d,%r15d
mov    %r13d,%r12d
shr    %r12d
cmp    $0x4,%r13d
jb     1835 <dfst+0x6e5>
cmp    $0x3,%r12d
mov    $0x2,%esi
cmovae %r12d,%esi
lea    -0x1(%rsi),%r8
mov    $0x1,%r11d
cmp    $0x2,%r8
jae    1660 <dfst+0x510>
mov    %r11,%rax
test   $0x1,%sil
jne    15c3 <dfst+0x473>
mov    %rcx,%rax
sub    %r11,%rax
lea    (%rax,%rcx,1),%rdx
mov    %rcx,%rdi
add    %r11,%rdi
movsd  (%rbx,%rdi,8),%xmm0
movsd  (%rbx,%rdx,8),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx,%r11,8)
movsd  (%rbx,%rdx,8),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx,%rax,8)
lea    0x1(%r11),%rax
not    %r11
add    %rsi,%r11
je     1420 <dfst+0x2d0>
lea    (%rax,%rcx,1),%rdx
add    $0x1,%rdx
mov    %rcx,%rdi
sub    %rax,%rdi
shl    $0x4,%rcx
lea    0x0(,%rax,8),%rbp
sub    %rbp,%rcx
nopl   0x0(%rax)
movsd  -0x8(%rbx,%rdx,8),%xmm0
movsd  (%rbx,%rcx,1),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx,%rax,8)
movsd  (%rbx,%rcx,1),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx,%rdi,8)
movsd  (%rbx,%rdx,8),%xmm0
movsd  -0x8(%rbx,%rcx,1),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rbx,%rax,8)
movsd  -0x8(%rbx,%rcx,1),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,-0x8(%rbx,%rdi,8)
add    $0x2,%rax
add    $0x2,%rdx
add    $0xfffffffffffffffe,%rdi
add    $0xfffffffffffffff0,%rcx
cmp    %rax,%rsi
jne    15f0 <dfst+0x4a0>
jmp    1420 <dfst+0x2d0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
lea    -0x2(%rsi),%rax
mov    $0x8,%edx
mul    %rdx
mov    0x48(%rsp),%rdx
lea    (%rdx,%rcx,8),%rbp
seto   %r9b
mov    %rbp,%rdx
sub    %rax,%rdx
cmp    %rbp,%rdx
ja     1584 <dfst+0x434>
test   %r9b,%r9b
jne    1584 <dfst+0x434>
mov    %rcx,%rdx
shl    $0x4,%rdx
mov    0x48(%rsp),%rdi
add    %rdx,%rdi
mov    %rdi,%rbp
sub    %rax,%rbp
cmp    %rdi,%rbp
ja     1584 <dfst+0x434>
test   %r9b,%r9b
jne    1584 <dfst+0x434>
mov    %r8,0x68(%rsp)
mov    %r15,0x40(%rsp)
lea    (%rcx,%rcx,1),%rdi
mov    0x58(%rsp),%r15
cmp    %rbx,%r15
setb   0x4(%rsp)
lea    (%rbx,%rsi,8),%r10
lea    (%r15,%rcx,8),%rax
mov    %rax,0x30(%rsp)
lea    (%rcx,%rsi,1),%rax
lea    (%rbx,%rax,8),%r9
mov    %r9,0x70(%rsp)
mov    %rcx,%rax
sub    %rsi,%rax
lea    (%r15,%rax,8),%r8
lea    (%rbx,%rcx,8),%rax
mov    %rdi,%rbp
mov    %rdi,0x60(%rsp)
sub    %rsi,%rdi
lea    (%r15,%rdi,8),%rbp
add    %rbx,%rdx
cmp    %r9,%r15
setb   0x7(%rsp)
cmp    %r10,0x30(%rsp)
setb   0x30(%rsp)
cmp    %rax,%r15
setb   %dil
cmp    %r10,%r8
setb   0x6(%rsp)
cmp    %rdx,%r15
setb   %r9b
cmp    %r10,%rbp
setb   0x5(%rsp)
cmp    0x70(%rsp),%r8
setb   %r15b
cmp    %rdx,%r8
setb   %r10b
cmp    %rax,%rbp
setb   %bpl
mov    0x30(%rsp),%dl
test   %dl,0x7(%rsp)
jne    1817 <dfst+0x6c7>
mov    %r15d,%r8d
and    0x6(%rsp),%dil
mov    0x40(%rsp),%r15
jne    1584 <dfst+0x434>
and    0x5(%rsp),%r9b
jne    1584 <dfst+0x434>
and    %r8b,0x4(%rsp)
jne    1584 <dfst+0x434>
and    %bpl,%r10b
jne    1584 <dfst+0x434>
mov    0x68(%rsp),%r10
mov    %r10,%r8
and    $0xfffffffffffffffe,%r8
mov    %r10,%r11
or     $0x1,%r11
mov    0x50(%rsp),%rdx
mov    0x60(%rsp),%rdi
lea    (%rdx,%rdi,8),%rdx
mov    %r8,%r9
neg    %r9
mov    $0x8,%edi
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movupd (%rdx,%rbp,8),%xmm0
shufpd $0x1,%xmm0,%xmm0
movupd (%rax,%rdi,1),%xmm1
addpd  %xmm1,%xmm0
movupd %xmm0,(%rbx,%rdi,1)
movupd (%rdx,%rbp,8),%xmm0
shufpd $0x1,%xmm1,%xmm1
subpd  %xmm1,%xmm0
movupd %xmm0,-0x10(%rax,%rbp,8)
add    $0xfffffffffffffffe,%rbp
add    $0x10,%rdi
cmp    %rbp,%r9
jne    17d0 <dfst+0x680>
cmp    %r8,%r10
je     1420 <dfst+0x2d0>
jmp    1584 <dfst+0x434>
mov    0x40(%rsp),%r15
jmp    1584 <dfst+0x434>
movsd  (%rbx),%xmm0
mov    $0x2,%eax
jmp    1845 <dfst+0x6f5>
add    %r15d,%r15d
mov    %r13d,%r12d
shr    %r12d
add    %r13d,%r12d
movsd  (%rbx,%r12,8),%xmm0
movsd  %xmm0,(%rbx)
movslq %r15d,%rax
movsd  %xmm0,(%r14,%rax,8)
movq   $0x0,(%r14)
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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