void bgsub_classify(float *image, int size, float *m, float *v,
                    unsigned char *bgsub) {
  int i;
  float val;
  for (i = 0; i < size; ++i) {
    bgsub[i] = (image[i] - m[i]) * (image[i] - m[i]) > v[i];
  }
}


x86 assembly (clang, -O0):
<bgsub_classify>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1191 <bgsub_classify+0x91>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm1
mulss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
ucomiss (%rax,%rcx,4),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
mov    %al,%dl
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    %dl,(%rax,%rcx,1)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    111e <bgsub_classify+0x1e>
pop    %rbp
ret

x86 assembly (clang, -O2):
<bgsub_classify>:
push   %rbp
push   %r14
push   %rbx
test   %esi,%esi
jle    124b <bgsub_classify+0x14b>
mov    %esi,%r9d
cmp    $0x8,%esi
jae    111b <bgsub_classify+0x1b>
xor    %esi,%esi
jmp    11e2 <bgsub_classify+0xe2>
lea    (%rdi,%r9,4),%rax
cmp    %r8,%rax
seta   %r14b
lea    (%r8,%r9,1),%rsi
cmp    %rdi,%rsi
seta   %bl
lea    (%rdx,%r9,4),%rax
cmp    %r8,%rax
seta   %al
lea    (%rcx,%r9,4),%rbp
cmp    %rdx,%rsi
seta   %r11b
cmp    %r8,%rbp
seta   %r10b
cmp    %rcx,%rsi
seta   %bpl
xor    %esi,%esi
test   %bl,%r14b
jne    11e2 <bgsub_classify+0xe2>
and    %r11b,%al
jne    11e2 <bgsub_classify+0xe2>
and    %bpl,%r10b
jne    11e2 <bgsub_classify+0xe2>
mov    %r9d,%esi
and    $0xfffffff8,%esi
xor    %eax,%eax
movdqa 0xe84(%rip),%xmm0
nopl   0x0(%rax)
movups (%rdi,%rax,4),%xmm1
movups 0x10(%rdi,%rax,4),%xmm2
movups (%rdx,%rax,4),%xmm3
subps  %xmm3,%xmm1
movups 0x10(%rdx,%rax,4),%xmm3
subps  %xmm3,%xmm2
mulps  %xmm1,%xmm1
mulps  %xmm2,%xmm2
movups (%rcx,%rax,4),%xmm3
cmpltps %xmm1,%xmm3
packssdw %xmm3,%xmm3
packsswb %xmm3,%xmm3
movups 0x10(%rcx,%rax,4),%xmm1
cmpltps %xmm2,%xmm1
packssdw %xmm1,%xmm1
packsswb %xmm1,%xmm1
pand   %xmm0,%xmm3
pand   %xmm0,%xmm1
movd   %xmm3,(%r8,%rax,1)
movd   %xmm1,0x4(%r8,%rax,1)
add    $0x8,%rax
cmp    %rax,%rsi
jne    1180 <bgsub_classify+0x80>
cmp    %r9,%rsi
je     124b <bgsub_classify+0x14b>
mov    %rsi,%rax
not    %rax
test   $0x1,%r9b
je     1209 <bgsub_classify+0x109>
movss  (%rdi,%rsi,4),%xmm0
subss  (%rdx,%rsi,4),%xmm0
mulss  %xmm0,%xmm0
ucomiss (%rcx,%rsi,4),%xmm0
seta   (%r8,%rsi,1)
or     $0x1,%rsi
add    %r9,%rax
je     124b <bgsub_classify+0x14b>
xchg   %ax,%ax
movss  (%rdi,%rsi,4),%xmm0
subss  (%rdx,%rsi,4),%xmm0
mulss  %xmm0,%xmm0
ucomiss (%rcx,%rsi,4),%xmm0
seta   (%r8,%rsi,1)
movss  0x4(%rdi,%rsi,4),%xmm0
subss  0x4(%rdx,%rsi,4),%xmm0
mulss  %xmm0,%xmm0
ucomiss 0x4(%rcx,%rsi,4),%xmm0
seta   0x1(%r8,%rsi,1)
add    $0x2,%rsi
cmp    %rsi,%r9
jne    1210 <bgsub_classify+0x110>
pop    %rbx
pop    %r14
pop    %rbp
ret

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

LLVM IR (-O2):
define void @bgsub_classify(float* nocapture %0, i32 %1, float* nocapture %2, float* nocapture %3, i8* nocapture writeonly %4) {
%6 = bitcast float* %0 to i8*
%7 = bitcast float* %2 to i8*
%8 = bitcast float* %3 to i8*
%9 = icmp sgt i32 %1, 0
br i1 %9, label %10, label %119
10: ; preds = %5
%11 = zext i32 %1 to i64
%12 = icmp ult i32 %1, 8
br i1 %12, label %70, label %13
13: ; preds = %10
%14 = getelementptr i8, i8* %4, i64 %11
%15 = getelementptr float, float* %0, i64 %11
%16 = bitcast float* %15 to i8*
%17 = getelementptr float, float* %2, i64 %11
%18 = bitcast float* %17 to i8*
%19 = getelementptr float, float* %3, i64 %11
%20 = bitcast float* %19 to i8*
%21 = icmp ugt i8* %16, %4
%22 = icmp ugt i8* %14, %6
%23 = and i1 %21, %22
%24 = icmp ugt i8* %18, %4
%25 = icmp ugt i8* %14, %7
%26 = and i1 %24, %25
%27 = or i1 %23, %26
%28 = icmp ugt i8* %20, %4
%29 = icmp ugt i8* %14, %8
%30 = and i1 %28, %29
%31 = or i1 %27, %30
br i1 %31, label %70, label %32
32: ; preds = %13
%33 = and i64 %11, 4294967288
br label %34
34: ; preds = %34, %32
%35 = phi i64 [ 0, %32 ], [ %66, %34 ]
%36 = getelementptr inbounds float, float* %0, i64 %35
%37 = bitcast float* %36 to <4 x float>*
%38 = load <4 x float>, <4 x float>* %37
%39 = getelementptr inbounds float, float* %36, i64 4
%40 = bitcast float* %39 to <4 x float>*
%41 = load <4 x float>, <4 x float>* %40
%42 = getelementptr inbounds float, float* %2, i64 %35
%43 = bitcast float* %42 to <4 x float>*
%44 = load <4 x float>, <4 x float>* %43
%45 = getelementptr inbounds float, float* %42, i64 4
%46 = bitcast float* %45 to <4 x float>*
%47 = load <4 x float>, <4 x float>* %46
%48 = fsub <4 x float> %38, %44
%49 = fsub <4 x float> %41, %47
%50 = fmul <4 x float> %48, %48
%51 = fmul <4 x float> %49, %49
%52 = getelementptr inbounds float, float* %3, i64 %35
%53 = bitcast float* %52 to <4 x float>*
%54 = load <4 x float>, <4 x float>* %53
%55 = getelementptr inbounds float, float* %52, i64 4
%56 = bitcast float* %55 to <4 x float>*
%57 = load <4 x float>, <4 x float>* %56
%58 = fcmp ogt <4 x float> %50, %54
%59 = fcmp ogt <4 x float> %51, %57
%60 = zext <4 x i1> %58 to <4 x i8>
%61 = zext <4 x i1> %59 to <4 x i8>
%62 = getelementptr inbounds i8, i8* %4, i64 %35
%63 = bitcast i8* %62 to <4 x i8>*
store <4 x i8> %60, <4 x i8>* %63
%64 = getelementptr inbounds i8, i8* %62, i64 4
%65 = bitcast i8* %64 to <4 x i8>*
store <4 x i8> %61, <4 x i8>* %65
%66 = add nuw i64 %35, 8
%67 = icmp eq i64 %66, %33
br i1 %67, label %68, label %34
68: ; preds = %34
%69 = icmp eq i64 %33, %11
br i1 %69, label %119, label %70
70: ; preds = %13, %10, %68
%71 = phi i64 [ 0, %13 ], [ 0, %10 ], [ %33, %68 ]
%72 = xor i64 %71, -1
%73 = and i64 %11, 1
%74 = icmp eq i64 %73, 0
br i1 %74, label %88, label %75
75: ; preds = %70
%76 = getelementptr inbounds float, float* %0, i64 %71
%77 = load float, float* %76
%78 = getelementptr inbounds float, float* %2, i64 %71
%79 = load float, float* %78
%80 = fsub float %77, %79
%81 = fmul float %80, %80
%82 = getelementptr inbounds float, float* %3, i64 %71
%83 = load float, float* %82
%84 = fcmp ogt float %81, %83
%85 = zext i1 %84 to i8
%86 = getelementptr inbounds i8, i8* %4, i64 %71
store i8 %85, i8* %86
%87 = or i64 %71, 1
br label %88
88: ; preds = %75, %70
%89 = phi i64 [ %71, %70 ], [ %87, %75 ]
%90 = sub nsw i64 0, %11
%91 = icmp eq i64 %72, %90
br i1 %91, label %119, label %92
92: ; preds = %88, %92
%93 = phi i64 [ %117, %92 ], [ %89, %88 ]
%94 = getelementptr inbounds float, float* %0, i64 %93
%95 = load float, float* %94
%96 = getelementptr inbounds float, float* %2, i64 %93
%97 = load float, float* %96
%98 = fsub float %95, %97
%99 = fmul float %98, %98
%100 = getelementptr inbounds float, float* %3, i64 %93
%101 = load float, float* %100
%102 = fcmp ogt float %99, %101
%103 = zext i1 %102 to i8
%104 = getelementptr inbounds i8, i8* %4, i64 %93
store i8 %103, i8* %104
%105 = add nuw nsw i64 %93, 1
%106 = getelementptr inbounds float, float* %0, i64 %105
%107 = load float, float* %106
%108 = getelementptr inbounds float, float* %2, i64 %105
%109 = load float, float* %108
%110 = fsub float %107, %109
%111 = fmul float %110, %110
%112 = getelementptr inbounds float, float* %3, i64 %105
%113 = load float, float* %112
%114 = fcmp ogt float %111, %113
%115 = zext i1 %114 to i8
%116 = getelementptr inbounds i8, i8* %4, i64 %105
store i8 %115, i8* %116
%117 = add nuw nsw i64 %93, 2
%118 = icmp eq i64 %117, %11
br i1 %118, label %119, label %92
119: ; preds = %88, %92, %68, %5
ret void
}