cor(s1, s2, span, p, it) float *s1, *s2, *p;
int span, it;
{
  int i;
  int i2 = 0;
  int i3 = 0;
  float fac = 1.0;
  float sum1 = 1.0e-15;
  float sum2 = 1.0e-15;
  float sum3 = 1.0e-15;
  float smin = 1.0e-10;
  for (i = 0; i < span; i = i + it) {
    sum1 = sum1 + s1[i] * s2[i];
    sum2 = sum2 + s1[i] * s1[i];
    sum3 = sum3 + s2[i] * s2[i];
  }
  if (sum2 == 0.0 || sum3 == 0.0)
    *p = -1.0;
  else {
    while (sum2 <= smin) {
      sum2 = sum2 / smin;
      i2++;
    }
    while (sum3 <= smin) {
      sum3 = sum3 / smin;
      i3++;
    }
    if (i2 || i3)
      fac = 1.0 / sqrt((i2 + i3) * smin);
    *p = (fac * sum1) / sqrt(sum2 * sum3);
  }
}


x86 assembly (clang, -O0):
<cor>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8d,-0x2c(%rbp)
movl   $0x0,-0x34(%rbp)
movl   $0x0,-0x38(%rbp)
movss  0xec7(%rip),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  0xeb6(%rip),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  0xea9(%rip),%xmm0
movss  %xmm0,-0x44(%rbp)
movss  0xe9c(%rip),%xmm0
movss  %xmm0,-0x48(%rbp)
movss  0xe8b(%rip),%xmm0
movss  %xmm0,-0x4c(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    121f <cor+0x10f>
movss  -0x40(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x44(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x48(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x30(%rbp),%eax
add    -0x2c(%rbp),%eax
mov    %eax,-0x30(%rbp)
jmp    1181 <cor+0x71>
movss  -0x44(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1240 <cor+0x130>
jp     1240 <cor+0x130>
jmp    125c <cor+0x14c>
movss  -0x48(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1271 <cor+0x161>
jp     1271 <cor+0x161>
mov    -0x28(%rbp),%rax
movss  0xda4(%rip),%xmm0
movss  %xmm0,(%rax)
jmp    135d <cor+0x24d>
jmp    1276 <cor+0x166>
movss  -0x44(%rbp),%xmm1
movss  -0x4c(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jb     12a6 <cor+0x196>
movss  -0x44(%rbp),%xmm0
divss  -0x4c(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1276 <cor+0x166>
jmp    12ab <cor+0x19b>
movss  -0x48(%rbp),%xmm1
movss  -0x4c(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jb     12db <cor+0x1cb>
movss  -0x48(%rbp),%xmm0
divss  -0x4c(%rbp),%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    12ab <cor+0x19b>
cmpl   $0x0,-0x34(%rbp)
jne    12ef <cor+0x1df>
cmpl   $0x0,-0x38(%rbp)
je     131f <cor+0x20f>
mov    -0x34(%rbp),%eax
add    -0x38(%rbp),%eax
cvtsi2ss %eax,%xmm0
mulss  -0x4c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
call   1030 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  0xcfe(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x3c(%rbp),%xmm0
mulss  -0x40(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
movss  -0x44(%rbp),%xmm0
mulss  -0x48(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
call   1030 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x58(%rbp),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x28(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x4(%rbp),%eax
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cor>:
push   %rbx
sub    $0x20,%rsp
mov    %rcx,%rbx
test   %edx,%edx
jle    11ab <cor+0x9b>
movslq %r8d,%rax
mov    %edx,%ecx
movaps 0xed4(%rip),%xmm2
movss  0xedc(%rip),%xmm4
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
movaps %xmm2,%xmm1
movss  (%rdi,%rdx,4),%xmm0
movss  (%rsi,%rdx,4),%xmm2
movaps %xmm0,%xmm3
mulss  %xmm2,%xmm3
addss  %xmm3,%xmm4
unpcklps %xmm2,%xmm0
mulps  %xmm0,%xmm0
movaps %xmm0,%xmm2
addps  %xmm1,%xmm2
add    %rax,%rdx
cmp    %rcx,%rdx
jl     1140 <cor+0x30>
addps  %xmm1,%xmm0
xorps  %xmm2,%xmm2
ucomiss %xmm2,%xmm0
movss  0xe97(%rip),%xmm1
jne    1185 <cor+0x75>
jnp    12a8 <cor+0x198>
movaps %xmm0,%xmm3
shufps $0x55,%xmm0,%xmm3
ucomiss %xmm2,%xmm3
jne    1197 <cor+0x87>
jnp    12a8 <cor+0x198>
movss  0xe79(%rip),%xmm1
ucomiss %xmm0,%xmm1
jae    11ba <cor+0xaa>
xor    %eax,%eax
movaps %xmm0,%xmm2
jmp    11e3 <cor+0xd3>
movaps 0xe4e(%rip),%xmm0
movss  0xe56(%rip),%xmm4
xor    %eax,%eax
movss  0xe54(%rip),%xmm1
movaps %xmm0,%xmm2
nopw   0x0(%rax,%rax,1)
divss  %xmm1,%xmm2
add    $0x1,%eax
ucomiss %xmm2,%xmm1
jae    11d0 <cor+0xc0>
shufps $0x55,%xmm0,%xmm0
movaps %xmm0,%xmm3
movss  0xe2d(%rip),%xmm0
xor    %ecx,%ecx
ucomiss %xmm3,%xmm0
jae    1200 <cor+0xf0>
test   %eax,%eax
jne    120c <cor+0xfc>
movss  0xe1e(%rip),%xmm0
jmp    1268 <cor+0x158>
divss  %xmm0,%xmm3
add    $0x1,%ecx
ucomiss %xmm3,%xmm0
jae    1200 <cor+0xf0>
add    %eax,%ecx
xorps  %xmm0,%xmm0
cvtsi2ss %ecx,%xmm0
mulss  0xdfb(%rip),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     1230 <cor+0x120>
sqrtsd %xmm0,%xmm0
jmp    1255 <cor+0x145>
movss  %xmm4,(%rsp)
movss  %xmm2,0xc(%rsp)
movaps %xmm3,0x10(%rsp)
call   1030 <sqrt@plt>
movaps 0x10(%rsp),%xmm3
movss  0xc(%rsp),%xmm2
movss  (%rsp),%xmm4
movsd  0xdc3(%rip),%xmm1
divsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
mulss  %xmm3,%xmm2
mulss  %xmm0,%xmm4
xorps  %xmm3,%xmm3
cvtss2sd %xmm4,%xmm3
xorps  %xmm0,%xmm0
cvtss2sd %xmm2,%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     128e <cor+0x17e>
sqrtsd %xmm0,%xmm0
jmp    129d <cor+0x18d>
movsd  %xmm3,(%rsp)
call   1030 <sqrt@plt>
movsd  (%rsp),%xmm3
divsd  %xmm0,%xmm3
xorps  %xmm1,%xmm1
cvtsd2ss %xmm3,%xmm1
movss  %xmm1,(%rbx)
add    $0x20,%rsp
pop    %rbx
ret

LLVM IR (-O0):
define i32 @cor(float* %0, float* %1, i32 %2, float* %3, i32 %4) {
%6 = alloca i32
%7 = alloca float*
%8 = alloca float*
%9 = alloca i32
%10 = alloca float*
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca float
%16 = alloca float
%17 = alloca float
%18 = alloca float
%19 = alloca float
store float* %0, float** %7
store float* %1, float** %8
store i32 %2, i32* %9
store float* %3, float** %10
store i32 %4, i32* %11
store i32 0, i32* %13
store i32 0, i32* %14
store float 1.000000e+00, float* %15
store float 0x3CD203AFA0000000, float* %16
store float 0x3CD203AFA0000000, float* %17
store float 0x3CD203AFA0000000, float* %18
store float 0x3DDB7CDFE0000000, float* %19
store i32 0, i32* %12
br label %20
20: ; preds = %61, %5
%21 = load i32, i32* %12
%22 = load i32, i32* %9
%23 = icmp slt i32 %21, %22
br i1 %23, label %24, label %65
24: ; preds = %20
%25 = load float, float* %16
%26 = load float*, float** %7
%27 = load i32, i32* %12
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds float, float* %26, i64 %28
%30 = load float, float* %29
%31 = load float*, float** %8
%32 = load i32, i32* %12
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds float, float* %31, i64 %33
%35 = load float, float* %34
%36 = call float @llvm.fmuladd.f32(float %30, float %35, float %25)
store float %36, float* %16
%37 = load float, float* %17
%38 = load float*, float** %7
%39 = load i32, i32* %12
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds float, float* %38, i64 %40
%42 = load float, float* %41
%43 = load float*, float** %7
%44 = load i32, i32* %12
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds float, float* %43, i64 %45
%47 = load float, float* %46
%48 = call float @llvm.fmuladd.f32(float %42, float %47, float %37)
store float %48, float* %17
%49 = load float, float* %18
%50 = load float*, float** %8
%51 = load i32, i32* %12
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds float, float* %50, i64 %52
%54 = load float, float* %53
%55 = load float*, float** %8
%56 = load i32, i32* %12
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds float, float* %55, i64 %57
%59 = load float, float* %58
%60 = call float @llvm.fmuladd.f32(float %54, float %59, float %49)
store float %60, float* %18
br label %61
61: ; preds = %24
%62 = load i32, i32* %12
%63 = load i32, i32* %11
%64 = add nsw i32 %62, %63
store i32 %64, i32* %12
br label %20
65: ; preds = %20
%66 = load float, float* %17
%67 = fpext float %66 to double
%68 = fcmp oeq double %67, 0.000000e+00
br i1 %68, label %73, label %69
69: ; preds = %65
%70 = load float, float* %18
%71 = fpext float %70 to double
%72 = fcmp oeq double %71, 0.000000e+00
br i1 %72, label %73, label %75
73: ; preds = %69, %65
%74 = load float*, float** %10
store float -1.000000e+00, float* %74
br label %127
75: ; preds = %69
br label %76
76: ; preds = %80, %75
%77 = load float, float* %17
%78 = load float, float* %19
%79 = fcmp ole float %77, %78
br i1 %79, label %80, label %86
80: ; preds = %76
%81 = load float, float* %17
%82 = load float, float* %19
%83 = fdiv float %81, %82
store float %83, float* %17
%84 = load i32, i32* %13
%85 = add nsw i32 %84, 1
store i32 %85, i32* %13
br label %76
86: ; preds = %76
br label %87
87: ; preds = %91, %86
%88 = load float, float* %18
%89 = load float, float* %19
%90 = fcmp ole float %88, %89
br i1 %90, label %91, label %97
91: ; preds = %87
%92 = load float, float* %18
%93 = load float, float* %19
%94 = fdiv float %92, %93
store float %94, float* %18
%95 = load i32, i32* %14
%96 = add nsw i32 %95, 1
store i32 %96, i32* %14
br label %87
97: ; preds = %87
%98 = load i32, i32* %13
%99 = icmp ne i32 %98, 0
br i1 %99, label %103, label %100
100: ; preds = %97
%101 = load i32, i32* %14
%102 = icmp ne i32 %101, 0
br i1 %102, label %103, label %114
103: ; preds = %100, %97
%104 = load i32, i32* %13
%105 = load i32, i32* %14
%106 = add nsw i32 %104, %105
%107 = sitofp i32 %106 to float
%108 = load float, float* %19
%109 = fmul float %107, %108
%110 = fpext float %109 to double
%111 = call double @sqrt(double %110)
%112 = fdiv double 1.000000e+00, %111
%113 = fptrunc double %112 to float
store float %113, float* %15
br label %114
114: ; preds = %103, %100
%115 = load float, float* %15
%116 = load float, float* %16
%117 = fmul float %115, %116
%118 = fpext float %117 to double
%119 = load float, float* %17
%120 = load float, float* %18
%121 = fmul float %119, %120
%122 = fpext float %121 to double
%123 = call double @sqrt(double %122)
%124 = fdiv double %118, %123
%125 = fptrunc double %124 to float
%126 = load float*, float** %10
store float %125, float* %126
br label %127
127: ; preds = %114, %73
%128 = load i32, i32* %6
ret i32 %128
}

LLVM IR (-O2):
define i32 @cor(float* nocapture %0, float* nocapture %1, i32 %2, float* nocapture writeonly %3, i32 %4) {
%6 = icmp sgt i32 %2, 0
br i1 %6, label %7, label %33
7: ; preds = %5
%8 = sext i32 %4 to i64
%9 = zext i32 %2 to i64
br label %10
10: ; preds = %7, %10
%11 = phi i64 [ 0, %7 ], [ %22, %10 ]
%12 = phi float [ 0x3CD203AFA0000000, %7 ], [ %18, %10 ]
%13 = phi <2 x float> [ <float 0x3CD203AFA0000000, float 0x3CD203AFA0000000>, %7 ], [ %21, %10 ]
%14 = getelementptr inbounds float, float* %0, i64 %11
%15 = load float, float* %14
%16 = getelementptr inbounds float, float* %1, i64 %11
%17 = load float, float* %16
%18 = tail call float @llvm.fmuladd.f32(float %15, float %17, float %12)
%19 = insertelement <2 x float> poison, float %15, i64 0
%20 = insertelement <2 x float> %19, float %17, i64 1
%21 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %20, <2 x float> %20, <2 x float> %13)
%22 = add i64 %11, %8
%23 = icmp slt i64 %22, %9
br i1 %23, label %10, label %24
24: ; preds = %10
%25 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %20, <2 x float> %20, <2 x float> %13)
%26 = extractelement <2 x float> %25, i64 0
%27 = fcmp oeq float %26, 0.000000e+00
%28 = extractelement <2 x float> %25, i64 1
%29 = fcmp oeq float %28, 0.000000e+00
%30 = select i1 %27, i1 true, i1 %29
br i1 %30, label %79, label %31
31: ; preds = %24
%32 = fcmp ugt float %26, 0x3DDB7CDFE0000000
br i1 %32, label %39, label %33
33: ; preds = %5, %31
%34 = phi float [ %18, %31 ], [ 0x3CD203AFA0000000, %5 ]
%35 = phi <2 x float> [ %25, %31 ], [ <float 0x3CD203AFA0000000, float 0x3CD203AFA0000000>, %5 ]
%36 = extractelement <2 x float> %35, i64 0
br label %45
37: ; preds = %45
%38 = extractelement <2 x float> %35, i64 1
br label %39
39: ; preds = %37, %31
%40 = phi float [ %18, %31 ], [ %34, %37 ]
%41 = phi float [ %28, %31 ], [ %38, %37 ]
%42 = phi i32 [ 0, %31 ], [ %49, %37 ]
%43 = phi float [ %26, %31 ], [ %48, %37 ]
%44 = fcmp ugt float %41, 0x3DDB7CDFE0000000
br i1 %44, label %57, label %51
45: ; preds = %33, %45
%46 = phi float [ %48, %45 ], [ %36, %33 ]
%47 = phi i32 [ %49, %45 ], [ 0, %33 ]
%48 = fdiv float %46, 0x3DDB7CDFE0000000
%49 = add nuw nsw i32 %47, 1
%50 = fcmp ugt float %48, 0x3DDB7CDFE0000000
br i1 %50, label %37, label %45
51: ; preds = %39, %51
%52 = phi float [ %54, %51 ], [ %41, %39 ]
%53 = phi i32 [ %55, %51 ], [ 0, %39 ]
%54 = fdiv float %52, 0x3DDB7CDFE0000000
%55 = add nuw nsw i32 %53, 1
%56 = fcmp ugt float %54, 0x3DDB7CDFE0000000
br i1 %56, label %59, label %51
57: ; preds = %39
%58 = icmp eq i32 %42, 0
br i1 %58, label %69, label %59
59: ; preds = %51, %57
%60 = phi float [ %41, %57 ], [ %54, %51 ]
%61 = phi i32 [ 0, %57 ], [ %55, %51 ]
%62 = add nuw nsw i32 %61, %42
%63 = sitofp i32 %62 to float
%64 = fmul float %63, 0x3DDB7CDFE0000000
%65 = fpext float %64 to double
%66 = tail call double @sqrt(double %65)
%67 = fdiv double 1.000000e+00, %66
%68 = fptrunc double %67 to float
br label %69
69: ; preds = %57, %59
%70 = phi float [ %60, %59 ], [ %41, %57 ]
%71 = phi float [ %68, %59 ], [ 1.000000e+00, %57 ]
%72 = fmul float %40, %71
%73 = fpext float %72 to double
%74 = fmul float %43, %70
%75 = fpext float %74 to double
%76 = tail call double @sqrt(double %75)
%77 = fdiv double %73, %76
%78 = fptrunc double %77 to float
br label %79
79: ; preds = %24, %69
%80 = phi float [ %78, %69 ], [ -1.000000e+00, %24 ]
store float %80, float* %3
ret i32 undef
}