void bitmap_set(unsigned long *map, int start, int nr) {
  unsigned long *p = map + ((start) / 64);
  const int size = start + nr;
  int bits_to_set = 64 - (start % 64);
  unsigned long mask_to_set = (~0UL << ((start) & (64 - 1)));
  while (nr - bits_to_set >= 0) {
    *p |= mask_to_set;
    nr -= bits_to_set;
    bits_to_set = 64;
    mask_to_set = ~0UL;
    p++;
  }
  if (nr) {
    mask_to_set &= (~0UL >> (-(size) & (64 - 1)));
    *p |= mask_to_set;
  }
}


x86 assembly (clang, -O0):
<bitmap_set>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0xc(%rbp),%eax
mov    $0x40,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x30(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x18(%rbp)
mov    -0xc(%rbp),%eax
add    -0x10(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0xc(%rbp),%eax
mov    $0x40,%ecx
cltd
idiv   %ecx
mov    $0x40,%eax
sub    %edx,%eax
mov    %eax,-0x20(%rbp)
mov    -0xc(%rbp),%eax
and    $0x3f,%eax
mov    %eax,%eax
mov    %eax,%ecx
mov    $0xffffffffffffffff,%rax
shl    %cl,%rax
mov    %rax,-0x28(%rbp)
mov    -0x10(%rbp),%eax
sub    -0x20(%rbp),%eax
cmp    $0x0,%eax
jl     11b3 <bitmap_set+0xb3>
mov    -0x28(%rbp),%rcx
mov    -0x18(%rbp),%rax
or     (%rax),%rcx
mov    %rcx,(%rax)
mov    -0x20(%rbp),%ecx
mov    -0x10(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x10(%rbp)
movl   $0x40,-0x20(%rbp)
movq   $0xffffffffffffffff,-0x28(%rbp)
ff
mov    -0x18(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x18(%rbp)
jmp    116b <bitmap_set+0x6b>
cmpl   $0x0,-0x10(%rbp)
je     11e9 <bitmap_set+0xe9>
xor    %eax,%eax
sub    -0x1c(%rbp),%eax
and    $0x3f,%eax
mov    %eax,%eax
mov    %eax,%ecx
mov    $0xffffffffffffffff,%rax
shr    %cl,%rax
and    -0x28(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x28(%rbp),%rcx
mov    -0x18(%rbp),%rax
or     (%rax),%rcx
mov    %rcx,(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<bitmap_set>:
push   %rbp
push   %r15
push   %r14
push   %rbx
mov    %esi,%ecx
lea    0x3f(%rcx),%esi
test   %ecx,%ecx
cmovns %ecx,%esi
mov    %esi,%eax
sar    $0x6,%eax
movslq %eax,%r9
lea    (%rdi,%r9,8),%rax
and    $0xffffffc0,%esi
mov    %ecx,%r8d
sub    %esi,%r8d
add    %edx,%r8d
mov    $0xffffffffffffffff,%rsi
shl    %cl,%rsi
mov    %r8d,%r15d
add    $0xffffffc0,%r15d
js     11a0 <bitmap_set+0xa0>
or     %rsi,(%rax)
add    $0x8,%rax
cmp    $0x40,%r15d
jl     12a7 <bitmap_set+0x1a7>
cmp    $0x7f,%r15d
mov    $0x7f,%esi
cmovl  %r15d,%esi
not    %esi
add    %esi,%r8d
cmp    $0xc0,%r8d
jb     1290 <bitmap_set+0x190>
shr    $0x6,%r8d
add    $0x1,%r8d
mov    %r8d,%r10d
and    $0xfffffffc,%r10d
mov    %r10d,%r11d
lea    -0x4(%r10),%rsi
mov    %rsi,%rbx
shr    $0x2,%rbx
add    $0x1,%rbx
mov    %ebx,%r14d
and    $0x7,%r14d
cmp    $0x1c,%rsi
jae    11b1 <bitmap_set+0xb1>
xor    %esi,%esi
jmp    124f <bitmap_set+0x14f>
mov    %edx,%r15d
test   %r15d,%r15d
jne    12b3 <bitmap_set+0x1b3>
jmp    12c0 <bitmap_set+0x1c0>
and    $0xfffffffffffffff8,%rbx
lea    (%rdi,%r9,8),%rbp
add    $0xf8,%rbp
xor    %esi,%esi
pcmpeqd %xmm0,%xmm0
cs nopw 0x0(%rax,%rax,1)
movdqu %xmm0,-0xf0(%rbp,%rsi,8)
ff ff
movdqu %xmm0,-0xe0(%rbp,%rsi,8)
ff ff
movdqu %xmm0,-0xd0(%rbp,%rsi,8)
ff ff
movdqu %xmm0,-0xc0(%rbp,%rsi,8)
ff ff
movdqu %xmm0,-0xb0(%rbp,%rsi,8)
ff ff
movdqu %xmm0,-0xa0(%rbp,%rsi,8)
ff ff
movdqu %xmm0,-0x90(%rbp,%rsi,8)
ff ff
movdqu %xmm0,-0x80(%rbp,%rsi,8)
movdqu %xmm0,-0x70(%rbp,%rsi,8)
movdqu %xmm0,-0x60(%rbp,%rsi,8)
movdqu %xmm0,-0x50(%rbp,%rsi,8)
movdqu %xmm0,-0x40(%rbp,%rsi,8)
movdqu %xmm0,-0x30(%rbp,%rsi,8)
movdqu %xmm0,-0x20(%rbp,%rsi,8)
movdqu %xmm0,-0x10(%rbp,%rsi,8)
movdqu %xmm0,0x0(%rbp,%rsi,8)
add    $0x20,%rsi
add    $0xfffffffffffffff8,%rbx
jne    11d0 <bitmap_set+0xd0>
shl    $0x6,%r11d
test   %r14,%r14
je     1284 <bitmap_set+0x184>
add    %r9,%rsi
lea    (%rdi,%rsi,8),%rsi
add    $0x18,%rsi
shl    $0x5,%r14
xor    %edi,%edi
pcmpeqd %xmm0,%xmm0
nopl   (%rax)
movdqu %xmm0,-0x10(%rsi,%rdi,1)
movdqu %xmm0,(%rsi,%rdi,1)
add    $0x20,%rdi
cmp    %rdi,%r14
jne    1270 <bitmap_set+0x170>
sub    %r11d,%r15d
lea    (%rax,%r10,8),%rax
cmp    %r8,%r10
je     12a7 <bitmap_set+0x1a7>
mov    %r15d,%esi
lea    -0x40(%rsi),%r15d
movq   $0xffffffffffffffff,(%rax)
add    $0x8,%rax
cmp    $0x7f,%esi
jg     1290 <bitmap_set+0x190>
mov    $0xffffffffffffffff,%rsi
test   %r15d,%r15d
je     12c0 <bitmap_set+0x1c0>
add    %edx,%ecx
neg    %cl
shl    %cl,%rsi
shr    %cl,%rsi
or     %rsi,(%rax)
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

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

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