void compositeSourceOut(unsigned char *dst, unsigned char *src, int width,
                        int height) {
  unsigned char *end = src + ((width * height) << 2);
  unsigned int *dst32 = (unsigned int *)dst;
  int d, s;
  float da_inv;
  while (src < end) {
    d = *(dst + 3);
    if (d == 0) {
      *dst32++ = *(unsigned int *)src;
    } else {
      s = *(src + 3);
      if (s == 0) {
        *dst32++ = 0;
      } else {
        da_inv = 1 - d / 255.0;
        *dst32++ = (int)(*src * da_inv) | (int)(*(src + 1) * da_inv) << 8 |
                   (int)(*(src + 2) * da_inv) << 16 | (int)(s * da_inv) << 24;
      }
    }
    src += 4;
    dst += 4;
  }
}


x86 assembly (clang, -O0):
<compositeSourceOut>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
shl    $0x2,%ecx
movslq %ecx,%rcx
add    %rcx,%rax
mov    %rax,-0x20(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x10(%rbp),%rax
cmp    -0x20(%rbp),%rax
jae    1254 <compositeSourceOut+0x154>
mov    -0x8(%rbp),%rax
movzbl 0x3(%rax),%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x0,-0x2c(%rbp)
jne    1171 <compositeSourceOut+0x71>
mov    -0x10(%rbp),%rax
mov    (%rax),%ecx
mov    -0x28(%rbp),%rax
mov    %rax,%rdx
add    $0x4,%rdx
mov    %rdx,-0x28(%rbp)
mov    %ecx,(%rax)
jmp    1237 <compositeSourceOut+0x137>
mov    -0x10(%rbp),%rax
movzbl 0x3(%rax),%eax
mov    %eax,-0x30(%rbp)
cmpl   $0x0,-0x30(%rbp)
jne    11a0 <compositeSourceOut+0xa0>
mov    -0x28(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x28(%rbp)
movl   $0x0,(%rax)
jmp    1232 <compositeSourceOut+0x132>
cvtsi2sdl -0x2c(%rbp),%xmm1
movsd  0xe5b(%rip),%xmm0
divsd  %xmm0,%xmm1
movsd  0xe47(%rip),%xmm0
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x10(%rbp),%rax
movzbl (%rax),%eax
cvtsi2ss %eax,%xmm0
mulss  -0x34(%rbp),%xmm0
cvttss2si %xmm0,%ecx
mov    -0x10(%rbp),%rax
movzbl 0x1(%rax),%eax
cvtsi2ss %eax,%xmm0
mulss  -0x34(%rbp),%xmm0
cvttss2si %xmm0,%eax
shl    $0x8,%eax
or     %eax,%ecx
mov    -0x10(%rbp),%rax
movzbl 0x2(%rax),%eax
cvtsi2ss %eax,%xmm0
mulss  -0x34(%rbp),%xmm0
cvttss2si %xmm0,%eax
shl    $0x10,%eax
or     %eax,%ecx
cvtsi2ssl -0x30(%rbp),%xmm0
mulss  -0x34(%rbp),%xmm0
cvttss2si %xmm0,%eax
shl    $0x18,%eax
or     %eax,%ecx
mov    -0x28(%rbp),%rax
mov    %rax,%rdx
add    $0x4,%rdx
mov    %rdx,-0x28(%rbp)
mov    %ecx,(%rax)
jmp    1237 <compositeSourceOut+0x137>
mov    -0x10(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x10(%rbp)
mov    -0x8(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x8(%rbp)
jmp    1132 <compositeSourceOut+0x32>
pop    %rbp
ret

x86 assembly (clang, -O2):
<compositeSourceOut>:
imul   %ecx,%edx
shl    $0x2,%edx
test   %edx,%edx
jle    11e7 <compositeSourceOut+0xe7>
movslq %edx,%r8
add    %rsi,%r8
movsd  0xee4(%rip),%xmm0
movsd  0xee4(%rip),%xmm1
pxor   %xmm2,%xmm2
jmp    1145 <compositeSourceOut+0x45>
nopw   0x0(%rax,%rax,1)
mov    (%rsi),%ecx
mov    %ecx,(%rdi)
add    $0x4,%rdi
add    $0x4,%rsi
cmp    %r8,%rsi
jae    11e7 <compositeSourceOut+0xe7>
movzbl 0x3(%rdi),%edx
test   %dl,%dl
je     1130 <compositeSourceOut+0x30>
movzbl 0x3(%rsi),%ecx
test   %cl,%cl
je     11e0 <compositeSourceOut+0xe0>
movzbl %dl,%edx
xorps  %xmm3,%xmm3
cvtsi2sd %edx,%xmm3
divsd  %xmm0,%xmm3
addsd  %xmm1,%xmm3
cvtsd2ss %xmm3,%xmm3
movzbl (%rsi),%edx
xorps  %xmm4,%xmm4
cvtsi2ss %edx,%xmm4
mulss  %xmm3,%xmm4
cvttss2si %xmm4,%edx
movzwl 0x1(%rsi),%eax
movd   %eax,%xmm4
punpcklbw %xmm2,%xmm4
punpcklwd %xmm2,%xmm4
movzbl %cl,%eax
xorps  %xmm5,%xmm5
cvtsi2ss %eax,%xmm5
cvtdq2ps %xmm4,%xmm4
mulss  %xmm3,%xmm5
shufps $0x0,%xmm3,%xmm3
mulps  %xmm4,%xmm3
cvttps2dq %xmm3,%xmm3
movdqa %xmm3,%xmm4
pslld  $0x10,%xmm4
cvttss2si %xmm5,%eax
pslld  $0x8,%xmm3
shl    $0x18,%eax
or     %edx,%eax
movd   %xmm3,%edx
or     %eax,%edx
pshufd $0x55,%xmm4,%xmm3
movd   %xmm3,%ecx
or     %edx,%ecx
jmp    1132 <compositeSourceOut+0x32>
nopw   0x0(%rax,%rax,1)
xor    %ecx,%ecx
jmp    1132 <compositeSourceOut+0x32>
ret

LLVM IR (-O0):
define void @compositeSourceOut(i8* %0, i8* %1, i32 %2, i32 %3) {
%5 = alloca i8*
%6 = alloca i8*
%7 = alloca i32
%8 = alloca i32
%9 = alloca i8*
%10 = alloca i32*
%11 = alloca i32
%12 = alloca i32
%13 = alloca float
store i8* %0, i8** %5
store i8* %1, i8** %6
store i32 %2, i32* %7
store i32 %3, i32* %8
%14 = load i8*, i8** %6
%15 = load i32, i32* %7
%16 = load i32, i32* %8
%17 = mul nsw i32 %15, %16
%18 = shl i32 %17, 2
%19 = sext i32 %18 to i64
%20 = getelementptr inbounds i8, i8* %14, i64 %19
store i8* %20, i8** %9
%21 = load i8*, i8** %5
%22 = bitcast i8* %21 to i32*
store i32* %22, i32** %10
br label %23
23: ; preds = %93, %4
%24 = load i8*, i8** %6
%25 = load i8*, i8** %9
%26 = icmp ult i8* %24, %25
br i1 %26, label %27, label %98
27: ; preds = %23
%28 = load i8*, i8** %5
%29 = getelementptr inbounds i8, i8* %28, i64 3
%30 = load i8, i8* %29
%31 = zext i8 %30 to i32
store i32 %31, i32* %11
%32 = load i32, i32* %11
%33 = icmp eq i32 %32, 0
br i1 %33, label %34, label %40
34: ; preds = %27
%35 = load i8*, i8** %6
%36 = bitcast i8* %35 to i32*
%37 = load i32, i32* %36
%38 = load i32*, i32** %10
%39 = getelementptr inbounds i32, i32* %38, i32 1
store i32* %39, i32** %10
store i32 %37, i32* %38
br label %93
40: ; preds = %27
%41 = load i8*, i8** %6
%42 = getelementptr inbounds i8, i8* %41, i64 3
%43 = load i8, i8* %42
%44 = zext i8 %43 to i32
store i32 %44, i32* %12
%45 = load i32, i32* %12
%46 = icmp eq i32 %45, 0
br i1 %46, label %47, label %50
47: ; preds = %40
%48 = load i32*, i32** %10
%49 = getelementptr inbounds i32, i32* %48, i32 1
store i32* %49, i32** %10
store i32 0, i32* %48
br label %92
50: ; preds = %40
%51 = load i32, i32* %11
%52 = sitofp i32 %51 to double
%53 = fdiv double %52, 2.550000e+02
%54 = fsub double 1.000000e+00, %53
%55 = fptrunc double %54 to float
store float %55, float* %13
%56 = load i8*, i8** %6
%57 = load i8, i8* %56
%58 = zext i8 %57 to i32
%59 = sitofp i32 %58 to float
%60 = load float, float* %13
%61 = fmul float %59, %60
%62 = fptosi float %61 to i32
%63 = load i8*, i8** %6
%64 = getelementptr inbounds i8, i8* %63, i64 1
%65 = load i8, i8* %64
%66 = zext i8 %65 to i32
%67 = sitofp i32 %66 to float
%68 = load float, float* %13
%69 = fmul float %67, %68
%70 = fptosi float %69 to i32
%71 = shl i32 %70, 8
%72 = or i32 %62, %71
%73 = load i8*, i8** %6
%74 = getelementptr inbounds i8, i8* %73, i64 2
%75 = load i8, i8* %74
%76 = zext i8 %75 to i32
%77 = sitofp i32 %76 to float
%78 = load float, float* %13
%79 = fmul float %77, %78
%80 = fptosi float %79 to i32
%81 = shl i32 %80, 16
%82 = or i32 %72, %81
%83 = load i32, i32* %12
%84 = sitofp i32 %83 to float
%85 = load float, float* %13
%86 = fmul float %84, %85
%87 = fptosi float %86 to i32
%88 = shl i32 %87, 24
%89 = or i32 %82, %88
%90 = load i32*, i32** %10
%91 = getelementptr inbounds i32, i32* %90, i32 1
store i32* %91, i32** %10
store i32 %89, i32* %90
br label %92
92: ; preds = %50, %47
br label %93
93: ; preds = %92, %34
%94 = load i8*, i8** %6
%95 = getelementptr inbounds i8, i8* %94, i64 4
store i8* %95, i8** %6
%96 = load i8*, i8** %5
%97 = getelementptr inbounds i8, i8* %96, i64 4
store i8* %97, i8** %5
br label %23
98: ; preds = %23
ret void
}

LLVM IR (-O2):
define void @compositeSourceOut(i8* nocapture %0, i8* %1, i32 %2, i32 %3) {
%5 = shl i32 %2, 2
%6 = mul i32 %5, %3
%7 = sext i32 %6 to i64
%8 = getelementptr inbounds i8, i8* %1, i64 %7
%9 = icmp sgt i32 %6, 0
br i1 %9, label %10, label %56
10: ; preds = %4, %51
%11 = phi i8* [ %53, %51 ], [ %0, %4 ]
%12 = phi i8* [ %54, %51 ], [ %1, %4 ]
%13 = bitcast i8* %11 to i32*
%14 = getelementptr inbounds i8, i8* %11, i64 3
%15 = load i8, i8* %14
%16 = icmp eq i8 %15, 0
br i1 %16, label %17, label %20
17: ; preds = %10
%18 = bitcast i8* %12 to i32*
%19 = load i32, i32* %18
br label %51
20: ; preds = %10
%21 = getelementptr inbounds i8, i8* %12, i64 3
%22 = load i8, i8* %21
%23 = icmp eq i8 %22, 0
br i1 %23, label %51, label %24
24: ; preds = %20
%25 = uitofp i8 %15 to double
%26 = fdiv double %25, 2.550000e+02
%27 = fsub double 1.000000e+00, %26
%28 = fptrunc double %27 to float
%29 = load i8, i8* %12
%30 = uitofp i8 %29 to float
%31 = fmul float %28, %30
%32 = fptosi float %31 to i32
%33 = getelementptr inbounds i8, i8* %12, i64 1
%34 = bitcast i8* %33 to <2 x i8>*
%35 = load <2 x i8>, <2 x i8>* %34
%36 = uitofp <2 x i8> %35 to <2 x float>
%37 = insertelement <2 x float> poison, float %28, i64 0
%38 = shufflevector <2 x float> %37, <2 x float> poison, <2 x i32> zeroinitializer
%39 = fmul <2 x float> %38, %36
%40 = fptosi <2 x float> %39 to <2 x i32>
%41 = shl <2 x i32> %40, <i32 8, i32 16>
%42 = uitofp i8 %22 to float
%43 = fmul float %28, %42
%44 = fptosi float %43 to i32
%45 = shl i32 %44, 24
%46 = or i32 %45, %32
%47 = extractelement <2 x i32> %41, i64 0
%48 = or i32 %46, %47
%49 = extractelement <2 x i32> %41, i64 1
%50 = or i32 %48, %49
br label %51
51: ; preds = %20, %24, %17
%52 = phi i32 [ %50, %24 ], [ %19, %17 ], [ 0, %20 ]
store i32 %52, i32* %13
%53 = getelementptr inbounds i8, i8* %11, i64 4
%54 = getelementptr inbounds i8, i8* %12, i64 4
%55 = icmp ult i8* %54, %8
br i1 %55, label %10, label %56
56: ; preds = %51, %4
ret void
}