void bitmap_or(unsigned char *dst, int dst_wrap, unsigned char *src,
               int src_wrap, int x, int y, int w, int h) {
  int i, j, bit, x1;
  dst += dst_wrap * y;
  for (i = 0; i < h; i++) {
    for (j = 0; j < w; j++) {
      bit = (src[j >> 3] >> (7 - (j & 7))) & 1;
      x1 = x + j;
      dst[x1 >> 3] |= bit << (7 - (x1 & 7));
    }
    dst += dst_wrap;
    src += src_wrap;
  }
}


x86 assembly (clang, -O0):
<bitmap_or>:
push   %rbp
mov    %rsp,%rbp
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    %r9d,-0x24(%rbp)
mov    -0xc(%rbp),%ecx
imul   -0x24(%rbp),%ecx
mov    -0x8(%rbp),%rax
movslq %ecx,%rcx
add    %rcx,%rax
mov    %rax,-0x8(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    11f0 <bitmap_or+0xf0>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    11c0 <bitmap_or+0xc0>
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sar    $0x3,%ecx
movslq %ecx,%rcx
movzbl (%rax,%rcx,1),%eax
mov    -0x2c(%rbp),%edx
and    $0x7,%edx
mov    $0x7,%ecx
sub    %edx,%ecx
sar    %cl,%eax
and    $0x1,%eax
mov    %eax,-0x30(%rbp)
mov    -0x20(%rbp),%eax
add    -0x2c(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x30(%rbp),%esi
mov    -0x34(%rbp),%eax
and    $0x7,%eax
mov    $0x7,%ecx
sub    %eax,%ecx
shl    %cl,%esi
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
sar    $0x3,%ecx
movslq %ecx,%rcx
movzbl (%rax,%rcx,1),%edx
or     %esi,%edx
mov    %dl,(%rax,%rcx,1)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    114f <bitmap_or+0x4f>
mov    -0xc(%rbp),%ecx
mov    -0x8(%rbp),%rax
movslq %ecx,%rcx
add    %rcx,%rax
mov    %rax,-0x8(%rbp)
mov    -0x1c(%rbp),%ecx
mov    -0x18(%rbp),%rax
movslq %ecx,%rcx
add    %rcx,%rax
mov    %rax,-0x18(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    113c <bitmap_or+0x3c>
pop    %rbp
ret

x86 assembly (clang, -O2):
<bitmap_or>:
push   %rbp
push   %r14
push   %rbx
mov    0x28(%rsp),%r10d
test   %r10d,%r10d
jle    1186 <bitmap_or+0x86>
mov    0x20(%rsp),%eax
imul   %esi,%r9d
movslq %r9d,%rbx
add    %rbx,%rdi
movslq %esi,%r9
movslq %ecx,%r11
xor    %r14d,%r14d
jmp    113f <bitmap_or+0x3f>
nopw   0x0(%rax,%rax,1)
add    %r9,%rdi
add    %r11,%rdx
add    $0x1,%r14d
cmp    %r10d,%r14d
je     1186 <bitmap_or+0x86>
test   %eax,%eax
jle    1130 <bitmap_or+0x30>
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nop
mov    %ebx,%ecx
shr    $0x3,%rcx
movzbl (%rdx,%rcx,1),%ebp
mov    %ebx,%ecx
not    %cl
and    $0x7,%cl
shr    %cl,%ebp
and    $0x1,%ebp
lea    (%r8,%rbx,1),%esi
mov    %esi,%ecx
not    %cl
and    $0x7,%cl
shl    %cl,%ebp
sar    $0x3,%esi
movslq %esi,%rcx
or     %bpl,(%rdi,%rcx,1)
add    $0x1,%ebx
cmp    %ebx,%eax
jne    1150 <bitmap_or+0x50>
jmp    1130 <bitmap_or+0x30>
pop    %rbx
pop    %r14
pop    %rbp
ret

LLVM IR (-O0):
define void @bitmap_or(i8* %0, i32 %1, i8* %2, i32 %3, i32 %4, i32 %5, i32 %6, i32 %7) {
%9 = alloca i8*
%10 = alloca i32
%11 = alloca i8*
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
store i8* %0, i8** %9
store i32 %1, i32* %10
store i8* %2, i8** %11
store i32 %3, i32* %12
store i32 %4, i32* %13
store i32 %5, i32* %14
store i32 %6, i32* %15
store i32 %7, i32* %16
%21 = load i32, i32* %10
%22 = load i32, i32* %14
%23 = mul nsw i32 %21, %22
%24 = load i8*, i8** %9
%25 = sext i32 %23 to i64
%26 = getelementptr inbounds i8, i8* %24, i64 %25
store i8* %26, i8** %9
store i32 0, i32* %17
br label %27
27: ; preds = %78, %8
%28 = load i32, i32* %17
%29 = load i32, i32* %16
%30 = icmp slt i32 %28, %29
br i1 %30, label %31, label %81
31: ; preds = %27
store i32 0, i32* %18
br label %32
32: ; preds = %66, %31
%33 = load i32, i32* %18
%34 = load i32, i32* %15
%35 = icmp slt i32 %33, %34
br i1 %35, label %36, label %69
36: ; preds = %32
%37 = load i8*, i8** %11
%38 = load i32, i32* %18
%39 = ashr i32 %38, 3
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds i8, i8* %37, i64 %40
%42 = load i8, i8* %41
%43 = zext i8 %42 to i32
%44 = load i32, i32* %18
%45 = and i32 %44, 7
%46 = sub nsw i32 7, %45
%47 = ashr i32 %43, %46
%48 = and i32 %47, 1
store i32 %48, i32* %19
%49 = load i32, i32* %13
%50 = load i32, i32* %18
%51 = add nsw i32 %49, %50
store i32 %51, i32* %20
%52 = load i32, i32* %19
%53 = load i32, i32* %20
%54 = and i32 %53, 7
%55 = sub nsw i32 7, %54
%56 = shl i32 %52, %55
%57 = load i8*, i8** %9
%58 = load i32, i32* %20
%59 = ashr i32 %58, 3
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds i8, i8* %57, i64 %60
%62 = load i8, i8* %61
%63 = zext i8 %62 to i32
%64 = or i32 %63, %56
%65 = trunc i32 %64 to i8
store i8 %65, i8* %61
br label %66
66: ; preds = %36
%67 = load i32, i32* %18
%68 = add nsw i32 %67, 1
store i32 %68, i32* %18
br label %32
69: ; preds = %32
%70 = load i32, i32* %10
%71 = load i8*, i8** %9
%72 = sext i32 %70 to i64
%73 = getelementptr inbounds i8, i8* %71, i64 %72
store i8* %73, i8** %9
%74 = load i32, i32* %12
%75 = load i8*, i8** %11
%76 = sext i32 %74 to i64
%77 = getelementptr inbounds i8, i8* %75, i64 %76
store i8* %77, i8** %11
br label %78
78: ; preds = %69
%79 = load i32, i32* %17
%80 = add nsw i32 %79, 1
store i32 %80, i32* %17
br label %27
81: ; preds = %27
ret void
}

LLVM IR (-O2):
define void @bitmap_or(i8* nocapture %0, i32 %1, i8* nocapture %2, i32 %3, i32 %4, i32 %5, i32 %6, i32 %7) {
%9 = icmp sgt i32 %7, 0
br i1 %9, label %10, label %49
10: ; preds = %8
%11 = mul nsw i32 %5, %1
%12 = sext i32 %11 to i64
%13 = getelementptr inbounds i8, i8* %0, i64 %12
%14 = icmp sgt i32 %6, 0
%15 = sext i32 %1 to i64
%16 = sext i32 %3 to i64
br label %17
17: ; preds = %10, %44
%18 = phi i8* [ %13, %10 ], [ %45, %44 ]
%19 = phi i8* [ %2, %10 ], [ %46, %44 ]
%20 = phi i32 [ 0, %10 ], [ %47, %44 ]
br i1 %14, label %21, label %44
21: ; preds = %17, %21
%22 = phi i32 [ %42, %21 ], [ 0, %17 ]
%23 = lshr i32 %22, 3
%24 = zext i32 %23 to i64
%25 = getelementptr inbounds i8, i8* %19, i64 %24
%26 = load i8, i8* %25
%27 = zext i8 %26 to i32
%28 = and i32 %22, 7
%29 = xor i32 %28, 7
%30 = lshr i32 %27, %29
%31 = and i32 %30, 1
%32 = add nsw i32 %22, %4
%33 = and i32 %32, 7
%34 = xor i32 %33, 7
%35 = shl nuw i32 %31, %34
%36 = ashr i32 %32, 3
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds i8, i8* %18, i64 %37
%39 = load i8, i8* %38
%40 = trunc i32 %35 to i8
%41 = or i8 %39, %40
store i8 %41, i8* %38
%42 = add nuw nsw i32 %22, 1
%43 = icmp eq i32 %42, %6
br i1 %43, label %44, label %21
44: ; preds = %21, %17
%45 = getelementptr inbounds i8, i8* %18, i64 %15
%46 = getelementptr inbounds i8, i8* %19, i64 %16
%47 = add nuw nsw i32 %20, 1
%48 = icmp eq i32 %47, %7
br i1 %48, label %49, label %17
49: ; preds = %44, %8
ret void
}