void demultiplyAlpha(int wide, int high, unsigned char *data) {
  int x, y;
  unsigned char *srcPtr;
  srcPtr = data;
  for (y = 0; y < high; y++) {
    srcPtr++;
    for (x = 0; x < 4 * wide; x += 4) {
      if (srcPtr[x + 3]) {
        srcPtr[x] = (srcPtr[x] * 255 + (srcPtr[x + 3] >> 1)) / srcPtr[x + 3];
        srcPtr[x + 1] =
            (srcPtr[x + 1] * 255 + (srcPtr[x + 3] >> 1)) / srcPtr[x + 3];
        srcPtr[x + 2] =
            (srcPtr[x + 2] * 255 + (srcPtr[x + 3] >> 1)) / srcPtr[x + 3];
      }
    }
    srcPtr += 4 * wide;
  }
}


x86 assembly (clang, -O0):
<demultiplyAlpha>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x20(%rbp)
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1288 <demultiplyAlpha+0x188>
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
mov    -0x4(%rbp),%ecx
shl    $0x2,%ecx
cmp    %ecx,%eax
jge    1266 <demultiplyAlpha+0x166>
mov    -0x20(%rbp),%rax
mov    -0x14(%rbp),%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
cmpb   $0x0,(%rax,%rcx,1)
je     1253 <demultiplyAlpha+0x153>
mov    -0x20(%rbp),%rax
movslq -0x14(%rbp),%rcx
movzbl (%rax,%rcx,1),%eax
imul   $0xff,%eax,%eax
mov    -0x20(%rbp),%rcx
mov    -0x14(%rbp),%edx
add    $0x3,%edx
movslq %edx,%rdx
movzbl (%rcx,%rdx,1),%ecx
sar    $0x1,%ecx
add    %ecx,%eax
mov    -0x20(%rbp),%rcx
mov    -0x14(%rbp),%edx
add    $0x3,%edx
movslq %edx,%rdx
movzbl (%rcx,%rdx,1),%ecx
cltd
idiv   %ecx
mov    %al,%dl
mov    -0x20(%rbp),%rax
movslq -0x14(%rbp),%rcx
mov    %dl,(%rax,%rcx,1)
mov    -0x20(%rbp),%rax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movzbl (%rax,%rcx,1),%eax
imul   $0xff,%eax,%eax
mov    -0x20(%rbp),%rcx
mov    -0x14(%rbp),%edx
add    $0x3,%edx
movslq %edx,%rdx
movzbl (%rcx,%rdx,1),%ecx
sar    $0x1,%ecx
add    %ecx,%eax
mov    -0x20(%rbp),%rcx
mov    -0x14(%rbp),%edx
add    $0x3,%edx
movslq %edx,%rdx
movzbl (%rcx,%rdx,1),%ecx
cltd
idiv   %ecx
mov    %al,%dl
mov    -0x20(%rbp),%rax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    %dl,(%rax,%rcx,1)
mov    -0x20(%rbp),%rax
mov    -0x14(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movzbl (%rax,%rcx,1),%eax
imul   $0xff,%eax,%eax
mov    -0x20(%rbp),%rcx
mov    -0x14(%rbp),%edx
add    $0x3,%edx
movslq %edx,%rdx
movzbl (%rcx,%rdx,1),%ecx
sar    $0x1,%ecx
add    %ecx,%eax
mov    -0x20(%rbp),%rcx
mov    -0x14(%rbp),%edx
add    $0x3,%edx
movslq %edx,%rdx
movzbl (%rcx,%rdx,1),%ecx
cltd
idiv   %ecx
mov    %al,%dl
mov    -0x20(%rbp),%rax
mov    -0x14(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
mov    %dl,(%rax,%rcx,1)
jmp    1258 <demultiplyAlpha+0x158>
mov    -0x14(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x14(%rbp)
jmp    113c <demultiplyAlpha+0x3c>
mov    -0x4(%rbp),%ecx
shl    $0x2,%ecx
mov    -0x20(%rbp),%rax
movslq %ecx,%rcx
add    %rcx,%rax
mov    %rax,-0x20(%rbp)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    111d <demultiplyAlpha+0x1d>
pop    %rbp
ret

x86 assembly (clang, -O2):
<demultiplyAlpha>:
push   %rbp
push   %rbx
test   %esi,%esi
jle    11a4 <demultiplyAlpha+0xa4>
mov    %rdx,%r11
lea    0x0(,%rdi,4),%eax
movslq %eax,%r8
xor    %r9d,%r9d
jmp    112f <demultiplyAlpha+0x2f>
nopl   0x0(%rax)
add    %r8,%r10
add    $0x1,%r9d
mov    %r10,%r11
cmp    %esi,%r9d
je     11a4 <demultiplyAlpha+0xa4>
lea    0x1(%r11),%r10
test   %edi,%edi
jle    1120 <demultiplyAlpha+0x20>
xor    %ecx,%ecx
jmp    1149 <demultiplyAlpha+0x49>
nopl   0x0(%rax,%rax,1)
add    $0x4,%rcx
cmp    %r8,%rcx
jge    1120 <demultiplyAlpha+0x20>
movzbl 0x4(%r11,%rcx,1),%ebx
test   %ebx,%ebx
je     1140 <demultiplyAlpha+0x40>
movzbl 0x1(%r11,%rcx,1),%edx
mov    %edx,%eax
shl    $0x8,%eax
sub    %edx,%eax
mov    %ebx,%ebp
shr    %ebp
add    %ebp,%eax
xor    %edx,%edx
div    %bx
mov    %al,0x1(%r11,%rcx,1)
movzbl 0x2(%r11,%rcx,1),%edx
mov    %edx,%eax
shl    $0x8,%eax
sub    %edx,%eax
add    %ebp,%eax
xor    %edx,%edx
div    %bx
mov    %al,0x2(%r11,%rcx,1)
movzbl 0x3(%r11,%rcx,1),%edx
mov    %edx,%eax
shl    $0x8,%eax
sub    %edx,%eax
add    %ebp,%eax
xor    %edx,%edx
div    %bx
mov    %al,0x3(%r11,%rcx,1)
jmp    1140 <demultiplyAlpha+0x40>
pop    %rbx
pop    %rbp
ret

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

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