double sokalmichener_distance_bool(const char *u, const char *v, int n) {
  int _i = 0;
  int ntt = 0, nft = 0, ntf = 0, nff = 0;
  for (_i = 0; _i < n; _i++) {
    ntt += (u[_i] && v[_i]);
    nff += (!u[_i] && !v[_i]);
    ntf += (u[_i] && !v[_i]);
    nft += (!u[_i] && v[_i]);
  }
  return (2.0 * (ntf + nft)) / (2.0 * (ntf + nft) + ntt + nff);
}


x86 assembly (clang, -O0):
<sokalmichener_distance_bool>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movl   $0x0,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
movl   $0x0,-0x28(%rbp)
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    123f <sokalmichener_distance_bool+0x13f>
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsbl (%rax,%rcx,1),%ecx
xor    %eax,%eax
cmp    $0x0,%ecx
mov    %al,-0x29(%rbp)
je     1174 <sokalmichener_distance_bool+0x74>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsbl (%rax,%rcx,1),%eax
cmp    $0x0,%eax
setne  %al
mov    %al,-0x29(%rbp)
mov    -0x29(%rbp),%al
and    $0x1,%al
movzbl %al,%eax
add    -0x1c(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x8(%rbp),%rcx
movslq -0x18(%rbp),%rdx
xor    %eax,%eax
cmpb   $0x0,(%rcx,%rdx,1)
mov    %al,-0x2a(%rbp)
jne    11ad <sokalmichener_distance_bool+0xad>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
cmpb   $0x0,(%rax,%rcx,1)
setne  %al
xor    $0xff,%al
mov    %al,-0x2a(%rbp)
mov    -0x2a(%rbp),%al
and    $0x1,%al
movzbl %al,%eax
add    -0x28(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsbl (%rax,%rcx,1),%ecx
xor    %eax,%eax
cmp    $0x0,%ecx
mov    %al,-0x2b(%rbp)
je     11e9 <sokalmichener_distance_bool+0xe9>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
cmpb   $0x0,(%rax,%rcx,1)
setne  %al
xor    $0xff,%al
mov    %al,-0x2b(%rbp)
mov    -0x2b(%rbp),%al
and    $0x1,%al
movzbl %al,%eax
add    -0x24(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x8(%rbp),%rcx
movslq -0x18(%rbp),%rdx
xor    %eax,%eax
cmpb   $0x0,(%rcx,%rdx,1)
mov    %al,-0x2c(%rbp)
jne    1223 <sokalmichener_distance_bool+0x123>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsbl (%rax,%rcx,1),%eax
cmp    $0x0,%eax
setne  %al
mov    %al,-0x2c(%rbp)
mov    -0x2c(%rbp),%al
and    $0x1,%al
movzbl %al,%eax
add    -0x20(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1139 <sokalmichener_distance_bool+0x39>
mov    -0x24(%rbp),%eax
add    -0x20(%rbp),%eax
cvtsi2sd %eax,%xmm1
movsd  0xdaf(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    -0x24(%rbp),%eax
add    -0x20(%rbp),%eax
cvtsi2sd %eax,%xmm1
cvtsi2sdl -0x1c(%rbp),%xmm2
addsd  %xmm1,%xmm1
addsd  %xmm2,%xmm1
cvtsi2sdl -0x28(%rbp),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<sokalmichener_distance_bool>:
test   %edx,%edx
jle    1120 <sokalmichener_distance_bool+0x20>
push   %rbp
push   %rbx
mov    %edx,%r8d
cmp    $0x4,%edx
jae    112f <sokalmichener_distance_bool+0x2f>
xor    %ebx,%ebx
xor    %r9d,%r9d
xor    %r11d,%r11d
xor    %eax,%eax
xor    %r10d,%r10d
jmp    12ab <sokalmichener_distance_bool+0x1ab>
xorpd  %xmm0,%xmm0
xor    %r10d,%r10d
xor    %r9d,%r9d
jmp    126f <sokalmichener_distance_bool+0x16f>
mov    %r8d,%ebx
and    $0xfffffffc,%ebx
pxor   %xmm9,%xmm9
xor    %eax,%eax
pcmpeqd %xmm8,%xmm8
movdqa 0xeb6(%rip),%xmm10
pxor   %xmm11,%xmm11
pxor   %xmm1,%xmm1
pxor   %xmm2,%xmm2
pxor   %xmm6,%xmm6
nopl   0x0(%rax,%rax,1)
mov    (%rdi,%rax,1),%edx
movd   %edx,%xmm7
mov    (%rsi,%rax,1),%ecx
or     %ecx,%edx
movd   %edx,%xmm4
pcmpeqb %xmm9,%xmm7
movdqa %xmm7,%xmm3
punpcklbw %xmm3,%xmm3
punpcklwd %xmm3,%xmm3
pxor   %xmm8,%xmm7
punpcklbw %xmm7,%xmm7
punpcklwd %xmm7,%xmm7
psrad  $0x18,%xmm7
movd   %ecx,%xmm0
pcmpeqb %xmm9,%xmm0
movdqa %xmm0,%xmm5
pxor   %xmm8,%xmm5
punpcklbw %xmm5,%xmm5
punpcklwd %xmm5,%xmm5
psrad  $0x18,%xmm5
punpcklbw %xmm0,%xmm0
punpcklwd %xmm0,%xmm0
pand   %xmm10,%xmm0
pand   %xmm7,%xmm0
paddd  %xmm0,%xmm1
pand   %xmm5,%xmm7
psubd  %xmm7,%xmm6
pcmpeqb %xmm9,%xmm4
punpcklbw %xmm4,%xmm4
punpcklwd %xmm4,%xmm4
pand   %xmm10,%xmm4
paddd  %xmm4,%xmm11
pand   %xmm10,%xmm3
pand   %xmm5,%xmm3
paddd  %xmm3,%xmm2
add    $0x4,%rax
cmp    %rax,%rbx
jne    1160 <sokalmichener_distance_bool+0x60>
pshufd $0xee,%xmm6,%xmm0
paddd  %xmm6,%xmm0
pshufd $0x55,%xmm0,%xmm3
paddd  %xmm0,%xmm3
movd   %xmm3,%r10d
pshufd $0xee,%xmm2,%xmm0
paddd  %xmm2,%xmm0
pshufd $0x55,%xmm0,%xmm2
paddd  %xmm0,%xmm2
movd   %xmm2,%eax
pshufd $0xee,%xmm1,%xmm0
paddd  %xmm1,%xmm0
pshufd $0x55,%xmm0,%xmm1
paddd  %xmm0,%xmm1
movd   %xmm1,%r11d
pshufd $0xee,%xmm11,%xmm0
paddd  %xmm11,%xmm0
pshufd $0x55,%xmm0,%xmm1
paddd  %xmm0,%xmm1
movd   %xmm1,%r9d
cmp    %r8,%rbx
jne    12ab <sokalmichener_distance_bool+0x1ab>
add    %eax,%r11d
xorps  %xmm0,%xmm0
cvtsi2sd %r11d,%xmm0
pop    %rbx
pop    %rbp
addsd  %xmm0,%xmm0
xorps  %xmm1,%xmm1
cvtsi2sd %r10d,%xmm1
addsd  %xmm0,%xmm1
xorps  %xmm2,%xmm2
cvtsi2sd %r9d,%xmm2
addsd  %xmm1,%xmm2
divsd  %xmm2,%xmm0
ret
xor    %ecx,%ecx
test   %dl,%dl
sete   %cl
setne  %dl
add    %ecx,%r9d
movzbl %dl,%ecx
add    %ecx,%eax
add    $0x1,%rbx
cmp    %rbx,%r8
je     1262 <sokalmichener_distance_bool+0x162>
cmpb   $0x0,(%rdi,%rbx,1)
movzbl (%rsi,%rbx,1),%edx
je     1290 <sokalmichener_distance_bool+0x190>
xor    %ecx,%ecx
xor    %ebp,%ebp
test   %dl,%dl
setne  %cl
sete   %bpl
add    %ecx,%r10d
add    %ebp,%r11d
xor    %edx,%edx
jmp    129d <sokalmichener_distance_bool+0x19d>

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

LLVM IR (-O2):
define double @sokalmichener_distance_bool(i8* nocapture %0, i8* nocapture %1, i32 %2) {
%4 = icmp sgt i32 %2, 0
br i1 %4, label %5, label %94
5: ; preds = %3
%6 = zext i32 %2 to i64
%7 = icmp ult i32 %2, 4
br i1 %7, label %49, label %8
8: ; preds = %5
%9 = and i64 %6, 4294967292
br label %10
10: ; preds = %10, %8
%11 = phi i64 [ 0, %8 ], [ %41, %10 ]
%12 = phi <4 x i32> [ zeroinitializer, %8 ], [ %37, %10 ]
%13 = phi <4 x i32> [ zeroinitializer, %8 ], [ %30, %10 ]
%14 = phi <4 x i32> [ zeroinitializer, %8 ], [ %40, %10 ]
%15 = phi <4 x i32> [ zeroinitializer, %8 ], [ %34, %10 ]
%16 = getelementptr inbounds i8, i8* %0, i64 %11
%17 = bitcast i8* %16 to <4 x i8>*
%18 = load <4 x i8>, <4 x i8>* %17
%19 = icmp eq <4 x i8> %18, zeroinitializer
%20 = getelementptr inbounds i8, i8* %1, i64 %11
%21 = bitcast i8* %20 to <4 x i8>*
%22 = load <4 x i8>, <4 x i8>* %21
%23 = icmp ne <4 x i8> %22, zeroinitializer
%24 = icmp eq <4 x i8> %22, zeroinitializer
%25 = icmp eq <4 x i8> %22, zeroinitializer
%26 = icmp ne <4 x i8> %22, zeroinitializer
%27 = xor <4 x i1> %19, <i1 true, i1 true, i1 true, i1 true>
%28 = select <4 x i1> %27, <4 x i1> %24, <4 x i1> zeroinitializer
%29 = zext <4 x i1> %28 to <4 x i32>
%30 = add <4 x i32> %13, %29
%31 = xor <4 x i1> %19, <i1 true, i1 true, i1 true, i1 true>
%32 = select <4 x i1> %31, <4 x i1> %23, <4 x i1> zeroinitializer
%33 = zext <4 x i1> %32 to <4 x i32>
%34 = add <4 x i32> %15, %33
%35 = select <4 x i1> %19, <4 x i1> %25, <4 x i1> zeroinitializer
%36 = zext <4 x i1> %35 to <4 x i32>
%37 = add <4 x i32> %12, %36
%38 = select <4 x i1> %19, <4 x i1> %26, <4 x i1> zeroinitializer
%39 = zext <4 x i1> %38 to <4 x i32>
%40 = add <4 x i32> %14, %39
%41 = add nuw i64 %11, 4
%42 = icmp eq i64 %41, %9
br i1 %42, label %43, label %10
43: ; preds = %10
%44 = call i32 @llvm.vector.reduce.add.v4i32(<4 x i32> %34)
%45 = call i32 @llvm.vector.reduce.add.v4i32(<4 x i32> %40)
%46 = call i32 @llvm.vector.reduce.add.v4i32(<4 x i32> %30)
%47 = call i32 @llvm.vector.reduce.add.v4i32(<4 x i32> %37)
%48 = icmp eq i64 %9, %6
br i1 %48, label %87, label %49
49: ; preds = %5, %43
%50 = phi i64 [ 0, %5 ], [ %9, %43 ]
%51 = phi i32 [ 0, %5 ], [ %47, %43 ]
%52 = phi i32 [ 0, %5 ], [ %46, %43 ]
%53 = phi i32 [ 0, %5 ], [ %45, %43 ]
%54 = phi i32 [ 0, %5 ], [ %44, %43 ]
br label %55
55: ; preds = %49, %78
%56 = phi i64 [ %85, %78 ], [ %50, %49 ]
%57 = phi i32 [ %81, %78 ], [ %51, %49 ]
%58 = phi i32 [ %79, %78 ], [ %52, %49 ]
%59 = phi i32 [ %84, %78 ], [ %53, %49 ]
%60 = phi i32 [ %80, %78 ], [ %54, %49 ]
%61 = getelementptr inbounds i8, i8* %0, i64 %56
%62 = load i8, i8* %61
%63 = icmp eq i8 %62, 0
%64 = getelementptr inbounds i8, i8* %1, i64 %56
%65 = load i8, i8* %64
br i1 %63, label %73, label %66
66: ; preds = %55
%67 = icmp ne i8 %65, 0
%68 = zext i1 %67 to i32
%69 = add nuw nsw i32 %60, %68
%70 = icmp eq i8 %65, 0
%71 = zext i1 %70 to i32
%72 = add nuw nsw i32 %58, %71
br label %78
73: ; preds = %55
%74 = icmp eq i8 %65, 0
%75 = zext i1 %74 to i32
%76 = add nuw nsw i32 %57, %75
%77 = icmp ne i8 %65, 0
br label %78
78: ; preds = %66, %73
%79 = phi i32 [ %72, %66 ], [ %58, %73 ]
%80 = phi i32 [ %69, %66 ], [ %60, %73 ]
%81 = phi i32 [ %57, %66 ], [ %76, %73 ]
%82 = phi i1 [ false, %66 ], [ %77, %73 ]
%83 = zext i1 %82 to i32
%84 = add nuw nsw i32 %59, %83
%85 = add nuw nsw i64 %56, 1
%86 = icmp eq i64 %85, %6
br i1 %86, label %87, label %55
87: ; preds = %78, %43
%88 = phi i32 [ %46, %43 ], [ %79, %78 ]
%89 = phi i32 [ %44, %43 ], [ %80, %78 ]
%90 = phi i32 [ %47, %43 ], [ %81, %78 ]
%91 = phi i32 [ %45, %43 ], [ %84, %78 ]
%92 = add nuw nsw i32 %88, %91
%93 = sitofp i32 %92 to double
br label %94
94: ; preds = %87, %3
%95 = phi i32 [ 0, %3 ], [ %89, %87 ]
%96 = phi i32 [ 0, %3 ], [ %90, %87 ]
%97 = phi double [ 0.000000e+00, %3 ], [ %93, %87 ]
%98 = fmul double %97, 2.000000e+00
%99 = sitofp i32 %95 to double
%100 = tail call double @llvm.fmuladd.f64(double %97, double 2.000000e+00, double %99)
%101 = sitofp i32 %96 to double
%102 = fadd double %100, %101
%103 = fdiv double %98, %102
ret double %103
}