float xCorrCoef(float *target, float *regressor, int subl) {
  int i;
  float ftmp1;
  float ftmp2;
  ftmp1 = 0.0f;
  ftmp2 = 0.0f;
  for (i = 0; i < subl; i++) {
    ftmp1 += target[i] * regressor[i];
    ftmp2 += regressor[i] * regressor[i];
  }
  if (ftmp1 > 0.0f)
    return (float)(ftmp1 * ftmp1 / ftmp2);
  return 0.0f;
}


x86 assembly (clang, -O0):
<xCorrCoef>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x24(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x28(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1198 <xCorrCoef+0x98>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x24(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x24(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x28(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x28(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1126 <xCorrCoef+0x26>
movss  -0x24(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jbe    11c2 <xCorrCoef+0xc2>
movss  -0x24(%rbp),%xmm0
mulss  -0x24(%rbp),%xmm0
divss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    11ca <xCorrCoef+0xca>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<xCorrCoef>:
test   %edx,%edx
jle    1112 <xCorrCoef+0x12>
mov    %edx,%eax
cmp    $0x1,%edx
jne    1116 <xCorrCoef+0x16>
xorps  %xmm1,%xmm1
xor    %ecx,%ecx
jmp    1162 <xCorrCoef+0x62>
xorps  %xmm0,%xmm0
ret
mov    %eax,%edx
and    $0xfffffffe,%edx
xorps  %xmm0,%xmm0
xor    %ecx,%ecx
movss  (%rdi,%rcx,4),%xmm2
movss  0x4(%rdi,%rcx,4),%xmm3
movss  (%rsi,%rcx,4),%xmm4
movss  0x4(%rsi,%rcx,4),%xmm1
movaps %xmm4,%xmm5
unpcklps %xmm2,%xmm5
shufps $0x0,%xmm4,%xmm4
mulps  %xmm5,%xmm4
addps  %xmm0,%xmm4
movaps %xmm1,%xmm0
unpcklps %xmm3,%xmm0
shufps $0x0,%xmm1,%xmm1
mulps  %xmm0,%xmm1
addps  %xmm4,%xmm1
add    $0x2,%rcx
movaps %xmm1,%xmm0
cmp    %rcx,%rdx
jne    1120 <xCorrCoef+0x20>
test   $0x1,%al
je     1180 <xCorrCoef+0x80>
movss  (%rdi,%rcx,4),%xmm0
movss  (%rsi,%rcx,4),%xmm2
movaps %xmm2,%xmm3
unpcklps %xmm0,%xmm3
shufps $0x0,%xmm2,%xmm2
mulps  %xmm3,%xmm2
addps  %xmm2,%xmm1
movaps %xmm1,%xmm2
shufps $0x55,%xmm1,%xmm2
xorps  %xmm0,%xmm0
ucomiss %xmm0,%xmm2
jbe    119a <xCorrCoef+0x9a>
mulss  %xmm2,%xmm2
divss  %xmm1,%xmm2
movaps %xmm2,%xmm0
ret

LLVM IR (-O0):
define float @xCorrCoef(float* %0, float* %1, i32 %2) {
%4 = alloca float
%5 = alloca float*
%6 = alloca float*
%7 = alloca i32
%8 = alloca i32
%9 = alloca float
%10 = alloca float
store float* %0, float** %5
store float* %1, float** %6
store i32 %2, i32* %7
store float 0.000000e+00, float* %9
store float 0.000000e+00, float* %10
store i32 0, i32* %8
br label %11
11: ; preds = %40, %3
%12 = load i32, i32* %8
%13 = load i32, i32* %7
%14 = icmp slt i32 %12, %13
br i1 %14, label %15, label %43
15: ; preds = %11
%16 = load float*, float** %5
%17 = load i32, i32* %8
%18 = sext i32 %17 to i64
%19 = getelementptr inbounds float, float* %16, i64 %18
%20 = load float, float* %19
%21 = load float*, float** %6
%22 = load i32, i32* %8
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds float, float* %21, i64 %23
%25 = load float, float* %24
%26 = load float, float* %9
%27 = call float @llvm.fmuladd.f32(float %20, float %25, float %26)
store float %27, float* %9
%28 = load float*, float** %6
%29 = load i32, i32* %8
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds float, float* %28, i64 %30
%32 = load float, float* %31
%33 = load float*, float** %6
%34 = load i32, i32* %8
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds float, float* %33, i64 %35
%37 = load float, float* %36
%38 = load float, float* %10
%39 = call float @llvm.fmuladd.f32(float %32, float %37, float %38)
store float %39, float* %10
br label %40
40: ; preds = %15
%41 = load i32, i32* %8
%42 = add nsw i32 %41, 1
store i32 %42, i32* %8
br label %11
43: ; preds = %11
%44 = load float, float* %9
%45 = fcmp ogt float %44, 0.000000e+00
br i1 %45, label %46, label %52
46: ; preds = %43
%47 = load float, float* %9
%48 = load float, float* %9
%49 = fmul float %47, %48
%50 = load float, float* %10
%51 = fdiv float %49, %50
store float %51, float* %4
br label %53
52: ; preds = %43
store float 0.000000e+00, float* %4
br label %53
53: ; preds = %52, %46
%54 = load float, float* %4
ret float %54
}

LLVM IR (-O2):
define float @xCorrCoef(float* nocapture %0, float* nocapture %1, i32 %2) {
%4 = icmp sgt i32 %2, 0
br i1 %4, label %5, label %57
5: ; preds = %3
%6 = zext i32 %2 to i64
%7 = and i64 %6, 1
%8 = icmp eq i32 %2, 1
br i1 %8, label %35, label %9
9: ; preds = %5
%10 = and i64 %6, 4294967294
br label %11
11: ; preds = %11, %9
%12 = phi i64 [ 0, %9 ], [ %32, %11 ]
%13 = phi <2 x float> [ zeroinitializer, %9 ], [ %31, %11 ]
%14 = phi i64 [ 0, %9 ], [ %33, %11 ]
%15 = getelementptr inbounds float, float* %0, i64 %12
%16 = load float, float* %15
%17 = getelementptr inbounds float, float* %1, i64 %12
%18 = load float, float* %17
%19 = insertelement <2 x float> poison, float %18, i64 0
%20 = insertelement <2 x float> %19, float %16, i64 1
%21 = shufflevector <2 x float> %19, <2 x float> poison, <2 x i32> zeroinitializer
%22 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %20, <2 x float> %21, <2 x float> %13)
%23 = or i64 %12, 1
%24 = getelementptr inbounds float, float* %0, i64 %23
%25 = load float, float* %24
%26 = getelementptr inbounds float, float* %1, i64 %23
%27 = load float, float* %26
%28 = insertelement <2 x float> poison, float %27, i64 0
%29 = insertelement <2 x float> %28, float %25, i64 1
%30 = shufflevector <2 x float> %28, <2 x float> poison, <2 x i32> zeroinitializer
%31 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %29, <2 x float> %30, <2 x float> %22)
%32 = add nuw nsw i64 %12, 2
%33 = add i64 %14, 2
%34 = icmp eq i64 %33, %10
br i1 %34, label %35, label %11
35: ; preds = %11, %5
%36 = phi <2 x float> [ undef, %5 ], [ %31, %11 ]
%37 = phi i64 [ 0, %5 ], [ %32, %11 ]
%38 = phi <2 x float> [ zeroinitializer, %5 ], [ %31, %11 ]
%39 = icmp eq i64 %7, 0
br i1 %39, label %49, label %40
40: ; preds = %35
%41 = getelementptr inbounds float, float* %0, i64 %37
%42 = load float, float* %41
%43 = getelementptr inbounds float, float* %1, i64 %37
%44 = load float, float* %43
%45 = insertelement <2 x float> poison, float %44, i64 0
%46 = insertelement <2 x float> %45, float %42, i64 1
%47 = shufflevector <2 x float> %45, <2 x float> poison, <2 x i32> zeroinitializer
%48 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %46, <2 x float> %47, <2 x float> %38)
br label %49
49: ; preds = %35, %40
%50 = phi <2 x float> [ %36, %35 ], [ %48, %40 ]
%51 = extractelement <2 x float> %50, i64 1
%52 = fcmp ogt float %51, 0.000000e+00
br i1 %52, label %53, label %57
53: ; preds = %49
%54 = fmul float %51, %51
%55 = extractelement <2 x float> %50, i64 0
%56 = fdiv float %54, %55
br label %57
57: ; preds = %3, %49, %53
%58 = phi float [ %56, %53 ], [ 0.000000e+00, %49 ], [ 0.000000e+00, %3 ]
ret float %58
}