double HammingDistance(char *v1, char *v2, int *filter, int length,
                       int numSelected) {
  int i;
  int d = 0;
  for (i = 0; i < length; i++)
    if (v1[i] != v2[i])
      d += filter[i];
  return ((double)d / numSelected);
}


x86 assembly (clang, -O0):
<HammingDistance>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
movl   $0x0,-0x28(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1175 <HammingDistance+0x75>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsbl (%rax,%rcx,1),%eax
mov    -0x10(%rbp),%rcx
movslq -0x24(%rbp),%rdx
movsbl (%rcx,%rdx,1),%ecx
cmp    %ecx,%eax
je     1162 <HammingDistance+0x62>
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
add    -0x28(%rbp),%eax
mov    %eax,-0x28(%rbp)
jmp    1167 <HammingDistance+0x67>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1125 <HammingDistance+0x25>
cvtsi2sdl -0x28(%rbp),%xmm0
cvtsi2sdl -0x20(%rbp),%xmm1
divsd  %xmm1,%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<HammingDistance>:
test   %ecx,%ecx
jle    1125 <HammingDistance+0x25>
mov    %ecx,%r9d
cmp    $0x1,%ecx
jne    1137 <HammingDistance+0x37>
xor    %ecx,%ecx
xor    %r11d,%r11d
test   $0x1,%r9b
je     1128 <HammingDistance+0x28>
mov    (%rdi,%rcx,1),%al
cmp    (%rsi,%rcx,1),%al
je     1128 <HammingDistance+0x28>
add    (%rdx,%rcx,4),%r11d
jmp    1128 <HammingDistance+0x28>
xor    %r11d,%r11d
cvtsi2sd %r11d,%xmm0
cvtsi2sd %r8d,%xmm1
divsd  %xmm1,%xmm0
ret
mov    %r9d,%r10d
and    $0xfffffffe,%r10d
xor    %ecx,%ecx
xor    %r11d,%r11d
jmp    1159 <HammingDistance+0x59>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x2,%rcx
cmp    %rcx,%r10
je     1111 <HammingDistance+0x11>
movzbl (%rdi,%rcx,1),%eax
cmp    (%rsi,%rcx,1),%al
je     1166 <HammingDistance+0x66>
add    (%rdx,%rcx,4),%r11d
movzbl 0x1(%rdi,%rcx,1),%eax
cmp    0x1(%rsi,%rcx,1),%al
je     1150 <HammingDistance+0x50>
add    0x4(%rdx,%rcx,4),%r11d
jmp    1150 <HammingDistance+0x50>

LLVM IR (-O0):
define double @HammingDistance(i8* %0, i8* %1, i32* %2, i32 %3, i32 %4) {
%6 = alloca i8*
%7 = alloca i8*
%8 = alloca i32*
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
store i8* %0, i8** %6
store i8* %1, i8** %7
store i32* %2, i32** %8
store i32 %3, i32* %9
store i32 %4, i32* %10
store i32 0, i32* %12
store i32 0, i32* %11
br label %13
13: ; preds = %40, %5
%14 = load i32, i32* %11
%15 = load i32, i32* %9
%16 = icmp slt i32 %14, %15
br i1 %16, label %17, label %43
17: ; preds = %13
%18 = load i8*, i8** %6
%19 = load i32, i32* %11
%20 = sext i32 %19 to i64
%21 = getelementptr inbounds i8, i8* %18, i64 %20
%22 = load i8, i8* %21
%23 = sext i8 %22 to i32
%24 = load i8*, i8** %7
%25 = load i32, i32* %11
%26 = sext i32 %25 to i64
%27 = getelementptr inbounds i8, i8* %24, i64 %26
%28 = load i8, i8* %27
%29 = sext i8 %28 to i32
%30 = icmp ne i32 %23, %29
br i1 %30, label %31, label %39
31: ; preds = %17
%32 = load i32*, i32** %8
%33 = load i32, i32* %11
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds i32, i32* %32, i64 %34
%36 = load i32, i32* %35
%37 = load i32, i32* %12
%38 = add nsw i32 %37, %36
store i32 %38, i32* %12
br label %39
39: ; preds = %31, %17
br label %40
40: ; preds = %39
%41 = load i32, i32* %11
%42 = add nsw i32 %41, 1
store i32 %42, i32* %11
br label %13
43: ; preds = %13
%44 = load i32, i32* %12
%45 = sitofp i32 %44 to double
%46 = load i32, i32* %10
%47 = sitofp i32 %46 to double
%48 = fdiv double %45, %47
ret double %48
}

LLVM IR (-O2):
define double @HammingDistance(i8* nocapture %0, i8* nocapture %1, i32* nocapture %2, i32 %3, i32 %4) {
%6 = icmp sgt i32 %3, 0
br i1 %6, label %7, label %58
7: ; preds = %5
%8 = zext i32 %3 to i64
%9 = and i64 %8, 1
%10 = icmp eq i32 %3, 1
br i1 %10, label %43, label %11
11: ; preds = %7
%12 = and i64 %8, 4294967294
br label %13
13: ; preds = %38, %11
%14 = phi i64 [ 0, %11 ], [ %40, %38 ]
%15 = phi i32 [ 0, %11 ], [ %39, %38 ]
%16 = phi i64 [ 0, %11 ], [ %41, %38 ]
%17 = getelementptr inbounds i8, i8* %0, i64 %14
%18 = load i8, i8* %17
%19 = getelementptr inbounds i8, i8* %1, i64 %14
%20 = load i8, i8* %19
%21 = icmp eq i8 %18, %20
br i1 %21, label %26, label %22
22: ; preds = %13
%23 = getelementptr inbounds i32, i32* %2, i64 %14
%24 = load i32, i32* %23
%25 = add nsw i32 %24, %15
br label %26
26: ; preds = %13, %22
%27 = phi i32 [ %25, %22 ], [ %15, %13 ]
%28 = or i64 %14, 1
%29 = getelementptr inbounds i8, i8* %0, i64 %28
%30 = load i8, i8* %29
%31 = getelementptr inbounds i8, i8* %1, i64 %28
%32 = load i8, i8* %31
%33 = icmp eq i8 %30, %32
br i1 %33, label %38, label %34
34: ; preds = %26
%35 = getelementptr inbounds i32, i32* %2, i64 %28
%36 = load i32, i32* %35
%37 = add nsw i32 %36, %27
br label %38
38: ; preds = %34, %26
%39 = phi i32 [ %37, %34 ], [ %27, %26 ]
%40 = add nuw nsw i64 %14, 2
%41 = add i64 %16, 2
%42 = icmp eq i64 %41, %12
br i1 %42, label %43, label %13
43: ; preds = %38, %7
%44 = phi i32 [ undef, %7 ], [ %39, %38 ]
%45 = phi i64 [ 0, %7 ], [ %40, %38 ]
%46 = phi i32 [ 0, %7 ], [ %39, %38 ]
%47 = icmp eq i64 %9, 0
br i1 %47, label %58, label %48
48: ; preds = %43
%49 = getelementptr inbounds i8, i8* %0, i64 %45
%50 = load i8, i8* %49
%51 = getelementptr inbounds i8, i8* %1, i64 %45
%52 = load i8, i8* %51
%53 = icmp eq i8 %50, %52
br i1 %53, label %58, label %54
54: ; preds = %48
%55 = getelementptr inbounds i32, i32* %2, i64 %45
%56 = load i32, i32* %55
%57 = add nsw i32 %56, %46
br label %58
58: ; preds = %43, %54, %48, %5
%59 = phi i32 [ 0, %5 ], [ %44, %43 ], [ %57, %54 ], [ %46, %48 ]
%60 = sitofp i32 %59 to double
%61 = sitofp i32 %4 to double
%62 = fdiv double %60, %61
ret double %62
}