void rdft(int n, int isgn, float *a, int *ip, float *w) {
  int j, nw, nc;
  float xi;
  void bitrv2(int n, int *ip, float *a), cftsub(int n, float *a, float *w),
      rftsub(int n, float *a, int nc, float *c);
  nw = ip[0];
  nc = ip[1];
  if (isgn < 0) {
    a[1] = 0.5 * (a[1] - a[0]);
    a[0] += a[1];
    for (j = 3; j <= n - 1; j += 2) {
      a[j] = -a[j];
    }
    if (n > 4) {
      rftsub(n, a, nc, w + nw);
      bitrv2(n, ip + 2, a);
    }
    cftsub(n, a, w);
    for (j = 1; j <= n - 1; j += 2) {
      a[j] = -a[j];
    }
  } else {
    if (n > 4) {
      bitrv2(n, ip + 2, a);
    }
    cftsub(n, a, w);
    if (n > 4) {
      rftsub(n, a, nc, w + nw);
    }
    xi = a[0] - a[1];
    a[0] += a[1];
    a[1] = xi;
  }
}


x86 assembly (clang, -O0):
<rdft>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    -0x18(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x28(%rbp)
mov    -0x18(%rbp),%rax
mov    0x4(%rax),%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x0,-0x8(%rbp)
jge    1295 <rdft+0x165>
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm0
mov    -0x10(%rbp),%rax
subss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe7c(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x10(%rbp),%rax
movss  %xmm0,0x4(%rax)
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm0
mov    -0x10(%rbp),%rax
addss  (%rax),%xmm0
movss  %xmm0,(%rax)
movl   $0x3,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     11f7 <rdft+0xc7>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x24(%rbp)
jmp    11b1 <rdft+0x81>
cmpl   $0x4,-0x4(%rbp)
jle    1233 <rdft+0x103>
mov    -0x4(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x2c(%rbp),%edx
mov    -0x20(%rbp),%rcx
movslq -0x28(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rcx
call   1050 <rftsub@plt>
mov    -0x4(%rbp),%edi
mov    -0x18(%rbp),%rsi
add    $0x8,%rsi
mov    -0x10(%rbp),%rdx
call   1040 <bitrv2@plt>
mov    -0x4(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x20(%rbp),%rdx
call   1030 <cftsub@plt>
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     1290 <rdft+0x160>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x24(%rbp)
jmp    124a <rdft+0x11a>
jmp    1324 <rdft+0x1f4>
cmpl   $0x4,-0x4(%rbp)
jle    12b3 <rdft+0x183>
mov    -0x4(%rbp),%edi
mov    -0x18(%rbp),%rsi
add    $0x8,%rsi
mov    -0x10(%rbp),%rdx
call   1040 <bitrv2@plt>
mov    -0x4(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x20(%rbp),%rdx
call   1030 <cftsub@plt>
cmpl   $0x4,-0x4(%rbp)
jle    12eb <rdft+0x1bb>
mov    -0x4(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x2c(%rbp),%edx
mov    -0x20(%rbp),%rcx
movslq -0x28(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rcx
call   1050 <rftsub@plt>
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x10(%rbp),%rax
subss  0x4(%rax),%xmm0
movss  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm0
mov    -0x10(%rbp),%rax
addss  (%rax),%xmm0
movss  %xmm0,(%rax)
movss  -0x30(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movss  %xmm0,0x4(%rax)
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<rdft>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    %r8,%r14
mov    %rcx,%r13
mov    %rdx,%rbx
mov    %edi,%r12d
movslq (%rcx),%r15
mov    0x4(%rcx),%ebp
test   %esi,%esi
js     1191 <rdft+0x61>
cmp    $0x5,%r12d
jl     121d <rdft+0xed>
add    $0x8,%r13
mov    %r12d,%edi
mov    %r13,%rsi
mov    %rbx,%rdx
call   1040 <bitrv2@plt>
mov    %r12d,%edi
mov    %rbx,%rsi
mov    %r14,%rdx
call   1030 <cftsub@plt>
lea    (%r14,%r15,4),%rcx
mov    %r12d,%edi
mov    %rbx,%rsi
mov    %ebp,%edx
call   1050 <rftsub@plt>
jmp    122b <rdft+0xfb>
mov    %ebp,%edx
movss  (%rbx),%xmm0
movss  0x4(%rbx),%xmm1
subss  %xmm0,%xmm1
mulss  0xe58(%rip),%xmm1
movss  %xmm1,0x4(%rbx)
addss  %xmm0,%xmm1
movss  %xmm1,(%rbx)
mov    %r12d,%ebp
cmp    $0x4,%r12d
jl     124a <rdft+0x11a>
mov    $0x3,%eax
movaps 0xe42(%rip),%xmm0
xchg   %ax,%ax
movss  (%rbx,%rax,4),%xmm1
xorps  %xmm0,%xmm1
movss  %xmm1,(%rbx,%rax,4)
add    $0x2,%rax
cmp    %rbp,%rax
jb     11d0 <rdft+0xa0>
cmp    $0x5,%r12d
jl     124a <rdft+0x11a>
lea    (%r14,%r15,4),%rcx
mov    %r12d,%edi
mov    %rbx,%rsi
call   1050 <rftsub@plt>
add    $0x8,%r13
mov    %r12d,%edi
mov    %r13,%rsi
mov    %rbx,%rdx
call   1040 <bitrv2@plt>
mov    %r12d,%edi
mov    %rbx,%rsi
mov    %r14,%rdx
call   1030 <cftsub@plt>
jmp    125e <rdft+0x12e>
mov    %r12d,%edi
mov    %rbx,%rsi
mov    %r14,%rdx
call   1030 <cftsub@plt>
movss  (%rbx),%xmm0
movss  0x4(%rbx),%xmm1
movaps %xmm0,%xmm2
subss  %xmm1,%xmm2
addss  %xmm1,%xmm0
movss  %xmm0,(%rbx)
movss  %xmm2,0x4(%rbx)
jmp    1286 <rdft+0x156>
mov    %r12d,%edi
mov    %rbx,%rsi
mov    %r14,%rdx
call   1030 <cftsub@plt>
cmp    $0x2,%r12d
jl     1286 <rdft+0x156>
mov    $0x1,%eax
movaps 0xda6(%rip),%xmm0
nopw   0x0(%rax,%rax,1)
movss  (%rbx,%rax,4),%xmm1
xorps  %xmm0,%xmm1
movss  %xmm1,(%rbx,%rax,4)
add    $0x2,%rax
cmp    %rbp,%rax
jb     1270 <rdft+0x140>
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @rdft(i32 %0, i32 %1, float* %2, i32* %3, float* %4) {
%6 = load i32, i32* %3
%7 = getelementptr inbounds i32, i32* %3, i64 1
%8 = load i32, i32* %7
%9 = icmp slt i32 %1, 0
br i1 %9, label %10, label %44
10: ; preds = %5
%11 = getelementptr inbounds float, float* %2, i64 1
%12 = load float, float* %11
%13 = load float, float* %2
%14 = fsub float %12, %13
%15 = fmul float %14, 5.000000e-01
store float %15, float* %11
%16 = fadd float %13, %15
store float %16, float* %2
%17 = icmp sgt i32 %0, 3
br i1 %17, label %18, label %33
18: ; preds = %10
%19 = zext i32 %0 to i64
br label %20
20: ; preds = %18, %20
%21 = phi i64 [ 3, %18 ], [ %25, %20 ]
%22 = getelementptr inbounds float, float* %2, i64 %21
%23 = load float, float* %22
%24 = fneg float %23
store float %24, float* %22
%25 = add nuw nsw i64 %21, 2
%26 = icmp ult i64 %25, %19
br i1 %26, label %20, label %27
27: ; preds = %20
%28 = icmp sgt i32 %0, 4
br i1 %28, label %29, label %33
29: ; preds = %27
%30 = sext i32 %6 to i64
%31 = getelementptr inbounds float, float* %4, i64 %30
tail call void @rftsub(i32 %0, float* %2, i32 %8, float* %31)
%32 = getelementptr inbounds i32, i32* %3, i64 2
tail call void @bitrv2(i32 %0, i32* %32, float* %2)
tail call void @cftsub(i32 %0, float* %2, float* %4)
br label %35
33: ; preds = %10, %27
tail call void @cftsub(i32 %0, float* %2, float* %4)
%34 = icmp sgt i32 %0, 1
br i1 %34, label %35, label %57
35: ; preds = %29, %33
%36 = zext i32 %0 to i64
br label %37
37: ; preds = %35, %37
%38 = phi i64 [ 1, %35 ], [ %42, %37 ]
%39 = getelementptr inbounds float, float* %2, i64 %38
%40 = load float, float* %39
%41 = fneg float %40
store float %41, float* %39
%42 = add nuw nsw i64 %38, 2
%43 = icmp ult i64 %42, %36
br i1 %43, label %37, label %57
44: ; preds = %5
%45 = icmp sgt i32 %0, 4
br i1 %45, label %46, label %50
46: ; preds = %44
%47 = getelementptr inbounds i32, i32* %3, i64 2
tail call void @bitrv2(i32 %0, i32* %47, float* %2)
tail call void @cftsub(i32 %0, float* %2, float* %4)
%48 = sext i32 %6 to i64
%49 = getelementptr inbounds float, float* %4, i64 %48
tail call void @rftsub(i32 %0, float* %2, i32 %8, float* %49)
br label %51
50: ; preds = %44
tail call void @cftsub(i32 %0, float* %2, float* %4)
br label %51
51: ; preds = %50, %46
%52 = load float, float* %2
%53 = getelementptr inbounds float, float* %2, i64 1
%54 = load float, float* %53
%55 = fsub float %52, %54
%56 = fadd float %52, %54
store float %56, float* %2
store float %55, float* %53
br label %57
57: ; preds = %37, %33, %51
ret void
}