void wt1(float a[], unsigned long n, int isign,
         void (*wtstep)(float[], unsigned long, int)) {
  unsigned long nn;
  if (n < 4)
    return;
  if (isign >= 0) {
    for (nn = n; nn >= 4; nn >>= 1)
      (*wtstep)(a, nn, isign);
  } else {
    for (nn = 4; nn <= n; nn <<= 1)
      (*wtstep)(a, nn, isign);
  }
}


x86 assembly (clang, -O0):
<wt1>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
cmpq   $0x4,-0x10(%rbp)
jae    1127 <wt1+0x27>
jmp    11a8 <wt1+0xa8>
cmpl   $0x0,-0x14(%rbp)
jl     116b <wt1+0x6b>
mov    -0x10(%rbp),%rax
mov    %rax,-0x28(%rbp)
cmpq   $0x4,-0x28(%rbp)
jb     1166 <wt1+0x66>
mov    -0x20(%rbp),%rax
mov    -0x8(%rbp),%rdi
mov    -0x28(%rbp),%rsi
mov    -0x14(%rbp),%edx
call   *%rax
mov    -0x28(%rbp),%rax
shr    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    1139 <wt1+0x39>
jmp    11a8 <wt1+0xa8>
movq   $0x4,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x10(%rbp),%rax
ja     11a3 <wt1+0xa3>
mov    -0x20(%rbp),%rax
mov    -0x8(%rbp),%rdi
mov    -0x28(%rbp),%rsi
mov    -0x14(%rbp),%edx
call   *%rax
mov    -0x28(%rbp),%rax
shl    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    1173 <wt1+0x73>
jmp    11a8 <wt1+0xa8>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<wt1>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
cmp    $0x4,%rsi
jb     1164 <wt1+0x64>
mov    %rcx,%r14
mov    %edx,%r15d
mov    %rsi,%rbp
mov    %rdi,%r12
test   %edx,%edx
js     114a <wt1+0x4a>
mov    %rbp,%rbx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %r12,%rdi
mov    %rbp,%rsi
mov    %r15d,%edx
call   *%r14
shr    %rbx
cmp    $0x7,%rbp
mov    %rbx,%rbp
ja     1130 <wt1+0x30>
jmp    1164 <wt1+0x64>
mov    $0x4,%ebx
nop
mov    %r12,%rdi
mov    %rbx,%rsi
mov    %r15d,%edx
call   *%r14
add    %rbx,%rbx
cmp    %rbp,%rbx
jbe    1150 <wt1+0x50>
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @wt1(float* %0, i64 %1, i32 %2, void (float*, i64, i32)* %3) {
%5 = alloca float*
%6 = alloca i64
%7 = alloca i32
%8 = alloca void (float*, i64, i32)*
%9 = alloca i64
store float* %0, float** %5
store i64 %1, i64* %6
store i32 %2, i32* %7
store void (float*, i64, i32)* %3, void (float*, i64, i32)** %8
%10 = load i64, i64* %6
%11 = icmp ult i64 %10, 4
br i1 %11, label %12, label %13
12: ; preds = %4
br label %44
13: ; preds = %4
%14 = load i32, i32* %7
%15 = icmp sge i32 %14, 0
br i1 %15, label %16, label %30
16: ; preds = %13
%17 = load i64, i64* %6
store i64 %17, i64* %9
br label %18
18: ; preds = %26, %16
%19 = load i64, i64* %9
%20 = icmp uge i64 %19, 4
br i1 %20, label %21, label %29
21: ; preds = %18
%22 = load void (float*, i64, i32)*, void (float*, i64, i32)** %8
%23 = load float*, float** %5
%24 = load i64, i64* %9
%25 = load i32, i32* %7
br label %26
26: ; preds = %21
%27 = load i64, i64* %9
%28 = lshr i64 %27, 1
store i64 %28, i64* %9
br label %18
29: ; preds = %18
br label %44
30: ; preds = %13
store i64 4, i64* %9
br label %31
31: ; preds = %40, %30
%32 = load i64, i64* %9
%33 = load i64, i64* %6
%34 = icmp ule i64 %32, %33
br i1 %34, label %35, label %43
35: ; preds = %31
%36 = load void (float*, i64, i32)*, void (float*, i64, i32)** %8
%37 = load float*, float** %5
%38 = load i64, i64* %9
%39 = load i32, i32* %7
br label %40
40: ; preds = %35
%41 = load i64, i64* %9
%42 = shl i64 %41, 1
store i64 %42, i64* %9
br label %31
43: ; preds = %31
br label %44
44: ; preds = %12, %43, %29
ret void
}

LLVM IR (-O2):
define void @wt1(float* %0, i64 %1, i32 %2, void (float*, i64, i32)* nocapture %3) {
%5 = icmp ult i64 %1, 4
br i1 %5, label %16, label %6
6: ; preds = %4
%7 = icmp sgt i32 %2, -1
br i1 %7, label %8, label %12
8: ; preds = %6, %8
%9 = phi i64 [ %10, %8 ], [ %1, %6 ]
tail call void %3(float* %0, i64 %9, i32 %2)
%10 = lshr i64 %9, 1
%11 = icmp ugt i64 %9, 7
br i1 %11, label %8, label %16
12: ; preds = %6, %12
%13 = phi i64 [ %14, %12 ], [ 4, %6 ]
tail call void %3(float* %0, i64 %13, i32 %2)
%14 = shl i64 %13, 1
%15 = icmp ugt i64 %14, %1
br i1 %15, label %16, label %12
16: ; preds = %12, %8, %4
ret void
}