#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 sin(double);
void sinft(float y[], int n) {
  void realft(float data[], unsigned long n, int isign);
  int j, n2 = n + 2;
  float sum, y1, y2;
  double theta, wi = 0.0, wr = 1.0, wpi, wpr, wtemp;
  theta = 3.14159265358979 / (double)n;
  wtemp = sin(0.5 * theta);
  wpr = -2.0 * wtemp * wtemp;
  wpi = sin(theta);
  y[1] = 0.0;
  for (j = 2; j <= (n >> 1) + 1; j++) {
    wr = (wtemp = wr) * wpr - wi * wpi + wr;
    wi = wi * wpr + wtemp * wpi + wi;
    y1 = wi * (y[j] + y[n2 - j]);
    y2 = 0.5 * (y[j] - y[n2 - j]);
    y[j] = y1 + y2;
    y[n2 - j] = y1 - y2;
  }
  realft(y, n, 1);
  y[1] *= 0.5;
  sum = y[2] = 0.0;
  for (j = 1; j <= n - 1; j += 2) {
    sum += y[j];
    y[j] = y[j + 1];
    y[j + 1] = sum;
  }
}


x86 assembly (clang, -O0):
<sinft>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    -0xc(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x14(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  0xed0(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
cvtsi2sdl -0xc(%rbp),%xmm1
movsd  0xeb6(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  0xe9d(%rip),%xmm0
mulsd  -0x28(%rbp),%xmm0
call   1040 <sin@plt>
movsd  %xmm0,-0x50(%rbp)
movsd  0xe7e(%rip),%xmm0
mulsd  -0x50(%rbp),%xmm0
mulsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x28(%rbp),%xmm0
call   1040 <sin@plt>
movsd  %xmm0,-0x40(%rbp)
mov    -0x8(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,0x4(%rax)
movl   $0x2,-0x10(%rbp)
mov    -0x10(%rbp),%eax
mov    -0xc(%rbp),%ecx
sar    $0x1,%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     12dd <sinft+0x1bd>
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x48(%rbp),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  -0x40(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    -0x10(%rbp),%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x1c(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    -0x10(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd79(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x1c(%rbp),%xmm0
addss  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x1c(%rbp),%xmm0
subss  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    -0x10(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    11b3 <sinft+0x93>
mov    -0x8(%rbp),%rdi
movslq -0xc(%rbp),%rsi
mov    $0x1,%edx
call   1030 <realft@plt>
mov    -0x8(%rbp),%rax
movss  0x4(%rax),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xd04(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x4(%rax)
mov    -0x8(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,0x8(%rax)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x18(%rbp)
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     1398 <sinft+0x278>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x18(%rbp),%xmm0
movss  %xmm0,-0x18(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x18(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x10(%rbp)
jmp    132c <sinft+0x20c>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<sinft>:
push   %rbp
push   %r14
push   %rbx
sub    $0x10,%rsp
mov    %esi,%r14d
mov    %rdi,%rbx
cvtsi2sd %esi,%xmm0
movsd  0xec6(%rip),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsp)
movsd  0xebc(%rip),%xmm0
mulsd  %xmm1,%xmm0
call   1040 <sin@plt>
movsd  %xmm0,(%rsp)
movsd  0x8(%rsp),%xmm0
call   1040 <sin@plt>
movl   $0x0,0x4(%rbx)
cmp    $0x2,%r14d
jl     1238 <sinft+0x118>
movsd  0xe92(%rip),%xmm1
movsd  (%rsp),%xmm2
mulsd  %xmm2,%xmm1
mulsd  %xmm2,%xmm1
mov    %r14d,%eax
shr    %eax
unpcklpd %xmm1,%xmm1
unpcklpd %xmm0,%xmm0
movsd  0xe78(%rip),%xmm3
xor    %ecx,%ecx
movapd 0xe75(%rip),%xmm8
movss  0xe7d(%rip),%xmm4
mov    %r14d,%edx
cs nopw 0x0(%rax,%rax,1)
movapd %xmm3,%xmm5
unpckhpd %xmm3,%xmm5
xorpd  %xmm8,%xmm5
unpcklpd %xmm3,%xmm5
mulpd  %xmm0,%xmm5
movapd %xmm3,%xmm6
mulpd  %xmm1,%xmm6
addpd  %xmm5,%xmm6
addpd  %xmm6,%xmm3
movss  0x8(%rbx,%rcx,4),%xmm5
movslq %edx,%rdx
movss  (%rbx,%rdx,4),%xmm6
movaps %xmm5,%xmm7
addss  %xmm6,%xmm7
cvtss2sd %xmm7,%xmm7
movapd %xmm3,%xmm2
unpckhpd %xmm3,%xmm2
mulsd  %xmm7,%xmm2
cvtsd2ss %xmm2,%xmm2
subss  %xmm6,%xmm5
mulss  %xmm4,%xmm5
movaps %xmm5,%xmm6
addss  %xmm2,%xmm6
movss  %xmm6,0x8(%rbx,%rcx,4)
subss  %xmm5,%xmm2
movss  %xmm2,(%rbx,%rdx,4)
add    $0xffffffff,%edx
add    $0x1,%rcx
cmp    %rcx,%rax
jne    11c0 <sinft+0xa0>
movslq %r14d,%rbp
mov    %rbx,%rdi
mov    %rbp,%rsi
mov    $0x1,%edx
call   1030 <realft@plt>
movss  0x4(%rbx),%xmm0
mulss  0xdd8(%rip),%xmm0
movss  %xmm0,0x4(%rbx)
movl   $0x0,0x8(%rbx)
cmp    $0x2,%r14d
jl     1303 <sinft+0x1e3>
xorpd  %xmm1,%xmm1
addss  %xmm1,%xmm0
movss  0x8(%rbx),%xmm1
movss  %xmm1,0x4(%rbx)
movss  %xmm0,0x8(%rbx)
cmp    $0x4,%r14d
jl     1303 <sinft+0x1e3>
add    $0xfffffffffffffffc,%rbp
mov    %rbp,%rcx
shr    %rcx
add    $0x1,%rcx
mov    $0x3,%eax
cmp    $0x2,%rbp
jb     12e8 <sinft+0x1c8>
mov    %rcx,%rdx
and    $0xfffffffffffffffe,%rdx
mov    $0x3,%eax
addss  (%rbx,%rax,4),%xmm0
movss  0x4(%rbx,%rax,4),%xmm1
movss  0xc(%rbx,%rax,4),%xmm2
movss  %xmm1,(%rbx,%rax,4)
movss  %xmm0,0x4(%rbx,%rax,4)
addss  0x8(%rbx,%rax,4),%xmm0
movss  %xmm2,0x8(%rbx,%rax,4)
movss  %xmm0,0xc(%rbx,%rax,4)
add    $0x4,%rax
add    $0xfffffffffffffffe,%rdx
jne    12b0 <sinft+0x190>
test   $0x1,%cl
je     1303 <sinft+0x1e3>
addss  (%rbx,%rax,4),%xmm0
movss  0x4(%rbx,%rax,4),%xmm1
movss  %xmm1,(%rbx,%rax,4)
movss  %xmm0,0x4(%rbx,%rax,4)
add    $0x10,%rsp
pop    %rbx
pop    %r14
pop    %rbp
ret

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

LLVM IR (-O2):
define void @sinft(float* %0, i32 %1) {
%3 = add nsw i32 %1, 2
%4 = sitofp i32 %1 to double
%5 = fdiv double 0x400921FB54442D11, %4
%6 = fmul double %5, 5.000000e-01
%7 = tail call double @sin(double %6)
%8 = tail call double @sin(double %5)
%9 = getelementptr inbounds float, float* %0, i64 1
store float 0.000000e+00, float* %9
%10 = icmp slt i32 %1, 2
br i1 %10, label %49, label %11
11: ; preds = %2
%12 = fmul double %7, -2.000000e+00
%13 = fmul double %7, %12
%14 = lshr i32 %1, 1
%15 = add nuw nsw i32 %14, 2
%16 = zext i32 %15 to i64
%17 = insertelement <2 x double> poison, double %13, i64 0
%18 = shufflevector <2 x double> %17, <2 x double> poison, <2 x i32> zeroinitializer
%19 = insertelement <2 x double> poison, double %8, i64 0
%20 = shufflevector <2 x double> %19, <2 x double> poison, <2 x i32> zeroinitializer
br label %21
21: ; preds = %11, %21
%22 = phi i64 [ 2, %11 ], [ %47, %21 ]
%23 = phi <2 x double> [ <double 1.000000e+00, double 0.000000e+00>, %11 ], [ %30, %21 ]
%24 = extractelement <2 x double> %23, i64 1
%25 = fneg double %24
%26 = insertelement <2 x double> poison, double %25, i64 0
%27 = shufflevector <2 x double> %26, <2 x double> %23, <2 x i32> <i32 0, i32 2>
%28 = fmul <2 x double> %20, %27
%29 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %23, <2 x double> %18, <2 x double> %28)
%30 = fadd <2 x double> %23, %29
%31 = getelementptr inbounds float, float* %0, i64 %22
%32 = load float, float* %31
%33 = trunc i64 %22 to i32
%34 = sub i32 %3, %33
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds float, float* %0, i64 %35
%37 = load float, float* %36
%38 = fadd float %32, %37
%39 = fpext float %38 to double
%40 = extractelement <2 x double> %30, i64 1
%41 = fmul double %40, %39
%42 = fptrunc double %41 to float
%43 = fsub float %32, %37
%44 = fmul float %43, 5.000000e-01
%45 = fadd float %44, %42
store float %45, float* %31
%46 = fsub float %42, %44
store float %46, float* %36
%47 = add nuw nsw i64 %22, 1
%48 = icmp eq i64 %47, %16
br i1 %48, label %49, label %21
49: ; preds = %21, %2
%50 = sext i32 %1 to i64
tail call void @realft(float* %0, i64 %50, i32 1)
%51 = load float, float* %9
%52 = fmul float %51, 5.000000e-01
store float %52, float* %9
%53 = getelementptr inbounds float, float* %0, i64 2
store float 0.000000e+00, float* %53
%54 = icmp sgt i32 %1, 1
br i1 %54, label %55, label %103
55: ; preds = %49
%56 = getelementptr inbounds float, float* %0, i64 1
%57 = fadd float %52, 0.000000e+00
%58 = getelementptr inbounds float, float* %0, i64 2
%59 = load float, float* %58
store float %59, float* %56
store float %57, float* %58
%60 = icmp sgt i32 %1, 3
br i1 %60, label %61, label %103
61: ; preds = %55
%62 = add nsw i64 %50, -4
%63 = lshr i64 %62, 1
%64 = add nuw i64 %63, 1
%65 = and i64 %64, 1
%66 = icmp ult i64 %62, 2
br i1 %66, label %91, label %67
67: ; preds = %61
%68 = and i64 %64, -2
br label %69
69: ; preds = %69, %67
%70 = phi i64 [ 3, %67 ], [ %88, %69 ]
%71 = phi float [ %57, %67 ], [ %84, %69 ]
%72 = phi i64 [ 0, %67 ], [ %89, %69 ]
%73 = getelementptr inbounds float, float* %0, i64 %70
%74 = load float, float* %73
%75 = getelementptr inbounds float, float* %0, i64 %70
%76 = fadd float %71, %74
%77 = add nuw nsw i64 %70, 1
%78 = getelementptr inbounds float, float* %0, i64 %77
%79 = load float, float* %78
store float %79, float* %75
store float %76, float* %78
%80 = add nuw nsw i64 %70, 2
%81 = getelementptr inbounds float, float* %0, i64 %80
%82 = load float, float* %81
%83 = getelementptr inbounds float, float* %0, i64 %80
%84 = fadd float %76, %82
%85 = add nuw nsw i64 %70, 3
%86 = getelementptr inbounds float, float* %0, i64 %85
%87 = load float, float* %86
store float %87, float* %83
store float %84, float* %86
%88 = add nuw nsw i64 %70, 4
%89 = add nuw i64 %72, 2
%90 = icmp eq i64 %89, %68
br i1 %90, label %91, label %69
91: ; preds = %69, %61
%92 = phi i64 [ 3, %61 ], [ %88, %69 ]
%93 = phi float [ %57, %61 ], [ %84, %69 ]
%94 = icmp eq i64 %65, 0
br i1 %94, label %103, label %95
95: ; preds = %91
%96 = getelementptr inbounds float, float* %0, i64 %92
%97 = load float, float* %96
%98 = getelementptr inbounds float, float* %0, i64 %92
%99 = fadd float %93, %97
%100 = add nuw nsw i64 %92, 1
%101 = getelementptr inbounds float, float* %0, i64 %100
%102 = load float, float* %101
store float %102, float* %98
store float %99, float* %101
br label %103
103: ; preds = %95, %91, %55, %49
ret void
}