void calculate_cluster_sizes(const double *Z, double *CS, int n) {
  int i, j, k;
  const double *row;
  for (k = 0; k < n - 1; k++) {
    row = Z + (k * 3);
    i = (int)row[0];
    j = (int)row[1];
    if (i >= n) {
      CS[k] = CS[i - n];
    } else {
      CS[k] = 1.0;
    }
    if (j >= n) {
      CS[k] = CS[k] + CS[j - n];
    } else {
      CS[k] = CS[k] + 1.0;
    }
  }
}


x86 assembly (clang, -O0):
<calculate_cluster_sizes>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    120b <calculate_cluster_sizes+0x10b>
mov    -0x8(%rbp),%rax
imul   $0x3,-0x20(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cvttsd2si (%rax),%eax
mov    %eax,-0x18(%rbp)
mov    -0x28(%rbp),%rax
cvttsd2si 0x8(%rax),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jl     1184 <calculate_cluster_sizes+0x84>
mov    -0x10(%rbp),%rax
mov    -0x18(%rbp),%ecx
sub    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1199 <calculate_cluster_sizes+0x99>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  0xe6c(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jl     11d6 <calculate_cluster_sizes+0xd6>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    -0x14(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    11f8 <calculate_cluster_sizes+0xf8>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  0xe1a(%rip),%xmm0
addsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    11fd <calculate_cluster_sizes+0xfd>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1116 <calculate_cluster_sizes+0x16>
pop    %rbp
ret

x86 assembly (clang, -O2):
<calculate_cluster_sizes>:
cmp    $0x2,%edx
jl     1164 <calculate_cluster_sizes+0x64>
lea    -0x1(%rdx),%r8d
add    $0x8,%rdi
xor    %ecx,%ecx
movsd  0xee9(%rip),%xmm0
jmp    1136 <calculate_cluster_sizes+0x36>
nopl   0x0(%rax)
addsd  %xmm2,%xmm1
movsd  %xmm1,(%rsi,%rcx,8)
add    $0x1,%rcx
add    $0x18,%rdi
cmp    %rcx,%r8
je     1164 <calculate_cluster_sizes+0x64>
cvttsd2si -0x8(%rdi),%eax
movapd %xmm0,%xmm1
sub    %edx,%eax
jl     114a <calculate_cluster_sizes+0x4a>
cltq
movsd  (%rsi,%rax,8),%xmm1
cvttsd2si (%rdi),%eax
movsd  %xmm1,(%rsi,%rcx,8)
movapd %xmm0,%xmm2
sub    %edx,%eax
jl     1120 <calculate_cluster_sizes+0x20>
cltq
movsd  (%rsi,%rax,8),%xmm2
jmp    1120 <calculate_cluster_sizes+0x20>
ret

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

LLVM IR (-O2):
define void @calculate_cluster_sizes(double* nocapture %0, double* nocapture %1, i32 %2) {
%4 = icmp sgt i32 %2, 1
br i1 %4, label %5, label %37
5: ; preds = %3
%6 = add nsw i32 %2, -1
%7 = zext i32 %6 to i64
br label %8
8: ; preds = %5, %32
%9 = phi i64 [ 0, %5 ], [ %35, %32 ]
%10 = mul nuw nsw i64 %9, 3
%11 = getelementptr inbounds double, double* %0, i64 %10
%12 = load double, double* %11
%13 = fptosi double %12 to i32
%14 = getelementptr inbounds double, double* %11, i64 1
%15 = load double, double* %14
%16 = fptosi double %15 to i32
%17 = icmp slt i32 %13, %2
br i1 %17, label %23, label %18
18: ; preds = %8
%19 = sub nsw i32 %13, %2
%20 = sext i32 %19 to i64
%21 = getelementptr inbounds double, double* %1, i64 %20
%22 = load double, double* %21
br label %23
23: ; preds = %8, %18
%24 = phi double [ %22, %18 ], [ 1.000000e+00, %8 ]
%25 = getelementptr inbounds double, double* %1, i64 %9
store double %24, double* %25
%26 = icmp slt i32 %16, %2
br i1 %26, label %32, label %27
27: ; preds = %23
%28 = sub nsw i32 %16, %2
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds double, double* %1, i64 %29
%31 = load double, double* %30
br label %32
32: ; preds = %23, %27
%33 = phi double [ %31, %27 ], [ 1.000000e+00, %23 ]
%34 = fadd double %24, %33
store double %34, double* %25
%35 = add nuw nsw i64 %9, 1
%36 = icmp eq i64 %35, %7
br i1 %36, label %37, label %8
37: ; preds = %32, %3
ret void
}