#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);
void cov2cor(double *S, double *P, int p) {
  double z;
  int i, j;
  for (i = 1; i < p; i++) {
    for (j = 0; j < i; j++) {
      z = S[i * p + j] / sqrt(S[i * p + i] * S[j * p + j]);
      P[i * p + j] = P[j * p + i] = z;
    }
  }
  for (i = 0; i < p; i++)
    P[i * p + i] = 1.0;
}


x86 assembly (clang, -O0):
<cov2cor>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11f8 <cov2cor+0xe8>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    11e5 <cov2cor+0xd5>
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
call   1030 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x30(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    113d <cov2cor+0x2d>
jmp    11ea <cov2cor+0xda>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    112a <cov2cor+0x1a>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1237 <cov2cor+0x127>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  0xddc(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11ff <cov2cor+0xef>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cov2cor>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %rsi,0x8(%rsp)
mov    %edx,%ebx
cmp    $0x2,%edx
jl     121a <cov2cor+0x10a>
mov    0x8(%rsp),%rax
lea    0x8(%rax),%r14
lea    (%rdi,%rbx,8),%r12
lea    0x0(,%rbx,8),%rbp
mov    %rbp,0x10(%rsp)
lea    (%rax,%rbx,8),%rcx
mov    $0x1,%r15d
xorpd  %xmm1,%xmm1
mov    %rdx,0x18(%rsp)
mov    %rdi,0x20(%rsp)
mov    %rbx,(%rsp)
jmp    1191 <cov2cor+0x81>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%r15
add    $0x8,%r14
mov    0x10(%rsp),%rax
add    %rax,%r12
add    %rax,%rcx
cmp    %rbx,%r15
mov    0x18(%rsp),%rdx
je     121a <cov2cor+0x10a>
mov    %r15d,%eax
imul   %edx,%eax
movslq %eax,%rdx
add    %r15,%rdx
xor    %r13d,%r13d
xor    %ebp,%ebp
mov    %rdx,0x28(%rsp)
jmp    11f8 <cov2cor+0xe8>
nopl   0x0(%rax)
mov    %rcx,%rbx
movsd  %xmm2,0x30(%rsp)
call   1030 <sqrt@plt>
movsd  0x30(%rsp),%xmm2
mov    0x28(%rsp),%rdx
xorpd  %xmm1,%xmm1
mov    0x20(%rsp),%rdi
mov    %rbx,%rcx
mov    (%rsp),%rbx
divsd  %xmm0,%xmm2
movsd  %xmm2,(%r14,%r13,8)
movsd  %xmm2,(%rcx,%rbp,8)
add    $0x1,%rbp
add    %rbx,%r13
cmp    %rbp,%r15
je     1170 <cov2cor+0x60>
movslq %r13d,%rax
movsd  (%r12,%rbp,8),%xmm2
movsd  (%rdi,%rdx,8),%xmm0
add    %rbp,%rax
mulsd  (%rdi,%rax,8),%xmm0
ucomisd %xmm1,%xmm0
jb     11b0 <cov2cor+0xa0>
sqrtsd %xmm0,%xmm0
jmp    11d9 <cov2cor+0xc9>
test   %edx,%edx
mov    0x8(%rsp),%r13
jle    1317 <cov2cor+0x207>
lea    -0x1(%rbx),%rax
mov    %ebx,%r10d
and    $0x3,%r10d
xor    %edi,%edi
cmp    $0x3,%rax
jae    1252 <cov2cor+0x142>
movabs $0x3ff0000000000000,%r15
test   %r10,%r10
jne    12f1 <cov2cor+0x1e1>
jmp    1317 <cov2cor+0x207>
lea    (%rdx,%rdx,1),%esi
lea    0x0(,%rdx,4),%r9d
lea    (%rdx,%rdx,2),%ebp
and    $0xfffffffc,%ebx
neg    %rbx
mov    %rbx,(%rsp)
mov    $0x2,%r11d
mov    $0x1,%r14d
mov    $0x3,%r8d
xor    %r12d,%r12d
mov    %edx,%ebx
xor    %ecx,%ecx
movabs $0x3ff0000000000000,%r15
nopl   (%rax)
movslq %esi,%rsi
movslq %ebx,%rbx
movslq %edi,%rdi
lea    (%rcx,%rdi,1),%rax
mov    %r15,0x0(%r13,%rax,8)
lea    (%r14,%rbx,1),%rax
mov    %r15,0x0(%r13,%rax,8)
lea    (%r11,%rsi,1),%rax
movslq %ebp,%rbp
mov    %r15,0x0(%r13,%rax,8)
lea    (%r8,%rbp,1),%rax
mov    %r15,0x0(%r13,%rax,8)
add    $0x4,%rcx
add    %r9d,%edi
add    $0x4,%r11
add    %r9d,%esi
add    $0x4,%r14
add    %r9d,%ebx
add    $0x4,%r8
add    %r9d,%ebp
add    $0xfffffffffffffffc,%r12
cmp    %r12,(%rsp)
jne    1290 <cov2cor+0x180>
neg    %r12
mov    %r12,%rdi
test   %r10,%r10
je     1317 <cov2cor+0x207>
mov    %edx,%eax
imul   %edi,%eax
cs nopw 0x0(%rax,%rax,1)
cltq
lea    (%rdi,%rax,1),%rcx
mov    %r15,0x0(%r13,%rcx,8)
add    $0x1,%rdi
add    %edx,%eax
add    $0xffffffffffffffff,%r10
jne    1300 <cov2cor+0x1f0>
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @cov2cor(double* nocapture %0, double* nocapture writeonly %1, i32 %2) {
%4 = icmp sgt i32 %2, 1
br i1 %4, label %5, label %16
5: ; preds = %3
%6 = zext i32 %2 to i64
%7 = zext i32 %2 to i64
br label %8
8: ; preds = %5, %46
%9 = phi i64 [ 1, %5 ], [ %47, %46 ]
%10 = trunc i64 %9 to i32
%11 = mul nsw i64 %9, %6
%12 = mul nsw i32 %10, %2
%13 = sext i32 %12 to i64
%14 = add nsw i64 %9, %13
%15 = getelementptr inbounds double, double* %0, i64 %14
br label %25
16: ; preds = %46, %3
%17 = icmp sgt i32 %2, 0
br i1 %17, label %18, label %92
18: ; preds = %16
%19 = zext i32 %2 to i64
%20 = add nsw i64 %19, -1
%21 = and i64 %19, 3
%22 = icmp ult i64 %20, 3
br i1 %22, label %78, label %23
23: ; preds = %18
%24 = and i64 %19, 4294967292
br label %49
25: ; preds = %8, %25
%26 = phi i64 [ 0, %8 ], [ %44, %25 ]
%27 = trunc i64 %26 to i32
%28 = add nuw nsw i64 %26, %11
%29 = getelementptr inbounds double, double* %0, i64 %28
%30 = load double, double* %29
%31 = load double, double* %15
%32 = mul nsw i64 %26, %6
%33 = mul nsw i32 %27, %2
%34 = sext i32 %33 to i64
%35 = add nsw i64 %26, %34
%36 = getelementptr inbounds double, double* %0, i64 %35
%37 = load double, double* %36
%38 = fmul double %31, %37
%39 = tail call double @sqrt(double %38)
%40 = fdiv double %30, %39
%41 = add nuw nsw i64 %32, %9
%42 = getelementptr inbounds double, double* %1, i64 %41
store double %40, double* %42
%43 = getelementptr inbounds double, double* %1, i64 %28
store double %40, double* %43
%44 = add nuw nsw i64 %26, 1
%45 = icmp eq i64 %44, %9
br i1 %45, label %46, label %25
46: ; preds = %25
%47 = add nuw nsw i64 %9, 1
%48 = icmp eq i64 %47, %7
br i1 %48, label %16, label %8
49: ; preds = %49, %23
%50 = phi i64 [ 0, %23 ], [ %75, %49 ]
%51 = phi i64 [ 0, %23 ], [ %76, %49 ]
%52 = trunc i64 %50 to i32
%53 = mul nsw i32 %52, %2
%54 = sext i32 %53 to i64
%55 = add nsw i64 %50, %54
%56 = getelementptr inbounds double, double* %1, i64 %55
store double 1.000000e+00, double* %56
%57 = or i64 %50, 1
%58 = trunc i64 %57 to i32
%59 = mul nsw i32 %58, %2
%60 = sext i32 %59 to i64
%61 = add nsw i64 %57, %60
%62 = getelementptr inbounds double, double* %1, i64 %61
store double 1.000000e+00, double* %62
%63 = or i64 %50, 2
%64 = trunc i64 %63 to i32
%65 = mul nsw i32 %64, %2
%66 = sext i32 %65 to i64
%67 = add nsw i64 %63, %66
%68 = getelementptr inbounds double, double* %1, i64 %67
store double 1.000000e+00, double* %68
%69 = or i64 %50, 3
%70 = trunc i64 %69 to i32
%71 = mul nsw i32 %70, %2
%72 = sext i32 %71 to i64
%73 = add nsw i64 %69, %72
%74 = getelementptr inbounds double, double* %1, i64 %73
store double 1.000000e+00, double* %74
%75 = add nuw nsw i64 %50, 4
%76 = add i64 %51, 4
%77 = icmp eq i64 %76, %24
br i1 %77, label %78, label %49
78: ; preds = %49, %18
%79 = phi i64 [ 0, %18 ], [ %75, %49 ]
%80 = icmp eq i64 %21, 0
br i1 %80, label %92, label %81
81: ; preds = %78, %81
%82 = phi i64 [ %89, %81 ], [ %79, %78 ]
%83 = phi i64 [ %90, %81 ], [ 0, %78 ]
%84 = trunc i64 %82 to i32
%85 = mul nsw i32 %84, %2
%86 = sext i32 %85 to i64
%87 = add nsw i64 %82, %86
%88 = getelementptr inbounds double, double* %1, i64 %87
store double 1.000000e+00, double* %88
%89 = add nuw nsw i64 %82, 1
%90 = add i64 %83, 1
%91 = icmp eq i64 %90, %21
br i1 %91, label %92, label %81
92: ; preds = %78, %81, %16
ret void
}