#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[order + 1];
  int i;
  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    $0x840,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %rsp,%rcx
mov    %rcx,-0x828(%rbp)
lea    0xf(,%rax,4),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x840(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x830(%rbp)
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>
mov    -0x840(%rbp),%rsi
lea    -0x820(%rbp),%rdi
mov    -0x14(%rbp),%edx
mov    -0x18(%rbp),%ecx
call   1060 <autocorrelate@plt>
mov    -0x840(%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,-0x834(%rbp)
mov    -0x834(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     1252 <find_aks+0x112>
mov    -0x840(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq -0x834(%rbp),%rdx
movss  (%rcx,%rdx,4),%xmm0
movslq -0x834(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax)
mov    -0x834(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x834(%rbp)
jmp    11f8 <find_aks+0xb8>
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jbe    127b <find_aks+0x13b>
mov    -0x20(%rbp),%rax
movss  0xd89(%rip),%xmm0
movss  %xmm0,(%rax)
mov    -0x828(%rbp),%rax
mov    %rax,%rsp
mov    %rbp,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<find_aks>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x818,%rsp
mov    %r8,%r15
mov    %rcx,-0x38(%rbp)
mov    %edx,%r12d
mov    %rsi,%r14
mov    %rdi,-0x30(%rbp)
mov    %rsp,-0x40(%rbp)
lea    0x1(%rcx),%r13d
mov    %rsp,%rbx
lea    0xf(,%r13,4),%rax
and    $0xfffffffffffffff0,%rax
sub    %rax,%rbx
mov    %rbx,%rsp
xor    %edi,%edi
cmp    $0x200,%edx
setl   %dil
call   1040 <assert@plt>
lea    -0x840(%rbp),%rsi
mov    -0x30(%rbp),%rdi
mov    %r12d,%edx
call   1030 <hanning_window@plt>
lea    -0x840(%rbp),%rdi
mov    %rbx,%rsi
mov    %r12d,%edx
mov    -0x38(%rbp),%r12
mov    %r12d,%ecx
call   1060 <autocorrelate@plt>
mov    %rbx,%rdi
mov    %r14,%rsi
mov    %r12d,%edx
call   1050 <levinson_durbin@plt>
movl   $0x0,(%r15)
test   %r12d,%r12d
js     129b <find_aks+0x15b>
lea    -0x1(%r13),%rcx
mov    %r13d,%eax
and    $0x3,%eax
cmp    $0x3,%rcx
jae    11f4 <find_aks+0xb4>
xorps  %xmm0,%xmm0
xor    %ecx,%ecx
jmp    125f <find_aks+0x11f>
and    $0xfffffffc,%r13d
xorps  %xmm0,%xmm0
xor    %ecx,%ecx
nopl   (%rax)
movss  (%r14,%rcx,4),%xmm1
mulss  (%rbx,%rcx,4),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%r15)
movss  0x4(%r14,%rcx,4),%xmm0
mulss  0x4(%rbx,%rcx,4),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%r15)
movss  0x8(%r14,%rcx,4),%xmm1
mulss  0x8(%rbx,%rcx,4),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%r15)
movss  0xc(%r14,%rcx,4),%xmm0
mulss  0xc(%rbx,%rcx,4),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%r15)
add    $0x4,%rcx
cmp    %rcx,%r13
jne    1200 <find_aks+0xc0>
test   %rax,%rax
je     128c <find_aks+0x14c>
lea    (%rbx,%rcx,4),%rdx
lea    (%r14,%rcx,4),%rcx
xor    %esi,%esi
xchg   %ax,%ax
movss  (%rcx,%rsi,4),%xmm1
mulss  (%rdx,%rsi,4),%xmm1
addss  %xmm1,%xmm0
movss  %xmm0,(%r15)
add    $0x1,%rsi
cmp    %rsi,%rax
jne    1270 <find_aks+0x130>
xorps  %xmm1,%xmm1
ucomiss %xmm0,%xmm1
jbe    129b <find_aks+0x15b>
movl   $0x2b8cbccc,(%r15)
mov    -0x40(%rbp),%rsp
lea    -0x28(%rbp),%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 i8*
%13 = alloca i64
%14 = 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
%15 = load i32, i32* %9
%16 = add nsw i32 %15, 1
%17 = zext i32 %16 to i64
%18 = call i8* @llvm.stacksave()
store i8* %18, i8** %12
%19 = alloca float, i64 %17
store i64 %17, i64* %13
%20 = load i32, i32* %8
%21 = icmp slt i32 %20, 512
%22 = zext i1 %21 to i32
%23 = call i32 @assert(i32 %22)
%24 = load float*, float** %6
%25 = getelementptr inbounds [512 x float], [512 x float]* %11, i64 0, i64 0
%26 = load i32, i32* %8
%27 = call i32 @hanning_window(float* %24, float* %25, i32 %26)
%28 = getelementptr inbounds [512 x float], [512 x float]* %11, i64 0, i64 0
%29 = load i32, i32* %8
%30 = load i32, i32* %9
%31 = call i32 @autocorrelate(float* %28, float* %19, i32 %29, i32 %30)
%32 = load float*, float** %7
%33 = load i32, i32* %9
%34 = call i32 @levinson_durbin(float* %19, float* %32, i32 %33)
%35 = load float*, float** %10
store float 0.000000e+00, float* %35
store i32 0, i32* %14
br label %36
36: ; preds = %53, %5
%37 = load i32, i32* %14
%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* %14
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds float, float* %41, i64 %43
%45 = load float, float* %44
%46 = load i32, i32* %14
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds float, float* %19, 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* %14
%55 = add nsw i32 %54, 1
store i32 %55, i32* %14
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
%64 = load i8*, i8** %12
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 = bitcast [512 x float]* %6 to i8*
%8 = add i32 %3, 1
%9 = zext i32 %8 to i64
%10 = tail call i8* @llvm.stacksave()
%11 = alloca float, i64 %9
%12 = icmp slt i32 %2, 512
%13 = zext i1 %12 to i32
%14 = tail call i32 @assert(i32 %13)
%15 = getelementptr inbounds [512 x float], [512 x float]* %6, i64 0, i64 0
%16 = call i32 @hanning_window(float* %0, float* %15, i32 %2)
%17 = call i32 @autocorrelate(float* %15, float* %11, i32 %2, i32 %3)
%18 = call i32 @levinson_durbin(float* %11, float* %1, i32 %3)
store float 0.000000e+00, float* %4
%19 = icmp slt i32 %3, 0
br i1 %19, label %77, label %20
20: ; preds = %5
%21 = add nsw i64 %9, -1
%22 = and i64 %9, 3
%23 = icmp ult i64 %21, 3
br i1 %23, label %56, label %24
24: ; preds = %20
%25 = and i64 %9, 4294967292
br label %26
26: ; preds = %26, %24
%27 = phi float [ 0.000000e+00, %24 ], [ %52, %26 ]
%28 = phi i64 [ 0, %24 ], [ %53, %26 ]
%29 = phi i64 [ 0, %24 ], [ %54, %26 ]
%30 = getelementptr inbounds float, float* %1, i64 %28
%31 = load float, float* %30
%32 = getelementptr inbounds float, float* %11, i64 %28
%33 = load float, float* %32
%34 = call float @llvm.fmuladd.f32(float %31, float %33, float %27)
store float %34, float* %4
%35 = or i64 %28, 1
%36 = getelementptr inbounds float, float* %1, i64 %35
%37 = load float, float* %36
%38 = getelementptr inbounds float, float* %11, i64 %35
%39 = load float, float* %38
%40 = call float @llvm.fmuladd.f32(float %37, float %39, float %34)
store float %40, float* %4
%41 = or i64 %28, 2
%42 = getelementptr inbounds float, float* %1, i64 %41
%43 = load float, float* %42
%44 = getelementptr inbounds float, float* %11, i64 %41
%45 = load float, float* %44
%46 = call float @llvm.fmuladd.f32(float %43, float %45, float %40)
store float %46, float* %4
%47 = or i64 %28, 3
%48 = getelementptr inbounds float, float* %1, i64 %47
%49 = load float, float* %48
%50 = getelementptr inbounds float, float* %11, i64 %47
%51 = load float, float* %50
%52 = call float @llvm.fmuladd.f32(float %49, float %51, float %46)
store float %52, float* %4
%53 = add nuw nsw i64 %28, 4
%54 = add i64 %29, 4
%55 = icmp eq i64 %54, %25
br i1 %55, label %56, label %26
56: ; preds = %26, %20
%57 = phi float [ undef, %20 ], [ %52, %26 ]
%58 = phi float [ 0.000000e+00, %20 ], [ %52, %26 ]
%59 = phi i64 [ 0, %20 ], [ %53, %26 ]
%60 = icmp eq i64 %22, 0
br i1 %60, label %73, label %61
61: ; preds = %56, %61
%62 = phi float [ %69, %61 ], [ %58, %56 ]
%63 = phi i64 [ %70, %61 ], [ %59, %56 ]
%64 = phi i64 [ %71, %61 ], [ 0, %56 ]
%65 = getelementptr inbounds float, float* %1, i64 %63
%66 = load float, float* %65
%67 = getelementptr inbounds float, float* %11, i64 %63
%68 = load float, float* %67
%69 = call float @llvm.fmuladd.f32(float %66, float %68, float %62)
store float %69, float* %4
%70 = add nuw nsw i64 %63, 1
%71 = add i64 %64, 1
%72 = icmp eq i64 %71, %22
br i1 %72, label %73, label %61
73: ; preds = %61, %56
%74 = phi float [ %57, %56 ], [ %69, %61 ]
%75 = fcmp olt float %74, 0.000000e+00
br i1 %75, label %76, label %77
76: ; preds = %73
store float 0x3D71979980000000, float* %4
br label %77
77: ; preds = %5, %76, %73
ret void
}