#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 */
float sqrt(float);
float v_inner(float *, float *, int);
float v_magsq(float *, int);
float find_pitch(float sig_in[], float *pcorr, int lower, int upper,
                 int length) {
  int i, cbegin, ipitch, even_flag;
  float corr, maxcorr;
  float c0_0, cT_T;
  ipitch = lower;
  maxcorr = 0.0;
  even_flag = 1;
  cbegin = -((length + upper) / 2);
  c0_0 = v_magsq(&sig_in[cbegin], length);
  cT_T = v_magsq(&sig_in[cbegin + upper], length);
  for (i = upper; i >= lower; i--) {
    corr = v_inner(&sig_in[cbegin], &sig_in[cbegin + i], length);
    if (corr > 0.01)
      corr = (corr * corr) / (c0_0 * cT_T);
    if (corr > maxcorr) {
      maxcorr = corr;
      ipitch = i;
    }
    if (even_flag) {
      even_flag = 0;
      c0_0 += (sig_in[cbegin + length] * sig_in[cbegin + length]);
      c0_0 -= (sig_in[cbegin] * sig_in[cbegin]);
      cbegin++;
    } else {
      even_flag = 1;
      cT_T += (sig_in[cbegin + i - 1] * sig_in[cbegin + i - 1]);
      cT_T -=
          (sig_in[cbegin + i - 1 + length] * sig_in[cbegin + i - 1 + length]);
    }
  }
  *pcorr = sqrt(maxcorr);
  return ((float)ipitch);
}


x86 assembly (clang, -O0):
<find_pitch>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x28(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x34(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x1c(%rbp),%eax
add    -0x18(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
xor    %eax,%eax
sub    %ecx,%eax
mov    %eax,-0x24(%rbp)
mov    -0x8(%rbp),%rdi
movslq -0x24(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
mov    -0x1c(%rbp),%esi
call   1050 <v_magsq@plt>
movss  %xmm0,-0x38(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0x24(%rbp),%eax
add    -0x18(%rbp),%eax
cltq
shl    $0x2,%rax
add    %rax,%rdi
mov    -0x1c(%rbp),%esi
call   1050 <v_magsq@plt>
movss  %xmm0,-0x3c(%rbp)
mov    -0x18(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x14(%rbp),%eax
jl     137d <find_pitch+0x24d>
mov    -0x8(%rbp),%rdi
movslq -0x24(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
mov    -0x8(%rbp),%rsi
mov    -0x24(%rbp),%eax
add    -0x20(%rbp),%eax
cltq
shl    $0x2,%rax
add    %rax,%rsi
mov    -0x1c(%rbp),%edx
call   1030 <v_inner@plt>
movss  %xmm0,-0x30(%rbp)
movss  -0x30(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xdfc(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    122b <find_pitch+0xfb>
movss  -0x30(%rbp),%xmm0
mulss  -0x30(%rbp),%xmm0
movss  -0x38(%rbp),%xmm1
mulss  -0x3c(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x30(%rbp),%xmm0
ucomiss -0x34(%rbp),%xmm0
jbe    124a <find_pitch+0x11a>
movss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x20(%rbp),%eax
mov    %eax,-0x28(%rbp)
cmpl   $0x0,-0x2c(%rbp)
je     12d8 <find_pitch+0x1a8>
movl   $0x0,-0x2c(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x38(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x38(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x38(%rbp),%xmm1
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    136a <find_pitch+0x23a>
movl   $0x1,-0x2c(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x20(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x20(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x3c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x3c(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x20(%rbp),%ecx
sub    $0x1,%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x20(%rbp),%ecx
sub    $0x1,%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x3c(%rbp),%xmm1
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x3c(%rbp)
jmp    136f <find_pitch+0x23f>
mov    -0x20(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x20(%rbp)
jmp    11b8 <find_pitch+0x88>
movss  -0x34(%rbp),%xmm0
call   1040 <sqrt@plt>
mov    -0x10(%rbp),%rax
movss  %xmm0,(%rax)
cvtsi2ssl -0x28(%rbp),%xmm0
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<find_pitch>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %r8d,%r15d
mov    %ecx,%ebp
mov    %edx,%r12d
mov    %rsi,0x30(%rsp)
mov    %rdi,%r13
lea    (%r15,%rbp,1),%eax
mov    %eax,%ebx
shr    $0x1f,%ebx
add    %eax,%ebx
sar    %ebx
mov    %ebx,%r14d
neg    %r14d
movslq %r14d,%rax
lea    (%rdi,%rax,4),%rdi
mov    %r8d,%esi
call   1050 <v_magsq@plt>
movss  %xmm0,0x10(%rsp)
mov    %ebp,%eax
sub    %ebx,%eax
cltq
mov    %r13,0x18(%rsp)
lea    0x0(,%rax,4),%rdi
add    %r13,%rdi
mov    %r15,0x20(%rsp)
mov    %r15d,%esi
call   1050 <v_magsq@plt>
mov    %r12d,%ecx
movss  %xmm0,0xc(%rsp)
cmp    %r12d,%ebp
jge    11c3 <find_pitch+0x93>
mov    %ecx,%r15d
xorps  %xmm0,%xmm0
ucomiss 0xe54(%rip),%xmm0
jb     12db <find_pitch+0x1ab>
sqrtss %xmm0,%xmm0
jmp    12e0 <find_pitch+0x1b0>
movslq %ebp,%r12
add    $0xffffffff,%ebp
mov    $0x1,%r13d
xorps  %xmm1,%xmm1
mov    %ecx,%r15d
mov    %ecx,0x14(%rsp)
jmp    123e <find_pitch+0x10e>
nopl   0x0(%rax,%rax,1)
lea    (%r14,%rbp,1),%eax
cltq
mov    0x18(%rsp),%rcx
movss  (%rcx,%rax,4),%xmm1
mulss  %xmm1,%xmm1
addss  0xc(%rsp),%xmm1
mov    0x20(%rsp),%rax
add    %r14d,%eax
add    %ebp,%eax
cltq
movss  (%rcx,%rax,4),%xmm2
mulss  %xmm2,%xmm2
subss  %xmm2,%xmm1
mov    $0x1,%r13d
movss  %xmm1,0xc(%rsp)
mov    0x14(%rsp),%ecx
maxss  %xmm3,%xmm0
lea    -0x1(%rbp),%eax
add    $0x1,%ebp
add    $0xffffffffffffffff,%r12
movaps %xmm0,%xmm1
cmp    %ecx,%ebp
mov    %eax,%ebp
jle    11ad <find_pitch+0x7d>
movss  %xmm1,0x2c(%rsp)
movslq %r14d,%rax
mov    0x18(%rsp),%rcx
lea    (%rcx,%rax,4),%rbx
add    %r12,%rax
lea    (%rcx,%rax,4),%rsi
mov    %rbx,%rdi
mov    0x20(%rsp),%rdx
call   1030 <v_inner@plt>
xorps  %xmm1,%xmm1
cvtss2sd %xmm0,%xmm1
ucomisd 0xd8d(%rip),%xmm1
jbe    1289 <find_pitch+0x159>
mulss  %xmm0,%xmm0
movss  0x10(%rsp),%xmm1
mulss  0xc(%rsp),%xmm1
divss  %xmm1,%xmm0
movss  0x2c(%rsp),%xmm3
ucomiss %xmm3,%xmm0
cmova  %r12d,%r15d
test   %r13d,%r13d
je     11e0 <find_pitch+0xb0>
mov    0x20(%rsp),%rax
add    %r14d,%eax
cltq
mov    0x18(%rsp),%rcx
movss  (%rcx,%rax,4),%xmm1
mulss  %xmm1,%xmm1
addss  0x10(%rsp),%xmm1
movss  (%rbx),%xmm2
mulss  %xmm2,%xmm2
subss  %xmm2,%xmm1
add    $0x1,%r14d
xor    %r13d,%r13d
movss  %xmm1,0x10(%rsp)
jmp    121f <find_pitch+0xef>
call   1040 <sqrt@plt>
mov    0x30(%rsp),%rax
movss  %xmm0,(%rax)
xorps  %xmm0,%xmm0
cvtsi2ss %r15d,%xmm0
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define float @find_pitch(float* %0, float* nocapture writeonly %1, i32 %2, i32 %3, i32 %4) {
%6 = add nsw i32 %4, %3
%7 = sdiv i32 %6, -2
%8 = sext i32 %7 to i64
%9 = getelementptr inbounds float, float* %0, i64 %8
%10 = tail call float @v_magsq(float* %9, i32 %4)
%11 = add nsw i32 %7, %3
%12 = sext i32 %11 to i64
%13 = getelementptr inbounds float, float* %0, i64 %12
%14 = tail call float @v_magsq(float* %13, i32 %4)
%15 = icmp slt i32 %3, %2
br i1 %15, label %71, label %16
16: ; preds = %5, %64
%17 = phi float [ %68, %64 ], [ %14, %5 ]
%18 = phi float [ %67, %64 ], [ %10, %5 ]
%19 = phi float [ %40, %64 ], [ 0.000000e+00, %5 ]
%20 = phi i32 [ %66, %64 ], [ 1, %5 ]
%21 = phi i32 [ %39, %64 ], [ %2, %5 ]
%22 = phi i32 [ %65, %64 ], [ %7, %5 ]
%23 = phi i32 [ %69, %64 ], [ %3, %5 ]
%24 = sext i32 %22 to i64
%25 = getelementptr inbounds float, float* %0, i64 %24
%26 = add nsw i32 %22, %23
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds float, float* %0, i64 %27
%29 = tail call float @v_inner(float* %25, float* %28, i32 %4)
%30 = fpext float %29 to double
%31 = fcmp ogt double %30, 1.000000e-02
br i1 %31, label %32, label %36
32: ; preds = %16
%33 = fmul float %29, %29
%34 = fmul float %18, %17
%35 = fdiv float %33, %34
br label %36
36: ; preds = %32, %16
%37 = phi float [ %35, %32 ], [ %29, %16 ]
%38 = fcmp ogt float %37, %19
%39 = select i1 %38, i32 %23, i32 %21
%40 = select i1 %38, float %37, float %19
%41 = icmp eq i32 %20, 0
br i1 %41, label %52, label %42
42: ; preds = %36
%43 = add nsw i32 %22, %4
%44 = sext i32 %43 to i64
%45 = getelementptr inbounds float, float* %0, i64 %44
%46 = load float, float* %45
%47 = tail call float @llvm.fmuladd.f32(float %46, float %46, float %18)
%48 = load float, float* %25
%49 = fneg float %48
%50 = tail call float @llvm.fmuladd.f32(float %49, float %48, float %47)
%51 = add nsw i32 %22, 1
br label %64
52: ; preds = %36
%53 = add nsw i32 %26, -1
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds float, float* %0, i64 %54
%56 = load float, float* %55
%57 = tail call float @llvm.fmuladd.f32(float %56, float %56, float %17)
%58 = add nsw i32 %53, %4
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds float, float* %0, i64 %59
%61 = load float, float* %60
%62 = fneg float %61
%63 = tail call float @llvm.fmuladd.f32(float %62, float %61, float %57)
br label %64
64: ; preds = %42, %52
%65 = phi i32 [ %51, %42 ], [ %22, %52 ]
%66 = phi i32 [ 0, %42 ], [ 1, %52 ]
%67 = phi float [ %50, %42 ], [ %18, %52 ]
%68 = phi float [ %17, %42 ], [ %63, %52 ]
%69 = add nsw i32 %23, -1
%70 = icmp sgt i32 %23, %2
br i1 %70, label %16, label %71
71: ; preds = %64, %5
%72 = phi i32 [ %2, %5 ], [ %39, %64 ]
%73 = phi float [ 0.000000e+00, %5 ], [ %40, %64 ]
%74 = tail call float @sqrt(float %73)
store float %74, float* %1
%75 = sitofp i32 %72 to float
ret float %75
}