void cftrec4(int n, double *a) {
  int cfttree(int n, int j, int k, double *a);
  void cftleaf(int n, int isplt, double *a);
  void cftmdl1(int n, double *a);
  int isplt, j, k, m;
  m = n;
  while (m > 512) {
    m >>= 2;
    cftmdl1(m, &a[n - m]);
  }
  cftleaf(m, 1, &a[n - m]);
  k = 0;
  for (j = n - m; j > 0; j -= m) {
    k++;
    isplt = cfttree(m, j, k, a);
    cftleaf(m, isplt, &a[j - m]);
  }
}


x86 assembly (clang, -O0):
<cftrec4>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    -0x4(%rbp),%eax
mov    %eax,-0x20(%rbp)
cmpl   $0x200,-0x20(%rbp)
jle    117b <cftrec4+0x4b>
mov    -0x20(%rbp),%eax
sar    $0x2,%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x4(%rbp),%eax
sub    -0x20(%rbp),%eax
cltq
shl    $0x3,%rax
add    %rax,%rsi
call   1040 <cftmdl1@plt>
jmp    1145 <cftrec4+0x15>
mov    -0x20(%rbp),%edi
mov    -0x10(%rbp),%rdx
mov    -0x4(%rbp),%eax
sub    -0x20(%rbp),%eax
cltq
shl    $0x3,%rax
add    %rax,%rdx
mov    $0x1,%esi
call   1030 <cftleaf@plt>
movl   $0x0,-0x1c(%rbp)
mov    -0x4(%rbp),%eax
sub    -0x20(%rbp),%eax
mov    %eax,-0x18(%rbp)
cmpl   $0x0,-0x18(%rbp)
jle    1201 <cftrec4+0xd1>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x20(%rbp),%edi
mov    -0x18(%rbp),%esi
mov    -0x1c(%rbp),%edx
mov    -0x10(%rbp),%rcx
call   1050 <cfttree@plt>
mov    %eax,-0x14(%rbp)
mov    -0x20(%rbp),%edi
mov    -0x14(%rbp),%esi
mov    -0x10(%rbp),%rdx
mov    -0x18(%rbp),%eax
sub    -0x20(%rbp),%eax
cltq
shl    $0x3,%rax
add    %rax,%rdx
call   1030 <cftleaf@plt>
mov    -0x20(%rbp),%ecx
mov    -0x18(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x18(%rbp)
jmp    11ab <cftrec4+0x7b>
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cftrec4>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    %rsi,%r14
mov    %edi,%ebx
mov    %edi,%r15d
cmp    $0x201,%edi
jl     1172 <cftrec4+0x42>
mov    %ebx,%ebp
mov    %ebx,%r15d
shr    $0x2,%r15d
mov    %ebx,%eax
sub    %r15d,%eax
cltq
lea    (%r14,%rax,8),%rsi
mov    %r15d,%edi
call   1040 <cftmdl1@plt>
cmp    $0x803,%ebp
mov    %r15d,%ebp
ja     1150 <cftrec4+0x20>
mov    %ebx,%ebp
sub    %r15d,%ebp
movslq %ebp,%r13
lea    (%r14,%r13,8),%rdx
mov    $0x1,%r12d
mov    %r15d,%edi
mov    $0x1,%esi
call   1030 <cftleaf@plt>
test   %r13d,%r13d
jle    11cf <cftrec4+0x9f>
sub    %r15d,%ebx
sub    %r15d,%ebx
nopl   0x0(%rax)
mov    %r15d,%edi
mov    %ebp,%esi
mov    %r12d,%edx
mov    %r14,%rcx
call   1050 <cfttree@plt>
sub    %r15d,%ebp
movslq %ebx,%rbx
lea    (%r14,%rbx,8),%rdx
mov    %r15d,%edi
mov    %eax,%esi
call   1030 <cftleaf@plt>
add    $0x1,%r12d
sub    %r15d,%ebx
test   %ebp,%ebp
jg     11a0 <cftrec4+0x70>
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @cftrec4(i32 %0, double* %1) {
%3 = icmp sgt i32 %0, 512
br i1 %3, label %4, label %11
4: ; preds = %2, %4
%5 = phi i32 [ %6, %4 ], [ %0, %2 ]
%6 = lshr i32 %5, 2
%7 = sub nsw i32 %0, %6
%8 = sext i32 %7 to i64
%9 = getelementptr inbounds double, double* %1, i64 %8
tail call void @cftmdl1(i32 %6, double* %9)
%10 = icmp ugt i32 %5, 2051
br i1 %10, label %4, label %11
11: ; preds = %4, %2
%12 = phi i32 [ %0, %2 ], [ %6, %4 ]
%13 = sub nsw i32 %0, %12
%14 = sext i32 %13 to i64
%15 = getelementptr inbounds double, double* %1, i64 %14
tail call void @cftleaf(i32 %12, i32 1, double* %15)
%16 = icmp sgt i32 %13, 0
br i1 %16, label %17, label %26
17: ; preds = %11, %17
%18 = phi i32 [ %20, %17 ], [ 0, %11 ]
%19 = phi i32 [ %22, %17 ], [ %13, %11 ]
%20 = add nuw nsw i32 %18, 1
%21 = tail call i32 @cfttree(i32 %12, i32 %19, i32 %20, double* %1)
%22 = sub nsw i32 %19, %12
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds double, double* %1, i64 %23
tail call void @cftleaf(i32 %12, i32 %21, double* %24)
%25 = icmp sgt i32 %22, 0
br i1 %25, label %17, label %26
26: ; preds = %17, %11
ret void
}