#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 */
double sqrt(double);
int correl(double x[], double y[], int n, double *mx, double *my, double *v11,
           double *v12, double *v22, double *r) {
  int i;
  *mx = *my = *v11 = *v12 = *v22 = 0.0;
  for (i = 0; i < n; i++) {
    *v11 += ((x[i] - *mx) * (x[i] - *mx)) * i / (i + 1.);
    *v22 += ((y[i] - *my) * (y[i] - *my)) * i / (i + 1.);
    *v12 += (x[i] - *mx) * (y[i] - *my) * i / (i + 1.);
    *mx = (*mx * i + x[i]) / (i + 1.);
    *my = (*my * i + y[i]) / (i + 1.);
  }
  *v11 /= n;
  *v22 /= n;
  *v12 /= n;
  if (*v11 > 0.0 && *v22 > 0.0)
    *r = *v12 / sqrt(*v11 * *v22);
  else
    *r = -9;
  return (0);
}


x86 assembly (clang, -O0):
<correl>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    0x18(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
mov    0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
mov    -0x30(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1319 <correl+0x209>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
subsd  (%rax),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
subsd  (%rax),%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x34(%rbp),%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x34(%rbp),%xmm1
movsd  0xe3f(%rip),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
subsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x28(%rbp),%rax
subsd  (%rax),%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x34(%rbp),%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x34(%rbp),%xmm1
movsd  0xde7(%rip),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    0x18(%rbp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
subsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x28(%rbp),%rax
subsd  (%rax),%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x34(%rbp),%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x34(%rbp),%xmm1
movsd  0xd8f(%rip),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    0x10(%rbp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
cvtsi2sdl -0x34(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsi2sdl -0x34(%rbp),%xmm1
movsd  0xd4c(%rip),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
cvtsi2sdl -0x34(%rbp),%xmm2
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsi2sdl -0x34(%rbp),%xmm1
movsd  0xd0d(%rip),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1179 <correl+0x69>
cvtsi2sdl -0x14(%rbp),%xmm1
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
cvtsi2sdl -0x14(%rbp),%xmm1
mov    0x18(%rbp),%rax
movsd  (%rax),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
cvtsi2sdl -0x14(%rbp),%xmm1
mov    0x10(%rbp),%rax
movsd  (%rax),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    13bd <correl+0x2ad>
mov    0x18(%rbp),%rax
movsd  (%rax),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    13bd <correl+0x2ad>
mov    0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
mov    0x18(%rbp),%rax
mulsd  (%rax),%xmm0
call   1030 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x40(%rbp),%xmm0
divsd  %xmm1,%xmm0
mov    0x20(%rbp),%rax
movsd  %xmm0,(%rax)
jmp    13cd <correl+0x2bd>
mov    0x20(%rbp),%rax
movsd  0xc37(%rip),%xmm0
movsd  %xmm0,(%rax)
xor    %eax,%eax
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<correl>:
push   %rbx
sub    $0x10,%rsp
mov    0x28(%rsp),%r10
mov    0x20(%rsp),%r11
movq   $0x0,(%r10)
movq   $0x0,(%r11)
movq   $0x0,(%r9)
movq   $0x0,(%r8)
movq   $0x0,(%rcx)
test   %edx,%edx
jle    1211 <correl+0x101>
mov    %edx,%eax
add    $0xffffffffffffffff,%rax
xorpd  %xmm1,%xmm1
xor    %ebx,%ebx
movsd  0xea2(%rip),%xmm0
xchg   %ax,%ax
movsd  (%rdi,%rbx,8),%xmm3
subsd  %xmm1,%xmm3
mulsd  %xmm3,%xmm3
xorps  %xmm1,%xmm1
cvtsi2sd %ebx,%xmm1
mulsd  %xmm1,%xmm3
movapd %xmm1,%xmm2
addsd  %xmm0,%xmm2
divsd  %xmm2,%xmm3
addsd  (%r9),%xmm3
movsd  %xmm3,(%r9)
movsd  (%rsi,%rbx,8),%xmm3
subsd  (%r8),%xmm3
mulsd  %xmm3,%xmm3
mulsd  %xmm1,%xmm3
divsd  %xmm2,%xmm3
addsd  (%r10),%xmm3
movsd  %xmm3,(%r10)
movsd  (%rdi,%rbx,8),%xmm3
subsd  (%rcx),%xmm3
movsd  (%rsi,%rbx,8),%xmm4
subsd  (%r8),%xmm4
mulsd  %xmm3,%xmm4
mulsd  %xmm1,%xmm4
divsd  %xmm2,%xmm4
addsd  (%r11),%xmm4
movsd  %xmm4,(%r11)
movsd  (%rcx),%xmm3
mulsd  %xmm1,%xmm3
addsd  (%rdi,%rbx,8),%xmm3
divsd  %xmm2,%xmm3
movsd  %xmm3,(%rcx)
mulsd  (%r8),%xmm1
addsd  (%rsi,%rbx,8),%xmm1
divsd  %xmm2,%xmm1
movsd  %xmm1,(%r8)
cmp    %rbx,%rax
je     1211 <correl+0x101>
movsd  (%rcx),%xmm1
add    $0x1,%rbx
jmp    1160 <correl+0x50>
mov    0x30(%rsp),%rbx
xorps  %xmm0,%xmm0
cvtsi2sd %edx,%xmm0
movsd  (%r9),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%r9)
movsd  (%r10),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%r10)
movsd  (%r11),%xmm4
divsd  %xmm0,%xmm4
movsd  %xmm4,(%r11)
movsd  (%r9),%xmm0
movsd  0xdb4(%rip),%xmm1
xorpd  %xmm3,%xmm3
ucomisd %xmm3,%xmm0
jbe    1296 <correl+0x186>
movsd  (%r10),%xmm2
ucomisd %xmm3,%xmm2
jbe    1296 <correl+0x186>
mulsd  %xmm2,%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     127d <correl+0x16d>
sqrtsd %xmm0,%xmm0
jmp    128e <correl+0x17e>
movsd  %xmm4,0x8(%rsp)
call   1030 <sqrt@plt>
movsd  0x8(%rsp),%xmm4
divsd  %xmm0,%xmm4
movapd %xmm4,%xmm1
movsd  %xmm1,(%rbx)
xor    %eax,%eax
add    $0x10,%rsp
pop    %rbx
ret

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

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