void adv_png_expand_1(unsigned width, unsigned height, unsigned char *ptr) {
  unsigned i, j;
  unsigned char *p8 = ptr + height * (width + 1) - 1;
  unsigned char *p1 = ptr + height * (width / 8 + 1) - 1;
  width /= 8;
  for (i = 0; i < height; ++i) {
    for (j = 0; j < width; ++j) {
      unsigned char v = *p1;
      *p8-- = v & 0x1;
      *p8-- = (v >> 1) & 0x1;
      *p8-- = (v >> 2) & 0x1;
      *p8-- = (v >> 3) & 0x1;
      *p8-- = (v >> 4) & 0x1;
      *p8-- = (v >> 5) & 0x1;
      *p8-- = (v >> 6) & 0x1;
      *p8-- = v >> 7;
      --p1;
    }
    --p8;
    --p1;
  }
}


x86 assembly (clang, -O0):
<adv_png_expand_1>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x8(%rbp),%ecx
mov    -0x4(%rbp),%edx
add    $0x1,%edx
imul   %edx,%ecx
mov    %ecx,%ecx
add    %rcx,%rax
add    $0xffffffffffffffff,%rax
mov    %rax,-0x20(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x8(%rbp),%ecx
mov    -0x4(%rbp),%edx
shr    $0x3,%edx
add    $0x1,%edx
imul   %edx,%ecx
mov    %ecx,%ecx
add    %rcx,%rax
add    $0xffffffffffffffff,%rax
mov    %rax,-0x28(%rbp)
mov    -0x4(%rbp),%eax
shr    $0x3,%eax
mov    %eax,-0x4(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x8(%rbp),%eax
jae    12a5 <adv_png_expand_1+0x1a5>
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x4(%rbp),%eax
jae    127f <adv_png_expand_1+0x17f>
mov    -0x28(%rbp),%rax
mov    (%rax),%al
mov    %al,-0x29(%rbp)
movzbl -0x29(%rbp),%eax
and    $0x1,%eax
mov    %al,%cl
mov    -0x20(%rbp),%rax
mov    %rax,%rdx
add    $0xffffffffffffffff,%rdx
mov    %rdx,-0x20(%rbp)
mov    %cl,(%rax)
movzbl -0x29(%rbp),%eax
sar    $0x1,%eax
and    $0x1,%eax
mov    %al,%cl
mov    -0x20(%rbp),%rax
mov    %rax,%rdx
add    $0xffffffffffffffff,%rdx
mov    %rdx,-0x20(%rbp)
mov    %cl,(%rax)
movzbl -0x29(%rbp),%eax
sar    $0x2,%eax
and    $0x1,%eax
mov    %al,%cl
mov    -0x20(%rbp),%rax
mov    %rax,%rdx
add    $0xffffffffffffffff,%rdx
mov    %rdx,-0x20(%rbp)
mov    %cl,(%rax)
movzbl -0x29(%rbp),%eax
sar    $0x3,%eax
and    $0x1,%eax
mov    %al,%cl
mov    -0x20(%rbp),%rax
mov    %rax,%rdx
add    $0xffffffffffffffff,%rdx
mov    %rdx,-0x20(%rbp)
mov    %cl,(%rax)
movzbl -0x29(%rbp),%eax
sar    $0x4,%eax
and    $0x1,%eax
mov    %al,%cl
mov    -0x20(%rbp),%rax
mov    %rax,%rdx
add    $0xffffffffffffffff,%rdx
mov    %rdx,-0x20(%rbp)
mov    %cl,(%rax)
movzbl -0x29(%rbp),%eax
sar    $0x5,%eax
and    $0x1,%eax
mov    %al,%cl
mov    -0x20(%rbp),%rax
mov    %rax,%rdx
add    $0xffffffffffffffff,%rdx
mov    %rdx,-0x20(%rbp)
mov    %cl,(%rax)
movzbl -0x29(%rbp),%eax
sar    $0x6,%eax
and    $0x1,%eax
mov    %al,%cl
mov    -0x20(%rbp),%rax
mov    %rax,%rdx
add    $0xffffffffffffffff,%rdx
mov    %rdx,-0x20(%rbp)
mov    %cl,(%rax)
movzbl -0x29(%rbp),%eax
sar    $0x7,%eax
mov    %al,%cl
mov    -0x20(%rbp),%rax
mov    %rax,%rdx
add    $0xffffffffffffffff,%rdx
mov    %rdx,-0x20(%rbp)
mov    %cl,(%rax)
mov    -0x28(%rbp),%rax
add    $0xffffffffffffffff,%rax
mov    %rax,-0x28(%rbp)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    116e <adv_png_expand_1+0x6e>
mov    -0x20(%rbp),%rax
add    $0xffffffffffffffff,%rax
mov    %rax,-0x20(%rbp)
mov    -0x28(%rbp),%rax
add    $0xffffffffffffffff,%rax
mov    %rax,-0x28(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    115b <adv_png_expand_1+0x5b>
pop    %rbp
ret

x86 assembly (clang, -O2):
<adv_png_expand_1>:
test   %esi,%esi
je     11bb <adv_png_expand_1+0xbb>
mov    %edi,%eax
shr    $0x3,%eax
lea    0x1(%rax),%r10d
imul   %esi,%r10d
add    %rdx,%r10
lea    0x1(%rdi),%ecx
imul   %esi,%ecx
add    %rdx,%rcx
cmp    $0x2,%eax
mov    $0x1,%r8d
cmovae %eax,%r8d
xor    %r9d,%r9d
jmp    1149 <adv_png_expand_1+0x49>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%r9d
cmp    %esi,%r9d
je     11bb <adv_png_expand_1+0xbb>
add    $0xffffffffffffffff,%rcx
add    $0xffffffffffffffff,%r10
cmp    $0x8,%edi
jb     1140 <adv_png_expand_1+0x40>
mov    %r8d,%r11d
nopl   0x0(%rax)
movzbl (%r10),%edx
mov    %edx,%eax
and    $0x1,%al
mov    %al,(%rcx)
mov    %edx,%eax
shr    %al
and    $0x1,%al
mov    %al,-0x1(%rcx)
mov    %edx,%eax
shr    $0x2,%al
and    $0x1,%al
mov    %al,-0x2(%rcx)
mov    %edx,%eax
shr    $0x3,%al
and    $0x1,%al
mov    %al,-0x3(%rcx)
mov    %edx,%eax
shr    $0x4,%al
and    $0x1,%al
mov    %al,-0x4(%rcx)
mov    %edx,%eax
shr    $0x5,%al
and    $0x1,%al
mov    %al,-0x5(%rcx)
mov    %edx,%eax
shr    $0x6,%al
and    $0x1,%al
mov    %al,-0x6(%rcx)
shr    $0x7,%dl
mov    %dl,-0x7(%rcx)
add    $0xfffffffffffffff8,%rcx
add    $0xffffffffffffffff,%r10
add    $0xffffffff,%r11d
jne    1160 <adv_png_expand_1+0x60>
jmp    1140 <adv_png_expand_1+0x40>
ret

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

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