#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 cov2(double *x, double *S, int n, int p) {
  int i, j, k, p2 = p * p, jpk;
  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 + n * j];
    mean[j] *= n1;
  }
  for (j = 0; j < p; j++) {
    for (k = 0; k <= j; k++) {
      i = 0;
      jpk = j + p * k;
      S[jpk] = (x[i + n * k] - mean[k]) * (x[i + n * j] - mean[j]);
      for (i = 1; i < n; i++)
        S[jpk] += (x[i + n * k] - mean[k]) * (x[i + n * j] - mean[j]);
      S[k + p * j] = S[jpk];
    }
  }
  for (i = 0; i < p2; i++)
    S[i] *= n2;
  if (mean) {
    free(mean);
    mean = NULL;
  };
  return 1;
}


x86 assembly (clang, -O0):
<cov2>:
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,-0x40(%rbp)
cvtsi2sdl -0x1c(%rbp),%xmm1
movsd  0xe9b(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
cvtsi2sdl -0x1c(%rbp),%xmm1
movsd  0xe85(%rip),%xmm0
subsd  %xmm0,%xmm1
movsd  0xe79(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x20(%rbp),%edi
mov    $0x8,%esi
call   1040 <calloc@plt>
mov    %rax,-0x40(%rbp)
cmp    $0x0,%rax
jne    11d4 <cov2+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 <cov2+0x2e4>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    1257 <cov2+0x127>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1232 <cov2+0x102>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x1c(%rbp),%edx
imul   -0x28(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%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 <cov2+0xbe>
movsd  -0x48(%rbp),%xmm0
mov    -0x40(%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 <cov2+0xab>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    13b9 <cov2+0x289>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x28(%rbp),%eax
jg     13a6 <cov2+0x276>
movl   $0x0,-0x24(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x20(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x34(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x1c(%rbp),%edx
imul   -0x2c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
movslq -0x2c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x1c(%rbp),%edx
imul   -0x28(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x40(%rbp),%rax
movslq -0x28(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1373 <cov2+0x243>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x1c(%rbp),%edx
imul   -0x2c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
movslq -0x2c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x1c(%rbp),%edx
imul   -0x28(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x40(%rbp),%rax
movslq -0x28(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm2
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%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    12f5 <cov2+0x1c5>
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    -0x20(%rbp),%edx
imul   -0x28(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1271 <cov2+0x141>
jmp    13ab <cov2+0x27b>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    125e <cov2+0x12e>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x30(%rbp),%eax
jge    13f1 <cov2+0x2c1>
movsd  -0x50(%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 <cov2+0x290>
cmpq   $0x0,-0x40(%rbp)
je     140d <cov2+0x2dd>
mov    -0x40(%rbp),%rdi
call   1030 <free@plt>
movq   $0x0,-0x40(%rbp)
movl   $0x1,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cov2>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %ecx,%r15d
mov    %edx,%r12d
mov    %rsi,%rbx
mov    %rdi,%rbp
mov    %ecx,%edi
mov    $0x8,%esi
call   1040 <calloc@plt>
mov    %rax,(%rsp)
test   %rax,%rax
je     141d <cov2+0x2ed>
mov    %r15d,%ecx
imul   %r15d,%ecx
mov    %ecx,0x18(%rsp)
cvtsi2sd %r12d,%xmm1
movsd  0xe85(%rip),%xmm0
test   %r15d,%r15d
jle    13f7 <cov2+0x2c7>
movapd %xmm0,%xmm2
divsd  %xmm1,%xmm2
movslq %r12d,%rcx
mov    %r15d,%edx
mov    %r12d,%r10d
lea    -0x1(%r10),%rax
mov    %rax,0x8(%rsp)
mov    %r10d,%r11d
and    $0x3,%r11d
and    $0xfffffffc,%r10d
lea    0x18(%rbp),%rsi
mov    %rcx,0x28(%rsp)
lea    0x0(,%rcx,8),%rax
mov    %rax,0x10(%rsp)
xor    %r9d,%r9d
mov    %rbp,%r8
jmp    11fc <cov2+0xcc>
nopw   0x0(%rax,%rax,1)
mov    (%rsp),%rax
movsd  (%rax,%r9,8),%xmm3
mulsd  %xmm2,%xmm3
mov    (%rsp),%rax
movsd  %xmm3,(%rax,%r9,8)
add    $0x1,%r9
mov    0x10(%rsp),%rax
add    %rax,%rsi
add    %rax,%r8
cmp    %rdx,%r9
je     1273 <cov2+0x143>
test   %r12d,%r12d
jle    11d0 <cov2+0xa0>
mov    (%rsp),%rax
movsd  (%rax,%r9,8),%xmm3
cmpq   $0x3,0x8(%rsp)
jae    1220 <cov2+0xf0>
xor    %edi,%edi
jmp    1250 <cov2+0x120>
nopw   0x0(%rax,%rax,1)
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
addsd  -0x18(%rsi,%rdi,8),%xmm3
addsd  -0x10(%rsi,%rdi,8),%xmm3
addsd  -0x8(%rsi,%rdi,8),%xmm3
addsd  (%rsi,%rdi,8),%xmm3
add    $0x4,%rdi
cmp    %rdi,%r10
jne    1230 <cov2+0x100>
test   %r11,%r11
je     11da <cov2+0xaa>
lea    (%r8,%rdi,8),%rcx
xor    %edi,%edi
nopl   0x0(%rax,%rax,1)
addsd  (%rcx,%rdi,8),%xmm3
add    $0x1,%rdi
cmp    %rdi,%r11
jne    1260 <cov2+0x130>
jmp    11da <cov2+0xaa>
test   %r15d,%r15d
jle    13f7 <cov2+0x2c7>
lea    0x8(%rbp),%rcx
mov    0x8(%rsp),%rdi
and    $0xfffffffffffffffe,%rdi
neg    %rdi
mov    $0x1,%eax
mov    %rax,0x20(%rsp)
xor    %r8d,%r8d
mov    %rcx,0x30(%rsp)
mov    %r12d,0x1c(%rsp)
jmp    12c8 <cov2+0x198>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%r8
addq   $0x1,0x20(%rsp)
add    0x10(%rsp),%rcx
cmp    %rdx,%r8
je     13f7 <cov2+0x2c7>
mov    %r8,%rsi
imul   0x28(%rsp),%rsi
mov    (%rsp),%rax
movsd  (%rax,%r8,8),%xmm2
mov    %r8,%rax
imul   %rdx,%rax
mov    0x30(%rsp),%r9
xor    %r14d,%r14d
jmp    130c <cov2+0x1dc>
nopl   0x0(%rax)
lea    (%r14,%rax,1),%rdx
movsd  %xmm3,(%rbx,%rdx,8)
add    $0x1,%r14
add    0x10(%rsp),%r9
cmp    0x20(%rsp),%r14
mov    %r15,%rdx
je     12b0 <cov2+0x180>
mov    %r14,%r10
mov    %rdx,%r15
imul   %rdx,%r10
add    %r8,%r10
mov    %r14,%r11
imul   0x28(%rsp),%r11
movsd  0x0(%rbp,%r11,8),%xmm5
mov    (%rsp),%rdx
movsd  (%rdx,%r14,8),%xmm4
subsd  %xmm4,%xmm5
movsd  0x0(%rbp,%rsi,8),%xmm3
subsd  %xmm2,%xmm3
mulsd  %xmm5,%xmm3
movsd  %xmm3,(%rbx,%r10,8)
cmp    $0x2,%r12d
jl     12f0 <cov2+0x1c0>
mov    $0x1,%r13d
je     13b8 <cov2+0x288>
mov    $0x1,%r13d
nop
movsd  -0x8(%r9,%r13,8),%xmm5
subsd  %xmm4,%xmm5
movsd  -0x8(%rcx,%r13,8),%xmm6
subsd  %xmm2,%xmm6
mulsd  %xmm5,%xmm6
addsd  %xmm3,%xmm6
movsd  %xmm6,(%rbx,%r10,8)
movsd  (%r9,%r13,8),%xmm5
subsd  %xmm4,%xmm5
movsd  (%rcx,%r13,8),%xmm3
subsd  %xmm2,%xmm3
mulsd  %xmm5,%xmm3
addsd  %xmm6,%xmm3
movsd  %xmm3,(%rbx,%r10,8)
lea    (%rdi,%r13,1),%r12
add    $0x2,%r12
add    $0x2,%r13
cmp    $0x1,%r12
jne    1360 <cov2+0x230>
testb  $0x1,0x8(%rsp)
mov    0x1c(%rsp),%r12d
je     12f0 <cov2+0x1c0>
add    %r13,%r11
movsd  0x0(%rbp,%r11,8),%xmm5
subsd  %xmm4,%xmm5
add    %rsi,%r13
movsd  0x0(%rbp,%r13,8),%xmm4
subsd  %xmm2,%xmm4
mulsd  %xmm5,%xmm4
addsd  %xmm4,%xmm3
movsd  %xmm3,(%rbx,%r10,8)
jmp    12f0 <cov2+0x1c0>
mov    0x18(%rsp),%edx
test   %edx,%edx
je     1507 <cov2+0x3d7>
addsd  0xbfd(%rip),%xmm1
divsd  %xmm1,%xmm0
mov    %edx,%ecx
cmp    $0x4,%edx
jae    1440 <cov2+0x310>
xor    %edx,%edx
jmp    14f0 <cov2+0x3c0>
mov    0x2bbc(%rip),%rax
mov    (%rax),%edi
movslq %r15d,%rdx
shl    $0x3,%rdx
lea    0xbdc(%rip),%rsi
call   1050 <fprintf@plt>
xor    %eax,%eax
jmp    1515 <cov2+0x3e5>
mov    %ecx,%edx
and    $0xfffffffc,%edx
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x4(%rdx),%rdi
mov    %rdi,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
test   %rdi,%rdi
je     1524 <cov2+0x3f4>
mov    %rsi,%rdi
and    $0xfffffffffffffffe,%rdi
xor    %ebp,%ebp
xchg   %ax,%ax
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,%rdi
jne    1470 <cov2+0x340>
test   $0x1,%sil
je     14dc <cov2+0x3ac>
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     1507 <cov2+0x3d7>
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    14f0 <cov2+0x3c0>
mov    (%rsp),%rdi
call   1030 <free@plt>
mov    $0x1,%eax
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %ebp,%ebp
test   $0x1,%sil
jne    14be <cov2+0x38e>
jmp    14dc <cov2+0x3ac>

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

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