void safeCPdir(int numCP, double *cpDir, double *corrSrfNrmls) {
  int i;
  for (i = 3; i <= 3 * (numCP - 2); i += 3) {
    if ((cpDir[i] * corrSrfNrmls[i] + cpDir[i + 1] * corrSrfNrmls[i + 1] +
         cpDir[i + 2] * corrSrfNrmls[i + 2]) < 0.96) {
      cpDir[i] = corrSrfNrmls[i];
      cpDir[i + 1] = corrSrfNrmls[i + 1];
      cpDir[i + 2] = corrSrfNrmls[i + 2];
    }
  }
  return;
}


x86 assembly (clang, -O0):
<safeCPdir>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x3,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x2,%ecx
imul   $0x3,%ecx,%ecx
cmp    %ecx,%eax
jg     1223 <safeCPdir+0x123>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  0xe5c(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    1210 <safeCPdir+0x110>
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1215 <safeCPdir+0x115>
mov    -0x1c(%rbp),%eax
add    $0x3,%eax
mov    %eax,-0x1c(%rbp)
jmp    1116 <safeCPdir+0x16>
pop    %rbp
ret

x86 assembly (clang, -O2):
<safeCPdir>:
lea    (%rdi,%rdi,2),%eax
add    $0xfffffffa,%eax
cmp    $0x3,%eax
jge    110c <safeCPdir+0xc>
ret
mov    %eax,%eax
mov    $0x3,%ecx
movsd  0xee5(%rip),%xmm0
jmp    1129 <safeCPdir+0x29>
nopl   (%rax)
add    $0x3,%rcx
cmp    %rax,%rcx
ja     110b <safeCPdir+0xb>
movsd  (%rdx,%rcx,8),%xmm1
movsd  0x8(%rsi,%rcx,8),%xmm2
movsd  0x10(%rsi,%rcx,8),%xmm3
mulsd  0x8(%rdx,%rcx,8),%xmm2
movsd  (%rsi,%rcx,8),%xmm4
mulsd  %xmm1,%xmm4
mulsd  0x10(%rdx,%rcx,8),%xmm3
addsd  %xmm2,%xmm4
addsd  %xmm4,%xmm3
ucomisd %xmm3,%xmm0
jbe    1120 <safeCPdir+0x20>
movsd  %xmm1,(%rsi,%rcx,8)
movsd  0x8(%rdx,%rcx,8),%xmm1
movsd  %xmm1,0x8(%rsi,%rcx,8)
movsd  0x10(%rdx,%rcx,8),%xmm1
movsd  %xmm1,0x10(%rsi,%rcx,8)
jmp    1120 <safeCPdir+0x20>

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

LLVM IR (-O2):
define void @safeCPdir(i32 %0, double* nocapture %1, double* nocapture %2) {
%4 = mul i32 %0, 3
%5 = add i32 %4, -6
%6 = icmp slt i32 %5, 3
br i1 %6, label %35, label %7
7: ; preds = %3
%8 = zext i32 %5 to i64
br label %9
9: ; preds = %7, %32
%10 = phi i64 [ 3, %7 ], [ %33, %32 ]
%11 = getelementptr inbounds double, double* %1, i64 %10
%12 = load double, double* %11
%13 = getelementptr inbounds double, double* %2, i64 %10
%14 = load double, double* %13
%15 = add nuw nsw i64 %10, 1
%16 = getelementptr inbounds double, double* %1, i64 %15
%17 = load double, double* %16
%18 = getelementptr inbounds double, double* %2, i64 %15
%19 = load double, double* %18
%20 = fmul double %17, %19
%21 = tail call double @llvm.fmuladd.f64(double %12, double %14, double %20)
%22 = add nuw nsw i64 %10, 2
%23 = getelementptr inbounds double, double* %1, i64 %22
%24 = load double, double* %23
%25 = getelementptr inbounds double, double* %2, i64 %22
%26 = load double, double* %25
%27 = tail call double @llvm.fmuladd.f64(double %24, double %26, double %21)
%28 = fcmp olt double %27, 0x3FEEB851EB851EB8
br i1 %28, label %29, label %32
29: ; preds = %9
store double %14, double* %11
%30 = load double, double* %18
store double %30, double* %16
%31 = load double, double* %25
store double %31, double* %23
br label %32
32: ; preds = %9, %29
%33 = add nuw nsw i64 %10, 3
%34 = icmp ugt i64 %33, %8
br i1 %34, label %35, label %9
35: ; preds = %32, %3
ret void
}