void realfft(n, real) int n;
double *real;
{
  double a, b;
  int i, j, k;
  fht(real, n);
  for (i = 1, j = n - 1, k = n / 2; i < k; i++, j--) {
    a = real[i];
    b = real[j];
    real[j] = (a - b) * 0.5;
    real[i] = (a + b) * 0.5;
  }
}


x86 assembly (clang, -O0):
<realfft>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    -0x10(%rbp),%rdi
mov    -0x4(%rbp),%esi
mov    $0x0,%al
call   1030 <fht@plt>
movl   $0x1,-0x24(%rbp)
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x28(%rbp)
mov    -0x4(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x2c(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    11d8 <realfft+0xc8>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x18(%rbp),%xmm0
subsd  -0x20(%rbp),%xmm0
movsd  0xe73(%rip),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x18(%rbp),%xmm0
addsd  -0x20(%rbp),%xmm0
movsd  0xe50(%rip),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x28(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x28(%rbp)
jmp    114b <realfft+0x3b>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<realfft>:
push   %r14
push   %rbx
push   %rax
mov    %rsi,%rbx
mov    %edi,%r14d
mov    %rsi,%rdi
mov    %r14d,%esi
xor    %eax,%eax
call   1030 <fht@plt>
cmp    $0x4,%r14d
jl     1217 <realfft+0x107>
mov    %r14d,%edx
shr    %r14d
lea    -0x1(%r14),%r8
mov    $0x1,%eax
cmp    $0x2,%r14d
jne    114e <realfft+0x3e>
mov    %rdx,%rdi
jmp    11df <realfft+0xcf>
mov    %r8,%rsi
and    $0xfffffffffffffffe,%rsi
neg    %rsi
mov    $0x1,%eax
movsd  0xe9b(%rip),%xmm0
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rbx,%rax,8),%xmm1
movsd  -0x8(%rbx,%rdx,8),%xmm2
movapd %xmm1,%xmm3
subsd  %xmm2,%xmm3
mulsd  %xmm0,%xmm3
movsd  %xmm3,-0x8(%rbx,%rdx,8)
addsd  %xmm2,%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx,%rax,8)
lea    -0x2(%rdx),%rdi
movsd  0x8(%rbx,%rax,8),%xmm1
movsd  -0x10(%rbx,%rdx,8),%xmm2
movapd %xmm1,%xmm3
subsd  %xmm2,%xmm3
mulsd  %xmm0,%xmm3
movsd  %xmm3,-0x10(%rbx,%rdx,8)
addsd  %xmm2,%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rbx,%rax,8)
lea    (%rsi,%rax,1),%rcx
add    $0x2,%rcx
add    $0x2,%rax
mov    %rdi,%rdx
cmp    $0x1,%rcx
jne    1170 <realfft+0x60>
test   $0x1,%r8b
je     1217 <realfft+0x107>
movsd  (%rbx,%rax,8),%xmm0
movsd  -0x8(%rbx,%rdi,8),%xmm1
movapd %xmm0,%xmm2
subsd  %xmm1,%xmm2
movsd  0xe00(%rip),%xmm3
mulsd  %xmm3,%xmm2
movsd  %xmm2,-0x8(%rbx,%rdi,8)
addsd  %xmm1,%xmm0
mulsd  %xmm3,%xmm0
movsd  %xmm0,(%rbx,%rax,8)
add    $0x8,%rsp
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define void @realfft(i32 %0, double* %1) {
%3 = alloca i32
%4 = alloca double*
%5 = alloca double
%6 = alloca double
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
store i32 %0, i32* %3
store double* %1, double** %4
%10 = load double*, double** %4
%11 = load i32, i32* %3
%12 = call i32 (double*, i32, ...) bitcast (i32 (...)* @fht to i32 (double*, i32, ...)*)(double* %10, i32 %11)
store i32 1, i32* %7
%13 = load i32, i32* %3
%14 = sub nsw i32 %13, 1
store i32 %14, i32* %8
%15 = load i32, i32* %3
%16 = sdiv i32 %15, 2
store i32 %16, i32* %9
br label %17
17: ; preds = %48, %2
%18 = load i32, i32* %7
%19 = load i32, i32* %9
%20 = icmp slt i32 %18, %19
br i1 %20, label %21, label %53
21: ; preds = %17
%22 = load double*, double** %4
%23 = load i32, i32* %7
%24 = sext i32 %23 to i64
%25 = getelementptr inbounds double, double* %22, i64 %24
%26 = load double, double* %25
store double %26, double* %5
%27 = load double*, double** %4
%28 = load i32, i32* %8
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds double, double* %27, i64 %29
%31 = load double, double* %30
store double %31, double* %6
%32 = load double, double* %5
%33 = load double, double* %6
%34 = fsub double %32, %33
%35 = fmul double %34, 5.000000e-01
%36 = load double*, double** %4
%37 = load i32, i32* %8
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds double, double* %36, i64 %38
store double %35, double* %39
%40 = load double, double* %5
%41 = load double, double* %6
%42 = fadd double %40, %41
%43 = fmul double %42, 5.000000e-01
%44 = load double*, double** %4
%45 = load i32, i32* %7
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds double, double* %44, i64 %46
store double %43, double* %47
br label %48
48: ; preds = %21
%49 = load i32, i32* %7
%50 = add nsw i32 %49, 1
store i32 %50, i32* %7
%51 = load i32, i32* %8
%52 = add nsw i32 %51, -1
store i32 %52, i32* %8
br label %17
53: ; preds = %17
ret void
}

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