#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
scalar_t__ cos(float);
scalar_t__ sin(float);
void rftbsub(int n, float *a) {
  int i, i0, j, k;
  float ec, w1r, w1i, wkr, wki, wdr, wdi, ss, xr, xi, yr, yi;
  ec = 2 * ((float)1.570796326794896619231321691639751442098584699687) / n;
  wkr = 0;
  wki = 0;
  wdi = (float)cos(ec);
  wdr = (float)sin(ec);
  wdi *= wdr;
  wdr *= wdr;
  w1r = 1 - 2 * wdr;
  w1i = 2 * wdi;
  ss = 2 * w1i;
  i = n >> 1;
  for (;;) {
    i0 = i - 4 * 64;
    if (i0 < 4) {
      i0 = 4;
    }
    for (j = i - 4; j >= i0; j -= 4) {
      k = n - j;
      xr = a[j + 2] - a[k - 2];
      xi = a[j + 3] + a[k - 1];
      yr = wdr * xr + wdi * xi;
      yi = wdr * xi - wdi * xr;
      a[j + 2] -= yr;
      a[j + 3] -= yi;
      a[k - 2] += yr;
      a[k - 1] -= yi;
      wkr += ss * wdi;
      wki += ss * (((float)0.5) - wdr);
      xr = a[j] - a[k];
      xi = a[j + 1] + a[k + 1];
      yr = wkr * xr + wki * xi;
      yi = wkr * xi - wki * xr;
      a[j] -= yr;
      a[j + 1] -= yi;
      a[k] += yr;
      a[k + 1] -= yi;
      wdr += ss * wki;
      wdi += ss * (((float)0.5) - wkr);
    }
    if (i0 == 4) {
      break;
    }
    wkr = ((float)0.5) * (float)sin(ec * i0);
    wki = ((float)0.5) * (float)cos(ec * i0);
    wdr = ((float)0.5) - (wkr * w1r - wki * w1i);
    wdi = wkr * w1i + wki * w1r;
    wkr = ((float)0.5) - wkr;
    i = i0;
  }
  xr = a[2] - a[n - 2];
  xi = a[3] + a[n - 1];
  yr = wdr * xr + wdi * xi;
  yi = wdr * xi - wdi * xr;
  a[2] -= yr;
  a[3] -= yi;
  a[n - 2] += yr;
  a[n - 1] -= yi;
}


x86 assembly (clang, -O0):
<rftbsub>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
cvtsi2ssl -0x4(%rbp),%xmm1
movss  0xecc(%rip),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,-0x24(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x24(%rbp),%xmm0
call   1040 <cos@plt>
cvtsi2ss %rax,%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x24(%rbp),%xmm0
call   1030 <sin@plt>
cvtsi2ss %rax,%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x38(%rbp),%xmm0
mulss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x38(%rbp),%xmm0
mulss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x38(%rbp),%xmm1
addss  %xmm1,%xmm1
movss  0xe58(%rip),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,-0x28(%rbp)
movss  0xe43(%rip),%xmm0
mulss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  0xe31(%rip),%xmm0
mulss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
mov    -0x4(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x14(%rbp)
mov    -0x14(%rbp),%eax
sub    $0x100,%eax
mov    %eax,-0x18(%rbp)
cmpl   $0x4,-0x18(%rbp)
jge    11fe <rftbsub+0xde>
movl   $0x4,-0x18(%rbp)
mov    -0x14(%rbp),%eax
sub    $0x4,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jl     14d0 <rftbsub+0x3b0>
mov    -0x4(%rbp),%eax
sub    -0x1c(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x38(%rbp),%xmm0
movss  -0x44(%rbp),%xmm2
movss  -0x3c(%rbp),%xmm1
mulss  -0x48(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x38(%rbp),%xmm0
movss  -0x48(%rbp),%xmm2
movss  -0x3c(%rbp),%xmm1
mulss  -0x44(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x4c(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movss  -0x50(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movss  -0x4c(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
movss  -0x50(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movss  -0x40(%rbp),%xmm0
movss  -0x3c(%rbp),%xmm2
movss  -0x30(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x40(%rbp),%xmm0
movss  0xcaa(%rip),%xmm2
subss  -0x38(%rbp),%xmm2
movss  -0x34(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x30(%rbp),%xmm0
movss  -0x44(%rbp),%xmm2
movss  -0x34(%rbp),%xmm1
mulss  -0x48(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x30(%rbp),%xmm0
movss  -0x48(%rbp),%xmm2
movss  -0x34(%rbp),%xmm1
mulss  -0x44(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x4c(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movss  -0x50(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movss  -0x4c(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
movss  -0x50(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movss  -0x40(%rbp),%xmm0
movss  -0x34(%rbp),%xmm2
movss  -0x38(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x40(%rbp),%xmm0
movss  0xb61(%rip),%xmm2
subss  -0x30(%rbp),%xmm2
movss  -0x3c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x3c(%rbp)
mov    -0x1c(%rbp),%eax
sub    $0x4,%eax
mov    %eax,-0x1c(%rbp)
jmp    1207 <rftbsub+0xe7>
cmpl   $0x4,-0x18(%rbp)
jne    14df <rftbsub+0x3bf>
jmp    15a9 <rftbsub+0x489>
movss  -0x24(%rbp),%xmm0
cvtsi2ssl -0x18(%rbp),%xmm1
mulss  %xmm1,%xmm0
call   1030 <sin@plt>
cvtsi2ss %rax,%xmm1
movss  0xb0d(%rip),%xmm0
mulss  %xmm1,%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x24(%rbp),%xmm0
cvtsi2ssl -0x18(%rbp),%xmm1
mulss  %xmm1,%xmm0
call   1040 <cos@plt>
cvtsi2ss %rax,%xmm1
movss  0xae4(%rip),%xmm0
mulss  %xmm1,%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x30(%rbp),%xmm1
movss  -0x28(%rbp),%xmm2
movss  -0x34(%rbp),%xmm0
mulss  -0x2c(%rbp),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movss  0xaaa(%rip),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x30(%rbp),%xmm0
movss  -0x2c(%rbp),%xmm2
movss  -0x34(%rbp),%xmm1
mulss  -0x28(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  0xa78(%rip),%xmm0
subss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
mov    -0x18(%rbp),%eax
mov    %eax,-0x14(%rbp)
jmp    11e2 <rftbsub+0xc2>
mov    -0x10(%rbp),%rax
movss  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x10(%rbp),%rax
movss  0xc(%rax),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x38(%rbp),%xmm0
movss  -0x44(%rbp),%xmm2
movss  -0x3c(%rbp),%xmm1
mulss  -0x48(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x38(%rbp),%xmm0
movss  -0x48(%rbp),%xmm2
movss  -0x3c(%rbp),%xmm1
mulss  -0x44(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x4c(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movss  0x8(%rax),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,0x8(%rax)
movss  -0x50(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movss  0xc(%rax),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,0xc(%rax)
movss  -0x4c(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
movss  -0x50(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<rftbsub>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x20,%rsp
mov    %edi,%ebp
cvtsi2ss %edi,%xmm1
mov    %rsi,%rbx
movss  0xec3(%rip),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,0xc(%rsp)
call   1040 <cos@plt>
xorps  %xmm0,%xmm0
cvtsi2ss %rax,%xmm0
movss  %xmm0,0x8(%rsp)
movss  0xc(%rsp),%xmm0
call   1030 <sin@plt>
movss  0x8(%rsp),%xmm9
cvtsi2ss %rax,%xmm8
mulss  %xmm8,%xmm9
mulss  %xmm8,%xmm8
movaps %xmm8,%xmm1
addss  %xmm8,%xmm1
movss  0xe78(%rip),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,0x18(%rsp)
movaps %xmm9,%xmm0
addss  %xmm9,%xmm0
movss  %xmm0,0x14(%rsp)
movaps %xmm0,%xmm10
addss  %xmm0,%xmm10
movslq %ebp,%r14
sar    %ebp
lea    (%rbx,%r14,4),%r15
add    $0x14,%r15
xorps  %xmm2,%xmm2
movss  0xe41(%rip),%xmm11
xorps  %xmm1,%xmm1
movss  %xmm10,0x10(%rsp)
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
cmp    $0x105,%ebp
mov    $0x104,%eax
cmovge %ebp,%eax
lea    -0x100(%rax),%r12d
add    $0xffffff04,%eax
cmp    %eax,%ebp
jl     133c <rftbsub+0x21c>
movslq %ebp,%rsi
lea    -0x4(%rsi),%rax
mov    %r12d,%ecx
shl    $0x2,%rsi
mov    %r15,%rdx
sub    %rsi,%rdx
cs nopw 0x0(%rax,%rax,1)
movss  0x8(%rbx,%rax,4),%xmm3
movss  0xc(%rbx,%rax,4),%xmm4
movaps %xmm3,%xmm5
subss  -0xc(%rdx),%xmm5
movss  -0x8(%rdx),%xmm6
addss  %xmm4,%xmm6
movaps %xmm9,%xmm7
mulss  %xmm6,%xmm7
movaps %xmm8,%xmm0
mulss  %xmm5,%xmm0
addss  %xmm7,%xmm0
mulss  %xmm9,%xmm5
mulss  %xmm8,%xmm6
subss  %xmm5,%xmm6
subss  %xmm0,%xmm3
movss  %xmm3,0x8(%rbx,%rax,4)
subss  %xmm6,%xmm4
movss  %xmm4,0xc(%rbx,%rax,4)
addss  -0xc(%rdx),%xmm0
movss  %xmm0,-0xc(%rdx)
movss  -0x8(%rdx),%xmm0
subss  %xmm6,%xmm0
movss  %xmm0,-0x8(%rdx)
movaps %xmm10,%xmm0
mulss  %xmm9,%xmm0
addss  %xmm0,%xmm2
movaps %xmm11,%xmm0
subss  %xmm8,%xmm0
mulss  %xmm10,%xmm0
addss  %xmm0,%xmm1
movss  (%rbx,%rax,4),%xmm0
movss  0x4(%rbx,%rax,4),%xmm3
movaps %xmm0,%xmm4
subss  -0x4(%rdx),%xmm4
movss  (%rdx),%xmm5
addss  %xmm3,%xmm5
movaps %xmm1,%xmm6
mulss  %xmm5,%xmm6
movaps %xmm2,%xmm7
mulss  %xmm4,%xmm7
addss  %xmm6,%xmm7
mulss  %xmm1,%xmm4
mulss  %xmm2,%xmm5
subss  %xmm4,%xmm5
subss  %xmm7,%xmm0
movss  %xmm0,(%rbx,%rax,4)
subss  %xmm5,%xmm3
movss  %xmm3,0x4(%rbx,%rax,4)
addss  -0x4(%rdx),%xmm7
movss  %xmm7,-0x4(%rdx)
movss  (%rdx),%xmm0
subss  %xmm5,%xmm0
movss  %xmm0,(%rdx)
movaps %xmm10,%xmm0
mulss  %xmm1,%xmm0
addss  %xmm0,%xmm8
movaps %xmm11,%xmm0
subss  %xmm2,%xmm0
mulss  %xmm10,%xmm0
addss  %xmm0,%xmm9
add    $0xfffffffffffffffc,%rax
add    $0x10,%rdx
cmp    %rcx,%rax
jge    1220 <rftbsub+0x100>
cmp    $0x4,%r12d
je     13ee <rftbsub+0x2ce>
xorps  %xmm0,%xmm0
cvtsi2ss %r12d,%xmm0
mulss  0xc(%rsp),%xmm0
movss  %xmm0,0x1c(%rsp)
call   1030 <sin@plt>
xorps  %xmm0,%xmm0
cvtsi2ss %rax,%xmm0
mulss  0xc99(%rip),%xmm0
movss  %xmm0,0x8(%rsp)
movss  0x1c(%rsp),%xmm0
call   1040 <cos@plt>
movss  0xc7f(%rip),%xmm11
movss  0x10(%rsp),%xmm10
movss  0x8(%rsp),%xmm9
xorps  %xmm1,%xmm1
cvtsi2ss %rax,%xmm1
mulss  %xmm11,%xmm1
movss  0x14(%rsp),%xmm4
movaps %xmm4,%xmm0
mulss  %xmm1,%xmm0
movaps %xmm9,%xmm2
movss  0x18(%rsp),%xmm3
mulss  %xmm3,%xmm2
subss  %xmm0,%xmm2
movaps %xmm11,%xmm8
subss  %xmm2,%xmm8
movaps %xmm3,%xmm0
mulss  %xmm1,%xmm0
movaps %xmm11,%xmm2
subss  %xmm9,%xmm2
mulss  %xmm4,%xmm9
addss  %xmm0,%xmm9
mov    %r12d,%ebp
jmp    11e0 <rftbsub+0xc0>
movss  0x8(%rbx),%xmm0
movss  0xc(%rbx),%xmm1
movaps %xmm0,%xmm2
subss  -0x8(%rbx,%r14,4),%xmm2
movss  -0x4(%rbx,%r14,4),%xmm3
addss  %xmm1,%xmm3
movaps %xmm9,%xmm4
mulss  %xmm3,%xmm4
movaps %xmm8,%xmm5
mulss  %xmm2,%xmm5
addss  %xmm4,%xmm5
mulss  %xmm9,%xmm2
mulss  %xmm8,%xmm3
subss  %xmm2,%xmm3
subss  %xmm5,%xmm0
movss  %xmm0,0x8(%rbx)
subss  %xmm3,%xmm1
movss  %xmm1,0xc(%rbx)
addss  -0x8(%rbx,%r14,4),%xmm5
movss  %xmm5,-0x8(%rbx,%r14,4)
movss  -0x4(%rbx,%r14,4),%xmm0
subss  %xmm3,%xmm0
movss  %xmm0,-0x4(%rbx,%r14,4)
add    $0x20,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @rftbsub(i32 %0, float* nocapture %1) {
%3 = sitofp i32 %0 to float
%4 = fdiv float 0x400921FB60000000, %3
%5 = tail call i64 @cos(float %4)
%6 = sitofp i64 %5 to float
%7 = tail call i64 @sin(float %4)
%8 = sitofp i64 %7 to float
%9 = fmul float %6, %8
%10 = fmul float %8, %8
%11 = tail call float @llvm.fmuladd.f32(float %10, float -2.000000e+00, float 1.000000e+00)
%12 = fmul float %9, 2.000000e+00
%13 = fmul float %12, 2.000000e+00
%14 = ashr i32 %0, 1
%15 = sext i32 %0 to i64
br label %16
16: ; preds = %100, %2
%17 = phi float [ 0.000000e+00, %2 ], [ %115, %100 ]
%18 = phi float [ 0.000000e+00, %2 ], [ %108, %100 ]
%19 = phi float [ %10, %2 ], [ %112, %100 ]
%20 = phi float [ %9, %2 ], [ %114, %100 ]
%21 = phi i32 [ %14, %2 ], [ %24, %100 ]
%22 = icmp sgt i32 %21, 260
%23 = select i1 %22, i32 %21, i32 260
%24 = add nsw i32 %23, -256
%25 = add nsw i32 %23, -252
%26 = icmp slt i32 %21, %25
br i1 %26, label %96, label %27
27: ; preds = %16
%28 = sext i32 %21 to i64
%29 = add nsw i64 %28, -4
%30 = zext i32 %24 to i64
br label %31
31: ; preds = %27, %31
%32 = phi i64 [ %28, %27 ], [ %38, %31 ]
%33 = phi i64 [ %29, %27 ], [ %94, %31 ]
%34 = phi float [ %20, %27 ], [ %93, %31 ]
%35 = phi float [ %19, %27 ], [ %91, %31 ]
%36 = phi float [ %18, %27 ], [ %67, %31 ]
%37 = phi float [ %17, %27 ], [ %65, %31 ]
%38 = add nsw i64 %32, -4
%39 = sub nsw i64 %15, %33
%40 = add nsw i64 %32, -2
%41 = getelementptr inbounds float, float* %1, i64 %40
%42 = load float, float* %41
%43 = add nsw i64 %39, -2
%44 = getelementptr inbounds float, float* %1, i64 %43
%45 = load float, float* %44
%46 = fsub float %42, %45
%47 = add nsw i64 %32, -1
%48 = getelementptr inbounds float, float* %1, i64 %47
%49 = load float, float* %48
%50 = add nsw i64 %39, -1
%51 = getelementptr inbounds float, float* %1, i64 %50
%52 = load float, float* %51
%53 = fadd float %49, %52
%54 = fmul float %34, %53
%55 = tail call float @llvm.fmuladd.f32(float %35, float %46, float %54)
%56 = fneg float %34
%57 = fmul float %46, %56
%58 = tail call float @llvm.fmuladd.f32(float %35, float %53, float %57)
%59 = fsub float %42, %55
store float %59, float* %41
%60 = fsub float %49, %58
store float %60, float* %48
%61 = load float, float* %44
%62 = fadd float %61, %55
store float %62, float* %44
%63 = load float, float* %51
%64 = fsub float %63, %58
store float %64, float* %51
%65 = tail call float @llvm.fmuladd.f32(float %13, float %34, float %37)
%66 = fsub float 5.000000e-01, %35
%67 = tail call float @llvm.fmuladd.f32(float %13, float %66, float %36)
%68 = getelementptr inbounds float, float* %1, i64 %33
%69 = load float, float* %68
%70 = getelementptr inbounds float, float* %1, i64 %39
%71 = load float, float* %70
%72 = fsub float %69, %71
%73 = add nsw i64 %32, -3
%74 = getelementptr inbounds float, float* %1, i64 %73
%75 = load float, float* %74
%76 = add nsw i64 %39, 1
%77 = getelementptr inbounds float, float* %1, i64 %76
%78 = load float, float* %77
%79 = fadd float %75, %78
%80 = fmul float %67, %79
%81 = tail call float @llvm.fmuladd.f32(float %65, float %72, float %80)
%82 = fneg float %67
%83 = fmul float %72, %82
%84 = tail call float @llvm.fmuladd.f32(float %65, float %79, float %83)
%85 = fsub float %69, %81
store float %85, float* %68
%86 = fsub float %75, %84
store float %86, float* %74
%87 = load float, float* %70
%88 = fadd float %87, %81
store float %88, float* %70
%89 = load float, float* %77
%90 = fsub float %89, %84
store float %90, float* %77
%91 = tail call float @llvm.fmuladd.f32(float %13, float %67, float %35)
%92 = fsub float 5.000000e-01, %65
%93 = tail call float @llvm.fmuladd.f32(float %13, float %92, float %34)
%94 = add nsw i64 %33, -4
%95 = icmp slt i64 %94, %30
br i1 %95, label %96, label %31
96: ; preds = %31, %16
%97 = phi float [ %19, %16 ], [ %91, %31 ]
%98 = phi float [ %20, %16 ], [ %93, %31 ]
%99 = icmp eq i32 %24, 4
br i1 %99, label %116, label %100
100: ; preds = %96
%101 = sitofp i32 %24 to float
%102 = fmul float %4, %101
%103 = tail call i64 @sin(float %102)
%104 = sitofp i64 %103 to float
%105 = fmul float %104, 5.000000e-01
%106 = tail call i64 @cos(float %102)
%107 = sitofp i64 %106 to float
%108 = fmul float %107, 5.000000e-01
%109 = fneg float %108
%110 = fmul float %12, %109
%111 = tail call float @llvm.fmuladd.f32(float %105, float %11, float %110)
%112 = fsub float 5.000000e-01, %111
%113 = fmul float %11, %108
%114 = tail call float @llvm.fmuladd.f32(float %105, float %12, float %113)
%115 = fsub float 5.000000e-01, %105
br label %16
116: ; preds = %96
%117 = getelementptr inbounds float, float* %1, i64 2
%118 = load float, float* %117
%119 = add nsw i32 %0, -2
%120 = sext i32 %119 to i64
%121 = getelementptr inbounds float, float* %1, i64 %120
%122 = load float, float* %121
%123 = fsub float %118, %122
%124 = getelementptr inbounds float, float* %1, i64 3
%125 = load float, float* %124
%126 = add nsw i32 %0, -1
%127 = sext i32 %126 to i64
%128 = getelementptr inbounds float, float* %1, i64 %127
%129 = load float, float* %128
%130 = fadd float %125, %129
%131 = fmul float %98, %130
%132 = tail call float @llvm.fmuladd.f32(float %97, float %123, float %131)
%133 = fneg float %98
%134 = fmul float %123, %133
%135 = tail call float @llvm.fmuladd.f32(float %97, float %130, float %134)
%136 = fsub float %118, %132
store float %136, float* %117
%137 = fsub float %125, %135
store float %137, float* %124
%138 = load float, float* %121
%139 = fadd float %138, %132
store float %139, float* %121
%140 = load float, float* %128
%141 = fsub float %140, %135
store float %141, float* %128
ret void
}