#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 corrc_precise(float *templates, float sum_square_template, float *data,
                    int n_samples_template, int normalize) {
  int i;
  float numerator = 0, sum_square_data = 0, denominator = 0, cc = 0, mean = 0,
        sample;
  if (normalize > 0) {
    for (i = 0; i < n_samples_template; i++) {
      mean += data[i];
    }
    mean /= n_samples_template;
  }
  for (i = 0; i < n_samples_template; i++) {
    sample = data[i] - mean;
    numerator += templates[i] * sample;
    sum_square_data += sample * sample;
  }
  denominator = sqrt(sum_square_template * sum_square_data);
  if (denominator > 0.000001f) {
    cc = numerator / denominator;
  }
  return cc;
}


x86 assembly (clang, -O0):
<corrc_precise>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
movss  %xmm0,-0xc(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
mov    %ecx,-0x20(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x28(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x2c(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x34(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x38(%rbp)
cmpl   $0x0,-0x20(%rbp)
jle    11a8 <corrc_precise+0x98>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1195 <corrc_precise+0x85>
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1164 <corrc_precise+0x54>
cvtsi2ssl -0x1c(%rbp),%xmm1
movss  -0x38(%rbp),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,-0x38(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1220 <corrc_precise+0x110>
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  -0x3c(%rbp),%xmm2
movss  -0x28(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x3c(%rbp),%xmm0
movss  -0x3c(%rbp),%xmm2
movss  -0x2c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11af <corrc_precise+0x9f>
movss  -0xc(%rbp),%xmm0
mulss  -0x2c(%rbp),%xmm0
call   1030 <sqrt@plt>
movss  %xmm0,-0x30(%rbp)
movss  -0x30(%rbp),%xmm0
movss  0xdbf(%rip),%xmm1
ucomiss %xmm1,%xmm0
jbe    1259 <corrc_precise+0x149>
movss  -0x28(%rbp),%xmm0
divss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x34(%rbp),%xmm0
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<corrc_precise>:
sub    $0x18,%rsp
mov    %edx,%r8d
test   %ecx,%ecx
jle    1137 <corrc_precise+0x27>
test   %edx,%edx
jle    114a <corrc_precise+0x3a>
lea    -0x1(%r8),%rax
mov    %r8d,%r9d
and    $0x7,%r9d
cmp    $0x7,%rax
jae    114f <corrc_precise+0x3f>
xorps  %xmm1,%xmm1
xor    %eax,%eax
jmp    1198 <corrc_precise+0x88>
xorps  %xmm1,%xmm1
test   %edx,%edx
jg     11ce <corrc_precise+0xbe>
xorps  %xmm6,%xmm6
jmp    1258 <corrc_precise+0x148>
xorps  %xmm1,%xmm1
jmp    11be <corrc_precise+0xae>
mov    %r8d,%ecx
and    $0xfffffff8,%ecx
xorps  %xmm1,%xmm1
xor    %eax,%eax
nopw   0x0(%rax,%rax,1)
addss  (%rsi,%rax,4),%xmm1
addss  0x4(%rsi,%rax,4),%xmm1
addss  0x8(%rsi,%rax,4),%xmm1
addss  0xc(%rsi,%rax,4),%xmm1
addss  0x10(%rsi,%rax,4),%xmm1
addss  0x14(%rsi,%rax,4),%xmm1
addss  0x18(%rsi,%rax,4),%xmm1
addss  0x1c(%rsi,%rax,4),%xmm1
add    $0x8,%rax
cmp    %rax,%rcx
jne    1160 <corrc_precise+0x50>
test   %r9,%r9
je     11be <corrc_precise+0xae>
lea    (%rsi,%rax,4),%rax
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
addss  (%rax,%rcx,4),%xmm1
add    $0x1,%rcx
cmp    %rcx,%r9
jne    11b0 <corrc_precise+0xa0>
cvtsi2ss %edx,%xmm2
divss  %xmm2,%xmm1
test   %edx,%edx
jle    1142 <corrc_precise+0x32>
cmp    $0x1,%edx
jne    11e0 <corrc_precise+0xd0>
xorps  %xmm2,%xmm2
xor    %ecx,%ecx
test   $0x1,%r8b
jne    123d <corrc_precise+0x12d>
jmp    1258 <corrc_precise+0x148>
mov    %r8d,%eax
and    $0xfffffffe,%eax
xorps  %xmm2,%xmm2
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
movss  (%rsi,%rcx,4),%xmm3
movss  0x4(%rsi,%rcx,4),%xmm6
subss  %xmm1,%xmm3
movss  (%rdi,%rcx,4),%xmm4
movss  0x4(%rdi,%rcx,4),%xmm5
unpcklps %xmm3,%xmm4
shufps $0x0,%xmm3,%xmm3
mulps  %xmm4,%xmm3
addps  %xmm2,%xmm3
subss  %xmm1,%xmm6
unpcklps %xmm6,%xmm5
shufps $0x0,%xmm6,%xmm6
mulps  %xmm5,%xmm6
movaps %xmm6,%xmm2
addps  %xmm3,%xmm2
add    $0x2,%rcx
cmp    %rcx,%rax
jne    11f0 <corrc_precise+0xe0>
addps  %xmm3,%xmm6
test   $0x1,%r8b
je     1258 <corrc_precise+0x148>
movss  (%rsi,%rcx,4),%xmm6
subss  %xmm1,%xmm6
movss  (%rdi,%rcx,4),%xmm1
unpcklps %xmm6,%xmm1
shufps $0x0,%xmm6,%xmm6
mulps  %xmm1,%xmm6
addps  %xmm2,%xmm6
movaps %xmm6,%xmm1
shufps $0x55,%xmm6,%xmm1
mulss  %xmm0,%xmm1
xorps  %xmm2,%xmm2
ucomiss %xmm2,%xmm1
jb     127d <corrc_precise+0x16d>
xorps  %xmm0,%xmm0
sqrtss %xmm1,%xmm0
ucomiss 0xd87(%rip),%xmm0
ja     1299 <corrc_precise+0x189>
jmp    12a0 <corrc_precise+0x190>
movaps %xmm1,%xmm0
movaps %xmm6,(%rsp)
call   1030 <sqrt@plt>
xorps  %xmm2,%xmm2
movaps (%rsp),%xmm6
ucomiss 0xd69(%rip),%xmm0
jbe    12a0 <corrc_precise+0x190>
divss  %xmm0,%xmm6
movaps %xmm6,%xmm2
movaps %xmm2,%xmm0
add    $0x18,%rsp
ret

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

LLVM IR (-O2):
define float @corrc_precise(float* nocapture %0, float %1, float* nocapture %2, i32 %3, i32 %4) {
%6 = icmp sgt i32 %4, 0
br i1 %6, label %7, label %73
7: ; preds = %5
%8 = icmp sgt i32 %3, 0
br i1 %8, label %9, label %69
9: ; preds = %7
%10 = zext i32 %3 to i64
%11 = add nsw i64 %10, -1
%12 = and i64 %10, 7
%13 = icmp ult i64 %11, 7
br i1 %13, label %54, label %14
14: ; preds = %9
%15 = and i64 %10, 4294967288
br label %16
16: ; preds = %16, %14
%17 = phi i64 [ 0, %14 ], [ %51, %16 ]
%18 = phi float [ 0.000000e+00, %14 ], [ %50, %16 ]
%19 = phi i64 [ 0, %14 ], [ %52, %16 ]
%20 = getelementptr inbounds float, float* %2, i64 %17
%21 = load float, float* %20
%22 = fadd float %18, %21
%23 = or i64 %17, 1
%24 = getelementptr inbounds float, float* %2, i64 %23
%25 = load float, float* %24
%26 = fadd float %22, %25
%27 = or i64 %17, 2
%28 = getelementptr inbounds float, float* %2, i64 %27
%29 = load float, float* %28
%30 = fadd float %26, %29
%31 = or i64 %17, 3
%32 = getelementptr inbounds float, float* %2, i64 %31
%33 = load float, float* %32
%34 = fadd float %30, %33
%35 = or i64 %17, 4
%36 = getelementptr inbounds float, float* %2, i64 %35
%37 = load float, float* %36
%38 = fadd float %34, %37
%39 = or i64 %17, 5
%40 = getelementptr inbounds float, float* %2, i64 %39
%41 = load float, float* %40
%42 = fadd float %38, %41
%43 = or i64 %17, 6
%44 = getelementptr inbounds float, float* %2, i64 %43
%45 = load float, float* %44
%46 = fadd float %42, %45
%47 = or i64 %17, 7
%48 = getelementptr inbounds float, float* %2, i64 %47
%49 = load float, float* %48
%50 = fadd float %46, %49
%51 = add nuw nsw i64 %17, 8
%52 = add i64 %19, 8
%53 = icmp eq i64 %52, %15
br i1 %53, label %54, label %16
54: ; preds = %16, %9
%55 = phi float [ undef, %9 ], [ %50, %16 ]
%56 = phi i64 [ 0, %9 ], [ %51, %16 ]
%57 = phi float [ 0.000000e+00, %9 ], [ %50, %16 ]
%58 = icmp eq i64 %12, 0
br i1 %58, label %69, label %59
59: ; preds = %54, %59
%60 = phi i64 [ %66, %59 ], [ %56, %54 ]
%61 = phi float [ %65, %59 ], [ %57, %54 ]
%62 = phi i64 [ %67, %59 ], [ 0, %54 ]
%63 = getelementptr inbounds float, float* %2, i64 %60
%64 = load float, float* %63
%65 = fadd float %61, %64
%66 = add nuw nsw i64 %60, 1
%67 = add i64 %62, 1
%68 = icmp eq i64 %67, %12
br i1 %68, label %69, label %59
69: ; preds = %54, %59, %7
%70 = phi float [ 0.000000e+00, %7 ], [ %55, %54 ], [ %65, %59 ]
%71 = sitofp i32 %3 to float
%72 = fdiv float %70, %71
br label %73
73: ; preds = %69, %5
%74 = phi float [ %72, %69 ], [ 0.000000e+00, %5 ]
%75 = icmp sgt i32 %3, 0
br i1 %75, label %76, label %128
76: ; preds = %73
%77 = zext i32 %3 to i64
%78 = and i64 %77, 1
%79 = icmp eq i32 %3, 1
br i1 %79, label %112, label %80
80: ; preds = %76
%81 = and i64 %77, 4294967294
br label %82
82: ; preds = %82, %80
%83 = phi i64 [ 0, %80 ], [ %107, %82 ]
%84 = phi <2 x float> [ zeroinitializer, %80 ], [ %106, %82 ]
%85 = phi i64 [ 0, %80 ], [ %108, %82 ]
%86 = getelementptr inbounds float, float* %2, i64 %83
%87 = load float, float* %86
%88 = fsub float %87, %74
%89 = getelementptr inbounds float, float* %0, i64 %83
%90 = load float, float* %89
%91 = insertelement <2 x float> poison, float %90, i64 0
%92 = insertelement <2 x float> %91, float %88, i64 1
%93 = insertelement <2 x float> poison, float %88, i64 0
%94 = shufflevector <2 x float> %93, <2 x float> poison, <2 x i32> zeroinitializer
%95 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %92, <2 x float> %94, <2 x float> %84)
%96 = or i64 %83, 1
%97 = getelementptr inbounds float, float* %2, i64 %96
%98 = load float, float* %97
%99 = fsub float %98, %74
%100 = getelementptr inbounds float, float* %0, i64 %96
%101 = load float, float* %100
%102 = insertelement <2 x float> poison, float %101, i64 0
%103 = insertelement <2 x float> %102, float %99, i64 1
%104 = insertelement <2 x float> poison, float %99, i64 0
%105 = shufflevector <2 x float> %104, <2 x float> poison, <2 x i32> zeroinitializer
%106 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %103, <2 x float> %105, <2 x float> %95)
%107 = add nuw nsw i64 %83, 2
%108 = add i64 %85, 2
%109 = icmp eq i64 %108, %81
br i1 %109, label %110, label %82
110: ; preds = %82
%111 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %103, <2 x float> %105, <2 x float> %95)
br label %112
112: ; preds = %110, %76
%113 = phi <2 x float> [ undef, %76 ], [ %111, %110 ]
%114 = phi i64 [ 0, %76 ], [ %107, %110 ]
%115 = phi <2 x float> [ zeroinitializer, %76 ], [ %106, %110 ]
%116 = icmp eq i64 %78, 0
br i1 %116, label %128, label %117
117: ; preds = %112
%118 = getelementptr inbounds float, float* %2, i64 %114
%119 = load float, float* %118
%120 = fsub float %119, %74
%121 = getelementptr inbounds float, float* %0, i64 %114
%122 = load float, float* %121
%123 = insertelement <2 x float> poison, float %122, i64 0
%124 = insertelement <2 x float> %123, float %120, i64 1
%125 = insertelement <2 x float> poison, float %120, i64 0
%126 = shufflevector <2 x float> %125, <2 x float> poison, <2 x i32> zeroinitializer
%127 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %124, <2 x float> %126, <2 x float> %115)
br label %128
128: ; preds = %117, %112, %73
%129 = phi <2 x float> [ zeroinitializer, %73 ], [ %113, %112 ], [ %127, %117 ]
%130 = extractelement <2 x float> %129, i64 1
%131 = fmul float %130, %1
%132 = tail call float @sqrt(float %131)
%133 = fcmp ogt float %132, 0x3EB0C6F7A0000000
%134 = extractelement <2 x float> %129, i64 0
%135 = fdiv float %134, %132
%136 = select i1 %133, float %135, float 0.000000e+00
ret float %136
}