#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 /*<<< orphan*/ assert(int);
int /*<<< orphan*/ autocorrelate(float *, float *, int, int);
int /*<<< orphan*/ hanning_window(float *, float *, int);
int /*<<< orphan*/ levinson_durbin(float *, float *, int);
void find_aks(float Sn[], float a[], int Nsam, int order, float *E) {
  float Wn[512];
  float R[20 + 1];
  int i;
  assert(order < 20);
  assert(Nsam < 512);
  hanning_window(Sn, Wn, Nsam);
  autocorrelate(Wn, R, Nsam, order);
  levinson_durbin(R, a, order);
  *E = 0.0;
  for (i = 0; i <= order; i++)
    *E += a[i] * R[i];
  if (*E < 0.0)
    *E = 1E-12;
}


x86 assembly (clang, -O0):
<find_aks>:
push   %rbp
mov    %rsp,%rbp
sub    $0x890,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
cmpl   $0x14,-0x18(%rbp)
setl   %al
and    $0x1,%al
movzbl %al,%edi
call   1040 <assert@plt>
cmpl   $0x200,-0x14(%rbp)
setl   %al
and    $0x1,%al
movzbl %al,%edi
call   1040 <assert@plt>
mov    -0x8(%rbp),%rdi
lea    -0x820(%rbp),%rsi
mov    -0x14(%rbp),%edx
call   1030 <hanning_window@plt>
lea    -0x820(%rbp),%rdi
lea    -0x880(%rbp),%rsi
mov    -0x14(%rbp),%edx
mov    -0x18(%rbp),%ecx
call   1060 <autocorrelate@plt>
lea    -0x880(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    -0x18(%rbp),%edx
call   1050 <levinson_durbin@plt>
mov    -0x20(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax)
movl   $0x0,-0x884(%rbp)
mov    -0x884(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     122d <find_aks+0xed>
mov    -0x10(%rbp),%rax
movslq -0x884(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movslq -0x884(%rbp),%rax
movss  -0x880(%rbp,%rax,4),%xmm2
ff ff
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax)
mov    -0x884(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x884(%rbp)
jmp    11d6 <find_aks+0x96>
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jbe    1256 <find_aks+0x116>
mov    -0x20(%rbp),%rax
movss  0xdae(%rip),%xmm0
movss  %xmm0,(%rax)
add    $0x890,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<find_aks>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x868,%rsp
mov    %r8,%rbx
mov    %ecx,%ebp
mov    %edx,%r15d
mov    %rsi,%r14
mov    %rdi,%r12
xor    %edi,%edi
cmp    $0x14,%ecx
setl   %dil
call   1040 <assert@plt>
xor    %edi,%edi
cmp    $0x200,%r15d
setl   %dil
call   1040 <assert@plt>
lea    0x60(%rsp),%r13
mov    %r12,%rdi
mov    %r13,%rsi
mov    %r15d,%edx
call   1030 <hanning_window@plt>
mov    %rsp,%r12
mov    %r13,%rdi
mov    %r12,%rsi
mov    %r15d,%edx
mov    %ebp,%ecx
call   1060 <autocorrelate@plt>
mov    %r12,%rdi
mov    %r14,%rsi
mov    %ebp,%edx
call   1050 <levinson_durbin@plt>
movl   $0x0,(%rbx)
test   %ebp,%ebp
js     1279 <find_aks+0x139>
mov    %ebp,%ecx
add    $0x1,%rcx
mov    %ecx,%eax
and    $0x3,%eax
cmp    $0x3,%ebp
jae    11d7 <find_aks+0x97>
xorps  %xmm0,%xmm0
xor    %edx,%edx
jmp    123b <find_aks+0xfb>
and    $0xfffffffc,%ecx
xorps  %xmm0,%xmm0
xor    %edx,%edx
nop
movss  (%r14,%rdx,4),%xmm1
mulss  (%rsp,%rdx,4),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%rbx)
movss  0x4(%r14,%rdx,4),%xmm0
mulss  0x4(%rsp,%rdx,4),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rbx)
movss  0x8(%r14,%rdx,4),%xmm1
mulss  0x8(%rsp,%rdx,4),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%rbx)
movss  0xc(%r14,%rdx,4),%xmm0
mulss  0xc(%rsp,%rdx,4),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rbx)
add    $0x4,%rdx
cmp    %rdx,%rcx
jne    11e0 <find_aks+0xa0>
test   %rax,%rax
je     126b <find_aks+0x12b>
lea    (%rsp,%rdx,4),%rcx
lea    (%r14,%rdx,4),%rdx
xor    %esi,%esi
nopw   0x0(%rax,%rax,1)
movss  (%rdx,%rsi,4),%xmm1
mulss  (%rcx,%rsi,4),%xmm1
addss  %xmm1,%xmm0
movss  %xmm0,(%rbx)
add    $0x1,%rsi
cmp    %rsi,%rax
jne    1250 <find_aks+0x110>
xorps  %xmm1,%xmm1
ucomiss %xmm0,%xmm1
jbe    1279 <find_aks+0x139>
movl   $0x2b8cbccc,(%rbx)
add    $0x868,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @find_aks(float* %0, float* %1, i32 %2, i32 %3, float* %4) {
%6 = alloca float*
%7 = alloca float*
%8 = alloca i32
%9 = alloca i32
%10 = alloca float*
%11 = alloca [512 x float]
%12 = alloca [21 x float]
%13 = alloca i32
store float* %0, float** %6
store float* %1, float** %7
store i32 %2, i32* %8
store i32 %3, i32* %9
store float* %4, float** %10
%14 = load i32, i32* %9
%15 = icmp slt i32 %14, 20
%16 = zext i1 %15 to i32
%17 = call i32 @assert(i32 %16)
%18 = load i32, i32* %8
%19 = icmp slt i32 %18, 512
%20 = zext i1 %19 to i32
%21 = call i32 @assert(i32 %20)
%22 = load float*, float** %6
%23 = getelementptr inbounds [512 x float], [512 x float]* %11, i64 0, i64 0
%24 = load i32, i32* %8
%25 = call i32 @hanning_window(float* %22, float* %23, i32 %24)
%26 = getelementptr inbounds [512 x float], [512 x float]* %11, i64 0, i64 0
%27 = getelementptr inbounds [21 x float], [21 x float]* %12, i64 0, i64 0
%28 = load i32, i32* %8
%29 = load i32, i32* %9
%30 = call i32 @autocorrelate(float* %26, float* %27, i32 %28, i32 %29)
%31 = getelementptr inbounds [21 x float], [21 x float]* %12, i64 0, i64 0
%32 = load float*, float** %7
%33 = load i32, i32* %9
%34 = call i32 @levinson_durbin(float* %31, float* %32, i32 %33)
%35 = load float*, float** %10
store float 0.000000e+00, float* %35
store i32 0, i32* %13
br label %36
36: ; preds = %53, %5
%37 = load i32, i32* %13
%38 = load i32, i32* %9
%39 = icmp sle i32 %37, %38
br i1 %39, label %40, label %56
40: ; preds = %36
%41 = load float*, float** %7
%42 = load i32, i32* %13
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds float, float* %41, i64 %43
%45 = load float, float* %44
%46 = load i32, i32* %13
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds [21 x float], [21 x float]* %12, i64 0, i64 %47
%49 = load float, float* %48
%50 = load float*, float** %10
%51 = load float, float* %50
%52 = call float @llvm.fmuladd.f32(float %45, float %49, float %51)
store float %52, float* %50
br label %53
53: ; preds = %40
%54 = load i32, i32* %13
%55 = add nsw i32 %54, 1
store i32 %55, i32* %13
br label %36
56: ; preds = %36
%57 = load float*, float** %10
%58 = load float, float* %57
%59 = fpext float %58 to double
%60 = fcmp olt double %59, 0.000000e+00
br i1 %60, label %61, label %63
61: ; preds = %56
%62 = load float*, float** %10
store float 0x3D71979980000000, float* %62
br label %63
63: ; preds = %61, %56
ret void
}

LLVM IR (-O2):
define void @find_aks(float* %0, float* %1, i32 %2, i32 %3, float* nocapture %4) {
%6 = alloca [512 x float]
%7 = alloca [21 x float]
%8 = bitcast [512 x float]* %6 to i8*
%9 = bitcast [21 x float]* %7 to i8*
%10 = icmp slt i32 %3, 20
%11 = zext i1 %10 to i32
%12 = tail call i32 @assert(i32 %11)
%13 = icmp slt i32 %2, 512
%14 = zext i1 %13 to i32
%15 = tail call i32 @assert(i32 %14)
%16 = getelementptr inbounds [512 x float], [512 x float]* %6, i64 0, i64 0
%17 = call i32 @hanning_window(float* %0, float* %16, i32 %2)
%18 = getelementptr inbounds [21 x float], [21 x float]* %7, i64 0, i64 0
%19 = call i32 @autocorrelate(float* %16, float* %18, i32 %2, i32 %3)
%20 = call i32 @levinson_durbin(float* %18, float* %1, i32 %3)
store float 0.000000e+00, float* %4
%21 = icmp slt i32 %3, 0
br i1 %21, label %81, label %22
22: ; preds = %5
%23 = add nuw i32 %3, 1
%24 = zext i32 %23 to i64
%25 = add nsw i64 %24, -1
%26 = and i64 %24, 3
%27 = icmp ult i64 %25, 3
br i1 %27, label %60, label %28
28: ; preds = %22
%29 = and i64 %24, 4294967292
br label %30
30: ; preds = %30, %28
%31 = phi float [ 0.000000e+00, %28 ], [ %56, %30 ]
%32 = phi i64 [ 0, %28 ], [ %57, %30 ]
%33 = phi i64 [ 0, %28 ], [ %58, %30 ]
%34 = getelementptr inbounds float, float* %1, i64 %32
%35 = load float, float* %34
%36 = getelementptr inbounds [21 x float], [21 x float]* %7, i64 0, i64 %32
%37 = load float, float* %36
%38 = call float @llvm.fmuladd.f32(float %35, float %37, float %31)
store float %38, float* %4
%39 = or i64 %32, 1
%40 = getelementptr inbounds float, float* %1, i64 %39
%41 = load float, float* %40
%42 = getelementptr inbounds [21 x float], [21 x float]* %7, i64 0, i64 %39
%43 = load float, float* %42
%44 = call float @llvm.fmuladd.f32(float %41, float %43, float %38)
store float %44, float* %4
%45 = or i64 %32, 2
%46 = getelementptr inbounds float, float* %1, i64 %45
%47 = load float, float* %46
%48 = getelementptr inbounds [21 x float], [21 x float]* %7, i64 0, i64 %45
%49 = load float, float* %48
%50 = call float @llvm.fmuladd.f32(float %47, float %49, float %44)
store float %50, float* %4
%51 = or i64 %32, 3
%52 = getelementptr inbounds float, float* %1, i64 %51
%53 = load float, float* %52
%54 = getelementptr inbounds [21 x float], [21 x float]* %7, i64 0, i64 %51
%55 = load float, float* %54
%56 = call float @llvm.fmuladd.f32(float %53, float %55, float %50)
store float %56, float* %4
%57 = add nuw nsw i64 %32, 4
%58 = add i64 %33, 4
%59 = icmp eq i64 %58, %29
br i1 %59, label %60, label %30
60: ; preds = %30, %22
%61 = phi float [ undef, %22 ], [ %56, %30 ]
%62 = phi float [ 0.000000e+00, %22 ], [ %56, %30 ]
%63 = phi i64 [ 0, %22 ], [ %57, %30 ]
%64 = icmp eq i64 %26, 0
br i1 %64, label %77, label %65
65: ; preds = %60, %65
%66 = phi float [ %73, %65 ], [ %62, %60 ]
%67 = phi i64 [ %74, %65 ], [ %63, %60 ]
%68 = phi i64 [ %75, %65 ], [ 0, %60 ]
%69 = getelementptr inbounds float, float* %1, i64 %67
%70 = load float, float* %69
%71 = getelementptr inbounds [21 x float], [21 x float]* %7, i64 0, i64 %67
%72 = load float, float* %71
%73 = call float @llvm.fmuladd.f32(float %70, float %72, float %66)
store float %73, float* %4
%74 = add nuw nsw i64 %67, 1
%75 = add i64 %68, 1
%76 = icmp eq i64 %75, %26
br i1 %76, label %77, label %65
77: ; preds = %65, %60
%78 = phi float [ %61, %60 ], [ %73, %65 ]
%79 = fcmp olt float %78, 0.000000e+00
br i1 %79, label %80, label %81
80: ; preds = %77
store float 0x3D71979980000000, float* %4
br label %81
81: ; preds = %5, %80, %77
ret void
}