#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 */
double cos(double);
int /*<<< orphan*/ factorize(size_t, size_t const *, int, size_t *);
double sin(double);
void rffti1(size_t n, double wa[], size_t ifac[]) {
  static const size_t ntryh[4] = {4, 2, 3, 5};
  static const double twopi = 6.28318530717958647692;
  size_t i, j, k, fi;
  double argh = twopi / n;
  size_t is = 0, l1 = 1;
  factorize(n, ntryh, 4, ifac);
  for (k = 1; k < ifac[1]; k++) {
    size_t ip = ifac[k + 1], ido = n / (l1 * ip);
    for (j = 1; j < ip; ++j) {
      double argld = j * l1 * argh;
      for (i = is, fi = 1; i <= ido + is - 3; i += 2, ++fi) {
        double arg = fi * argld;
        wa[i] = cos(arg);
        wa[i + 1] = sin(arg);
      }
      is += ido;
    }
    l1 *= ip;
  }
}


x86 assembly (clang, -O0):
<rffti1>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movaps 0xec5(%rip),%xmm0
movq   -0x8(%rbp),%xmm1
punpckldq %xmm0,%xmm1
movapd 0xec4(%rip),%xmm0
subpd  %xmm0,%xmm1
movaps %xmm1,%xmm0
unpckhpd %xmm1,%xmm1
addsd  %xmm0,%xmm1
movsd  0xe8d(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movq   $0x0,-0x48(%rbp)
movq   $0x1,-0x50(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0x18(%rbp),%rcx
lea    0xe95(%rip),%rsi
mov    $0x4,%edx
call   1050 <factorize@plt>
movq   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%rax
mov    -0x18(%rbp),%rcx
cmp    0x8(%rcx),%rax
jae    131a <rffti1+0x1ea>
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    0x8(%rax,%rcx,8),%rax
mov    %rax,-0x58(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%rcx
imul   -0x58(%rbp),%rcx
xor    %edx,%edx
div    %rcx
mov    %rax,-0x60(%rbp)
movq   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x58(%rbp),%rax
jae    12fc <rffti1+0x1cc>
mov    -0x28(%rbp),%rax
mov    -0x50(%rbp),%rcx
imul   %rcx,%rax
movq   %rax,%xmm0
movaps 0xdfc(%rip),%xmm1
punpckldq %xmm1,%xmm0
movapd 0xe00(%rip),%xmm1
subpd  %xmm1,%xmm0
movaps %xmm0,%xmm1
unpckhpd %xmm0,%xmm0
addsd  %xmm1,%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x48(%rbp),%rax
mov    %rax,-0x20(%rbp)
movq   $0x1,-0x38(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x60(%rbp),%rcx
add    -0x48(%rbp),%rcx
sub    $0x3,%rcx
cmp    %rcx,%rax
ja     12df <rffti1+0x1af>
movaps 0xda7(%rip),%xmm1
movq   -0x38(%rbp),%xmm0
punpckldq %xmm1,%xmm0
movapd 0xda6(%rip),%xmm1
subpd  %xmm1,%xmm0
movaps %xmm0,%xmm1
unpckhpd %xmm0,%xmm0
addsd  %xmm1,%xmm0
mulsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
call   1040 <cos@plt>
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
call   1030 <sin@plt>
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
movsd  %xmm0,0x8(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
add    $0x2,%rax
mov    %rax,-0x20(%rbp)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    1249 <rffti1+0x119>
mov    -0x60(%rbp),%rax
add    -0x48(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    11ee <rffti1+0xbe>
mov    -0x58(%rbp),%rax
imul   -0x50(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    11ad <rffti1+0x7d>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<rffti1>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %rdx,%r14
mov    %rsi,%rbx
mov    %rdi,%rbp
lea    0xee2(%rip),%rsi
mov    $0x4,%edx
mov    %r14,%rcx
call   1050 <factorize@plt>
mov    0x8(%r14),%r8
cmp    $0x2,%r8
jb     1312 <rffti1+0x1e2>
movq   %rbp,%xmm0
punpckldq 0xe8a(%rip),%xmm0
subpd  0xe92(%rip),%xmm0
movapd %xmm0,%xmm1
unpckhpd %xmm0,%xmm1
addsd  %xmm0,%xmm1
movsd  0xe8e(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,0x20(%rsp)
mov    $0x1,%esi
xor    %r12d,%r12d
movdqa 0xe54(%rip),%xmm1
movapd 0xe5c(%rip),%xmm2
mov    $0x1,%ecx
mov    %r14,0x8(%rsp)
mov    %rbp,(%rsp)
jmp    11f0 <rffti1+0xc0>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    0x8(%rsp),%r14
mov    0x8(%r14),%r8
mov    (%rsp),%rbp
mov    0x18(%rsp),%rsi
mov    0x10(%rsp),%rcx
cmp    %r8,%rsi
jae    1312 <rffti1+0x1e2>
mov    %rcx,%rax
mov    0x8(%r14,%rsi,8),%rdi
mov    %rdi,%rcx
mov    %rax,0x38(%rsp)
imul   %rax,%rcx
mov    %rbp,%rax
or     %rcx,%rax
shr    $0x20,%rax
je     1220 <rffti1+0xf0>
mov    %rbp,%rax
xor    %edx,%edx
div    %rcx
jmp    1226 <rffti1+0xf6>
nopw   0x0(%rax,%rax,1)
mov    %ebp,%eax
xor    %edx,%edx
div    %ecx
add    $0x1,%rsi
cmp    $0x2,%rdi
jb     11e7 <rffti1+0xb7>
mov    %rcx,0x10(%rsp)
mov    %rsi,0x18(%rsp)
mov    $0x1,%r13d
mov    %rdi,0x30(%rsp)
mov    %rax,0x28(%rsp)
jmp    126a <rffti1+0x13a>
nopl   0x0(%rax)
mov    0x28(%rsp),%rax
add    %rax,%r12
add    $0x1,%r13
mov    0x30(%rsp),%rdi
cmp    %rdi,%r13
je     11d0 <rffti1+0xa0>
mov    %r13,%rcx
imul   0x38(%rsp),%rcx
movq   %rcx,%xmm0
punpckldq %xmm1,%xmm0
subpd  %xmm2,%xmm0
lea    (%r12,%rax,1),%r14
add    $0xfffffffffffffffd,%r14
cmp    %r14,%r12
ja     1250 <rffti1+0x120>
movapd %xmm0,%xmm3
unpckhpd %xmm0,%xmm3
addsd  %xmm0,%xmm3
mulsd  0x20(%rsp),%xmm3
movapd %xmm3,0x40(%rsp)
mov    $0x1,%r15d
mov    %r12,%rbp
xchg   %ax,%ax
movq   %r15,%xmm3
punpckldq %xmm1,%xmm3
subpd  %xmm2,%xmm3
movapd %xmm3,%xmm0
unpckhpd %xmm3,%xmm0
addsd  %xmm3,%xmm0
mulsd  0x40(%rsp),%xmm0
movapd %xmm0,0x50(%rsp)
call   1040 <cos@plt>
movsd  %xmm0,(%rbx,%rbp,8)
movapd 0x50(%rsp),%xmm0
call   1030 <sin@plt>
movapd 0xd1e(%rip),%xmm2
movdqa 0xd06(%rip),%xmm1
movsd  %xmm0,0x8(%rbx,%rbp,8)
add    $0x2,%rbp
add    $0x1,%r15
cmp    %r14,%rbp
jbe    12b0 <rffti1+0x180>
jmp    1250 <rffti1+0x120>
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@rffti1.ntryh = internal constant [4 x i64] [i64 4, i64 2, i64 3, i64 5]
@rffti1.twopi = internal constant double 0x401921FB54442D18
define void @rffti1(i64 %0, double* %1, i64* %2) {
%4 = alloca i64
%5 = alloca double*
%6 = alloca i64*
%7 = alloca i64
%8 = alloca i64
%9 = alloca i64
%10 = alloca i64
%11 = alloca double
%12 = alloca i64
%13 = alloca i64
%14 = alloca i64
%15 = alloca i64
%16 = alloca double
%17 = alloca double
store i64 %0, i64* %4
store double* %1, double** %5
store i64* %2, i64** %6
%18 = load i64, i64* %4
%19 = uitofp i64 %18 to double
%20 = fdiv double 0x401921FB54442D18, %19
store double %20, double* %11
store i64 0, i64* %12
store i64 1, i64* %13
%21 = load i64, i64* %4
%22 = load i64*, i64** %6
%23 = call i32 @factorize(i64 %21, i64* getelementptr inbounds ([4 x i64], [4 x i64]* @rffti1.ntryh, i64 0, i64 0), i32 4, i64* %22)
store i64 1, i64* %9
br label %24
24: ; preds = %92, %3
%25 = load i64, i64* %9
%26 = load i64*, i64** %6
%27 = getelementptr inbounds i64, i64* %26, i64 1
%28 = load i64, i64* %27
%29 = icmp ult i64 %25, %28
br i1 %29, label %30, label %95
30: ; preds = %24
%31 = load i64*, i64** %6
%32 = load i64, i64* %9
%33 = add i64 %32, 1
%34 = getelementptr inbounds i64, i64* %31, i64 %33
%35 = load i64, i64* %34
store i64 %35, i64* %14
%36 = load i64, i64* %4
%37 = load i64, i64* %13
%38 = load i64, i64* %14
%39 = mul i64 %37, %38
%40 = udiv i64 %36, %39
store i64 %40, i64* %15
store i64 1, i64* %8
br label %41
41: ; preds = %85, %30
%42 = load i64, i64* %8
%43 = load i64, i64* %14
%44 = icmp ult i64 %42, %43
br i1 %44, label %45, label %88
45: ; preds = %41
%46 = load i64, i64* %8
%47 = load i64, i64* %13
%48 = mul i64 %46, %47
%49 = uitofp i64 %48 to double
%50 = load double, double* %11
%51 = fmul double %49, %50
store double %51, double* %16
%52 = load i64, i64* %12
store i64 %52, i64* %7
store i64 1, i64* %10
br label %53
53: ; preds = %76, %45
%54 = load i64, i64* %7
%55 = load i64, i64* %15
%56 = load i64, i64* %12
%57 = add i64 %55, %56
%58 = sub i64 %57, 3
%59 = icmp ule i64 %54, %58
br i1 %59, label %60, label %81
60: ; preds = %53
%61 = load i64, i64* %10
%62 = uitofp i64 %61 to double
%63 = load double, double* %16
%64 = fmul double %62, %63
store double %64, double* %17
%65 = load double, double* %17
%66 = call double @cos(double %65)
%67 = load double*, double** %5
%68 = load i64, i64* %7
%69 = getelementptr inbounds double, double* %67, i64 %68
store double %66, double* %69
%70 = load double, double* %17
%71 = call double @sin(double %70)
%72 = load double*, double** %5
%73 = load i64, i64* %7
%74 = add i64 %73, 1
%75 = getelementptr inbounds double, double* %72, i64 %74
store double %71, double* %75
br label %76
76: ; preds = %60
%77 = load i64, i64* %7
%78 = add i64 %77, 2
store i64 %78, i64* %7
%79 = load i64, i64* %10
%80 = add i64 %79, 1
store i64 %80, i64* %10
br label %53
81: ; preds = %53
%82 = load i64, i64* %15
%83 = load i64, i64* %12
%84 = add i64 %83, %82
store i64 %84, i64* %12
br label %85
85: ; preds = %81
%86 = load i64, i64* %8
%87 = add i64 %86, 1
store i64 %87, i64* %8
br label %41
88: ; preds = %41
%89 = load i64, i64* %14
%90 = load i64, i64* %13
%91 = mul i64 %90, %89
store i64 %91, i64* %13
br label %92
92: ; preds = %88
%93 = load i64, i64* %9
%94 = add i64 %93, 1
store i64 %94, i64* %9
br label %24
95: ; preds = %24
ret void
}

LLVM IR (-O2):
@rffti1.ntryh = internal constant [4 x i64] [i64 4, i64 2, i64 3, i64 5]
define void @rffti1(i64 %0, double* nocapture writeonly %1, i64* %2) {
%4 = uitofp i64 %0 to double
%5 = fdiv double 0x401921FB54442D18, %4
%6 = tail call i32 @factorize(i64 %0, i64* getelementptr inbounds ([4 x i64], [4 x i64]* @rffti1.ntryh, i64 0, i64 0), i32 4, i64* %2)
%7 = getelementptr inbounds i64, i64* %2, i64 1
%8 = load i64, i64* %7
%9 = icmp ugt i64 %8, 1
br i1 %9, label %16, label %52
10: ; preds = %49
%11 = load i64, i64* %7
br label %12
12: ; preds = %10, %16
%13 = phi i64 [ %17, %16 ], [ %11, %10 ]
%14 = phi i64 [ %20, %16 ], [ %33, %10 ]
%15 = icmp ult i64 %21, %13
br i1 %15, label %16, label %52
16: ; preds = %3, %12
%17 = phi i64 [ %13, %12 ], [ %8, %3 ]
%18 = phi i64 [ %21, %12 ], [ 1, %3 ]
%19 = phi i64 [ %24, %12 ], [ 1, %3 ]
%20 = phi i64 [ %14, %12 ], [ 0, %3 ]
%21 = add nuw i64 %18, 1
%22 = getelementptr inbounds i64, i64* %2, i64 %21
%23 = load i64, i64* %22
%24 = mul i64 %23, %19
%25 = udiv i64 %0, %24
%26 = icmp ugt i64 %23, 1
br i1 %26, label %27, label %12
27: ; preds = %16, %49
%28 = phi i64 [ %50, %49 ], [ 1, %16 ]
%29 = phi i64 [ %33, %49 ], [ %20, %16 ]
%30 = mul i64 %28, %19
%31 = uitofp i64 %30 to double
%32 = fmul double %5, %31
%33 = add i64 %29, %25
%34 = add i64 %33, -3
%35 = icmp ugt i64 %29, %34
br i1 %35, label %49, label %36
36: ; preds = %27, %36
%37 = phi i64 [ %46, %36 ], [ %29, %27 ]
%38 = phi i64 [ %47, %36 ], [ 1, %27 ]
%39 = uitofp i64 %38 to double
%40 = fmul double %32, %39
%41 = tail call double @cos(double %40)
%42 = getelementptr inbounds double, double* %1, i64 %37
store double %41, double* %42
%43 = tail call double @sin(double %40)
%44 = add i64 %37, 1
%45 = getelementptr inbounds double, double* %1, i64 %44
store double %43, double* %45
%46 = add i64 %37, 2
%47 = add i64 %38, 1
%48 = icmp ugt i64 %46, %34
br i1 %48, label %49, label %36
49: ; preds = %36, %27
%50 = add nuw i64 %28, 1
%51 = icmp eq i64 %50, %23
br i1 %51, label %10, label %27
52: ; preds = %12, %3
ret void
}