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 + 2, nw, w);
      rftfsub(m, a, nc, w + nw);
    } else if (m == 4) {
      cftfsub(m, a, ip + 2, 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 + 2, nw, w);
        rftfsub(m, t, nc, w + nw);
      } else if (m == 4) {
        cftfsub(m, t, ip + 2, 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    173f <dfst+0x5ef>
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    13cb <dfst+0x27b>
mov    -0x38(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x20(%rbp),%rdx
add    $0x8,%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    13f5 <dfst+0x2a5>
cmpl   $0x4,-0x38(%rbp)
jne    13f0 <dfst+0x2a0>
mov    -0x38(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x20(%rbp),%rdx
add    $0x8,%rdx
mov    -0x40(%rbp),%ecx
mov    -0x28(%rbp),%r8
call   1070 <cftfsub@plt>
jmp    13f5 <dfst+0x2a5>
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     14d2 <dfst+0x382>
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    143b <dfst+0x2eb>
movl   $0x2,-0x34(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x38(%rbp)
cmpl   $0x2,-0x38(%rbp)
jl     172a <dfst+0x5da>
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    154f <dfst+0x3ff>
mov    -0x38(%rbp),%edi
mov    -0x18(%rbp),%rsi
mov    -0x20(%rbp),%rdx
add    $0x8,%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    1579 <dfst+0x429>
cmpl   $0x4,-0x38(%rbp)
jne    1574 <dfst+0x424>
mov    -0x38(%rbp),%edi
mov    -0x18(%rbp),%rsi
mov    -0x20(%rbp),%rdx
add    $0x8,%rdx
mov    -0x40(%rbp),%ecx
mov    -0x28(%rbp),%r8
call   1070 <cftfsub@plt>
jmp    1579 <dfst+0x429>
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    1667 <dfst+0x517>
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    15c8 <dfst+0x478>
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    1705 <dfst+0x5b5>
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    1680 <dfst+0x530>
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    14df <dfst+0x38f>
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,0x10(%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,0x8(%rsp)
mov    %r15d,%edi
shr    %edi
mov    %r15d,%r12d
shr    $0x2,%r12d
mov    %edi,%r8d
mov    %r15,0x18(%rsp)
cmp    $0x8,%r15d
jb     12a6 <dfst+0x156>
mov    0x18(%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    0x18(%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    0x10(%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     1341 <dfst+0x1f1>
mov    0x8(%rsp),%rax
lea    0x8(%rax),%rdx
mov    %ebp,%edi
mov    %r14,%rsi
mov    0x20(%rsp),%rcx
mov    0x10(%rsp),%r8
call   1070 <cftfsub@plt>
mov    %ebp,%edi
mov    %r14,%rsi
mov    0x28(%rsp),%rdx
mov    %r13,%rcx
call   1040 <rftfsub@plt>
jmp    1366 <dfst+0x216>
cmp    $0x4,%ebp
mov    0x20(%rsp),%rcx
jne    1366 <dfst+0x216>
mov    0x8(%rsp),%rax
lea    0x8(%rax),%rdx
mov    $0x4,%edi
mov    %r14,%rsi
mov    0x10(%rsp),%r8
call   1070 <cftfsub@plt>
mov    0x18(%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 0xc56(%rip),%xmm2
nopw   0x0(%rax,%rax,1)
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,0x18(%rsp)
jbe    1821 <dfst+0x6d1>
addq   $0x8,0x8(%rsp)
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>
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    0x8(%rsp),%rdx
mov    0x20(%rsp),%rcx
mov    0x10(%rsp),%r8
call   1070 <cftfsub@plt>
mov    %r13d,%edi
mov    %rbx,%rsi
mov    %ebp,%edx
mov    %r12,%rcx
call   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    0x8(%rsp),%rdx
mov    0x20(%rsp),%rcx
mov    0x10(%rsp),%r8
call   1070 <cftfsub@plt>
movapd 0xb3e(%rip),%xmm1
movsd  0x8(%rbx),%xmm0
subsd  (%rbx),%xmm0
mov    0x18(%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 %472
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 %206
192: ; preds = %147
%193 = load i32, i32* %14
%194 = load double*, double** %7
%195 = load i32*, i32** %9
%196 = getelementptr inbounds i32, i32* %195, i64 2
%197 = load i32, i32* %16
%198 = load double*, double** %10
%199 = load i32, i32* %14
%200 = load double*, double** %7
%201 = load i32, i32* %17
%202 = load double*, double** %10
%203 = load i32, i32* %16
%204 = sext i32 %203 to i64
%205 = getelementptr inbounds double, double* %202, i64 %204
br label %217
206: ; preds = %147
%207 = load i32, i32* %14
%208 = icmp eq i32 %207, 4
br i1 %208, label %209, label %216
209: ; preds = %206
%210 = load i32, i32* %14
%211 = load double*, double** %7
%212 = load i32*, i32** %9
%213 = getelementptr inbounds i32, i32* %212, i64 2
%214 = load i32, i32* %16
%215 = load double*, double** %10
br label %216
216: ; preds = %209, %206
br label %217
217: ; preds = %216, %192
%218 = load double*, double** %7
%219 = getelementptr inbounds double, double* %218, i64 1
%220 = load double, double* %219
%221 = load double*, double** %7
%222 = getelementptr inbounds double, double* %221, i64 0
%223 = load double, double* %222
%224 = fsub double %220, %223
%225 = load double*, double** %7
%226 = load i32, i32* %6
%227 = sub nsw i32 %226, 1
%228 = sext i32 %227 to i64
%229 = getelementptr inbounds double, double* %225, i64 %228
store double %224, double* %229
%230 = load double*, double** %7
%231 = getelementptr inbounds double, double* %230, i64 0
%232 = load double, double* %231
%233 = load double*, double** %7
%234 = getelementptr inbounds double, double* %233, i64 1
%235 = load double, double* %234
%236 = fadd double %232, %235
%237 = load double*, double** %7
%238 = getelementptr inbounds double, double* %237, i64 1
store double %236, double* %238
%239 = load i32, i32* %14
%240 = sub nsw i32 %239, 2
store i32 %240, i32* %11
br label %241
241: ; preds = %282, %217
%242 = load i32, i32* %11
%243 = icmp sge i32 %242, 2
br i1 %243, label %244, label %285
244: ; preds = %241
%245 = load double*, double** %7
%246 = load i32, i32* %11
%247 = sext i32 %246 to i64
%248 = getelementptr inbounds double, double* %245, i64 %247
%249 = load double, double* %248
%250 = load double*, double** %7
%251 = load i32, i32* %11
%252 = add nsw i32 %251, 1
%253 = sext i32 %252 to i64
%254 = getelementptr inbounds double, double* %250, i64 %253
%255 = load double, double* %254
%256 = fsub double %249, %255
%257 = load double*, double** %7
%258 = load i32, i32* %11
%259 = mul nsw i32 2, %258
%260 = add nsw i32 %259, 1
%261 = sext i32 %260 to i64
%262 = getelementptr inbounds double, double* %257, i64 %261
store double %256, double* %262
%263 = load double*, double** %7
%264 = load i32, i32* %11
%265 = sext i32 %264 to i64
%266 = getelementptr inbounds double, double* %263, i64 %265
%267 = load double, double* %266
%268 = fneg double %267
%269 = load double*, double** %7
%270 = load i32, i32* %11
%271 = add nsw i32 %270, 1
%272 = sext i32 %271 to i64
%273 = getelementptr inbounds double, double* %269, i64 %272
%274 = load double, double* %273
%275 = fsub double %268, %274
%276 = load double*, double** %7
%277 = load i32, i32* %11
%278 = mul nsw i32 2, %277
%279 = sub nsw i32 %278, 1
%280 = sext i32 %279 to i64
%281 = getelementptr inbounds double, double* %276, i64 %280
store double %275, double* %281
br label %282
282: ; preds = %244
%283 = load i32, i32* %11
%284 = sub nsw i32 %283, 2
store i32 %284, i32* %11
br label %241
285: ; preds = %241
store i32 2, i32* %13
%286 = load i32, i32* %15
store i32 %286, i32* %14
br label %287
287: ; preds = %453, %285
%288 = load i32, i32* %14
%289 = icmp sge i32 %288, 2
br i1 %289, label %290, label %464
290: ; preds = %287
%291 = load i32, i32* %14
%292 = load double*, double** %8
%293 = load i32, i32* %17
%294 = load double*, double** %10
%295 = load i32, i32* %16
%296 = sext i32 %295 to i64
%297 = getelementptr inbounds double, double* %294, i64 %296
%298 = load i32, i32* %14
%299 = icmp sgt i32 %298, 4
br i1 %299, label %300, label %314
300: ; preds = %290
%301 = load i32, i32* %14
%302 = load double*, double** %8
%303 = load i32*, i32** %9
%304 = getelementptr inbounds i32, i32* %303, i64 2
%305 = load i32, i32* %16
%306 = load double*, double** %10
%307 = load i32, i32* %14
%308 = load double*, double** %8
%309 = load i32, i32* %17
%310 = load double*, double** %10
%311 = load i32, i32* %16
%312 = sext i32 %311 to i64
%313 = getelementptr inbounds double, double* %310, i64 %312
br label %325
314: ; preds = %290
%315 = load i32, i32* %14
%316 = icmp eq i32 %315, 4
br i1 %316, label %317, label %324
317: ; preds = %314
%318 = load i32, i32* %14
%319 = load double*, double** %8
%320 = load i32*, i32** %9
%321 = getelementptr inbounds i32, i32* %320, i64 2
%322 = load i32, i32* %16
%323 = load double*, double** %10
br label %324
324: ; preds = %317, %314
br label %325
325: ; preds = %324, %300
%326 = load double*, double** %8
%327 = getelementptr inbounds double, double* %326, i64 1
%328 = load double, double* %327
%329 = load double*, double** %8
%330 = getelementptr inbounds double, double* %329, i64 0
%331 = load double, double* %330
%332 = fsub double %328, %331
%333 = load double*, double** %7
%334 = load i32, i32* %6
%335 = load i32, i32* %13
%336 = sub nsw i32 %334, %335
%337 = sext i32 %336 to i64
%338 = getelementptr inbounds double, double* %333, i64 %337
store double %332, double* %338
%339 = load double*, double** %8
%340 = getelementptr inbounds double, double* %339, i64 0
%341 = load double, double* %340
%342 = load double*, double** %8
%343 = getelementptr inbounds double, double* %342, i64 1
%344 = load double, double* %343
%345 = fadd double %341, %344
%346 = load double*, double** %7
%347 = load i32, i32* %13
%348 = sext i32 %347 to i64
%349 = getelementptr inbounds double, double* %346, i64 %348
store double %345, double* %349
store i32 0, i32* %12
store i32 2, i32* %11
br label %350
350: ; preds = %396, %325
%351 = load i32, i32* %11
%352 = load i32, i32* %14
%353 = icmp slt i32 %351, %352
br i1 %353, label %354, label %399
354: ; preds = %350
%355 = load i32, i32* %13
%356 = shl i32 %355, 2
%357 = load i32, i32* %12
%358 = add nsw i32 %357, %356
store i32 %358, i32* %12
%359 = load double*, double** %8
%360 = load i32, i32* %11
%361 = sext i32 %360 to i64
%362 = getelementptr inbounds double, double* %359, i64 %361
%363 = load double, double* %362
%364 = fneg double %363
%365 = load double*, double** %8
%366 = load i32, i32* %11
%367 = add nsw i32 %366, 1
%368 = sext i32 %367 to i64
%369 = getelementptr inbounds double, double* %365, i64 %368
%370 = load double, double* %369
%371 = fsub double %364, %370
%372 = load double*, double** %7
%373 = load i32, i32* %12
%374 = load i32, i32* %13
%375 = sub nsw i32 %373, %374
%376 = sext i32 %375 to i64
%377 = getelementptr inbounds double, double* %372, i64 %376
store double %371, double* %377
%378 = load double*, double** %8
%379 = load i32, i32* %11
%380 = sext i32 %379 to i64
%381 = getelementptr inbounds double, double* %378, i64 %380
%382 = load double, double* %381
%383 = load double*, double** %8
%384 = load i32, i32* %11
%385 = add nsw i32 %384, 1
%386 = sext i32 %385 to i64
%387 = getelementptr inbounds double, double* %383, i64 %386
%388 = load double, double* %387
%389 = fsub double %382, %388
%390 = load double*, double** %7
%391 = load i32, i32* %12
%392 = load i32, i32* %13
%393 = add nsw i32 %391, %392
%394 = sext i32 %393 to i64
%395 = getelementptr inbounds double, double* %390, i64 %394
store double %389, double* %395
br label %396
396: ; preds = %354
%397 = load i32, i32* %11
%398 = add nsw i32 %397, 2
store i32 %398, i32* %11
br label %350
399: ; preds = %350
%400 = load i32, i32* %13
%401 = shl i32 %400, 1
store i32 %401, i32* %13
%402 = load i32, i32* %14
%403 = ashr i32 %402, 1
store i32 %403, i32* %15
store i32 1, i32* %11
br label %404
404: ; preds = %450, %399
%405 = load i32, i32* %11
%406 = load i32, i32* %15
%407 = icmp slt i32 %405, %406
br i1 %407, label %408, label %453
408: ; preds = %404
%409 = load i32, i32* %14
%410 = load i32, i32* %11
%411 = sub nsw i32 %409, %410
store i32 %411, i32* %12
%412 = load double*, double** %8
%413 = load i32, i32* %14
%414 = load i32, i32* %12
%415 = add nsw i32 %413, %414
%416 = sext i32 %415 to i64
%417 = getelementptr inbounds double, double* %412, i64 %416
%418 = load double, double* %417
%419 = load double*, double** %8
%420 = load i32, i32* %14
%421 = load i32, i32* %11
%422 = add nsw i32 %420, %421
%423 = sext i32 %422 to i64
%424 = getelementptr inbounds double, double* %419, i64 %423
%425 = load double, double* %424
%426 = fadd double %418, %425
%427 = load double*, double** %8
%428 = load i32, i32* %11
%429 = sext i32 %428 to i64
%430 = getelementptr inbounds double, double* %427, i64 %429
store double %426, double* %430
%431 = load double*, double** %8
%432 = load i32, i32* %14
%433 = load i32, i32* %12
%434 = add nsw i32 %432, %433
%435 = sext i32 %434 to i64
%436 = getelementptr inbounds double, double* %431, i64 %435
%437 = load double, double* %436
%438 = load double*, double** %8
%439 = load i32, i32* %14
%440 = load i32, i32* %11
%441 = add nsw i32 %439, %440
%442 = sext i32 %441 to i64
%443 = getelementptr inbounds double, double* %438, i64 %442
%444 = load double, double* %443
%445 = fsub double %437, %444
%446 = load double*, double** %8
%447 = load i32, i32* %12
%448 = sext i32 %447 to i64
%449 = getelementptr inbounds double, double* %446, i64 %448
store double %445, double* %449
br label %450
450: ; preds = %408
%451 = load i32, i32* %11
%452 = add nsw i32 %451, 1
store i32 %452, i32* %11
br label %404
453: ; preds = %404
%454 = load double*, double** %8
%455 = load i32, i32* %14
%456 = load i32, i32* %15
%457 = add nsw i32 %455, %456
%458 = sext i32 %457 to i64
%459 = getelementptr inbounds double, double* %454, i64 %458
%460 = load double, double* %459
%461 = load double*, double** %8
%462 = getelementptr inbounds double, double* %461, i64 0
store double %460, double* %462
%463 = load i32, i32* %15
store i32 %463, i32* %14
br label %287
464: ; preds = %287
%465 = load double*, double** %8
%466 = getelementptr inbounds double, double* %465, i64 0
%467 = load double, double* %466
%468 = load double*, double** %7
%469 = load i32, i32* %13
%470 = sext i32 %469 to i64
%471 = getelementptr inbounds double, double* %468, i64 %470
store double %467, double* %471
br label %472
472: ; preds = %464, %52
%473 = load double*, double** %7
%474 = getelementptr inbounds double, double* %473, i64 0
store double 0.000000e+00, double* %474
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 %340
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 %76
74: ; preds = %56
%75 = getelementptr inbounds i32, i32* %3, i64 2
tail call void @cftfsub(i32 %25, double* %1, i32* %75, i32 %12, double* %4)
tail call void @rftfsub(i32 %25, double* %1, i32 %22, double* %72)
br label %80
76: ; preds = %56
%77 = icmp eq i32 %25, 4
br i1 %77, label %78, label %80
78: ; preds = %76
%79 = getelementptr inbounds i32, i32* %3, i64 2
tail call void @cftfsub(i32 4, double* %1, i32* %79, i32 %12, double* %4)
br label %80
80: ; preds = %76, %78, %74
%81 = getelementptr inbounds double, double* %1, i64 1
%82 = load double, double* %81
%83 = load double, double* %1
%84 = fsub double %82, %83
%85 = add nsw i32 %0, -1
%86 = zext i32 %85 to i64
%87 = getelementptr inbounds double, double* %1, i64 %86
store double %84, double* %87
%88 = load double, double* %1
%89 = load double, double* %81
%90 = fadd double %88, %89
store double %90, double* %81
br i1 %27, label %106, label %92
91: ; preds = %106
br i1 %27, label %94, label %92
92: ; preds = %80, %91
%93 = load double, double* %2
br label %336
94: ; preds = %91
%95 = getelementptr inbounds i32, i32* %3, i64 2
%96 = getelementptr inbounds 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 = getelementptr double, double* %2, i64 1
%101 = getelementptr double, double* %2, i64 1
%102 = getelementptr double, double* %2, i64 1
%103 = icmp ult double* %100, %2
%104 = getelementptr inbounds double, double* %2, i64 -1
%105 = getelementptr inbounds double, double* %2, i64 -1
br label %128
106: ; preds = %80, %106
%107 = phi i32 [ %108, %106 ], [ %25, %80 ]
%108 = add nsw i32 %107, -2
%109 = zext i32 %108 to i64
%110 = getelementptr inbounds double, double* %1, i64 %109
%111 = load double, double* %110
%112 = add nsw i32 %107, -1
%113 = zext i32 %112 to i64
%114 = getelementptr inbounds double, double* %1, i64 %113
%115 = load double, double* %114
%116 = fsub double %111, %115
%117 = shl nuw nsw i32 %108, 1
%118 = or i32 %117, 1
%119 = zext i32 %118 to i64
%120 = getelementptr inbounds double, double* %1, i64 %119
store double %116, double* %120
%121 = fneg double %111
%122 = load double, double* %114
%123 = fsub double %121, %122
%124 = add nsw i32 %117, -1
%125 = zext i32 %124 to i64
%126 = getelementptr inbounds double, double* %1, i64 %125
store double %123, double* %126
%127 = icmp sgt i32 %107, 5
br i1 %127, label %106, label %91
128: ; preds = %94, %320
%129 = phi i32 [ 2, %94 ], [ %177, %320 ]
%130 = phi i32 [ %26, %94 ], [ %178, %320 ]
tail call void @dstsub(i32 %130, double* %2, i32 %22, double* %72)
%131 = icmp ugt i32 %130, 4
br i1 %131, label %132, label %133
132: ; preds = %128
tail call void @cftfsub(i32 %130, double* %2, i32* %95, i32 %12, double* %4)
tail call void @rftfsub(i32 %130, double* %2, i32 %22, double* %72)
br label %136
133: ; preds = %128
%134 = icmp eq i32 %130, 4
br i1 %134, label %135, label %136
135: ; preds = %133
tail call void @cftfsub(i32 4, double* %2, i32* %95, i32 %12, double* %4)
br label %136
136: ; preds = %133, %135, %132
%137 = load double, double* %96
%138 = load double, double* %2
%139 = fsub double %137, %138
%140 = sub nsw i32 %0, %129
%141 = sext i32 %140 to i64
%142 = getelementptr inbounds double, double* %1, i64 %141
store double %139, double* %142
%143 = load double, double* %2
%144 = load double, double* %96
%145 = fadd double %143, %144
%146 = sext i32 %129 to i64
%147 = getelementptr inbounds double, double* %1, i64 %146
store double %145, double* %147
%148 = icmp ugt i32 %130, 2
br i1 %148, label %152, label %149
149: ; preds = %136
%150 = shl i32 %129, 1
%151 = lshr i32 %130, 1
br label %325
152: ; preds = %136
%153 = shl i32 %129, 2
%154 = zext i32 %130 to i64
%155 = sext i32 %153 to i64
br label %156
156: ; preds = %152, %156
%157 = phi i64 [ 0, %152 ], [ %159, %156 ]
%158 = phi i64 [ 2, %152 ], [ %174, %156 ]
%159 = add nsw i64 %157, %155
%160 = getelementptr inbounds double, double* %2, i64 %158
%161 = load double, double* %160
%162 = fneg double %161
%163 = or i64 %158, 1
%164 = getelementptr inbounds double, double* %2, i64 %163
%165 = load double, double* %164
%166 = fsub double %162, %165
%167 = sub nsw i64 %159, %146
%168 = getelementptr inbounds double, double* %1, i64 %167
store double %166, double* %168
%169 = load double, double* %160
%170 = load double, double* %164
%171 = fsub double %169, %170
%172 = add nsw i64 %159, %146
%173 = getelementptr inbounds double, double* %1, i64 %172
store double %171, double* %173
%174 = add nuw nsw i64 %158, 2
%175 = icmp ult i64 %174, %154
br i1 %175, label %156, label %176
176: ; preds = %156
%177 = shl i32 %129, 1
%178 = lshr i32 %130, 1
%179 = icmp ugt i32 %130, 3
br i1 %179, label %180, label %325
180: ; preds = %176
%181 = zext i32 %130 to i64
%182 = call i32 @llvm.umax.i32(i32 %178, i32 2)
%183 = zext i32 %182 to i64
%184 = add nsw i64 %183, -1
%185 = icmp ult i64 %184, 2
br i1 %185, label %268, label %186
186: ; preds = %180
%187 = add nsw i64 %183, -2
%188 = getelementptr double, double* %97, i64 %154
%189 = call { i64, i1 } @llvm.umul.with.overflow.i64(i64 %187, i64 8)
%190 = extractvalue { i64, i1 } %189, 0
%191 = extractvalue { i64, i1 } %189, 1
%192 = bitcast double* %188 to i8*
%193 = sub i64 0, %190
%194 = getelementptr i8, i8* %192, i64 %193
%195 = icmp ugt i8* %194, %192
%196 = or i1 %195, %191
%197 = shl nuw nsw i64 %154, 1
%198 = getelementptr double, double* %98, i64 %197
%199 = call { i64, i1 } @llvm.umul.with.overflow.i64(i64 %187, i64 8)
%200 = extractvalue { i64, i1 } %199, 0
%201 = extractvalue { i64, i1 } %199, 1
%202 = bitcast double* %198 to i8*
%203 = sub i64 0, %200
%204 = getelementptr i8, i8* %202, i64 %203
%205 = icmp ugt i8* %204, %202
%206 = or i1 %205, %201
%207 = or i1 %196, %206
br i1 %207, label %268, label %208
208: ; preds = %186
%209 = getelementptr double, double* %2, i64 %183
%210 = getelementptr double, double* %100, i64 %154
%211 = add nuw nsw i64 %154, %183
%212 = getelementptr double, double* %2, i64 %211
%213 = sub nsw i64 %154, %183
%214 = getelementptr double, double* %101, i64 %213
%215 = getelementptr double, double* %2, i64 %154
%216 = shl nuw nsw i64 %154, 1
%217 = sub nsw i64 %216, %183
%218 = getelementptr double, double* %102, i64 %217
%219 = getelementptr double, double* %2, i64 %216
%220 = icmp ult double* %99, %212
%221 = icmp ult double* %210, %209
%222 = and i1 %220, %221
%223 = icmp ult double* %99, %215
%224 = icmp ult double* %214, %209
%225 = and i1 %223, %224
%226 = or i1 %222, %225
%227 = icmp ult double* %99, %219
%228 = icmp ult double* %218, %209
%229 = and i1 %227, %228
%230 = or i1 %226, %229
%231 = icmp ult double* %214, %212
%232 = and i1 %103, %231
%233 = or i1 %230, %232
%234 = icmp ult double* %214, %219
%235 = icmp ult double* %218, %215
%236 = and i1 %234, %235
%237 = or i1 %233, %236
br i1 %237, label %268, label %238
238: ; preds = %208
%239 = and i64 %184, -2
%240 = or i64 %184, 1
br label %241
241: ; preds = %241, %238
%242 = phi i64 [ 0, %238 ], [ %264, %241 ]
%243 = or i64 %242, 1
%244 = sub nsw i64 %181, %243
%245 = add nsw i64 %244, %181
%246 = getelementptr inbounds double, double* %104, i64 %245
%247 = bitcast double* %246 to <2 x double>*
%248 = load <2 x double>, <2 x double>* %247
%249 = shufflevector <2 x double> %248, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%250 = add nuw nsw i64 %243, %181
%251 = getelementptr inbounds double, double* %2, i64 %250
%252 = bitcast double* %251 to <2 x double>*
%253 = load <2 x double>, <2 x double>* %252
%254 = fadd <2 x double> %249, %253
%255 = getelementptr inbounds double, double* %2, i64 %243
%256 = bitcast double* %255 to <2 x double>*
store <2 x double> %254, <2 x double>* %256
%257 = bitcast double* %246 to <2 x double>*
%258 = load <2 x double>, <2 x double>* %257
%259 = shufflevector <2 x double> %258, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%260 = fsub <2 x double> %259, %253
%261 = shufflevector <2 x double> %260, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%262 = getelementptr inbounds double, double* %105, i64 %244
%263 = bitcast double* %262 to <2 x double>*
store <2 x double> %261, <2 x double>* %263
%264 = add nuw i64 %242, 2
%265 = icmp eq i64 %264, %239
br i1 %265, label %266, label %241
266: ; preds = %241
%267 = icmp eq i64 %184, %239
br i1 %267, label %320, label %268
268: ; preds = %208, %186, %180, %266
%269 = phi i64 [ 1, %208 ], [ 1, %186 ], [ 1, %180 ], [ %240, %266 ]
%270 = xor i64 %269, -1
%271 = and i64 %183, 1
%272 = icmp eq i64 %271, 0
br i1 %272, label %273, label %287
273: ; preds = %268
%274 = sub nsw i64 %181, %269
%275 = add nsw i64 %274, %181
%276 = getelementptr inbounds double, double* %2, i64 %275
%277 = load double, double* %276
%278 = add nuw nsw i64 %269, %181
%279 = getelementptr inbounds double, double* %2, i64 %278
%280 = load double, double* %279
%281 = fadd double %277, %280
%282 = getelementptr inbounds double, double* %2, i64 %269
store double %281, double* %282
%283 = load double, double* %276
%284 = fsub double %283, %280
%285 = getelementptr inbounds double, double* %2, i64 %274
store double %284, double* %285
%286 = add nuw nsw i64 %269, 1
br label %287
287: ; preds = %273, %268
%288 = phi i64 [ %269, %268 ], [ %286, %273 ]
%289 = sub nsw i64 0, %183
%290 = icmp eq i64 %270, %289
br i1 %290, label %320, label %291
291: ; preds = %287, %291
%292 = phi i64 [ %318, %291 ], [ %288, %287 ]
%293 = sub nsw i64 %181, %292
%294 = add nsw i64 %293, %181
%295 = getelementptr inbounds double, double* %2, i64 %294
%296 = load double, double* %295
%297 = add nuw nsw i64 %292, %181
%298 = getelementptr inbounds double, double* %2, i64 %297
%299 = load double, double* %298
%300 = fadd double %296, %299
%301 = getelementptr inbounds double, double* %2, i64 %292
store double %300, double* %301
%302 = load double, double* %295
%303 = fsub double %302, %299
%304 = getelementptr inbounds double, double* %2, i64 %293
store double %303, double* %304
%305 = add nuw nsw i64 %292, 1
%306 = sub nsw i64 %181, %305
%307 = add nsw i64 %306, %181
%308 = getelementptr inbounds double, double* %2, i64 %307
%309 = load double, double* %308
%310 = add nuw nsw i64 %305, %181
%311 = getelementptr inbounds double, double* %2, i64 %310
%312 = load double, double* %311
%313 = fadd double %309, %312
%314 = getelementptr inbounds double, double* %2, i64 %305
store double %313, double* %314
%315 = load double, double* %308
%316 = fsub double %315, %312
%317 = getelementptr inbounds double, double* %2, i64 %306
store double %316, double* %317
%318 = add nuw nsw i64 %292, 2
%319 = icmp eq i64 %318, %183
br i1 %319, label %320, label %291
320: ; preds = %287, %291, %266
%321 = add nuw nsw i32 %178, %130
%322 = zext i32 %321 to i64
%323 = getelementptr inbounds double, double* %2, i64 %322
%324 = load double, double* %323
store double %324, double* %2
br i1 %179, label %128, label %332
325: ; preds = %176, %149
%326 = phi i32 [ %151, %149 ], [ %178, %176 ]
%327 = phi i32 [ %150, %149 ], [ %177, %176 ]
%328 = add nuw nsw i32 %326, %130
%329 = zext i32 %328 to i64
%330 = getelementptr inbounds double, double* %2, i64 %329
%331 = load double, double* %330
store double %331, double* %2
br label %332
332: ; preds = %320, %325
%333 = phi double [ %331, %325 ], [ %324, %320 ]
%334 = phi i32 [ %327, %325 ], [ %177, %320 ]
%335 = sext i32 %334 to i64
br label %336
336: ; preds = %92, %332
%337 = phi double [ %93, %92 ], [ %333, %332 ]
%338 = phi i64 [ 2, %92 ], [ %335, %332 ]
%339 = getelementptr inbounds double, double* %1, i64 %338
store double %337, double* %339
br label %340
340: ; preds = %336, %21
store double 0.000000e+00, double* %1
ret void
}