void getUnweightedCP(int ncp, int kcp, double *cpw, double *cp) {
  int i;
  for (i = 0; i < kcp * ncp; i++) {
    cp[3 * i] = cpw[4 * i] / cpw[4 * i + 3];
    cp[3 * i + 1] = cpw[4 * i + 1] / cpw[4 * i + 3];
    cp[3 * i + 2] = cpw[4 * i + 2] / cpw[4 * i + 3];
  }
}


x86 assembly (clang, -O0):
<getUnweightedCP>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x8(%rbp),%ecx
imul   -0x4(%rbp),%ecx
cmp    %ecx,%eax
jge    11ea <getUnweightedCP+0xea>
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
shl    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
shl    $0x2,%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
imul   $0x3,-0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
shl    $0x2,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
shl    $0x2,%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
imul   $0x3,-0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
shl    $0x2,%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
shl    $0x2,%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
imul   $0x3,-0x1c(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1119 <getUnweightedCP+0x19>
pop    %rbp
ret

x86 assembly (clang, -O2):
<getUnweightedCP>:
imul   %esi,%edi
test   %edi,%edi
jle    1227 <getUnweightedCP+0x127>
mov    %edi,%eax
cmp    $0x3,%edi
jb     112e <getUnweightedCP+0x2e>
mov    %rax,%rsi
shl    $0x5,%rsi
add    %rdx,%rsi
cmp    %rcx,%rsi
jbe    1136 <getUnweightedCP+0x36>
lea    (%rax,%rax,2),%rsi
lea    (%rcx,%rsi,8),%rsi
cmp    %rdx,%rsi
jbe    1136 <getUnweightedCP+0x36>
xor    %r8d,%r8d
jmp    11cd <getUnweightedCP+0xcd>
mov    %eax,%esi
and    $0x1,%esi
test   %rsi,%rsi
mov    $0x2,%edi
cmovne %rsi,%rdi
mov    %rax,%r8
sub    %rdi,%r8
mov    %r8,%r9
mov    %rcx,%rsi
mov    %rdx,%rdi
cs nopw 0x0(%rax,%rax,1)
movupd 0x30(%rdi),%xmm0
movupd (%rdi),%xmm1
movupd 0x10(%rdi),%xmm2
movupd 0x18(%rdi),%xmm3
movupd 0x20(%rdi),%xmm4
movapd %xmm1,%xmm5
unpckhpd %xmm4,%xmm5
movapd %xmm2,%xmm6
unpcklpd %xmm0,%xmm6
unpckhpd %xmm0,%xmm2
divpd  %xmm2,%xmm6
movhpd 0x38(%rdi),%xmm3
unpcklpd %xmm4,%xmm1
divpd  %xmm3,%xmm1
divpd  %xmm3,%xmm5
movapd %xmm1,%xmm0
unpcklpd %xmm5,%xmm0
unpckhpd %xmm6,%xmm5
movsd  %xmm6,%xmm1
movupd %xmm1,0x10(%rsi)
movupd %xmm5,0x20(%rsi)
movupd %xmm0,(%rsi)
add    $0x40,%rdi
add    $0x30,%rsi
add    $0xfffffffffffffffe,%r9
jne    1160 <getUnweightedCP+0x60>
lea    (%r8,%r8,2),%rsi
lea    (%rcx,%rsi,8),%rcx
add    $0x10,%rcx
sub    %r8,%rax
shl    $0x5,%r8
add    %r8,%rdx
add    $0x18,%rdx
nopw   0x0(%rax,%rax,1)
movsd  -0x18(%rdx),%xmm0
divsd  (%rdx),%xmm0
movsd  %xmm0,-0x10(%rcx)
movsd  -0x10(%rdx),%xmm0
divsd  (%rdx),%xmm0
movsd  %xmm0,-0x8(%rcx)
movsd  -0x8(%rdx),%xmm0
divsd  (%rdx),%xmm0
movsd  %xmm0,(%rcx)
add    $0x18,%rcx
add    $0x20,%rdx
add    $0xffffffffffffffff,%rax
jne    11f0 <getUnweightedCP+0xf0>
ret

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

LLVM IR (-O2):
define void @getUnweightedCP(i32 %0, i32 %1, double* nocapture %2, double* nocapture writeonly %3) {
%5 = mul i32 %1, %0
%6 = icmp sgt i32 %5, 0
br i1 %6, label %7, label %78
7: ; preds = %4
%8 = zext i32 %5 to i64
%9 = icmp ult i32 %5, 3
br i1 %9, label %49, label %10
10: ; preds = %7
%11 = mul nuw nsw i64 %8, 3
%12 = getelementptr double, double* %3, i64 %11
%13 = shl nuw nsw i64 %8, 2
%14 = getelementptr double, double* %2, i64 %13
%15 = icmp ugt double* %14, %3
%16 = icmp ugt double* %12, %2
%17 = and i1 %15, %16
br i1 %17, label %49, label %18
18: ; preds = %10
%19 = and i64 %8, 1
%20 = icmp eq i64 %19, 0
%21 = select i1 %20, i64 2, i64 %19
%22 = sub nsw i64 %8, %21
br label %23
23: ; preds = %23, %18
%24 = phi i64 [ 0, %18 ], [ %47, %23 ]
%25 = shl nsw i64 %24, 2
%26 = getelementptr inbounds double, double* %2, i64 %25
%27 = bitcast double* %26 to <8 x double>*
%28 = load <8 x double>, <8 x double>* %27
%29 = shufflevector <8 x double> %28, <8 x double> poison, <2 x i32> <i32 1, i32 5>
%30 = shufflevector <8 x double> %28, <8 x double> poison, <2 x i32> <i32 2, i32 6>
%31 = shufflevector <8 x double> %28, <8 x double> poison, <2 x i32> <i32 3, i32 7>
%32 = or i64 %25, 3
%33 = getelementptr inbounds double, double* %2, i64 %32
%34 = bitcast double* %33 to <8 x double>*
%35 = load <8 x double>, <8 x double>* %34
%36 = fdiv <8 x double> %28, %35
%37 = shufflevector <8 x double> %36, <8 x double> poison, <2 x i32> <i32 0, i32 4>
%38 = mul nuw nsw i64 %24, 3
%39 = shufflevector <8 x double> %35, <8 x double> poison, <2 x i32> <i32 0, i32 4>
%40 = fdiv <2 x double> %29, %39
%41 = fdiv <2 x double> %30, %31
%42 = getelementptr inbounds double, double* %3, i64 %38
%43 = bitcast double* %42 to <6 x double>*
%44 = shufflevector <2 x double> %37, <2 x double> %40, <4 x i32> <i32 0, i32 1, i32 2, i32 3>
%45 = shufflevector <2 x double> %41, <2 x double> poison, <4 x i32> <i32 0, i32 1, i32 undef, i32 undef>
%46 = shufflevector <4 x double> %44, <4 x double> %45, <6 x i32> <i32 0, i32 2, i32 4, i32 1, i32 3, i32 5>
store <6 x double> %46, <6 x double>* %43
%47 = add nuw i64 %24, 2
%48 = icmp eq i64 %47, %22
br i1 %48, label %49, label %23
49: ; preds = %23, %10, %7
%50 = phi i64 [ 0, %10 ], [ 0, %7 ], [ %22, %23 ]
br label %51
51: ; preds = %49, %51
%52 = phi i64 [ %76, %51 ], [ %50, %49 ]
%53 = shl nsw i64 %52, 2
%54 = getelementptr inbounds double, double* %2, i64 %53
%55 = load double, double* %54
%56 = or i64 %53, 3
%57 = getelementptr inbounds double, double* %2, i64 %56
%58 = load double, double* %57
%59 = fdiv double %55, %58
%60 = mul nuw nsw i64 %52, 3
%61 = getelementptr inbounds double, double* %3, i64 %60
store double %59, double* %61
%62 = or i64 %53, 1
%63 = getelementptr inbounds double, double* %2, i64 %62
%64 = load double, double* %63
%65 = load double, double* %57
%66 = fdiv double %64, %65
%67 = add nuw nsw i64 %60, 1
%68 = getelementptr inbounds double, double* %3, i64 %67
store double %66, double* %68
%69 = or i64 %53, 2
%70 = getelementptr inbounds double, double* %2, i64 %69
%71 = load double, double* %70
%72 = load double, double* %57
%73 = fdiv double %71, %72
%74 = add nuw nsw i64 %60, 2
%75 = getelementptr inbounds double, double* %3, i64 %74
store double %73, double* %75
%76 = add nuw nsw i64 %52, 1
%77 = icmp eq i64 %76, %8
br i1 %77, label %78, label %51
78: ; preds = %51, %4
ret void
}