#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int M_PI;
float cos(float);
int decompose(int, int *, int const *);
float sin(float);
void cffti1_ps(int n, float *wa, int *ifac) {
  static const int ntryh[] = {5, 3, 4, 2, 0};
  int k1, j, ii;
  int nf = decompose(n, ifac, ntryh);
  float argh = (float)((2 * M_PI) / n);
  int i = 1;
  int l1 = 1;
  for (k1 = 1; k1 <= nf; k1++) {
    int ip = ifac[k1 + 1];
    int ld = 0;
    int l2 = l1 * ip;
    int ido = n / l2;
    int idot = ido + ido + 2;
    int ipm = ip - 1;
    for (j = 1; j <= ipm; j++) {
      float argld;
      int i1 = i, fi = 0;
      wa[i - 1] = 1;
      wa[i] = 0;
      ld += l1;
      argld = (float)ld * argh;
      for (ii = 4; ii <= idot; ii += 2) {
        i += 2;
        fi += 1;
        wa[i - 1] = cos((float)fi * argld);
        wa[i] = sin((float)fi * argld);
      }
      if (ip > 5) {
        wa[i1 - 1] = wa[i - 1];
        wa[i1] = wa[i];
      }
    }
    l1 = l2;
  }
}


x86 assembly (clang, -O0):
<cffti1_ps>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x4(%rbp),%edi
mov    -0x18(%rbp),%rsi
lea    0xebf(%rip),%rdx
call   1050 <decompose@plt>
mov    %eax,-0x28(%rbp)
mov    0x2e88(%rip),%rax
mov    (%rax),%eax
shl    $0x1,%eax
cltd
idivl  -0x4(%rbp)
cvtsi2ss %eax,%xmm0
movss  %xmm0,-0x2c(%rbp)
movl   $0x1,-0x30(%rbp)
movl   $0x1,-0x34(%rbp)
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x28(%rbp),%eax
jg     1317 <cffti1_ps+0x1e7>
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x38(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x34(%rbp),%eax
imul   -0x38(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x4(%rbp),%eax
cltd
idivl  -0x40(%rbp)
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
add    -0x44(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x48(%rbp)
mov    -0x38(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x4c(%rbp)
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x4c(%rbp),%eax
jg     1303 <cffti1_ps+0x1d3>
mov    -0x30(%rbp),%eax
mov    %eax,-0x54(%rbp)
movl   $0x0,-0x58(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  0xdf5(%rip),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x3c(%rbp)
cvtsi2ssl -0x3c(%rbp),%xmm0
mulss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x50(%rbp)
movl   $0x4,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x48(%rbp),%eax
jg     12a8 <cffti1_ps+0x178>
mov    -0x30(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x30(%rbp)
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
cvtsi2ssl -0x58(%rbp),%xmm0
mulss  -0x50(%rbp),%xmm0
call   1040 <cos@plt>
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
cvtsi2ssl -0x58(%rbp),%xmm0
mulss  -0x50(%rbp),%xmm0
call   1030 <sin@plt>
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x24(%rbp)
jmp    123f <cffti1_ps+0x10f>
cmpl   $0x5,-0x38(%rbp)
jle    12f0 <cffti1_ps+0x1c0>
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x54(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x54(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    12f5 <cffti1_ps+0x1c5>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    11dd <cffti1_ps+0xad>
mov    -0x40(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1187 <cffti1_ps+0x57>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cffti1_ps>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %rdx,%rax
mov    %rsi,%rbx
mov    %edi,%ebp
lea    0xeb3(%rip),%rdx
mov    %rax,0x30(%rsp)
mov    %rax,%rsi
call   1050 <decompose@plt>
mov    %eax,%ecx
mov    0x2e85(%rip),%rax
mov    (%rax),%eax
add    %eax,%eax
cltd
mov    %ebp,0x10(%rsp)
idiv   %ebp
mov    %rcx,0x8(%rsp)
test   %ecx,%ecx
jle    12e9 <cffti1_ps+0x1b9>
cvtsi2ss %eax,%xmm0
movss  %xmm0,0x18(%rsp)
mov    0x8(%rsp),%rax
add    $0x1,%eax
mov    %rax,0x8(%rsp)
mov    $0x1,%ebp
mov    $0x1,%ecx
mov    $0x1,%eax
mov    %rbx,0x40(%rsp)
jmp    11c4 <cffti1_ps+0x94>
nopl   0x0(%rax,%rax,1)
mov    0x38(%rsp),%rcx
cmp    0x8(%rsp),%rcx
mov    0x14(%rsp),%eax
je     12e9 <cffti1_ps+0x1b9>
mov    %eax,%edi
mov    0x30(%rsp),%rax
mov    0x4(%rax,%rcx,4),%r8d
add    $0x1,%rcx
mov    %r8d,%esi
imul   %edi,%esi
mov    0x10(%rsp),%eax
cltd
idiv   %esi
cmp    $0x2,%r8d
mov    %rcx,0x38(%rsp)
mov    %esi,0x14(%rsp)
jl     11b0 <cffti1_ps+0x80>
lea    (%rax,%rax,1),%r14d
xor    %ecx,%ecx
mov    $0x1,%r13d
mov    %edi,0x20(%rsp)
mov    %r8d,0x1c(%rsp)
mov    %rax,0x48(%rsp)
jmp    1226 <cffti1_ps+0xf6>
nopl   0x0(%rax)
add    $0x1,%r13d
cmp    %edx,%r13d
mov    0x20(%rsp),%edi
mov    0x48(%rsp),%rax
mov    0x24(%rsp),%ecx
je     11b0 <cffti1_ps+0x80>
movslq %ebp,%rdx
add    %edi,%ecx
movl   $0x3f800000,-0x4(%rbx,%rdx,4)
3f
movl   $0x0,(%rbx,%rdx,4)
test   %eax,%eax
mov    %ecx,0x24(%rsp)
mov    %rdx,0x50(%rsp)
jle    12b1 <cffti1_ps+0x181>
xorps  %xmm0,%xmm0
cvtsi2ss %ecx,%xmm0
mulss  0x18(%rsp),%xmm0
movss  %xmm0,0x28(%rsp)
lea    (%rbx,%rdx,4),%r15
mov    $0x2,%ebx
mov    $0x1,%r12d
nopl   0x0(%rax)
xorps  %xmm0,%xmm0
cvtsi2ss %r12d,%xmm0
mulss  0x28(%rsp),%xmm0
movss  %xmm0,0x2c(%rsp)
call   1040 <cos@plt>
movss  %xmm0,-0x4(%r15,%rbx,4)
movss  0x2c(%rsp),%xmm0
call   1030 <sin@plt>
movss  %xmm0,(%r15,%rbx,4)
add    $0x2,%rbx
add    $0x2,%ebp
add    $0x1,%r12d
cmp    %r14d,%ebx
jle    1270 <cffti1_ps+0x140>
mov    0x1c(%rsp),%edx
cmp    $0x6,%edx
mov    0x40(%rsp),%rbx
jb     1210 <cffti1_ps+0xe0>
mov    0x50(%rsp),%rsi
lea    -0x1(%rsi),%rax
movslq %ebp,%rcx
movss  -0x4(%rbx,%rcx,4),%xmm0
movss  %xmm0,(%rbx,%rax,4)
movss  (%rbx,%rcx,4),%xmm0
movss  %xmm0,(%rbx,%rsi,4)
jmp    1210 <cffti1_ps+0xe0>
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
@cffti1_ps.ntryh = internal constant [5 x i32] [i32 5, i32 3, i32 4, i32 2, i32 0]
@M_PI = global i32 0
define void @cffti1_ps(i32 %0, float* nocapture %1, i32* %2) {
%4 = tail call i32 @decompose(i32 %0, i32* %2, i32* getelementptr inbounds ([5 x i32], [5 x i32]* @cffti1_ps.ntryh, i64 0, i64 0))
%5 = load i32, i32* @M_PI
%6 = shl nsw i32 %5, 1
%7 = sdiv i32 %6, %0
%8 = sitofp i32 %7 to float
%9 = icmp slt i32 %4, 1
br i1 %9, label %72, label %10
10: ; preds = %3
%11 = add nuw i32 %4, 1
%12 = zext i32 %11 to i64
br label %16
13: ; preds = %69, %16
%14 = phi i32 [ %18, %16 ], [ %60, %69 ]
%15 = icmp eq i64 %20, %12
br i1 %15, label %72, label %16
16: ; preds = %10, %13
%17 = phi i64 [ 1, %10 ], [ %20, %13 ]
%18 = phi i32 [ 1, %10 ], [ %14, %13 ]
%19 = phi i32 [ 1, %10 ], [ %23, %13 ]
%20 = add nuw nsw i64 %17, 1
%21 = getelementptr inbounds i32, i32* %2, i64 %20
%22 = load i32, i32* %21
%23 = mul nsw i32 %22, %19
%24 = sdiv i32 %0, %23
%25 = shl nsw i32 %24, 1
%26 = icmp sgt i32 %22, 1
br i1 %26, label %27, label %13
27: ; preds = %16
%28 = icmp slt i32 %24, 1
%29 = icmp ugt i32 %22, 5
br label %30
30: ; preds = %27, %69
%31 = phi i32 [ 1, %27 ], [ %70, %69 ]
%32 = phi i32 [ %18, %27 ], [ %60, %69 ]
%33 = phi i32 [ 0, %27 ], [ %39, %69 ]
%34 = add nsw i32 %32, -1
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds float, float* %1, i64 %35
store float 1.000000e+00, float* %36
%37 = sext i32 %32 to i64
%38 = getelementptr inbounds float, float* %1, i64 %37
store float 0.000000e+00, float* %38
%39 = add nsw i32 %33, %19
%40 = sitofp i32 %39 to float
%41 = fmul float %8, %40
br i1 %28, label %59, label %42
42: ; preds = %30, %42
%43 = phi i64 [ %46, %42 ], [ %37, %30 ]
%44 = phi i32 [ %47, %42 ], [ 0, %30 ]
%45 = phi i32 [ %55, %42 ], [ 4, %30 ]
%46 = add nsw i64 %43, 2
%47 = add nuw nsw i32 %44, 1
%48 = sitofp i32 %47 to float
%49 = fmul float %41, %48
%50 = tail call float @cos(float %49)
%51 = add nsw i64 %43, 1
%52 = getelementptr inbounds float, float* %1, i64 %51
store float %50, float* %52
%53 = tail call float @sin(float %49)
%54 = getelementptr inbounds float, float* %1, i64 %46
store float %53, float* %54
%55 = add nuw nsw i32 %45, 2
%56 = icmp sgt i32 %45, %25
br i1 %56, label %57, label %42
57: ; preds = %42
%58 = trunc i64 %46 to i32
br label %59
59: ; preds = %57, %30
%60 = phi i32 [ %32, %30 ], [ %58, %57 ]
br i1 %29, label %61, label %69
61: ; preds = %59
%62 = add nsw i32 %60, -1
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds float, float* %1, i64 %63
%65 = load float, float* %64
store float %65, float* %36
%66 = sext i32 %60 to i64
%67 = getelementptr inbounds float, float* %1, i64 %66
%68 = load float, float* %67
store float %68, float* %38
br label %69
69: ; preds = %61, %59
%70 = add nuw nsw i32 %31, 1
%71 = icmp eq i32 %70, %22
br i1 %71, label %13, label %30
72: ; preds = %13, %3
ret void
}