void dfst(int n, float *a, float *t, int *ip, float *w) {
  void makewt(int nw, int *ip, float *w);
  void makect(int nc, int *ip, float *c);
  void cftfsub(int n, float *a, int *ip, int nw, float *w);
  void rftfsub(int n, float *a, int nc, float *c);
  void dstsub(int n, float *a, int nc, float *c);
  int j, k, l, m, mh, nw, nc;
  float 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    $0x60,%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    $0x2,%rax
add    %rax,%rdx
call   1060 <makect@plt>
cmpl   $0x2,-0x4(%rbp)
jle    171b <dfst+0x5cb>
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
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x4c(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    -0x30(%rbp),%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    -0x30(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x54(%rbp)
movss  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x50(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x4c(%rbp),%xmm0
addss  -0x54(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x4c(%rbp),%xmm0
subss  -0x54(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1200 <dfst+0xb0>
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x38(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x44(%rbp),%edx
mov    -0x28(%rbp),%rcx
movslq -0x40(%rbp),%rax
shl    $0x2,%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    $0x2,%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
movss  0x4(%rax),%xmm0
mov    -0x10(%rbp),%rax
subss  (%rax),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x10(%rbp),%rax
addss  0x4(%rax),%xmm0
mov    -0x10(%rbp),%rax
movss  %xmm0,0x4(%rax)
mov    -0x38(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x2,-0x2c(%rbp)
jl     14c0 <dfst+0x370>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
shl    $0x1,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
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     1706 <dfst+0x5b6>
mov    -0x38(%rbp),%edi
mov    -0x18(%rbp),%rsi
mov    -0x44(%rbp),%edx
mov    -0x28(%rbp),%rcx
movslq -0x40(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rcx
call   1050 <dstsub@plt>
cmpl   $0x4,-0x38(%rbp)
jle    1539 <dfst+0x3e9>
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    $0x2,%rax
add    %rax,%rcx
call   1040 <rftfsub@plt>
jmp    155f <dfst+0x40f>
cmpl   $0x4,-0x38(%rbp)
jne    155a <dfst+0x40a>
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    155f <dfst+0x40f>
mov    -0x18(%rbp),%rax
movss  0x4(%rax),%xmm0
mov    -0x18(%rbp),%rax
subss  (%rax),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x18(%rbp),%rax
addss  0x4(%rax),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movl   $0x0,-0x30(%rbp)
movl   $0x2,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x38(%rbp),%eax
jge    1643 <dfst+0x4f3>
mov    -0x34(%rbp),%eax
shl    $0x2,%eax
add    -0x30(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x2c(%rbp)
jmp    15ae <dfst+0x45e>
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    16e1 <dfst+0x591>
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
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    165c <dfst+0x50c>
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x38(%rbp)
jmp    14cd <dfst+0x37d>
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax)
add    $0x60,%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,0x30(%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,4),%rdx
mov    %rdi,0x30(%rsp)
mov    %r13,%rsi
mov    %rcx,%r12
call   1060 <makect@plt>
mov    %r12,%rcx
cmp    $0x3,%r15d
jl     1808 <dfst+0x6b8>
mov    %rcx,0x48(%rsp)
mov    %rbp,0x28(%rsp)
mov    %r13,0x20(%rsp)
mov    %r15d,%edi
shr    %edi
mov    %r15d,%r12d
shr    $0x2,%r12d
mov    %edi,%r8d
mov    %r15,0x18(%rsp)
cmp    $0x8,%r15d
jb     12a3 <dfst+0x153>
mov    0x18(%rsp),%eax
cmp    $0x3,%r12d
mov    $0x2,%ebp
cmovae %r12d,%ebp
lea    (%r14,%rax,4),%r11
mov    $0x1,%edx
mov    $0x1,%esi
sub    %r8,%rsi
lea    (%r14,%r8,4),%rcx
add    $0xfffffffffffffffc,%rcx
lea    (%rbx,%r8,4),%r9
add    $0xfffffffffffffffc,%r9
lea    -0x4(,%rbp,4),%r10
ff
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  (%r14,%rdx,4),%xmm0
movss  -0x4(%r11,%rax,1),%xmm1
movaps %xmm0,%xmm2
addss  %xmm1,%xmm2
subss  %xmm1,%xmm0
movss  (%rcx,%rax,1),%xmm1
movss  (%r11,%rsi,4),%xmm3
movaps %xmm1,%xmm4
addss  %xmm3,%xmm4
subss  %xmm3,%xmm1
movss  %xmm2,(%r14,%rdx,4)
movss  %xmm4,(%rcx,%rax,1)
movaps %xmm0,%xmm2
addss  %xmm1,%xmm2
movss  %xmm2,(%rbx,%rdx,4)
subss  %xmm1,%xmm0
movss  %xmm0,(%r9,%rax,1)
add    $0x1,%rdx
add    $0xfffffffffffffffc,%rax
add    $0x1,%rsi
mov    %r10,%rbp
add    %rax,%rbp
jne    1240 <dfst+0xf0>
mov    %r12d,%eax
movss  (%r14,%rax,4),%xmm0
mov    0x18(%rsp),%r15
mov    %r15d,%ecx
sub    %r12d,%ecx
movslq %ecx,%rcx
subss  (%r14,%rcx,4),%xmm0
movss  %xmm0,(%rbx)
movss  (%r14,%rcx,4),%xmm0
addss  (%r14,%rax,4),%xmm0
movss  %xmm0,(%r14,%rax,4)
movss  (%r14,%r8,4),%xmm0
movss  %xmm0,(%r14)
mov    0x48(%rsp),%rax
mov    0x10(%rsp),%rcx
lea    (%rcx,%rax,4),%r13
mov    %r14,%rsi
mov    0x30(%rsp),%rdx
mov    %r13,%rcx
mov    %rdi,%rbp
call   1050 <dstsub@plt>
cmp    $0xa,%r15d
mov    %r13,0x48(%rsp)
jl     133a <dfst+0x1ea>
mov    %ebp,%edi
mov    %r14,%rsi
mov    0x20(%rsp),%rdx
mov    0x28(%rsp),%rcx
mov    0x10(%rsp),%r8
call   1070 <cftfsub@plt>
mov    %ebp,%edi
mov    %r14,%rsi
mov    0x30(%rsp),%rdx
mov    %r13,%rcx
call   1040 <rftfsub@plt>
jmp    135b <dfst+0x20b>
cmp    $0x4,%ebp
mov    0x28(%rsp),%rcx
jne    135b <dfst+0x20b>
mov    $0x4,%edi
mov    %r14,%rsi
mov    0x20(%rsp),%rdx
mov    0x10(%rsp),%r8
call   1070 <cftfsub@plt>
mov    0x18(%rsp),%rcx
lea    -0x1(%rcx),%eax
movss  0x4(%r14),%xmm0
subss  (%r14),%xmm0
movss  %xmm0,(%r14,%rax,4)
movss  (%r14),%xmm0
addss  0x4(%r14),%xmm0
movss  %xmm0,0x4(%r14)
cmp    $0x8,%ecx
jb     17de <dfst+0x68e>
mov    %ebp,%eax
add    $0x2,%eax
lea    -0x5(%rbp,%rbp,1),%ecx
movaps 0xc62(%rip),%xmm2
xchg   %ax,%ax
lea    -0x4(%rax),%edx
movss  (%r14,%rdx,4),%xmm0
lea    -0x3(%rax),%edx
movaps %xmm0,%xmm1
subss  (%r14,%rdx,4),%xmm1
lea    0x2(%rcx),%esi
movss  %xmm1,(%r14,%rsi,4)
xorps  %xmm2,%xmm0
subss  (%r14,%rdx,4),%xmm0
mov    %ecx,%edx
movss  %xmm0,(%r14,%rdx,4)
add    $0xfffffffe,%eax
add    $0xfffffffc,%ecx
cmp    $0x5,%eax
jg     13a0 <dfst+0x250>
cmpl   $0x7,0x18(%rsp)
jbe    17de <dfst+0x68e>
lea    -0x4(%rbx),%rax
mov    %rax,0x50(%rsp)
lea    0x4(%rbx),%rax
mov    %rax,0x60(%rsp)
lea    -0x10(%rbx),%rax
mov    %rax,0x58(%rsp)
mov    $0x2,%r15d
jmp    1427 <dfst+0x2d7>
nopl   0x0(%rax,%rax,1)
lea    (%r12,%r13,1),%eax
movss  (%rbx,%rax,4),%xmm0
movss  %xmm0,(%rbx)
cmp    $0x3,%r13d
jbe    17ff <dfst+0x6af>
mov    %r12d,%r13d
mov    %r12d,%edi
mov    %rbx,%rsi
mov    0x30(%rsp),%rbp
mov    %ebp,%edx
mov    0x48(%rsp),%r12
mov    %r12,%rcx
call   1050 <dstsub@plt>
cmp    $0x5,%r13d
jb     1480 <dfst+0x330>
mov    %r13d,%edi
mov    %rbx,%rsi
mov    0x20(%rsp),%rdx
mov    0x28(%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    14a9 <dfst+0x359>
cs nopw 0x0(%rax,%rax,1)
cmp    $0x4,%r13d
movaps 0xb75(%rip),%xmm1
jne    14b0 <dfst+0x360>
mov    $0x4,%edi
mov    %rbx,%rsi
mov    0x20(%rsp),%rdx
mov    0x28(%rsp),%rcx
mov    0x10(%rsp),%r8
call   1070 <cftfsub@plt>
movaps 0xb50(%rip),%xmm1
movss  0x4(%rbx),%xmm0
subss  (%rbx),%xmm0
mov    0x18(%rsp),%rax
sub    %r15d,%eax
cltq
movss  %xmm0,(%r14,%rax,4)
movss  (%rbx),%xmm0
addss  0x4(%rbx),%xmm0
movslq %r15d,%rax
movss  %xmm0,(%r14,%rax,4)
cmp    $0x2,%r13d
jbe    17e9 <dfst+0x699>
lea    0x0(,%r15,4),%ecx
mov    %r13d,%esi
movslq %ecx,%rbp
lea    (%r14,%rax,4),%rcx
shl    $0x2,%rax
mov    %r14,%rdx
sub    %rax,%rdx
mov    $0x2,%edi
mov    %rbp,%rax
nopl   0x0(%rax)
movss  (%rbx,%rdi,4),%xmm0
xorps  %xmm1,%xmm0
subss  0x4(%rbx,%rdi,4),%xmm0
movss  %xmm0,(%rdx,%rax,4)
movss  (%rbx,%rdi,4),%xmm0
subss  0x4(%rbx,%rdi,4),%xmm0
movss  %xmm0,(%rcx,%rax,4)
add    $0x2,%rdi
add    %rbp,%rax
cmp    %rsi,%rdi
jb     1510 <dfst+0x3c0>
add    %r15d,%r15d
mov    %r13d,%r12d
shr    %r12d
cmp    $0x4,%r13d
jb     17f2 <dfst+0x6a2>
cmp    $0x3,%r12d
mov    $0x2,%ecx
cmovae %r12d,%ecx
lea    -0x1(%rcx),%r8
mov    $0x1,%r11d
cmp    $0x4,%r8
jae    1650 <dfst+0x500>
mov    %r11,%rax
test   $0x1,%cl
jne    15b1 <dfst+0x461>
mov    %rsi,%rax
sub    %r11,%rax
lea    (%rax,%rsi,1),%rdx
mov    %rsi,%rdi
add    %r11,%rdi
movss  (%rbx,%rdi,4),%xmm0
movss  (%rbx,%rdx,4),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%rbx,%r11,4)
movss  (%rbx,%rdx,4),%xmm1
subss  %xmm0,%xmm1
movss  %xmm1,(%rbx,%rax,4)
lea    0x1(%r11),%rax
not    %r11
add    %rcx,%r11
je     1410 <dfst+0x2c0>
lea    (%rax,%rsi,1),%rdx
add    $0x1,%rdx
mov    %rsi,%rdi
sub    %rax,%rdi
lea    0x0(,%rsi,8),%rsi
lea    0x0(,%rax,4),%rbp
sub    %rbp,%rsi
xchg   %ax,%ax
movss  -0x4(%rbx,%rdx,4),%xmm0
movss  (%rbx,%rsi,1),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%rbx,%rax,4)
movss  (%rbx,%rsi,1),%xmm1
subss  %xmm0,%xmm1
movss  %xmm1,(%rbx,%rdi,4)
movss  (%rbx,%rdx,4),%xmm0
movss  -0x4(%rbx,%rsi,1),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,0x4(%rbx,%rax,4)
movss  -0x4(%rbx,%rsi,1),%xmm1
subss  %xmm0,%xmm1
movss  %xmm1,-0x4(%rbx,%rdi,4)
add    $0x2,%rax
add    $0x2,%rdx
add    $0xfffffffffffffffe,%rdi
add    $0xfffffffffffffff8,%rsi
cmp    %rax,%rcx
jne    15e0 <dfst+0x490>
jmp    1410 <dfst+0x2c0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
lea    -0x2(%rcx),%rax
mov    $0x4,%edx
mul    %rdx
mov    0x50(%rsp),%rdx
lea    (%rdx,%rsi,4),%rbp
seto   %dl
mov    %rbp,%rdi
sub    %rax,%rdi
cmp    %rbp,%rdi
ja     1573 <dfst+0x423>
test   %dl,%dl
jne    1573 <dfst+0x423>
mov    0x50(%rsp),%rdi
lea    (%rdi,%rsi,8),%rbp
mov    %rbp,%rdi
sub    %rax,%rdi
cmp    %rbp,%rdi
ja     1573 <dfst+0x423>
test   %dl,%dl
jne    1573 <dfst+0x423>
mov    %r8,0x70(%rsp)
lea    (%rsi,%rsi,1),%rbp
mov    0x60(%rsp),%rdx
cmp    %rbx,%rdx
setb   0xc(%rsp)
lea    (%rbx,%rcx,4),%r10
lea    (%rdx,%rsi,4),%rax
mov    %rax,0x40(%rsp)
lea    (%rsi,%rcx,1),%rax
lea    (%rbx,%rax,4),%r8
mov    %r8,0x38(%rsp)
mov    %rsi,%rax
sub    %rcx,%rax
lea    (%rdx,%rax,4),%r9
lea    (%rbx,%rsi,4),%rax
mov    %rbp,%rdi
mov    %rbp,0x68(%rsp)
sub    %rcx,%rbp
lea    (%rdx,%rbp,4),%rdi
lea    (%rbx,%rsi,8),%rbp
cmp    %r8,%rdx
setb   0xf(%rsp)
cmp    %r10,0x40(%rsp)
setb   0x40(%rsp)
cmp    %rax,%rdx
setb   %r8b
cmp    %r10,%r9
setb   0xe(%rsp)
cmp    %rbp,%rdx
setb   0xd(%rsp)
cmp    %r10,%rdi
setb   %r10b
cmp    0x38(%rsp),%r9
setb   0x38(%rsp)
cmp    %rbp,%r9
setb   %r9b
cmp    %rax,%rdi
setb   %dil
mov    0x40(%rsp),%dl
test   %dl,0xf(%rsp)
jne    1573 <dfst+0x423>
and    0xe(%rsp),%r8b
jne    1573 <dfst+0x423>
and    %r10b,0xd(%rsp)
jne    1573 <dfst+0x423>
mov    0x38(%rsp),%dl
and    %dl,0xc(%rsp)
jne    1573 <dfst+0x423>
and    %dil,%r9b
jne    1573 <dfst+0x423>
mov    0x70(%rsp),%r8
mov    %r8,%r10
and    $0xfffffffffffffffc,%r10
lea    0x1(%r10),%r11
mov    0x58(%rsp),%rdx
mov    0x68(%rsp),%rdi
lea    (%rdx,%rdi,4),%rdx
mov    %r10,%r9
neg    %r9
mov    $0x4,%edi
xor    %ebp,%ebp
xchg   %ax,%ax
movups (%rdx,%rbp,4),%xmm0
shufps $0x1b,%xmm0,%xmm0
movups (%rax,%rdi,1),%xmm1
addps  %xmm1,%xmm0
movups %xmm0,(%rbx,%rdi,1)
movups (%rdx,%rbp,4),%xmm0
shufps $0x1b,%xmm1,%xmm1
subps  %xmm1,%xmm0
movups %xmm0,-0x10(%rax,%rbp,4)
add    $0xfffffffffffffffc,%rbp
add    $0x10,%rdi
cmp    %rbp,%r9
jne    17a0 <dfst+0x650>
cmp    %r10,%r8
je     1410 <dfst+0x2c0>
jmp    1573 <dfst+0x423>
movss  (%rbx),%xmm0
mov    $0x2,%eax
jmp    1802 <dfst+0x6b2>
add    %r15d,%r15d
mov    %r13d,%r12d
shr    %r12d
add    %r13d,%r12d
movss  (%rbx,%r12,4),%xmm0
movss  %xmm0,(%rbx)
movslq %r15d,%rax
movss  %xmm0,(%r14,%rax,4)
movl   $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, float* %1, float* %2, i32* %3, float* %4) {
%6 = alloca i32
%7 = alloca float*
%8 = alloca float*
%9 = alloca i32*
%10 = alloca float*
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca float
%19 = alloca float
%20 = alloca float
%21 = alloca float
store i32 %0, i32* %6
store float* %1, float** %7
store float* %2, float** %8
store i32* %3, i32** %9
store float* %4, float** %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 float*, float** %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 float*, float** %10
%49 = load i32, i32* %16
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds float, float* %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 float*, float** %7
%69 = load i32, i32* %11
%70 = sext i32 %69 to i64
%71 = getelementptr inbounds float, float* %68, i64 %70
%72 = load float, float* %71
%73 = load float*, float** %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 float, float* %73, i64 %77
%79 = load float, float* %78
%80 = fadd float %72, %79
store float %80, float* %18
%81 = load float*, float** %7
%82 = load i32, i32* %11
%83 = sext i32 %82 to i64
%84 = getelementptr inbounds float, float* %81, i64 %83
%85 = load float, float* %84
%86 = load float*, float** %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 float, float* %86, i64 %90
%92 = load float, float* %91
%93 = fsub float %85, %92
store float %93, float* %19
%94 = load float*, float** %7
%95 = load i32, i32* %12
%96 = sext i32 %95 to i64
%97 = getelementptr inbounds float, float* %94, i64 %96
%98 = load float, float* %97
%99 = load float*, float** %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 float, float* %99, i64 %103
%105 = load float, float* %104
%106 = fadd float %98, %105
store float %106, float* %20
%107 = load float*, float** %7
%108 = load i32, i32* %12
%109 = sext i32 %108 to i64
%110 = getelementptr inbounds float, float* %107, i64 %109
%111 = load float, float* %110
%112 = load float*, float** %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 float, float* %112, i64 %116
%118 = load float, float* %117
%119 = fsub float %111, %118
store float %119, float* %21
%120 = load float, float* %18
%121 = load float*, float** %7
%122 = load i32, i32* %11
%123 = sext i32 %122 to i64
%124 = getelementptr inbounds float, float* %121, i64 %123
store float %120, float* %124
%125 = load float, float* %20
%126 = load float*, float** %7
%127 = load i32, i32* %12
%128 = sext i32 %127 to i64
%129 = getelementptr inbounds float, float* %126, i64 %128
store float %125, float* %129
%130 = load float, float* %19
%131 = load float, float* %21
%132 = fadd float %130, %131
%133 = load float*, float** %8
%134 = load i32, i32* %11
%135 = sext i32 %134 to i64
%136 = getelementptr inbounds float, float* %133, i64 %135
store float %132, float* %136
%137 = load float, float* %19
%138 = load float, float* %21
%139 = fsub float %137, %138
%140 = load float*, float** %8
%141 = load i32, i32* %12
%142 = sext i32 %141 to i64
%143 = getelementptr inbounds float, float* %140, i64 %142
store float %139, float* %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 float*, float** %7
%149 = load i32, i32* %15
%150 = sext i32 %149 to i64
%151 = getelementptr inbounds float, float* %148, i64 %150
%152 = load float, float* %151
%153 = load float*, float** %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 float, float* %153, i64 %157
%159 = load float, float* %158
%160 = fsub float %152, %159
%161 = load float*, float** %8
%162 = getelementptr inbounds float, float* %161, i64 0
store float %160, float* %162
%163 = load float*, float** %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 float, float* %163, i64 %167
%169 = load float, float* %168
%170 = load float*, float** %7
%171 = load i32, i32* %15
%172 = sext i32 %171 to i64
%173 = getelementptr inbounds float, float* %170, i64 %172
%174 = load float, float* %173
%175 = fadd float %174, %169
store float %175, float* %173
%176 = load float*, float** %7
%177 = load i32, i32* %14
%178 = sext i32 %177 to i64
%179 = getelementptr inbounds float, float* %176, i64 %178
%180 = load float, float* %179
%181 = load float*, float** %7
%182 = getelementptr inbounds float, float* %181, i64 0
store float %180, float* %182
%183 = load i32, i32* %14
%184 = load float*, float** %7
%185 = load i32, i32* %17
%186 = load float*, float** %10
%187 = load i32, i32* %16
%188 = sext i32 %187 to i64
%189 = getelementptr inbounds float, float* %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 float*, float** %7
%195 = load i32*, i32** %9
%196 = load i32, i32* %16
%197 = load float*, float** %10
%198 = load i32, i32* %14
%199 = load float*, float** %7
%200 = load i32, i32* %17
%201 = load float*, float** %10
%202 = load i32, i32* %16
%203 = sext i32 %202 to i64
%204 = getelementptr inbounds float, float* %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 float*, float** %7
%211 = load i32*, i32** %9
%212 = load i32, i32* %16
%213 = load float*, float** %10
br label %214
214: ; preds = %208, %205
br label %215
215: ; preds = %214, %192
%216 = load float*, float** %7
%217 = getelementptr inbounds float, float* %216, i64 1
%218 = load float, float* %217
%219 = load float*, float** %7
%220 = getelementptr inbounds float, float* %219, i64 0
%221 = load float, float* %220
%222 = fsub float %218, %221
%223 = load float*, float** %7
%224 = load i32, i32* %6
%225 = sub nsw i32 %224, 1
%226 = sext i32 %225 to i64
%227 = getelementptr inbounds float, float* %223, i64 %226
store float %222, float* %227
%228 = load float*, float** %7
%229 = getelementptr inbounds float, float* %228, i64 0
%230 = load float, float* %229
%231 = load float*, float** %7
%232 = getelementptr inbounds float, float* %231, i64 1
%233 = load float, float* %232
%234 = fadd float %230, %233
%235 = load float*, float** %7
%236 = getelementptr inbounds float, float* %235, i64 1
store float %234, float* %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 float*, float** %7
%244 = load i32, i32* %11
%245 = sext i32 %244 to i64
%246 = getelementptr inbounds float, float* %243, i64 %245
%247 = load float, float* %246
%248 = load float*, float** %7
%249 = load i32, i32* %11
%250 = add nsw i32 %249, 1
%251 = sext i32 %250 to i64
%252 = getelementptr inbounds float, float* %248, i64 %251
%253 = load float, float* %252
%254 = fsub float %247, %253
%255 = load float*, float** %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 float, float* %255, i64 %259
store float %254, float* %260
%261 = load float*, float** %7
%262 = load i32, i32* %11
%263 = sext i32 %262 to i64
%264 = getelementptr inbounds float, float* %261, i64 %263
%265 = load float, float* %264
%266 = fneg float %265
%267 = load float*, float** %7
%268 = load i32, i32* %11
%269 = add nsw i32 %268, 1
%270 = sext i32 %269 to i64
%271 = getelementptr inbounds float, float* %267, i64 %270
%272 = load float, float* %271
%273 = fsub float %266, %272
%274 = load float*, float** %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 float, float* %274, i64 %278
store float %273, float* %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 float*, float** %8
%291 = load i32, i32* %17
%292 = load float*, float** %10
%293 = load i32, i32* %16
%294 = sext i32 %293 to i64
%295 = getelementptr inbounds float, float* %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 float*, float** %8
%301 = load i32*, i32** %9
%302 = load i32, i32* %16
%303 = load float*, float** %10
%304 = load i32, i32* %14
%305 = load float*, float** %8
%306 = load i32, i32* %17
%307 = load float*, float** %10
%308 = load i32, i32* %16
%309 = sext i32 %308 to i64
%310 = getelementptr inbounds float, float* %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 float*, float** %8
%317 = load i32*, i32** %9
%318 = load i32, i32* %16
%319 = load float*, float** %10
br label %320
320: ; preds = %314, %311
br label %321
321: ; preds = %320, %298
%322 = load float*, float** %8
%323 = getelementptr inbounds float, float* %322, i64 1
%324 = load float, float* %323
%325 = load float*, float** %8
%326 = getelementptr inbounds float, float* %325, i64 0
%327 = load float, float* %326
%328 = fsub float %324, %327
%329 = load float*, float** %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 float, float* %329, i64 %333
store float %328, float* %334
%335 = load float*, float** %8
%336 = getelementptr inbounds float, float* %335, i64 0
%337 = load float, float* %336
%338 = load float*, float** %8
%339 = getelementptr inbounds float, float* %338, i64 1
%340 = load float, float* %339
%341 = fadd float %337, %340
%342 = load float*, float** %7
%343 = load i32, i32* %13
%344 = sext i32 %343 to i64
%345 = getelementptr inbounds float, float* %342, i64 %344
store float %341, float* %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 float*, float** %8
%356 = load i32, i32* %11
%357 = sext i32 %356 to i64
%358 = getelementptr inbounds float, float* %355, i64 %357
%359 = load float, float* %358
%360 = fneg float %359
%361 = load float*, float** %8
%362 = load i32, i32* %11
%363 = add nsw i32 %362, 1
%364 = sext i32 %363 to i64
%365 = getelementptr inbounds float, float* %361, i64 %364
%366 = load float, float* %365
%367 = fsub float %360, %366
%368 = load float*, float** %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 float, float* %368, i64 %372
store float %367, float* %373
%374 = load float*, float** %8
%375 = load i32, i32* %11
%376 = sext i32 %375 to i64
%377 = getelementptr inbounds float, float* %374, i64 %376
%378 = load float, float* %377
%379 = load float*, float** %8
%380 = load i32, i32* %11
%381 = add nsw i32 %380, 1
%382 = sext i32 %381 to i64
%383 = getelementptr inbounds float, float* %379, i64 %382
%384 = load float, float* %383
%385 = fsub float %378, %384
%386 = load float*, float** %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 float, float* %386, i64 %390
store float %385, float* %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 float*, float** %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 float, float* %408, i64 %412
%414 = load float, float* %413
%415 = load float*, float** %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 float, float* %415, i64 %419
%421 = load float, float* %420
%422 = fadd float %414, %421
%423 = load float*, float** %8
%424 = load i32, i32* %11
%425 = sext i32 %424 to i64
%426 = getelementptr inbounds float, float* %423, i64 %425
store float %422, float* %426
%427 = load float*, float** %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 float, float* %427, i64 %431
%433 = load float, float* %432
%434 = load float*, float** %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 float, float* %434, i64 %438
%440 = load float, float* %439
%441 = fsub float %433, %440
%442 = load float*, float** %8
%443 = load i32, i32* %12
%444 = sext i32 %443 to i64
%445 = getelementptr inbounds float, float* %442, i64 %444
store float %441, float* %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 float*, float** %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 float, float* %450, i64 %454
%456 = load float, float* %455
%457 = load float*, float** %8
%458 = getelementptr inbounds float, float* %457, i64 0
store float %456, float* %458
%459 = load i32, i32* %15
store i32 %459, i32* %14
br label %285
460: ; preds = %285
%461 = load float*, float** %8
%462 = getelementptr inbounds float, float* %461, i64 0
%463 = load float, float* %462
%464 = load float*, float** %7
%465 = load i32, i32* %13
%466 = sext i32 %465 to i64
%467 = getelementptr inbounds float, float* %464, i64 %466
store float %463, float* %467
br label %468
468: ; preds = %460, %52
%469 = load float*, float** %7
%470 = getelementptr inbounds float, float* %469, i64 0
store float 0.000000e+00, float* %470
ret void
}

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