#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 *calloc(int, int);
int /*<<< orphan*/ fprintf(int /*<<< orphan*/, char *, long);
int /*<<< orphan*/ free(double *);
int /*<<< orphan*/ stderr;
int cov(double *x, double *S, int n, int p) {
  int i, j, k, p2 = p * p;
  double *mean = NULL;
  double n1 = 1.0 / n, n2 = 1.0 / (n - 1.0);
  if (!(mean = calloc(p, sizeof(double)))) {
    fprintf(stderr, "can't calloc %ld bytes\n", (long)p * sizeof(double));
    return 0;
  };
  for (j = 0; j < p; j++) {
    for (i = 0; i < n; i++)
      mean[j] += x[i * p + j];
    mean[j] *= n1;
  }
  for (j = 0; j < p; j++) {
    for (k = 0; k <= j; k++) {
      i = 0;
      S[j * p + k] = (x[i * p + k] - mean[k]) * (x[i * p + j] - mean[j]);
      for (i = 1; i < n; i++)
        S[j * p + k] += (x[i * p + k] - mean[k]) * (x[i * p + j] - mean[j]);
      S[k * p + j] = S[j * p + k];
    }
  }
  for (i = 0; i < p2; i++)
    S[i] *= n2;
  if (mean) {
    free(mean);
    mean = NULL;
  };
  return 1;
}


x86 assembly (clang, -O0):
<cov>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
mov    %ecx,-0x20(%rbp)
mov    -0x20(%rbp),%eax
imul   -0x20(%rbp),%eax
mov    %eax,-0x30(%rbp)
movq   $0x0,-0x38(%rbp)
cvtsi2sdl -0x1c(%rbp),%xmm1
movsd  0xe9b(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
cvtsi2sdl -0x1c(%rbp),%xmm1
movsd  0xe85(%rip),%xmm0
subsd  %xmm0,%xmm1
movsd  0xe79(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x20(%rbp),%edi
mov    $0x8,%esi
call   1040 <calloc@plt>
mov    %rax,-0x38(%rbp)
cmp    $0x0,%rax
jne    11d4 <cov+0xa4>
mov    0x2e2e(%rip),%rax
mov    (%rax),%edi
movslq -0x20(%rbp),%rdx
shl    $0x3,%rdx
lea    0xe45(%rip),%rsi
call   1050 <fprintf@plt>
movl   $0x0,-0x4(%rbp)
jmp    1414 <cov+0x2e4>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    1255 <cov+0x125>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1230 <cov+0x100>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x20(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x28(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11ee <cov+0xbe>
movsd  -0x40(%rbp),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x28(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    11db <cov+0xab>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    13b9 <cov+0x289>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x28(%rbp),%eax
jg     13a6 <cov+0x276>
movl   $0x0,-0x24(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x20(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x2c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x20(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x38(%rbp),%rax
movslq -0x28(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x20(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    136c <cov+0x23c>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x20(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x2c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x20(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x38(%rbp),%rax
movslq -0x28(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm2
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x20(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    12e9 <cov+0x1b9>
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x20(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x20(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    126f <cov+0x13f>
jmp    13ab <cov+0x27b>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    125c <cov+0x12c>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x30(%rbp),%eax
jge    13f1 <cov+0x2c1>
movsd  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    13c0 <cov+0x290>
cmpq   $0x0,-0x38(%rbp)
je     140d <cov+0x2dd>
mov    -0x38(%rbp),%rdi
call   1030 <free@plt>
movq   $0x0,-0x38(%rbp)
movl   $0x1,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cov>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %ecx,%r15d
mov    %edx,%r14d
mov    %rsi,%rbx
mov    %rdi,%rbp
mov    %ecx,%edi
mov    $0x8,%esi
call   1040 <calloc@plt>
test   %rax,%rax
je     1485 <cov+0x355>
mov    %r15d,%ecx
imul   %r15d,%ecx
mov    %ecx,0x18(%rsp)
mov    %rax,%rcx
cvtsi2sd %r14d,%xmm1
movsd  0xe86(%rip),%xmm0
test   %r15d,%r15d
jle    145c <cov+0x32c>
mov    %rax,%r13
movapd %xmm0,%xmm2
divsd  %xmm1,%xmm2
mov    %r15d,%ecx
mov    %r14d,%eax
mov    %r14d,%r9d
lea    -0x1(%r9),%rdx
mov    %rdx,0x10(%rsp)
mov    %r9d,%r8d
and    $0x3,%r8d
and    $0xfffffffc,%r9d
mov    %rcx,%rdi
shl    $0x4,%rdi
mov    %rcx,%r14
shl    $0x5,%r14
mov    %rcx,0x8(%rsp)
lea    0x0(,%rcx,8),%rdx
mov    %r13,%rcx
lea    (%rdx,%rdx,2),%r12
xor    %r11d,%r11d
mov    %rbp,%r10
mov    %rdx,0x20(%rsp)
jmp    1201 <cov+0xd1>
nopw   0x0(%rax,%rax,1)
mulsd  %xmm2,%xmm3
mov    %r13,%rcx
movsd  %xmm3,0x0(%r13,%r11,8)
add    $0x1,%r11
add    $0x8,%r10
cmp    0x8(%rsp),%r11
je     1292 <cov+0x162>
movsd  (%rcx,%r11,8),%xmm3
test   %eax,%eax
jle    11e0 <cov+0xb0>
cmpq   $0x3,0x10(%rsp)
jae    1220 <cov+0xf0>
xor    %esi,%esi
jmp    1250 <cov+0x120>
nopw   0x0(%rax,%rax,1)
mov    %r10,%rcx
xor    %esi,%esi
mov    0x8(%rsp),%rdx
nopw   0x0(%rax,%rax,1)
addsd  (%rcx),%xmm3
addsd  (%rcx,%rdx,8),%xmm3
addsd  (%rcx,%rdi,1),%xmm3
addsd  (%rcx,%r12,1),%xmm3
add    $0x4,%rsi
add    %r14,%rcx
cmp    %rsi,%r9
jne    1230 <cov+0x100>
test   %r8,%r8
mov    0x20(%rsp),%rdx
je     11e0 <cov+0xb0>
imul   0x8(%rsp),%rsi
add    %r11,%rsi
lea    0x0(,%rsi,8),%rsi
add    %rbp,%rsi
mov    %r8,%rcx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
addsd  (%rsi),%xmm3
add    %rdx,%rsi
add    $0xffffffffffffffff,%rcx
jne    1280 <cov+0x150>
jmp    11e0 <cov+0xb0>
test   %r15d,%r15d
mov    0x8(%rsp),%r8
jle    145c <cov+0x32c>
mov    0x10(%rsp),%rsi
and    $0xfffffffffffffffe,%rsi
neg    %rsi
mov    $0x1,%r10d
xor    %r11d,%r11d
mov    %rdi,0x30(%rsp)
mov    %rdx,0x28(%rsp)
mov    %rcx,0x38(%rsp)
mov    %eax,0x1c(%rsp)
jmp    12f5 <cov+0x1c5>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r11
add    $0x1,%r10
addq   $0x8,0x28(%rsp)
addq   $0x8,0x30(%rsp)
cmp    %r8,%r11
mov    0x20(%rsp),%rdx
mov    %r13,%rcx
je     145c <cov+0x32c>
mov    %rcx,%r13
movsd  (%rcx,%r11,8),%xmm2
mov    %r11,%rcx
imul   %r8,%rcx
mov    %rcx,0x50(%rsp)
mov    %rdi,%r14
xor    %r9d,%r9d
mov    %r10,0x48(%rsp)
mov    %r11,0x40(%rsp)
jmp    1340 <cov+0x210>
nopl   0x0(%rax)
mov    %r9,%rcx
imul   %r8,%rcx
add    %r11,%rcx
movsd  %xmm3,(%rbx,%rcx,8)
add    $0x1,%r9
add    $0x8,%rdx
add    $0x8,%r14
cmp    %r10,%r9
je     12d0 <cov+0x1a0>
mov    0x50(%rsp),%rcx
lea    (%r9,%rcx,1),%r15
movsd  0x0(%rbp,%r9,8),%xmm5
mov    %r13,%rcx
movsd  0x0(%r13,%r9,8),%xmm4
subsd  %xmm4,%xmm5
movsd  0x0(%rbp,%r11,8),%xmm3
subsd  %xmm2,%xmm3
mulsd  %xmm5,%xmm3
movsd  %xmm3,(%rbx,%r15,8)
cmp    $0x2,%eax
jl     1320 <cov+0x1f0>
mov    $0x1,%r12d
je     1406 <cov+0x2d6>
mov    $0x1,%r12d
mov    %r14,%r8
mov    0x30(%rsp),%r11
mov    %rdx,%r10
mov    0x28(%rsp),%r13
nopw   0x0(%rax,%rax,1)
movsd  0x0(%rbp,%r10,1),%xmm5
subsd  %xmm4,%xmm5
movsd  0x0(%rbp,%r13,1),%xmm6
subsd  %xmm2,%xmm6
mulsd  %xmm5,%xmm6
addsd  %xmm3,%xmm6
movsd  %xmm6,(%rbx,%r15,8)
movsd  0x0(%rbp,%r8,1),%xmm5
subsd  %xmm4,%xmm5
movsd  0x0(%rbp,%r11,1),%xmm3
subsd  %xmm2,%xmm3
mulsd  %xmm5,%xmm3
addsd  %xmm6,%xmm3
movsd  %xmm3,(%rbx,%r15,8)
add    %rdi,%r13
add    %rdi,%r10
add    %rdi,%r11
add    %rdi,%r8
lea    (%rsi,%r12,1),%rcx
add    $0x2,%rcx
add    $0x2,%r12
cmp    $0x1,%rcx
jne    13a0 <cov+0x270>
testb  $0x1,0x10(%rsp)
mov    0x38(%rsp),%r13
mov    0x1c(%rsp),%eax
mov    0x8(%rsp),%r8
mov    0x48(%rsp),%r10
mov    0x40(%rsp),%r11
je     1320 <cov+0x1f0>
imul   %r8,%r12
lea    (%r12,%r9,1),%rcx
movsd  0x0(%rbp,%rcx,8),%xmm5
subsd  %xmm4,%xmm5
add    %r11,%r12
movsd  0x0(%rbp,%r12,8),%xmm4
subsd  %xmm2,%xmm4
mulsd  %xmm5,%xmm4
addsd  %xmm4,%xmm3
movsd  %xmm3,(%rbx,%r15,8)
jmp    1320 <cov+0x1f0>
mov    %rcx,%rdi
mov    0x18(%rsp),%edx
test   %edx,%edx
je     1577 <cov+0x447>
addsd  0xb95(%rip),%xmm1
divsd  %xmm1,%xmm0
mov    %edx,%ecx
cmp    $0x4,%edx
jae    14a8 <cov+0x378>
xor    %edx,%edx
jmp    1560 <cov+0x430>
mov    0x2b54(%rip),%rax
mov    (%rax),%edi
movslq %r15d,%rdx
shl    $0x3,%rdx
lea    0xb74(%rip),%rsi
call   1050 <fprintf@plt>
xor    %eax,%eax
jmp    1581 <cov+0x451>
mov    %ecx,%edx
and    $0xfffffffc,%edx
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x4(%rdx),%rax
mov    %rax,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
test   %rax,%rax
je     1590 <cov+0x460>
mov    %rsi,%rax
and    $0xfffffffffffffffe,%rax
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
movupd (%rbx,%rbp,8),%xmm2
movupd 0x10(%rbx,%rbp,8),%xmm3
movupd 0x20(%rbx,%rbp,8),%xmm4
movupd 0x30(%rbx,%rbp,8),%xmm5
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,(%rbx,%rbp,8)
movupd %xmm3,0x10(%rbx,%rbp,8)
mulpd  %xmm1,%xmm4
mulpd  %xmm1,%xmm5
movupd %xmm4,0x20(%rbx,%rbp,8)
movupd %xmm5,0x30(%rbx,%rbp,8)
add    $0x8,%rbp
add    $0xfffffffffffffffe,%rax
jne    14e0 <cov+0x3b0>
test   $0x1,%sil
je     154c <cov+0x41c>
movupd (%rbx,%rbp,8),%xmm2
movupd 0x10(%rbx,%rbp,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,(%rbx,%rbp,8)
movupd %xmm3,0x10(%rbx,%rbp,8)
cmp    %rcx,%rdx
je     1577 <cov+0x447>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%rbx,%rdx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx,%rdx,8)
add    $0x1,%rdx
cmp    %rdx,%rcx
jne    1560 <cov+0x430>
call   1030 <free@plt>
mov    $0x1,%eax
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %ebp,%ebp
test   $0x1,%sil
jne    152e <cov+0x3fe>
jmp    154c <cov+0x41c>

LLVM IR (-O0):
@stderr = global i32 0
@.str = [24 x i8] c"can't calloc %ld bytes\0A\00"
define i32 @cov(double* %0, double* %1, i32 %2, i32 %3) {
%5 = alloca i32
%6 = alloca double*
%7 = alloca double*
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca double*
%15 = alloca double
%16 = alloca double
store double* %0, double** %6
store double* %1, double** %7
store i32 %2, i32* %8
store i32 %3, i32* %9
%17 = load i32, i32* %9
%18 = load i32, i32* %9
%19 = mul nsw i32 %17, %18
store i32 %19, i32* %13
store double* null, double** %14
%20 = load i32, i32* %8
%21 = sitofp i32 %20 to double
%22 = fdiv double 1.000000e+00, %21
store double %22, double* %15
%23 = load i32, i32* %8
%24 = sitofp i32 %23 to double
%25 = fsub double %24, 1.000000e+00
%26 = fdiv double 1.000000e+00, %25
store double %26, double* %16
%27 = load i32, i32* %9
%28 = call double* @calloc(i32 %27, i32 8)
store double* %28, double** %14
%29 = icmp ne double* %28, null
br i1 %29, label %36, label %30
30: ; preds = %4
%31 = load i32, i32* @stderr
%32 = load i32, i32* %9
%33 = sext i32 %32 to i64
%34 = mul i64 %33, 8
%35 = call i32 @fprintf(i32 %31, i8* getelementptr inbounds ([24 x i8], [24 x i8]* @.str, i64 0, i64 0), i64 %34)
store i32 0, i32* %5
br label %222
36: ; preds = %4
store i32 0, i32* %11
br label %37
37: ; preds = %73, %36
%38 = load i32, i32* %11
%39 = load i32, i32* %9
%40 = icmp slt i32 %38, %39
br i1 %40, label %41, label %76
41: ; preds = %37
store i32 0, i32* %10
br label %42
42: ; preds = %62, %41
%43 = load i32, i32* %10
%44 = load i32, i32* %8
%45 = icmp slt i32 %43, %44
br i1 %45, label %46, label %65
46: ; preds = %42
%47 = load double*, double** %6
%48 = load i32, i32* %10
%49 = load i32, i32* %9
%50 = mul nsw i32 %48, %49
%51 = load i32, i32* %11
%52 = add nsw i32 %50, %51
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds double, double* %47, i64 %53
%55 = load double, double* %54
%56 = load double*, double** %14
%57 = load i32, i32* %11
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds double, double* %56, i64 %58
%60 = load double, double* %59
%61 = fadd double %60, %55
store double %61, double* %59
br label %62
62: ; preds = %46
%63 = load i32, i32* %10
%64 = add nsw i32 %63, 1
store i32 %64, i32* %10
br label %42
65: ; preds = %42
%66 = load double, double* %15
%67 = load double*, double** %14
%68 = load i32, i32* %11
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds double, double* %67, i64 %69
%71 = load double, double* %70
%72 = fmul double %71, %66
store double %72, double* %70
br label %73
73: ; preds = %65
%74 = load i32, i32* %11
%75 = add nsw i32 %74, 1
store i32 %75, i32* %11
br label %37
76: ; preds = %37
store i32 0, i32* %11
br label %77
77: ; preds = %196, %76
%78 = load i32, i32* %11
%79 = load i32, i32* %9
%80 = icmp slt i32 %78, %79
br i1 %80, label %81, label %199
81: ; preds = %77
store i32 0, i32* %12
br label %82
82: ; preds = %192, %81
%83 = load i32, i32* %12
%84 = load i32, i32* %11
%85 = icmp sle i32 %83, %84
br i1 %85, label %86, label %195
86: ; preds = %82
store i32 0, i32* %10
%87 = load double*, double** %6
%88 = load i32, i32* %10
%89 = load i32, i32* %9
%90 = mul nsw i32 %88, %89
%91 = load i32, i32* %12
%92 = add nsw i32 %90, %91
%93 = sext i32 %92 to i64
%94 = getelementptr inbounds double, double* %87, i64 %93
%95 = load double, double* %94
%96 = load double*, double** %14
%97 = load i32, i32* %12
%98 = sext i32 %97 to i64
%99 = getelementptr inbounds double, double* %96, i64 %98
%100 = load double, double* %99
%101 = fsub double %95, %100
%102 = load double*, double** %6
%103 = load i32, i32* %10
%104 = load i32, i32* %9
%105 = mul nsw i32 %103, %104
%106 = load i32, i32* %11
%107 = add nsw i32 %105, %106
%108 = sext i32 %107 to i64
%109 = getelementptr inbounds double, double* %102, i64 %108
%110 = load double, double* %109
%111 = load double*, double** %14
%112 = load i32, i32* %11
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds double, double* %111, i64 %113
%115 = load double, double* %114
%116 = fsub double %110, %115
%117 = fmul double %101, %116
%118 = load double*, double** %7
%119 = load i32, i32* %11
%120 = load i32, i32* %9
%121 = mul nsw i32 %119, %120
%122 = load i32, i32* %12
%123 = add nsw i32 %121, %122
%124 = sext i32 %123 to i64
%125 = getelementptr inbounds double, double* %118, i64 %124
store double %117, double* %125
store i32 1, i32* %10
br label %126
126: ; preds = %171, %86
%127 = load i32, i32* %10
%128 = load i32, i32* %8
%129 = icmp slt i32 %127, %128
br i1 %129, label %130, label %174
130: ; preds = %126
%131 = load double*, double** %6
%132 = load i32, i32* %10
%133 = load i32, i32* %9
%134 = mul nsw i32 %132, %133
%135 = load i32, i32* %12
%136 = add nsw i32 %134, %135
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds double, double* %131, i64 %137
%139 = load double, double* %138
%140 = load double*, double** %14
%141 = load i32, i32* %12
%142 = sext i32 %141 to i64
%143 = getelementptr inbounds double, double* %140, i64 %142
%144 = load double, double* %143
%145 = fsub double %139, %144
%146 = load double*, double** %6
%147 = load i32, i32* %10
%148 = load i32, i32* %9
%149 = mul nsw i32 %147, %148
%150 = load i32, i32* %11
%151 = add nsw i32 %149, %150
%152 = sext i32 %151 to i64
%153 = getelementptr inbounds double, double* %146, i64 %152
%154 = load double, double* %153
%155 = load double*, double** %14
%156 = load i32, i32* %11
%157 = sext i32 %156 to i64
%158 = getelementptr inbounds double, double* %155, i64 %157
%159 = load double, double* %158
%160 = fsub double %154, %159
%161 = load double*, double** %7
%162 = load i32, i32* %11
%163 = load i32, i32* %9
%164 = mul nsw i32 %162, %163
%165 = load i32, i32* %12
%166 = add nsw i32 %164, %165
%167 = sext i32 %166 to i64
%168 = getelementptr inbounds double, double* %161, i64 %167
%169 = load double, double* %168
%170 = call double @llvm.fmuladd.f64(double %145, double %160, double %169)
store double %170, double* %168
br label %171
171: ; preds = %130
%172 = load i32, i32* %10
%173 = add nsw i32 %172, 1
store i32 %173, i32* %10
br label %126
174: ; preds = %126
%175 = load double*, double** %7
%176 = load i32, i32* %11
%177 = load i32, i32* %9
%178 = mul nsw i32 %176, %177
%179 = load i32, i32* %12
%180 = add nsw i32 %178, %179
%181 = sext i32 %180 to i64
%182 = getelementptr inbounds double, double* %175, i64 %181
%183 = load double, double* %182
%184 = load double*, double** %7
%185 = load i32, i32* %12
%186 = load i32, i32* %9
%187 = mul nsw i32 %185, %186
%188 = load i32, i32* %11
%189 = add nsw i32 %187, %188
%190 = sext i32 %189 to i64
%191 = getelementptr inbounds double, double* %184, i64 %190
store double %183, double* %191
br label %192
192: ; preds = %174
%193 = load i32, i32* %12
%194 = add nsw i32 %193, 1
store i32 %194, i32* %12
br label %82
195: ; preds = %82
br label %196
196: ; preds = %195
%197 = load i32, i32* %11
%198 = add nsw i32 %197, 1
store i32 %198, i32* %11
br label %77
199: ; preds = %77
store i32 0, i32* %10
br label %200
200: ; preds = %212, %199
%201 = load i32, i32* %10
%202 = load i32, i32* %13
%203 = icmp slt i32 %201, %202
br i1 %203, label %204, label %215
204: ; preds = %200
%205 = load double, double* %16
%206 = load double*, double** %7
%207 = load i32, i32* %10
%208 = sext i32 %207 to i64
%209 = getelementptr inbounds double, double* %206, i64 %208
%210 = load double, double* %209
%211 = fmul double %210, %205
store double %211, double* %209
br label %212
212: ; preds = %204
%213 = load i32, i32* %10
%214 = add nsw i32 %213, 1
store i32 %214, i32* %10
br label %200
215: ; preds = %200
%216 = load double*, double** %14
%217 = icmp ne double* %216, null
br i1 %217, label %218, label %221
218: ; preds = %215
%219 = load double*, double** %14
%220 = call i32 @free(double* %219)
store double* null, double** %14
br label %221
221: ; preds = %218, %215
store i32 1, i32* %5
br label %222
222: ; preds = %221, %30
%223 = load i32, i32* %5
ret i32 %223
}

LLVM IR (-O2):
@stderr = global i32 0
@.str = [24 x i8] c"can't calloc %ld bytes\0A\00"
define i32 @cov(double* nocapture %0, double* nocapture %1, i32 %2, i32 %3) {
%5 = mul i32 %3, %3
%6 = sitofp i32 %2 to double
%7 = fdiv double 1.000000e+00, %6
%8 = fadd double %6, -1.000000e+00
%9 = fdiv double 1.000000e+00, %8
%10 = tail call double* @calloc(i32 %3, i32 8)
%11 = icmp eq double* %10, null
br i1 %11, label %24, label %12
12: ; preds = %4
%13 = icmp sgt i32 %3, 0
br i1 %13, label %14, label %102
14: ; preds = %12
%15 = icmp sgt i32 %2, 0
%16 = zext i32 %3 to i64
%17 = zext i32 %3 to i64
%18 = zext i32 %2 to i64
%19 = add nsw i64 %18, -1
%20 = and i64 %18, 3
%21 = icmp ult i64 %19, 3
%22 = and i64 %18, 4294967292
%23 = icmp eq i64 %20, 0
br label %29
24: ; preds = %4
%25 = load i32, i32* @stderr
%26 = sext i32 %3 to i64
%27 = shl nsw i64 %26, 3
%28 = tail call i32 @fprintf(i32 %25, i8* getelementptr inbounds ([24 x i8], [24 x i8]* @.str, i64 0, i64 0), i64 %27)
br label %241
29: ; preds = %14, %89
%30 = phi i64 [ 0, %14 ], [ %93, %89 ]
%31 = getelementptr inbounds double, double* %10, i64 %30
%32 = load double, double* %31
br i1 %15, label %33, label %89
33: ; preds = %29
br i1 %21, label %73, label %43
34: ; preds = %89
br i1 %13, label %35, label %102
35: ; preds = %34
%36 = icmp sgt i32 %2, 1
%37 = zext i32 %3 to i64
%38 = zext i32 %3 to i64
%39 = and i64 %19, 1
%40 = icmp eq i32 %2, 2
%41 = and i64 %19, -2
%42 = icmp eq i64 %39, 0
br label %95
43: ; preds = %33, %43
%44 = phi i64 [ %70, %43 ], [ 0, %33 ]
%45 = phi double [ %69, %43 ], [ %32, %33 ]
%46 = phi i64 [ %71, %43 ], [ 0, %33 ]
%47 = mul nsw i64 %44, %16
%48 = add nuw nsw i64 %47, %30
%49 = getelementptr inbounds double, double* %0, i64 %48
%50 = load double, double* %49
%51 = fadd double %50, %45
%52 = or i64 %44, 1
%53 = mul nsw i64 %52, %16
%54 = add nuw nsw i64 %53, %30
%55 = getelementptr inbounds double, double* %0, i64 %54
%56 = load double, double* %55
%57 = fadd double %56, %51
%58 = or i64 %44, 2
%59 = mul nsw i64 %58, %16
%60 = add nuw nsw i64 %59, %30
%61 = getelementptr inbounds double, double* %0, i64 %60
%62 = load double, double* %61
%63 = fadd double %62, %57
%64 = or i64 %44, 3
%65 = mul nsw i64 %64, %16
%66 = add nuw nsw i64 %65, %30
%67 = getelementptr inbounds double, double* %0, i64 %66
%68 = load double, double* %67
%69 = fadd double %68, %63
%70 = add nuw nsw i64 %44, 4
%71 = add i64 %46, 4
%72 = icmp eq i64 %71, %22
br i1 %72, label %73, label %43
73: ; preds = %43, %33
%74 = phi double [ undef, %33 ], [ %69, %43 ]
%75 = phi i64 [ 0, %33 ], [ %70, %43 ]
%76 = phi double [ %32, %33 ], [ %69, %43 ]
br i1 %23, label %89, label %77
77: ; preds = %73, %77
%78 = phi i64 [ %86, %77 ], [ %75, %73 ]
%79 = phi double [ %85, %77 ], [ %76, %73 ]
%80 = phi i64 [ %87, %77 ], [ 0, %73 ]
%81 = mul nsw i64 %78, %16
%82 = add nuw nsw i64 %81, %30
%83 = getelementptr inbounds double, double* %0, i64 %82
%84 = load double, double* %83
%85 = fadd double %84, %79
%86 = add nuw nsw i64 %78, 1
%87 = add i64 %80, 1
%88 = icmp eq i64 %87, %20
br i1 %88, label %89, label %77
89: ; preds = %73, %77, %29
%90 = phi double [ %32, %29 ], [ %74, %73 ], [ %85, %77 ]
%91 = getelementptr inbounds double, double* %10, i64 %30
%92 = fmul double %7, %90
store double %92, double* %91
%93 = add nuw nsw i64 %30, 1
%94 = icmp eq i64 %93, %17
br i1 %94, label %34, label %29
95: ; preds = %35, %228
%96 = phi i64 [ 0, %35 ], [ %229, %228 ]
%97 = phi i64 [ 1, %35 ], [ %230, %228 ]
%98 = getelementptr inbounds double, double* %0, i64 %96
%99 = getelementptr inbounds double, double* %10, i64 %96
%100 = load double, double* %99
%101 = mul nsw i64 %96, %37
br label %165
102: ; preds = %228, %12, %34
%103 = icmp eq i32 %5, 0
br i1 %103, label %239, label %104
104: ; preds = %102
%105 = zext i32 %5 to i64
%106 = icmp ult i32 %5, 4
br i1 %106, label %163, label %107
107: ; preds = %104
%108 = and i64 %105, 4294967292
%109 = insertelement <2 x double> poison, double %9, i64 0
%110 = shufflevector <2 x double> %109, <2 x double> poison, <2 x i32> zeroinitializer
%111 = insertelement <2 x double> poison, double %9, i64 0
%112 = shufflevector <2 x double> %111, <2 x double> poison, <2 x i32> zeroinitializer
%113 = add nsw i64 %108, -4
%114 = lshr exact i64 %113, 2
%115 = add nuw nsw i64 %114, 1
%116 = and i64 %115, 1
%117 = icmp eq i64 %113, 0
br i1 %117, label %147, label %118
118: ; preds = %107
%119 = and i64 %115, 9223372036854775806
br label %120
120: ; preds = %120, %118
%121 = phi i64 [ 0, %118 ], [ %144, %120 ]
%122 = phi i64 [ 0, %118 ], [ %145, %120 ]
%123 = getelementptr inbounds double, double* %1, i64 %121
%124 = bitcast double* %123 to <2 x double>*
%125 = load <2 x double>, <2 x double>* %124
%126 = getelementptr inbounds double, double* %123, i64 2
%127 = bitcast double* %126 to <2 x double>*
%128 = load <2 x double>, <2 x double>* %127
%129 = fmul <2 x double> %110, %125
%130 = fmul <2 x double> %112, %128
%131 = bitcast double* %123 to <2 x double>*
store <2 x double> %129, <2 x double>* %131
%132 = bitcast double* %126 to <2 x double>*
store <2 x double> %130, <2 x double>* %132
%133 = or i64 %121, 4
%134 = getelementptr inbounds double, double* %1, i64 %133
%135 = bitcast double* %134 to <2 x double>*
%136 = load <2 x double>, <2 x double>* %135
%137 = getelementptr inbounds double, double* %134, i64 2
%138 = bitcast double* %137 to <2 x double>*
%139 = load <2 x double>, <2 x double>* %138
%140 = fmul <2 x double> %110, %136
%141 = fmul <2 x double> %112, %139
%142 = bitcast double* %134 to <2 x double>*
store <2 x double> %140, <2 x double>* %142
%143 = bitcast double* %137 to <2 x double>*
store <2 x double> %141, <2 x double>* %143
%144 = add nuw i64 %121, 8
%145 = add i64 %122, 2
%146 = icmp eq i64 %145, %119
br i1 %146, label %147, label %120
147: ; preds = %120, %107
%148 = phi i64 [ 0, %107 ], [ %144, %120 ]
%149 = icmp eq i64 %116, 0
br i1 %149, label %161, label %150
150: ; preds = %147
%151 = getelementptr inbounds double, double* %1, i64 %148
%152 = bitcast double* %151 to <2 x double>*
%153 = load <2 x double>, <2 x double>* %152
%154 = getelementptr inbounds double, double* %151, i64 2
%155 = bitcast double* %154 to <2 x double>*
%156 = load <2 x double>, <2 x double>* %155
%157 = fmul <2 x double> %110, %153
%158 = fmul <2 x double> %112, %156
%159 = bitcast double* %151 to <2 x double>*
store <2 x double> %157, <2 x double>* %159
%160 = bitcast double* %154 to <2 x double>*
store <2 x double> %158, <2 x double>* %160
br label %161
161: ; preds = %147, %150
%162 = icmp eq i64 %108, %105
br i1 %162, label %239, label %163
163: ; preds = %104, %161
%164 = phi i64 [ 0, %104 ], [ %108, %161 ]
br label %232
165: ; preds = %95, %221
%166 = phi i64 [ 0, %95 ], [ %226, %221 ]
%167 = getelementptr inbounds double, double* %0, i64 %166
%168 = load double, double* %167
%169 = getelementptr inbounds double, double* %10, i64 %166
%170 = load double, double* %169
%171 = fsub double %168, %170
%172 = load double, double* %98
%173 = fsub double %172, %100
%174 = fmul double %171, %173
%175 = add nuw nsw i64 %166, %101
%176 = getelementptr inbounds double, double* %1, i64 %175
store double %174, double* %176
br i1 %36, label %177, label %221
177: ; preds = %165
br i1 %40, label %206, label %178
178: ; preds = %177, %178
%179 = phi double [ %202, %178 ], [ %174, %177 ]
%180 = phi i64 [ %203, %178 ], [ 1, %177 ]
%181 = phi i64 [ %204, %178 ], [ 0, %177 ]
%182 = mul nsw i64 %180, %37
%183 = add nuw nsw i64 %182, %166
%184 = getelementptr inbounds double, double* %0, i64 %183
%185 = load double, double* %184
%186 = fsub double %185, %170
%187 = add nuw nsw i64 %182, %96
%188 = getelementptr inbounds double, double* %0, i64 %187
%189 = load double, double* %188
%190 = fsub double %189, %100
%191 = tail call double @llvm.fmuladd.f64(double %186, double %190, double %179)
store double %191, double* %176
%192 = add nuw nsw i64 %180, 1
%193 = mul nsw i64 %192, %37
%194 = add nuw nsw i64 %193, %166
%195 = getelementptr inbounds double, double* %0, i64 %194
%196 = load double, double* %195
%197 = fsub double %196, %170
%198 = add nuw nsw i64 %193, %96
%199 = getelementptr inbounds double, double* %0, i64 %198
%200 = load double, double* %199
%201 = fsub double %200, %100
%202 = tail call double @llvm.fmuladd.f64(double %197, double %201, double %191)
store double %202, double* %176
%203 = add nuw nsw i64 %180, 2
%204 = add i64 %181, 2
%205 = icmp eq i64 %204, %41
br i1 %205, label %206, label %178
206: ; preds = %178, %177
%207 = phi double [ undef, %177 ], [ %202, %178 ]
%208 = phi double [ %174, %177 ], [ %202, %178 ]
%209 = phi i64 [ 1, %177 ], [ %203, %178 ]
br i1 %42, label %221, label %210
210: ; preds = %206
%211 = mul nsw i64 %209, %37
%212 = add nuw nsw i64 %211, %166
%213 = getelementptr inbounds double, double* %0, i64 %212
%214 = load double, double* %213
%215 = fsub double %214, %170
%216 = add nuw nsw i64 %211, %96
%217 = getelementptr inbounds double, double* %0, i64 %216
%218 = load double, double* %217
%219 = fsub double %218, %100
%220 = tail call double @llvm.fmuladd.f64(double %215, double %219, double %208)
store double %220, double* %176
br label %221
221: ; preds = %210, %206, %165
%222 = phi double [ %174, %165 ], [ %207, %206 ], [ %220, %210 ]
%223 = mul nsw i64 %166, %37
%224 = add nuw nsw i64 %223, %96
%225 = getelementptr inbounds double, double* %1, i64 %224
store double %222, double* %225
%226 = add nuw nsw i64 %166, 1
%227 = icmp eq i64 %226, %97
br i1 %227, label %228, label %165
228: ; preds = %221
%229 = add nuw nsw i64 %96, 1
%230 = add nuw nsw i64 %97, 1
%231 = icmp eq i64 %229, %38
br i1 %231, label %102, label %95
232: ; preds = %163, %232
%233 = phi i64 [ %237, %232 ], [ %164, %163 ]
%234 = getelementptr inbounds double, double* %1, i64 %233
%235 = load double, double* %234
%236 = fmul double %9, %235
store double %236, double* %234
%237 = add nuw nsw i64 %233, 1
%238 = icmp eq i64 %237, %105
br i1 %238, label %239, label %232
239: ; preds = %232, %161, %102
%240 = tail call i32 @free(double* %10)
br label %241
241: ; preds = %239, %24
%242 = phi i32 [ 1, %239 ], [ 0, %24 ]
ret i32 %242
}