void CPU_BCV(int *histogram, float *BCVtable, int size) {
  int i;
  for (i = 0; i < 256; i++) {
    BCVtable[i] = 0;
  }
  int wsumtotal = 0;
  for (i = 0; i < 256; i++) {
    wsumtotal = wsumtotal + i * histogram[i];
  }
  int sumb = 0;
  int sumf = size;
  int wsumb = 0;
  int wsumf = wsumtotal;
  float wb;
  float wf;
  float meanb;
  float meanf;
  for (i = 0; i < 256; i++) {
    wsumb = wsumb + i * histogram[i];
    wsumf = wsumtotal - wsumb;
    wb = sumb / (float)size;
    wf = sumf / (float)size;
    meanb = wsumb / (float)sumb;
    meanf = wsumf / (float)sumf;
    if (sumf == 0) {
      break;
    }
    BCVtable[i] = wb * wf * (meanb - meanf) * (meanb - meanf);
    if (sumb == 0) {
      BCVtable[i] = 0;
    }
    sumb = sumb + histogram[i];
    sumf = size - sumb;
  }
}


x86 assembly (clang, -O0):
<CPU_BCV>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movl   $0x0,-0x18(%rbp)
cmpl   $0x100,-0x18(%rbp)
jge    1141 <CPU_BCV+0x41>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1116 <CPU_BCV+0x16>
movl   $0x0,-0x1c(%rbp)
movl   $0x0,-0x18(%rbp)
cmpl   $0x100,-0x18(%rbp)
jge    1181 <CPU_BCV+0x81>
mov    -0x1c(%rbp),%eax
mov    -0x18(%rbp),%ecx
mov    -0x8(%rbp),%rdx
movslq -0x18(%rbp),%rsi
imul   (%rdx,%rsi,4),%ecx
add    %ecx,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    114f <CPU_BCV+0x4f>
movl   $0x0,-0x20(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x24(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x1c(%rbp),%eax
mov    %eax,-0x2c(%rbp)
movl   $0x0,-0x18(%rbp)
cmpl   $0x100,-0x18(%rbp)
jge    129f <CPU_BCV+0x19f>
mov    -0x28(%rbp),%eax
mov    -0x18(%rbp),%ecx
mov    -0x8(%rbp),%rdx
movslq -0x18(%rbp),%rsi
imul   (%rdx,%rsi,4),%ecx
add    %ecx,%eax
mov    %eax,-0x28(%rbp)
mov    -0x1c(%rbp),%eax
sub    -0x28(%rbp),%eax
mov    %eax,-0x2c(%rbp)
cvtsi2ssl -0x20(%rbp),%xmm0
cvtsi2ssl -0x14(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x30(%rbp)
cvtsi2ssl -0x24(%rbp),%xmm0
cvtsi2ssl -0x14(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x34(%rbp)
cvtsi2ssl -0x28(%rbp),%xmm0
cvtsi2ssl -0x20(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x38(%rbp)
cvtsi2ssl -0x2c(%rbp),%xmm0
cvtsi2ssl -0x24(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x3c(%rbp)
cmpl   $0x0,-0x24(%rbp)
jne    122a <CPU_BCV+0x12a>
jmp    129f <CPU_BCV+0x19f>
movss  -0x30(%rbp),%xmm0
mulss  -0x34(%rbp),%xmm0
movss  -0x38(%rbp),%xmm1
subss  -0x3c(%rbp),%xmm1
mulss  %xmm1,%xmm0
movss  -0x38(%rbp),%xmm1
subss  -0x3c(%rbp),%xmm1
mulss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
cmpl   $0x0,-0x20(%rbp)
jne    1277 <CPU_BCV+0x177>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%rbp),%eax
mov    -0x8(%rbp),%rcx
movslq -0x18(%rbp),%rdx
add    (%rcx,%rdx,4),%eax
mov    %eax,-0x20(%rbp)
mov    -0x14(%rbp),%eax
sub    -0x20(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    11a2 <CPU_BCV+0xa2>
pop    %rbp
ret

x86 assembly (clang, -O2):
<CPU_BCV>:
push   %rbp
push   %r15
push   %r14
push   %rbx
push   %rax
mov    %edx,%r15d
mov    %rsi,%r14
mov    %rdi,%rbx
xor    %ebp,%ebp
mov    $0x400,%edx
mov    %rsi,%rdi
xor    %esi,%esi
call   1030 <memset@plt>
movdqa 0xec7(%rip),%xmm4
pxor   %xmm5,%xmm5
movdqa 0xeca(%rip),%xmm8
movdqa 0xed1(%rip),%xmm9
movdqa 0xed8(%rip),%xmm10
movdqa 0xedf(%rip),%xmm11
pxor   %xmm6,%xmm6
cs nopw 0x0(%rax,%rax,1)
nop
movdqa %xmm4,%xmm0
paddd  %xmm8,%xmm0
movdqu (%rbx,%rbp,4),%xmm1
movdqu 0x10(%rbx,%rbp,4),%xmm2
movdqu 0x20(%rbx,%rbp,4),%xmm3
movdqu 0x30(%rbx,%rbp,4),%xmm13
pshufd $0xf5,%xmm1,%xmm12
pmuludq %xmm4,%xmm1
pshufd $0xe8,%xmm1,%xmm1
pshufd $0xf5,%xmm4,%xmm7
pmuludq %xmm12,%xmm7
pshufd $0xe8,%xmm7,%xmm7
punpckldq %xmm7,%xmm1
paddd  %xmm5,%xmm1
pshufd $0xf5,%xmm2,%xmm5
pmuludq %xmm0,%xmm2
pshufd $0xe8,%xmm2,%xmm2
pshufd $0xf5,%xmm0,%xmm0
pmuludq %xmm5,%xmm0
pshufd $0xe8,%xmm0,%xmm0
punpckldq %xmm0,%xmm2
paddd  %xmm6,%xmm2
movdqa %xmm4,%xmm0
paddd  %xmm9,%xmm0
movdqa %xmm4,%xmm7
paddd  %xmm10,%xmm7
pshufd $0xf5,%xmm3,%xmm6
pmuludq %xmm0,%xmm3
pshufd $0xe8,%xmm3,%xmm5
pshufd $0xf5,%xmm0,%xmm0
pmuludq %xmm6,%xmm0
pshufd $0xe8,%xmm0,%xmm0
punpckldq %xmm0,%xmm5
paddd  %xmm1,%xmm5
pshufd $0xf5,%xmm13,%xmm0
pmuludq %xmm7,%xmm13
pshufd $0xe8,%xmm13,%xmm6
pshufd $0xf5,%xmm7,%xmm1
pmuludq %xmm0,%xmm1
pshufd $0xe8,%xmm1,%xmm0
punpckldq %xmm0,%xmm6
paddd  %xmm2,%xmm6
add    $0x10,%rbp
paddd  %xmm11,%xmm4
cmp    $0x100,%rbp
jne    1170 <CPU_BCV+0x60>
paddd  %xmm5,%xmm6
pshufd $0xee,%xmm6,%xmm0
paddd  %xmm6,%xmm0
pshufd $0x55,%xmm0,%xmm1
paddd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsi2ss %r15d,%xmm0
movd   %xmm1,%r8d
xor    %ecx,%ecx
xor    %edx,%edx
mov    %r15d,%ebp
xor    %esi,%esi
jmp    129a <CPU_BCV+0x18a>
nopl   0x0(%rax)
movd   %xmm1,(%r14,%rcx,4)
add    %edi,%edx
mov    %r15d,%ebp
sub    %edx,%ebp
add    $0x1,%rcx
cmp    $0x100,%rcx
je     12f3 <CPU_BCV+0x1e3>
test   %ebp,%ebp
je     12f3 <CPU_BCV+0x1e3>
mov    (%rbx,%rcx,4),%edi
mov    %ecx,%eax
imul   %edi,%eax
add    %eax,%esi
pxor   %xmm1,%xmm1
test   %edx,%edx
je     1280 <CPU_BCV+0x170>
xorps  %xmm1,%xmm1
cvtsi2ss %edx,%xmm1
xorps  %xmm2,%xmm2
cvtsi2ss %ebp,%xmm2
mov    %r8d,%eax
sub    %esi,%eax
xorps  %xmm3,%xmm3
cvtsi2ss %eax,%xmm3
xorps  %xmm4,%xmm4
cvtsi2ss %esi,%xmm4
divss  %xmm2,%xmm3
divss  %xmm1,%xmm4
subss  %xmm3,%xmm4
divss  %xmm0,%xmm2
divss  %xmm0,%xmm1
mulss  %xmm2,%xmm1
mulss  %xmm4,%xmm1
mulss  %xmm4,%xmm1
jmp    1280 <CPU_BCV+0x170>
add    $0x8,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @CPU_BCV(i32* nocapture %0, float* nocapture writeonly %1, i32 %2) {
%4 = bitcast float* %1 to i8*
br label %5
5: ; preds = %5, %3
%6 = phi i64 [ 0, %3 ], [ %34, %5 ]
%7 = phi <4 x i32> [ zeroinitializer, %3 ], [ %32, %5 ]
%8 = phi <4 x i32> [ zeroinitializer, %3 ], [ %33, %5 ]
%9 = phi <4 x i32> [ <i32 0, i32 1, i32 2, i32 3>, %3 ], [ %35, %5 ]
%10 = add <4 x i32> %9, <i32 4, i32 4, i32 4, i32 4>
%11 = getelementptr inbounds i32, i32* %0, i64 %6
%12 = bitcast i32* %11 to <4 x i32>*
%13 = load <4 x i32>, <4 x i32>* %12
%14 = getelementptr inbounds i32, i32* %11, i64 4
%15 = bitcast i32* %14 to <4 x i32>*
%16 = load <4 x i32>, <4 x i32>* %15
%17 = mul nsw <4 x i32> %13, %9
%18 = mul nsw <4 x i32> %16, %10
%19 = add <4 x i32> %17, %7
%20 = add <4 x i32> %18, %8
%21 = or i64 %6, 8
%22 = add <4 x i32> %9, <i32 8, i32 8, i32 8, i32 8>
%23 = add <4 x i32> %9, <i32 12, i32 12, i32 12, i32 12>
%24 = getelementptr inbounds i32, i32* %0, i64 %21
%25 = bitcast i32* %24 to <4 x i32>*
%26 = load <4 x i32>, <4 x i32>* %25
%27 = getelementptr inbounds i32, i32* %24, i64 4
%28 = bitcast i32* %27 to <4 x i32>*
%29 = load <4 x i32>, <4 x i32>* %28
%30 = mul nsw <4 x i32> %26, %22
%31 = mul nsw <4 x i32> %29, %23
%32 = add <4 x i32> %30, %19
%33 = add <4 x i32> %31, %20
%34 = add nuw nsw i64 %6, 16
%35 = add <4 x i32> %9, <i32 16, i32 16, i32 16, i32 16>
%36 = icmp eq i64 %34, 256
br i1 %36, label %37, label %5
37: ; preds = %5
%38 = add <4 x i32> %33, %32
%39 = call i32 @llvm.vector.reduce.add.v4i32(<4 x i32> %38)
%40 = sitofp i32 %2 to float
br label %41
41: ; preds = %37, %54
%42 = phi i64 [ 0, %37 ], [ %71, %54 ]
%43 = phi i32 [ 0, %37 ], [ %69, %54 ]
%44 = phi i32 [ %2, %37 ], [ %70, %54 ]
%45 = phi i32 [ 0, %37 ], [ %50, %54 ]
%46 = getelementptr inbounds i32, i32* %0, i64 %42
%47 = load i32, i32* %46
%48 = trunc i64 %42 to i32
%49 = mul nsw i32 %47, %48
%50 = add nsw i32 %49, %45
%51 = sitofp i32 %43 to float
%52 = sitofp i32 %44 to float
%53 = icmp eq i32 %44, 0
br i1 %53, label %73, label %54
54: ; preds = %41
%55 = sub nsw i32 %39, %50
%56 = sitofp i32 %55 to float
%57 = fdiv float %56, %52
%58 = sitofp i32 %50 to float
%59 = fdiv float %58, %51
%60 = fdiv float %52, %40
%61 = fdiv float %51, %40
%62 = fmul float %60, %61
%63 = fsub float %59, %57
%64 = fmul float %62, %63
%65 = fmul float %63, %64
%66 = getelementptr inbounds float, float* %1, i64 %42
%67 = icmp eq i32 %43, 0
%68 = select i1 %67, float 0.000000e+00, float %65
store float %68, float* %66
%69 = add nsw i32 %47, %43
%70 = sub nsw i32 %2, %69
%71 = add nuw nsw i64 %42, 1
%72 = icmp eq i64 %71, 256
br i1 %72, label %73, label %41
73: ; preds = %41, %54
ret void
}