unsigned long find_next_zero_bit(const unsigned long *addr, unsigned long size,
                                 unsigned long offset) {
  const unsigned long *p = addr + ((offset) / (8 * sizeof(long)));
  unsigned long result = offset & ~((8 * sizeof(long)) - 1);
  unsigned long tmp;
  if (offset >= size)
    return size;
  size -= result;
  offset %= (8 * sizeof(long));
  if (offset) {
    tmp = *(p++);
    tmp |= ~0UL >> ((8 * sizeof(long)) - offset);
    if (size < (8 * sizeof(long)))
      goto found_first;
    if (~tmp)
      goto found_middle;
    size -= (8 * sizeof(long));
    result += (8 * sizeof(long));
  }
  while (size & ~((8 * sizeof(long)) - 1)) {
    if (~(tmp = *(p++)))
      goto found_middle;
    result += (8 * sizeof(long));
    size -= (8 * sizeof(long));
  }
  if (!size)
    return result;
  tmp = *p;
found_first:
  tmp |= ~0UL << size;
  if (tmp == ~0UL)
    return result + size;
found_middle:
  return result + (~(tmp) ? __builtin_ffsl(~(tmp)) - 1 : 0);
}


x86 assembly (clang, -O0):
<find_next_zero_bit>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
shr    $0x6,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x28(%rbp)
mov    -0x20(%rbp),%rax
and    $0xffffffffffffffc0,%rax
mov    %rax,-0x30(%rbp)
mov    -0x20(%rbp),%rax
cmp    -0x18(%rbp),%rax
jb     114e <find_next_zero_bit+0x4e>
mov    -0x18(%rbp),%rax
mov    %rax,-0x8(%rbp)
jmp    12f5 <find_next_zero_bit+0x1f5>
mov    -0x30(%rbp),%rcx
mov    -0x18(%rbp),%rax
sub    %rcx,%rax
mov    %rax,-0x18(%rbp)
mov    -0x20(%rbp),%rax
and    $0x3f,%rax
mov    %rax,-0x20(%rbp)
cmpq   $0x0,-0x20(%rbp)
je     11e4 <find_next_zero_bit+0xe4>
mov    -0x28(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x28(%rbp)
mov    (%rax),%rax
mov    %rax,-0x38(%rbp)
mov    $0x40,%ecx
sub    -0x20(%rbp),%rcx
mov    $0xffffffffffffffff,%rax
shr    %cl,%rax
or     -0x38(%rbp),%rax
mov    %rax,-0x38(%rbp)
cmpq   $0x40,-0x18(%rbp)
jae    11b5 <find_next_zero_bit+0xb5>
jmp    1264 <find_next_zero_bit+0x164>
mov    -0x38(%rbp),%rax
xor    $0xffffffffffffffff,%rax
cmp    $0x0,%rax
je     11cc <find_next_zero_bit+0xcc>
jmp    129b <find_next_zero_bit+0x19b>
mov    -0x18(%rbp),%rax
sub    $0x40,%rax
mov    %rax,-0x18(%rbp)
mov    -0x30(%rbp),%rax
add    $0x40,%rax
mov    %rax,-0x30(%rbp)
jmp    11e9 <find_next_zero_bit+0xe9>
mov    -0x18(%rbp),%rax
and    $0xffffffffffffffc0,%rax
cmp    $0x0,%rax
je     1241 <find_next_zero_bit+0x141>
mov    -0x28(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x28(%rbp)
mov    (%rax),%rax
mov    %rax,-0x38(%rbp)
xor    $0xffffffffffffffff,%rax
cmp    $0x0,%rax
je     1224 <find_next_zero_bit+0x124>
jmp    129b <find_next_zero_bit+0x19b>
mov    -0x30(%rbp),%rax
add    $0x40,%rax
mov    %rax,-0x30(%rbp)
mov    -0x18(%rbp),%rax
sub    $0x40,%rax
mov    %rax,-0x18(%rbp)
jmp    11e9 <find_next_zero_bit+0xe9>
cmpq   $0x0,-0x18(%rbp)
jne    1259 <find_next_zero_bit+0x159>
mov    -0x30(%rbp),%rax
mov    %rax,-0x8(%rbp)
jmp    12f5 <find_next_zero_bit+0x1f5>
mov    -0x28(%rbp),%rax
mov    (%rax),%rax
mov    %rax,-0x38(%rbp)
mov    -0x18(%rbp),%rcx
mov    $0xffffffffffffffff,%rax
shl    %cl,%rax
or     -0x38(%rbp),%rax
mov    %rax,-0x38(%rbp)
cmpq   $0xffffffffffffffff,-0x38(%rbp)
jne    1296 <find_next_zero_bit+0x196>
mov    -0x30(%rbp),%rax
add    -0x18(%rbp),%rax
mov    %rax,-0x8(%rbp)
jmp    12f5 <find_next_zero_bit+0x1f5>
jmp    129b <find_next_zero_bit+0x19b>
mov    -0x30(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x38(%rbp),%rax
xor    $0xffffffffffffffff,%rax
cmp    $0x0,%rax
je     12da <find_next_zero_bit+0x1da>
mov    -0x38(%rbp),%rdx
xor    $0xffffffffffffffff,%rdx
bsf    %rdx,%rax
add    $0x1,%rax
xor    %ecx,%ecx
cmp    $0x0,%rdx
cmove  %rcx,%rax
sub    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    12e4 <find_next_zero_bit+0x1e4>
xor    %eax,%eax
mov    %eax,-0x44(%rbp)
jmp    12e4 <find_next_zero_bit+0x1e4>
mov    -0x40(%rbp),%rax
mov    -0x44(%rbp),%ecx
movslq %ecx,%rcx
add    %rcx,%rax
mov    %rax,-0x8(%rbp)
mov    -0x8(%rbp),%rax
pop    %rbp
ret

x86 assembly (clang, -O2):
<find_next_zero_bit>:
mov    %rsi,%rax
cmp    %rsi,%rdx
jae    11a5 <find_next_zero_bit+0xa5>
mov    %rdx,%rcx
shr    $0x6,%rcx
lea    (%rdi,%rcx,8),%rdi
mov    %rdx,%r8
and    $0xffffffffffffffc0,%r8
sub    %r8,%rax
and    $0x3f,%rdx
je     1150 <find_next_zero_bit+0x50>
neg    %dl
mov    $0xffffffffffffffff,%rsi
mov    %edx,%ecx
shr    %cl,%rsi
or     (%rdi),%rsi
cmp    $0x40,%rax
jb     1183 <find_next_zero_bit+0x83>
cmp    $0xffffffffffffffff,%rsi
jne    119b <find_next_zero_bit+0x9b>
add    $0x8,%rdi
add    $0xffffffffffffffc0,%rax
add    $0x40,%r8
cmp    $0x40,%rax
jb     117b <find_next_zero_bit+0x7b>
cs nopw 0x0(%rax,%rax,1)
mov    (%rdi),%rsi
cmp    $0xffffffffffffffff,%rsi
jne    119b <find_next_zero_bit+0x9b>
add    $0x8,%rdi
add    $0x40,%r8
add    $0xffffffffffffffc0,%rax
cmp    $0x3f,%rax
ja     1160 <find_next_zero_bit+0x60>
test   %rax,%rax
je     11a6 <find_next_zero_bit+0xa6>
mov    (%rdi),%rsi
mov    %rsi,%rdx
mov    $0xffffffffffffffff,%rsi
mov    %eax,%ecx
shl    %cl,%rsi
or     %rdx,%rsi
cmp    $0xffffffffffffffff,%rsi
je     11a2 <find_next_zero_bit+0xa2>
not    %rsi
bsf    %rsi,%rax
add    %r8,%rax
ret
mov    %r8,%rax
ret

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

LLVM IR (-O2):
define i64 @find_next_zero_bit(i64* nocapture %0, i64 %1, i64 %2) {
%4 = lshr i64 %2, 6
%5 = getelementptr inbounds i64, i64* %0, i64 %4
%6 = and i64 %2, -64
%7 = icmp ult i64 %2, %1
br i1 %7, label %8, label %62
8: ; preds = %3
%9 = sub i64 %1, %6
%10 = and i64 %2, 63
%11 = icmp eq i64 %10, 0
br i1 %11, label %24, label %12
12: ; preds = %8
%13 = getelementptr inbounds i64, i64* %5, i64 1
%14 = load i64, i64* %5
%15 = sub nuw nsw i64 64, %10
%16 = lshr i64 -1, %15
%17 = or i64 %14, %16
%18 = icmp ult i64 %9, 64
br i1 %18, label %47, label %19
19: ; preds = %12
%20 = icmp eq i64 %17, -1
br i1 %20, label %21, label %56
21: ; preds = %19
%22 = add i64 %9, -64
%23 = add i64 %6, 64
br label %24
24: ; preds = %21, %8
%25 = phi i64 [ %22, %21 ], [ %9, %8 ]
%26 = phi i64* [ %13, %21 ], [ %5, %8 ]
%27 = phi i64 [ %23, %21 ], [ %6, %8 ]
%28 = icmp ult i64 %25, 64
br i1 %28, label %40, label %29
29: ; preds = %24, %35
%30 = phi i64 [ %37, %35 ], [ %27, %24 ]
%31 = phi i64* [ %36, %35 ], [ %26, %24 ]
%32 = phi i64 [ %38, %35 ], [ %25, %24 ]
%33 = load i64, i64* %31
%34 = icmp eq i64 %33, -1
br i1 %34, label %35, label %56
35: ; preds = %29
%36 = getelementptr inbounds i64, i64* %31, i64 1
%37 = add i64 %30, 64
%38 = add i64 %32, -64
%39 = icmp ult i64 %38, 64
br i1 %39, label %40, label %29
40: ; preds = %35, %24
%41 = phi i64 [ %25, %24 ], [ %38, %35 ]
%42 = phi i64* [ %26, %24 ], [ %36, %35 ]
%43 = phi i64 [ %27, %24 ], [ %37, %35 ]
%44 = icmp eq i64 %41, 0
br i1 %44, label %62, label %45
45: ; preds = %40
%46 = load i64, i64* %42
br label %47
47: ; preds = %12, %45
%48 = phi i64 [ %9, %12 ], [ %41, %45 ]
%49 = phi i64 [ %6, %12 ], [ %43, %45 ]
%50 = phi i64 [ %17, %12 ], [ %46, %45 ]
%51 = shl nsw i64 -1, %48
%52 = or i64 %50, %51
%53 = icmp eq i64 %52, -1
br i1 %53, label %54, label %56
54: ; preds = %47
%55 = add i64 %49, %48
br label %62
56: ; preds = %29, %47, %19
%57 = phi i64 [ %49, %47 ], [ %6, %19 ], [ %30, %29 ]
%58 = phi i64 [ %52, %47 ], [ %17, %19 ], [ %33, %29 ]
%59 = xor i64 %58, -1
%60 = tail call i64 @llvm.cttz.i64(i64 %59, i1 true)
%61 = add i64 %60, %57
br label %62
62: ; preds = %40, %3, %56, %54
%63 = phi i64 [ %55, %54 ], [ %61, %56 ], [ %1, %3 ], [ %43, %40 ]
ret i64 %63
}