int fft_real_radix_top(int npoint, int f_first) {
  int i, m;
  int fft_n_thresh = 128 * 17 / 20;
  int f_small_prohibit;
  static int trm[5] = {16, 5, 8, 3, 4};
  f_small_prohibit = (f_first == 1 && npoint > fft_n_thresh) ? 1 : 0;
  for (i = 0; i < 5; i++) {
    if (npoint % trm[i] == 0) {
      m = npoint / trm[i];
      if (m == 1)
        return trm[i];
      if (m % 8 == 0)
        return 8;
      if (m % 4 == 0)
        return 4;
      if (m % 7 == 0)
        return 7;
      if (m % 5 == 0)
        return 5;
      if (f_small_prohibit == 1) {
        if (m % 2 == 0)
          return 4;
        if (m % 3 == 0)
          return 4;
      } else {
        if (m % 2 == 0)
          return 2;
        if (m % 3 == 0)
          return 3;
      }
      return 0;
    }
  }
  return 0;
}


x86 assembly (clang, -O0):
<fft_real_radix_top>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movl   $0x6c,-0x18(%rbp)
xor    %eax,%eax
cmpl   $0x1,-0xc(%rbp)
mov    %al,-0x1d(%rbp)
jne    112c <fft_real_radix_top+0x2c>
mov    -0x8(%rbp),%eax
cmp    -0x18(%rbp),%eax
setg   %al
mov    %al,-0x1d(%rbp)
mov    -0x1d(%rbp),%dl
xor    %eax,%eax
mov    $0x1,%ecx
test   $0x1,%dl
cmovne %ecx,%eax
mov    %eax,-0x1c(%rbp)
movl   $0x0,-0x10(%rbp)
cmpl   $0x5,-0x10(%rbp)
jge    12d3 <fft_real_radix_top+0x1d3>
mov    -0x8(%rbp),%eax
movslq -0x10(%rbp),%rsi
lea    0x2ed2(%rip),%rcx
cltd
idivl  (%rcx,%rsi,4)
cmp    $0x0,%edx
jne    12c0 <fft_real_radix_top+0x1c0>
mov    -0x8(%rbp),%eax
movslq -0x10(%rbp),%rsi
lea    0x2eb7(%rip),%rcx
cltd
idivl  (%rcx,%rsi,4)
mov    %eax,-0x14(%rbp)
cmpl   $0x1,-0x14(%rbp)
jne    11a0 <fft_real_radix_top+0xa0>
movslq -0x10(%rbp),%rcx
lea    0x2e9b(%rip),%rax
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x4(%rbp)
jmp    12da <fft_real_radix_top+0x1da>
mov    -0x14(%rbp),%eax
mov    $0x8,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
jne    11c0 <fft_real_radix_top+0xc0>
movl   $0x8,-0x4(%rbp)
jmp    12da <fft_real_radix_top+0x1da>
mov    -0x14(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
jne    11e0 <fft_real_radix_top+0xe0>
movl   $0x4,-0x4(%rbp)
jmp    12da <fft_real_radix_top+0x1da>
mov    -0x14(%rbp),%eax
mov    $0x7,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
jne    1200 <fft_real_radix_top+0x100>
movl   $0x7,-0x4(%rbp)
jmp    12da <fft_real_radix_top+0x1da>
mov    -0x14(%rbp),%eax
mov    $0x5,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
jne    1220 <fft_real_radix_top+0x120>
movl   $0x5,-0x4(%rbp)
jmp    12da <fft_real_radix_top+0x1da>
cmpl   $0x1,-0x1c(%rbp)
jne    126f <fft_real_radix_top+0x16f>
mov    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
jne    124a <fft_real_radix_top+0x14a>
movl   $0x4,-0x4(%rbp)
jmp    12da <fft_real_radix_top+0x1da>
mov    -0x14(%rbp),%eax
mov    $0x3,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
jne    126a <fft_real_radix_top+0x16a>
movl   $0x4,-0x4(%rbp)
jmp    12da <fft_real_radix_top+0x1da>
jmp    12b4 <fft_real_radix_top+0x1b4>
mov    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
jne    128f <fft_real_radix_top+0x18f>
movl   $0x2,-0x4(%rbp)
jmp    12da <fft_real_radix_top+0x1da>
mov    -0x14(%rbp),%eax
mov    $0x3,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
jne    12af <fft_real_radix_top+0x1af>
movl   $0x3,-0x4(%rbp)
jmp    12da <fft_real_radix_top+0x1da>
jmp    12b4 <fft_real_radix_top+0x1b4>
movl   $0x0,-0x4(%rbp)
jmp    12da <fft_real_radix_top+0x1da>
jmp    12c5 <fft_real_radix_top+0x1c5>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1146 <fft_real_radix_top+0x46>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<fft_real_radix_top>:
cmp    $0x1,%esi
sete   %sil
cmp    $0x6d,%edi
setge  %r8b
test   $0xf,%dil
je     1149 <fft_real_radix_top+0x49>
imul   $0xcccccccd,%edi,%eax
add    $0x19999999,%eax
cmp    $0x33333333,%eax
jb     1150 <fft_real_radix_top+0x50>
test   $0x7,%dil
je     1157 <fft_real_radix_top+0x57>
imul   $0xaaaaaaab,%edi,%eax
add    $0x2aaaaaaa,%eax
cmp    $0x55555555,%eax
jb     115e <fft_real_radix_top+0x5e>
xor    %ecx,%ecx
test   $0x3,%dil
je     1165 <fft_real_radix_top+0x65>
mov    %ecx,%eax
ret
mov    $0x10,%ecx
jmp    116a <fft_real_radix_top+0x6a>
mov    $0x5,%ecx
jmp    116a <fft_real_radix_top+0x6a>
mov    $0x8,%ecx
jmp    116a <fft_real_radix_top+0x6a>
mov    $0x3,%ecx
jmp    116a <fft_real_radix_top+0x6a>
mov    $0x4,%ecx
mov    %edi,%eax
cltd
idiv   %ecx
cmp    $0x1,%eax
je     1146 <fft_real_radix_top+0x46>
test   $0x7,%al
je     11d5 <fft_real_radix_top+0xd5>
mov    $0x4,%ecx
test   $0x3,%al
je     1146 <fft_real_radix_top+0x46>
imul   $0xb6db6db7,%eax,%edx
add    $0x12492492,%edx
cmp    $0x24924925,%edx
jb     11db <fft_real_radix_top+0xdb>
imul   $0xcccccccd,%eax,%edx
add    $0x19999999,%edx
cmp    $0x33333333,%edx
jb     11e1 <fft_real_radix_top+0xe1>
and    %sil,%r8b
je     11e7 <fft_real_radix_top+0xe7>
test   $0x1,%al
je     1146 <fft_real_radix_top+0x46>
movslq %eax,%rdx
imul   $0x55555556,%rdx,%rdx
mov    %rdx,%rsi
shr    $0x3f,%rsi
shr    $0x20,%rdx
add    %esi,%edx
lea    (%rdx,%rdx,2),%edx
sub    %edx,%eax
jne    1206 <fft_real_radix_top+0x106>
jmp    1146 <fft_real_radix_top+0x46>
mov    $0x8,%eax
ret
mov    $0x7,%eax
ret
mov    $0x5,%eax
ret
mov    $0x2,%ecx
test   $0x1,%al
je     1146 <fft_real_radix_top+0x46>
imul   $0xaaaaaaab,%eax,%eax
add    $0x2aaaaaaa,%eax
cmp    $0x55555555,%eax
jb     1209 <fft_real_radix_top+0x109>
xor    %eax,%eax
ret
mov    $0x3,%ecx
mov    %ecx,%eax
ret

LLVM IR (-O0):
@fft_real_radix_top.trm = internal global [5 x i32] [i32 16, i32 5, i32 8, i32 3, i32 4]
define i32 @fft_real_radix_top(i32 %0, i32 %1) {
%3 = alloca i32
%4 = alloca i32
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
store i32 %0, i32* %4
store i32 %1, i32* %5
store i32 108, i32* %8
%10 = load i32, i32* %5
%11 = icmp eq i32 %10, 1
br i1 %11, label %12, label %16
12: ; preds = %2
%13 = load i32, i32* %4
%14 = load i32, i32* %8
%15 = icmp sgt i32 %13, %14
br label %16
16: ; preds = %12, %2
%17 = phi i1 [ false, %2 ], [ %15, %12 ]
%18 = zext i1 %17 to i64
%19 = select i1 %17, i32 1, i32 0
store i32 %19, i32* %9
store i32 0, i32* %6
br label %20
20: ; preds = %92, %16
%21 = load i32, i32* %6
%22 = icmp slt i32 %21, 5
br i1 %22, label %23, label %95
23: ; preds = %20
%24 = load i32, i32* %4
%25 = load i32, i32* %6
%26 = sext i32 %25 to i64
%27 = getelementptr inbounds [5 x i32], [5 x i32]* @fft_real_radix_top.trm, i64 0, i64 %26
%28 = load i32, i32* %27
%29 = srem i32 %24, %28
%30 = icmp eq i32 %29, 0
br i1 %30, label %31, label %91
31: ; preds = %23
%32 = load i32, i32* %4
%33 = load i32, i32* %6
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds [5 x i32], [5 x i32]* @fft_real_radix_top.trm, i64 0, i64 %34
%36 = load i32, i32* %35
%37 = sdiv i32 %32, %36
store i32 %37, i32* %7
%38 = load i32, i32* %7
%39 = icmp eq i32 %38, 1
br i1 %39, label %40, label %45
40: ; preds = %31
%41 = load i32, i32* %6
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds [5 x i32], [5 x i32]* @fft_real_radix_top.trm, i64 0, i64 %42
%44 = load i32, i32* %43
store i32 %44, i32* %3
br label %96
45: ; preds = %31
%46 = load i32, i32* %7
%47 = srem i32 %46, 8
%48 = icmp eq i32 %47, 0
br i1 %48, label %49, label %50
49: ; preds = %45
store i32 8, i32* %3
br label %96
50: ; preds = %45
%51 = load i32, i32* %7
%52 = srem i32 %51, 4
%53 = icmp eq i32 %52, 0
br i1 %53, label %54, label %55
54: ; preds = %50
store i32 4, i32* %3
br label %96
55: ; preds = %50
%56 = load i32, i32* %7
%57 = srem i32 %56, 7
%58 = icmp eq i32 %57, 0
br i1 %58, label %59, label %60
59: ; preds = %55
store i32 7, i32* %3
br label %96
60: ; preds = %55
%61 = load i32, i32* %7
%62 = srem i32 %61, 5
%63 = icmp eq i32 %62, 0
br i1 %63, label %64, label %65
64: ; preds = %60
store i32 5, i32* %3
br label %96
65: ; preds = %60
%66 = load i32, i32* %9
%67 = icmp eq i32 %66, 1
br i1 %67, label %68, label %79
68: ; preds = %65
%69 = load i32, i32* %7
%70 = srem i32 %69, 2
%71 = icmp eq i32 %70, 0
br i1 %71, label %72, label %73
72: ; preds = %68
store i32 4, i32* %3
br label %96
73: ; preds = %68
%74 = load i32, i32* %7
%75 = srem i32 %74, 3
%76 = icmp eq i32 %75, 0
br i1 %76, label %77, label %78
77: ; preds = %73
store i32 4, i32* %3
br label %96
78: ; preds = %73
br label %90
79: ; preds = %65
%80 = load i32, i32* %7
%81 = srem i32 %80, 2
%82 = icmp eq i32 %81, 0
br i1 %82, label %83, label %84
83: ; preds = %79
store i32 2, i32* %3
br label %96
84: ; preds = %79
%85 = load i32, i32* %7
%86 = srem i32 %85, 3
%87 = icmp eq i32 %86, 0
br i1 %87, label %88, label %89
88: ; preds = %84
store i32 3, i32* %3
br label %96
89: ; preds = %84
br label %90
90: ; preds = %89, %78
store i32 0, i32* %3
br label %96
91: ; preds = %23
br label %92
92: ; preds = %91
%93 = load i32, i32* %6
%94 = add nsw i32 %93, 1
store i32 %94, i32* %6
br label %20
95: ; preds = %20
store i32 0, i32* %3
br label %96
96: ; preds = %95, %90, %88, %83, %77, %72, %64, %59, %54, %49, %40
%97 = load i32, i32* %3
ret i32 %97
}

LLVM IR (-O2):
define i32 @fft_real_radix_top(i32 %0, i32 %1) {
%3 = icmp eq i32 %1, 1
%4 = icmp sgt i32 %0, 108
%5 = and i1 %4, %3
%6 = and i32 %0, 15
%7 = icmp eq i32 %6, 0
br i1 %7, label %20, label %8
8: ; preds = %2
%9 = srem i32 %0, 5
%10 = icmp eq i32 %9, 0
br i1 %10, label %20, label %11
11: ; preds = %8
%12 = and i32 %0, 7
%13 = icmp eq i32 %12, 0
br i1 %13, label %20, label %14
14: ; preds = %11
%15 = srem i32 %0, 3
%16 = icmp eq i32 %15, 0
br i1 %16, label %20, label %17
17: ; preds = %14
%18 = and i32 %0, 3
%19 = icmp eq i32 %18, 0
br i1 %19, label %20, label %48
20: ; preds = %17, %14, %11, %8, %2
%21 = phi i32 [ 16, %2 ], [ 5, %8 ], [ 8, %11 ], [ 3, %14 ], [ 4, %17 ]
%22 = sdiv i32 %0, %21
%23 = icmp eq i32 %22, 1
br i1 %23, label %48, label %24
24: ; preds = %20
%25 = and i32 %22, 7
%26 = icmp eq i32 %25, 0
br i1 %26, label %48, label %27
27: ; preds = %24
%28 = and i32 %22, 3
%29 = icmp eq i32 %28, 0
br i1 %29, label %48, label %30
30: ; preds = %27
%31 = srem i32 %22, 7
%32 = icmp eq i32 %31, 0
br i1 %32, label %48, label %33
33: ; preds = %30
%34 = srem i32 %22, 5
%35 = icmp eq i32 %34, 0
br i1 %35, label %48, label %36
36: ; preds = %33
%37 = and i32 %22, 1
%38 = icmp eq i32 %37, 0
br i1 %5, label %39, label %43
39: ; preds = %36
%40 = srem i32 %22, 3
%41 = icmp eq i32 %40, 0
%42 = or i1 %38, %41
br i1 %42, label %48, label %47
43: ; preds = %36
br i1 %38, label %48, label %44
44: ; preds = %43
%45 = srem i32 %22, 3
%46 = icmp eq i32 %45, 0
br i1 %46, label %48, label %47
47: ; preds = %39, %44
br label %48
48: ; preds = %17, %44, %43, %39, %33, %30, %27, %24, %20, %47
%49 = phi i32 [ 0, %47 ], [ %21, %20 ], [ 8, %24 ], [ 4, %27 ], [ 7, %30 ], [ 5, %33 ], [ 4, %39 ], [ 2, %43 ], [ 3, %44 ], [ 0, %17 ]
ret i32 %49
}