#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 */
typedef scalar_t__ uint8_t;

/* Variables and functions */
float fabsf(float);
int sqrtf(float);
float pearsonr(uint8_t *x, uint8_t *y, int len) {
  float mx = 0.0f, my = 0.0f;
  float xmi, ymi;
  float dx = 0.0f, dy = 0.0f;
  float r_num = 0.0f, r_den;
  int i;
  for (i = 0; i < len; i++) {
    mx += x[i];
    my += y[i];
  }
  mx /= len;
  my /= len;
  for (i = 0; i < len; i++) {
    xmi = (float)x[i] - mx;
    ymi = (float)y[i] - my;
    r_num += xmi * ymi;
    dx += xmi * xmi;
    dy += ymi * ymi;
  }
  r_num *= len;
  r_den = sqrtf(dx * dy) * len;
  return fabsf(r_num / r_den);
}


x86 assembly (clang, -O0):
<pearsonr>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x18(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x1c(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x28(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x2c(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x30(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    119c <pearsonr+0x8c>
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
cvtsi2ssq (%rax,%rcx,8),%xmm0
addss  -0x18(%rbp),%xmm0
movss  %xmm0,-0x18(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
cvtsi2ssq (%rax,%rcx,8),%xmm0
addss  -0x1c(%rbp),%xmm0
movss  %xmm0,-0x1c(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1152 <pearsonr+0x42>
cvtsi2ssl -0x14(%rbp),%xmm1
movss  -0x18(%rbp),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,-0x18(%rbp)
cvtsi2ssl -0x14(%rbp),%xmm1
movss  -0x1c(%rbp),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,-0x1c(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1267 <pearsonr+0x157>
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
cvtsi2ssq (%rax,%rcx,8),%xmm0
subss  -0x18(%rbp),%xmm0
movss  %xmm0,-0x20(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
cvtsi2ssq (%rax,%rcx,8),%xmm0
subss  -0x1c(%rbp),%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x20(%rbp),%xmm0
movss  -0x24(%rbp),%xmm2
movss  -0x30(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x20(%rbp),%xmm0
movss  -0x20(%rbp),%xmm2
movss  -0x28(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x24(%rbp),%xmm0
movss  -0x24(%rbp),%xmm2
movss  -0x2c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    11c9 <pearsonr+0xb9>
cvtsi2ssl -0x14(%rbp),%xmm0
mulss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x28(%rbp),%xmm0
mulss  -0x2c(%rbp),%xmm0
call   1030 <sqrtf@plt>
imul   -0x14(%rbp),%eax
cvtsi2ss %eax,%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x30(%rbp),%xmm0
divss  -0x34(%rbp),%xmm0
movaps 0xd5d(%rip),%xmm1
pand   %xmm1,%xmm0
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<pearsonr>:
push   %rbx
sub    $0x10,%rsp
mov    %edx,%ebx
test   %edx,%edx
jle    1129 <pearsonr+0x19>
mov    %ebx,%eax
cmp    $0x1,%ebx
jne    1132 <pearsonr+0x22>
xorps  %xmm0,%xmm0
xor    %ecx,%ecx
jmp    117b <pearsonr+0x6b>
cvtsi2ss %ebx,%xmm5
jmp    1203 <pearsonr+0xf3>
mov    %eax,%edx
and    $0xfffffffe,%edx
xorps  %xmm0,%xmm0
xor    %ecx,%ecx
nopl   0x0(%rax)
xorps  %xmm1,%xmm1
cvtsi2ssq (%rdi,%rcx,8),%xmm1
xorps  %xmm2,%xmm2
cvtsi2ssq (%rsi,%rcx,8),%xmm2
unpcklps %xmm1,%xmm2
xorps  %xmm1,%xmm1
cvtsi2ssq 0x8(%rdi,%rcx,8),%xmm1
addps  %xmm0,%xmm2
xorps  %xmm0,%xmm0
cvtsi2ssq 0x8(%rsi,%rcx,8),%xmm0
unpcklps %xmm1,%xmm0
addps  %xmm2,%xmm0
add    $0x2,%rcx
cmp    %rcx,%rdx
jne    1140 <pearsonr+0x30>
test   $0x1,%al
je     1197 <pearsonr+0x87>
xorps  %xmm1,%xmm1
cvtsi2ssq (%rdi,%rcx,8),%xmm1
xorps  %xmm2,%xmm2
cvtsi2ssq (%rsi,%rcx,8),%xmm2
unpcklps %xmm1,%xmm2
addps  %xmm2,%xmm0
cvtsi2ss %ebx,%xmm5
test   %ebx,%ebx
jle    1203 <pearsonr+0xf3>
movaps %xmm5,%xmm1
shufps $0x0,%xmm5,%xmm1
divps  %xmm1,%xmm0
xorps  %xmm2,%xmm2
xorps  %xmm1,%xmm1
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
xorps  %xmm3,%xmm3
cvtsi2ssq (%rdi,%rcx,8),%xmm3
xorps  %xmm4,%xmm4
cvtsi2ssq (%rsi,%rcx,8),%xmm4
unpcklps %xmm3,%xmm4
subps  %xmm0,%xmm4
movaps %xmm4,%xmm3
shufps $0x55,%xmm4,%xmm3
mulss  %xmm4,%xmm3
addss  %xmm3,%xmm1
mulps  %xmm4,%xmm4
addps  %xmm4,%xmm2
add    $0x1,%rcx
cmp    %rcx,%rax
jne    11c0 <pearsonr+0xb0>
movaps %xmm2,%xmm0
shufps $0x55,%xmm2,%xmm0
mulss  %xmm2,%xmm0
jmp    1209 <pearsonr+0xf9>
xorps  %xmm1,%xmm1
xorps  %xmm0,%xmm0
mulss  %xmm1,%xmm5
movaps %xmm5,(%rsp)
call   1030 <sqrtf@plt>
imul   %ebx,%eax
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
movaps (%rsp),%xmm1
divss  %xmm0,%xmm1
andps  0xdd1(%rip),%xmm1
movaps %xmm1,%xmm0
add    $0x10,%rsp
pop    %rbx
ret

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

LLVM IR (-O2):
define float @pearsonr(i64* nocapture %0, i64* nocapture %1, i32 %2) {
%4 = icmp sgt i32 %2, 0
br i1 %4, label %7, label %5
5: ; preds = %3
%6 = sitofp i32 %2 to float
br label %81
7: ; preds = %3
%8 = zext i32 %2 to i64
%9 = and i64 %8, 1
%10 = icmp eq i32 %2, 1
br i1 %10, label %37, label %11
11: ; preds = %7
%12 = and i64 %8, 4294967294
br label %13
13: ; preds = %13, %11
%14 = phi i64 [ 0, %11 ], [ %34, %13 ]
%15 = phi <2 x float> [ zeroinitializer, %11 ], [ %33, %13 ]
%16 = phi i64 [ 0, %11 ], [ %35, %13 ]
%17 = getelementptr inbounds i64, i64* %0, i64 %14
%18 = load i64, i64* %17
%19 = getelementptr inbounds i64, i64* %1, i64 %14
%20 = load i64, i64* %19
%21 = insertelement <2 x i64> poison, i64 %20, i64 0
%22 = insertelement <2 x i64> %21, i64 %18, i64 1
%23 = sitofp <2 x i64> %22 to <2 x float>
%24 = fadd <2 x float> %15, %23
%25 = or i64 %14, 1
%26 = getelementptr inbounds i64, i64* %0, i64 %25
%27 = load i64, i64* %26
%28 = getelementptr inbounds i64, i64* %1, i64 %25
%29 = load i64, i64* %28
%30 = insertelement <2 x i64> poison, i64 %29, i64 0
%31 = insertelement <2 x i64> %30, i64 %27, i64 1
%32 = sitofp <2 x i64> %31 to <2 x float>
%33 = fadd <2 x float> %24, %32
%34 = add nuw nsw i64 %14, 2
%35 = add i64 %16, 2
%36 = icmp eq i64 %35, %12
br i1 %36, label %37, label %13
37: ; preds = %13, %7
%38 = phi <2 x float> [ undef, %7 ], [ %33, %13 ]
%39 = phi i64 [ 0, %7 ], [ %34, %13 ]
%40 = phi <2 x float> [ zeroinitializer, %7 ], [ %33, %13 ]
%41 = icmp eq i64 %9, 0
br i1 %41, label %51, label %42
42: ; preds = %37
%43 = getelementptr inbounds i64, i64* %0, i64 %39
%44 = load i64, i64* %43
%45 = getelementptr inbounds i64, i64* %1, i64 %39
%46 = load i64, i64* %45
%47 = insertelement <2 x i64> poison, i64 %46, i64 0
%48 = insertelement <2 x i64> %47, i64 %44, i64 1
%49 = sitofp <2 x i64> %48 to <2 x float>
%50 = fadd <2 x float> %40, %49
br label %51
51: ; preds = %37, %42
%52 = phi <2 x float> [ %38, %37 ], [ %50, %42 ]
%53 = sitofp i32 %2 to float
%54 = insertelement <2 x float> poison, float %53, i64 0
%55 = shufflevector <2 x float> %54, <2 x float> poison, <2 x i32> zeroinitializer
%56 = fdiv <2 x float> %52, %55
br i1 %4, label %57, label %81
57: ; preds = %51
%58 = zext i32 %2 to i64
br label %59
59: ; preds = %57, %59
%60 = phi i64 [ 0, %57 ], [ %75, %59 ]
%61 = phi float [ 0.000000e+00, %57 ], [ %73, %59 ]
%62 = phi <2 x float> [ zeroinitializer, %57 ], [ %74, %59 ]
%63 = getelementptr inbounds i64, i64* %0, i64 %60
%64 = load i64, i64* %63
%65 = getelementptr inbounds i64, i64* %1, i64 %60
%66 = load i64, i64* %65
%67 = insertelement <2 x i64> poison, i64 %66, i64 0
%68 = insertelement <2 x i64> %67, i64 %64, i64 1
%69 = sitofp <2 x i64> %68 to <2 x float>
%70 = fsub <2 x float> %69, %56
%71 = extractelement <2 x float> %70, i64 0
%72 = extractelement <2 x float> %70, i64 1
%73 = tail call float @llvm.fmuladd.f32(float %72, float %71, float %61)
%74 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %70, <2 x float> %70, <2 x float> %62)
%75 = add nuw nsw i64 %60, 1
%76 = icmp eq i64 %75, %58
br i1 %76, label %77, label %59
77: ; preds = %59
%78 = shufflevector <2 x float> %74, <2 x float> poison, <2 x i32> <i32 1, i32 undef>
%79 = fmul <2 x float> %78, %74
%80 = extractelement <2 x float> %79, i64 0
br label %81
81: ; preds = %5, %77, %51
%82 = phi float [ %53, %51 ], [ %53, %77 ], [ %6, %5 ]
%83 = phi float [ 0.000000e+00, %51 ], [ %73, %77 ], [ 0.000000e+00, %5 ]
%84 = phi float [ 0.000000e+00, %51 ], [ %80, %77 ], [ 0.000000e+00, %5 ]
%85 = fmul float %83, %82
%86 = tail call i32 @sqrtf(float %84)
%87 = mul nsw i32 %86, %2
%88 = sitofp i32 %87 to float
%89 = fdiv float %85, %88
%90 = tail call float @llvm.fabs.f32(float %89)
ret float %90
}