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


x86 assembly (clang, -O0):
<rogerstanimoto_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 <rogerstanimoto_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 <rogerstanimoto_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),%rax
movslq -0x18(%rbp),%rcx
movsbl (%rax,%rcx,1),%ecx
xor    %eax,%eax
cmp    $0x0,%ecx
mov    %al,-0x2a(%rbp)
je     11b0 <rogerstanimoto_distance_bool+0xb0>
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),%rcx
movslq -0x18(%rbp),%rdx
xor    %eax,%eax
cmpb   $0x0,(%rcx,%rdx,1)
mov    %al,-0x2b(%rbp)
jne    11ea <rogerstanimoto_distance_bool+0xea>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsbl (%rax,%rcx,1),%eax
cmp    $0x0,%eax
setne  %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 <rogerstanimoto_distance_bool+0x123>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
cmpb   $0x0,(%rax,%rcx,1)
setne  %al
xor    $0xff,%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 <rogerstanimoto_distance_bool+0x39>
mov    -0x28(%rbp),%eax
add    -0x24(%rbp),%eax
cvtsi2sd %eax,%xmm1
movsd  0xdaf(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x1c(%rbp),%xmm2
cvtsi2sdl -0x20(%rbp),%xmm1
addsd  %xmm1,%xmm2
mov    -0x28(%rbp),%eax
add    -0x24(%rbp),%eax
cvtsi2sd %eax,%xmm1
addsd  %xmm1,%xmm1
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<rogerstanimoto_distance_bool>:
test   %edx,%edx
jle    111f <rogerstanimoto_distance_bool+0x1f>
push   %rbx
mov    %edx,%r8d
cmp    $0x4,%edx
jae    112f <rogerstanimoto_distance_bool+0x2f>
xor    %ecx,%ecx
xor    %r11d,%r11d
xor    %r9d,%r9d
xor    %eax,%eax
xor    %r10d,%r10d
jmp    12a6 <rogerstanimoto_distance_bool+0x1a6>
xorpd  %xmm1,%xmm1
xor    %r10d,%r10d
xorpd  %xmm0,%xmm0
jmp    1273 <rogerstanimoto_distance_bool+0x173>
mov    %r8d,%ecx
and    $0xfffffffc,%ecx
pxor   %xmm9,%xmm9
xor    %eax,%eax
pcmpeqd %xmm8,%xmm8
movdqa 0xeb6(%rip),%xmm10
pxor   %xmm11,%xmm11
pxor   %xmm12,%xmm12
pxor   %xmm2,%xmm2
pxor   %xmm6,%xmm6
nopl   0x0(%rax)
mov    (%rdi,%rax,1),%edx
movd   %edx,%xmm4
mov    (%rsi,%rax,1),%ebx
or     %ebx,%edx
movd   %edx,%xmm7
pcmpeqb %xmm9,%xmm4
punpcklbw %xmm4,%xmm3
pxor   %xmm8,%xmm4
punpcklbw %xmm4,%xmm4
punpcklwd %xmm4,%xmm4
punpcklwd %xmm3,%xmm1
psrad  $0x18,%xmm1
movd   %ebx,%xmm0
pcmpeqb %xmm9,%xmm0
movdqa %xmm0,%xmm5
pxor   %xmm8,%xmm5
punpcklbw %xmm5,%xmm5
punpcklwd %xmm5,%xmm5
pand   %xmm10,%xmm5
punpcklbw %xmm0,%xmm0
punpcklwd %xmm0,%xmm0
movdqa %xmm1,%xmm3
pand   %xmm5,%xmm3
paddd  %xmm3,%xmm12
pandn  %xmm5,%xmm1
paddd  %xmm1,%xmm6
pand   %xmm10,%xmm0
pand   %xmm4,%xmm0
paddd  %xmm0,%xmm11
pcmpeqb %xmm9,%xmm7
punpcklbw %xmm7,%xmm7
punpcklwd %xmm7,%xmm7
pand   %xmm10,%xmm7
paddd  %xmm7,%xmm2
add    $0x4,%rax
cmp    %rax,%rcx
jne    1160 <rogerstanimoto_distance_bool+0x60>
pshufd $0xee,%xmm6,%xmm0
paddd  %xmm6,%xmm0
pshufd $0x55,%xmm0,%xmm1
paddd  %xmm0,%xmm1
movd   %xmm1,%r10d
pshufd $0xee,%xmm2,%xmm0
paddd  %xmm2,%xmm0
pshufd $0x55,%xmm0,%xmm1
paddd  %xmm0,%xmm1
movd   %xmm1,%eax
pshufd $0xee,%xmm12,%xmm0
paddd  %xmm12,%xmm0
pshufd $0x55,%xmm0,%xmm1
paddd  %xmm0,%xmm1
movd   %xmm1,%r9d
pshufd $0xee,%xmm11,%xmm0
paddd  %xmm11,%xmm0
pshufd $0x55,%xmm0,%xmm1
paddd  %xmm0,%xmm1
movd   %xmm1,%r11d
cmp    %r8,%rcx
jne    12a6 <rogerstanimoto_distance_bool+0x1a6>
add    %r11d,%r9d
xorps  %xmm0,%xmm0
cvtsi2sd %r9d,%xmm0
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
pop    %rbx
xorps  %xmm2,%xmm2
cvtsi2sd %r10d,%xmm2
addsd  %xmm0,%xmm0
addsd  %xmm1,%xmm2
addsd  %xmm0,%xmm2
divsd  %xmm2,%xmm0
ret
nopl   0x0(%rax)
add    %ebx,%r9d
test   %dl,%dl
sete   %dl
movzbl %dl,%edx
add    %edx,%eax
add    $0x1,%rcx
cmp    %rcx,%r8
je     1260 <rogerstanimoto_distance_bool+0x160>
movzbl (%rsi,%rcx,1),%edx
xor    %ebx,%ebx
test   %dl,%dl
setne  %bl
cmpb   $0x0,(%rdi,%rcx,1)
je     1290 <rogerstanimoto_distance_bool+0x190>
add    %ebx,%r10d
cmp    $0x1,%dl
adc    $0x0,%r11d
xor    %edx,%edx
jmp    1298 <rogerstanimoto_distance_bool+0x198>

LLVM IR (-O0):
define double @rogerstanimoto_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 = sext i8 %41 to i32
%43 = icmp ne i32 %42, 0
br i1 %43, label %44, label %52
44: ; preds = %32
%45 = load i8*, i8** %5
%46 = load i32, i32* %7
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds i8, i8* %45, i64 %47
%49 = load i8, i8* %48
%50 = icmp ne i8 %49, 0
%51 = xor i1 %50, true
br label %52
52: ; preds = %44, %32
%53 = phi i1 [ false, %32 ], [ %51, %44 ]
%54 = zext i1 %53 to i32
%55 = load i32, i32* %11
%56 = add nsw i32 %55, %54
store i32 %56, i32* %11
%57 = load i8*, i8** %4
%58 = load i32, i32* %7
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds i8, i8* %57, i64 %59
%61 = load i8, i8* %60
%62 = icmp ne i8 %61, 0
br i1 %62, label %71, label %63
63: ; preds = %52
%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 = sext i8 %68 to i32
%70 = icmp ne i32 %69, 0
br label %71
71: ; preds = %63, %52
%72 = phi i1 [ false, %52 ], [ %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 = icmp ne i8 %87, 0
%89 = xor i1 %88, true
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* %11
%100 = load i32, i32* %10
%101 = add nsw i32 %99, %100
%102 = sitofp i32 %101 to double
%103 = fmul double 2.000000e+00, %102
%104 = load i32, i32* %8
%105 = sitofp i32 %104 to double
%106 = load i32, i32* %9
%107 = sitofp i32 %106 to double
%108 = fadd double %105, %107
%109 = load i32, i32* %11
%110 = load i32, i32* %10
%111 = add nsw i32 %109, %110
%112 = sitofp i32 %111 to double
%113 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %112, double %108)
%114 = fdiv double %103, %113
ret double %114
}

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