unsigned short blend_hardlight16(unsigned short color1, unsigned short color2) {
  unsigned r1 = (color1 >> 11), r2 = (color2 >> 11);
  unsigned g1 = ((color1 & 0x7E0) >> 5), g2 = ((color2 & 0x7E0) >> 5);
  unsigned b1 = (color1 & 0x1F), b2 = (color2 & 0x1F);
  return (((((r1) < (0x10) ? ((((r1) << 1)) * ((r2)) / ((0x1F)))
                           : (((((((r1) - (0x10)) << 1)) ^ ((0x1F))) *
                               (((r2)) ^ ((0x1F))) / ((0x1F))) ^
                              ((0x1F)))))
           << 11) |
          ((((g1) < (0x20) ? ((((g1) << 1)) * ((g2)) / ((0x3F)))
                           : (((((((g1) - (0x20)) << 1)) ^ ((0x3F))) *
                               (((g2)) ^ ((0x3F))) / ((0x3F))) ^
                              ((0x3F)))))
           << 5) |
          ((b1) < (0x10) ? ((((b1) << 1)) * ((b2)) / ((0x1F)))
                         : (((((((b1) - (0x10)) << 1)) ^ ((0x1F))) *
                             (((b2)) ^ ((0x1F))) / ((0x1F))) ^
                            ((0x1F)))));
}


x86 assembly (clang, -O0):
<blend_hardlight16>:
push   %rbp
mov    %rsp,%rbp
mov    %si,%ax
mov    %di,%cx
mov    %cx,-0x2(%rbp)
mov    %ax,-0x4(%rbp)
movzwl -0x2(%rbp),%eax
sar    $0xb,%eax
mov    %eax,-0x8(%rbp)
movzwl -0x4(%rbp),%eax
sar    $0xb,%eax
mov    %eax,-0xc(%rbp)
movzwl -0x2(%rbp),%eax
and    $0x7e0,%eax
sar    $0x5,%eax
mov    %eax,-0x10(%rbp)
movzwl -0x4(%rbp),%eax
and    $0x7e0,%eax
sar    $0x5,%eax
mov    %eax,-0x14(%rbp)
movzwl -0x2(%rbp),%eax
and    $0x1f,%eax
mov    %eax,-0x18(%rbp)
movzwl -0x4(%rbp),%eax
and    $0x1f,%eax
mov    %eax,-0x1c(%rbp)
cmpl   $0x10,-0x8(%rbp)
jae    117d <blend_hardlight16+0x7d>
mov    -0x8(%rbp),%eax
shl    $0x1,%eax
imul   -0xc(%rbp),%eax
mov    $0x1f,%ecx
xor    %edx,%edx
div    %ecx
mov    %eax,-0x20(%rbp)
jmp    11a1 <blend_hardlight16+0xa1>
mov    -0x8(%rbp),%eax
sub    $0x10,%eax
shl    $0x1,%eax
xor    $0x1f,%eax
mov    -0xc(%rbp),%ecx
xor    $0x1f,%ecx
imul   %ecx,%eax
mov    $0x1f,%ecx
xor    %edx,%edx
div    %ecx
xor    $0x1f,%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
shl    $0xb,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x20,-0x10(%rbp)
jae    11cf <blend_hardlight16+0xcf>
mov    -0x10(%rbp),%eax
shl    $0x1,%eax
imul   -0x14(%rbp),%eax
mov    $0x3f,%ecx
xor    %edx,%edx
div    %ecx
mov    %eax,-0x28(%rbp)
jmp    11f3 <blend_hardlight16+0xf3>
mov    -0x10(%rbp),%eax
sub    $0x20,%eax
shl    $0x1,%eax
xor    $0x3f,%eax
mov    -0x14(%rbp),%ecx
xor    $0x3f,%ecx
imul   %ecx,%eax
mov    $0x3f,%ecx
xor    %edx,%edx
div    %ecx
xor    $0x3f,%eax
mov    %eax,-0x28(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x28(%rbp),%ecx
shl    $0x5,%ecx
or     %ecx,%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x10,-0x18(%rbp)
jae    1226 <blend_hardlight16+0x126>
mov    -0x18(%rbp),%eax
shl    $0x1,%eax
imul   -0x1c(%rbp),%eax
mov    $0x1f,%ecx
xor    %edx,%edx
div    %ecx
mov    %eax,-0x30(%rbp)
jmp    124a <blend_hardlight16+0x14a>
mov    -0x18(%rbp),%eax
sub    $0x10,%eax
shl    $0x1,%eax
xor    $0x1f,%eax
mov    -0x1c(%rbp),%ecx
xor    $0x1f,%ecx
imul   %ecx,%eax
mov    $0x1f,%ecx
xor    %edx,%edx
div    %ecx
xor    $0x1f,%eax
mov    %eax,-0x30(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x30(%rbp),%ecx
or     %ecx,%eax
movzwl %ax,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<blend_hardlight16>:
mov    %esi,%ecx
shr    $0xb,%ecx
mov    %edi,%r8d
shr    $0x5,%r8d
and    $0x3f,%r8d
mov    %esi,%eax
shr    $0x5,%eax
mov    %edi,%edx
shr    $0xa,%edx
and    $0xfffffffe,%edx
test   %di,%di
js     1140 <blend_hardlight16+0x40>
imul   %ecx,%edx
imul   $0x843,%edx,%r9d
shr    $0x10,%r9d
sub    %r9d,%edx
movzwl %dx,%ecx
shr    %ecx
add    %r9d,%ecx
shr    $0x4,%ecx
jmp    1163 <blend_hardlight16+0x63>
add    $0xffffffe0,%edx
xor    $0x1f,%edx
xor    $0x1f,%ecx
imul   %edx,%ecx
imul   $0x8421085,%rcx,%rdx
shr    $0x20,%rdx
sub    %edx,%ecx
shr    %ecx
add    %edx,%ecx
shr    $0x4,%ecx
xor    $0x1f,%ecx
and    $0x3f,%eax
and    $0x1f,%edi
lea    (%r8,%r8,1),%edx
cmp    $0x1f,%r8d
ja     1191 <blend_hardlight16+0x91>
imul   %eax,%edx
imul   $0x411,%edx,%r8d
shr    $0x10,%r8d
sub    %r8d,%edx
movzwl %dx,%eax
shr    %eax
add    %r8d,%eax
shr    $0x5,%eax
jmp    11b4 <blend_hardlight16+0xb4>
add    $0xffffffc0,%edx
xor    $0x3f,%edx
xor    $0x3f,%eax
imul   %edx,%eax
imul   $0x4104105,%rax,%rdx
shr    $0x20,%rdx
sub    %edx,%eax
shr    %eax
add    %edx,%eax
shr    $0x5,%eax
xor    $0x3f,%eax
and    $0x1f,%esi
lea    (%rdi,%rdi,1),%edx
cmp    $0xf,%edi
ja     11d9 <blend_hardlight16+0xd9>
imul   %esi,%edx
imul   $0x843,%edx,%edi
shr    $0x10,%edi
sub    %edi,%edx
movzwl %dx,%esi
shr    %esi
add    %edi,%esi
shr    $0x4,%esi
jmp    11fc <blend_hardlight16+0xfc>
add    $0xffffffe0,%edx
xor    $0x1f,%edx
xor    $0x1f,%esi
imul   %edx,%esi
imul   $0x8421085,%rsi,%rdx
shr    $0x20,%rdx
sub    %edx,%esi
shr    %esi
add    %edx,%esi
shr    $0x4,%esi
xor    $0x1f,%esi
shl    $0xb,%ecx
or     %esi,%ecx
shl    $0x5,%eax
or     %ecx,%eax
ret

LLVM IR (-O0):
define zeroext i16 @blend_hardlight16(i16 zeroext %0, i16 zeroext %1) {
%3 = alloca i16
%4 = alloca i16
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
store i16 %0, i16* %3
store i16 %1, i16* %4
%11 = load i16, i16* %3
%12 = zext i16 %11 to i32
%13 = ashr i32 %12, 11
store i32 %13, i32* %5
%14 = load i16, i16* %4
%15 = zext i16 %14 to i32
%16 = ashr i32 %15, 11
store i32 %16, i32* %6
%17 = load i16, i16* %3
%18 = zext i16 %17 to i32
%19 = and i32 %18, 2016
%20 = ashr i32 %19, 5
store i32 %20, i32* %7
%21 = load i16, i16* %4
%22 = zext i16 %21 to i32
%23 = and i32 %22, 2016
%24 = ashr i32 %23, 5
store i32 %24, i32* %8
%25 = load i16, i16* %3
%26 = zext i16 %25 to i32
%27 = and i32 %26, 31
store i32 %27, i32* %9
%28 = load i16, i16* %4
%29 = zext i16 %28 to i32
%30 = and i32 %29, 31
store i32 %30, i32* %10
%31 = load i32, i32* %5
%32 = icmp ult i32 %31, 16
br i1 %32, label %33, label %39
33: ; preds = %2
%34 = load i32, i32* %5
%35 = shl i32 %34, 1
%36 = load i32, i32* %6
%37 = mul i32 %35, %36
%38 = udiv i32 %37, 31
br label %49
39: ; preds = %2
%40 = load i32, i32* %5
%41 = sub i32 %40, 16
%42 = shl i32 %41, 1
%43 = xor i32 %42, 31
%44 = load i32, i32* %6
%45 = xor i32 %44, 31
%46 = mul i32 %43, %45
%47 = udiv i32 %46, 31
%48 = xor i32 %47, 31
br label %49
49: ; preds = %39, %33
%50 = phi i32 [ %38, %33 ], [ %48, %39 ]
%51 = shl i32 %50, 11
%52 = load i32, i32* %7
%53 = icmp ult i32 %52, 32
br i1 %53, label %54, label %60
54: ; preds = %49
%55 = load i32, i32* %7
%56 = shl i32 %55, 1
%57 = load i32, i32* %8
%58 = mul i32 %56, %57
%59 = udiv i32 %58, 63
br label %70
60: ; preds = %49
%61 = load i32, i32* %7
%62 = sub i32 %61, 32
%63 = shl i32 %62, 1
%64 = xor i32 %63, 63
%65 = load i32, i32* %8
%66 = xor i32 %65, 63
%67 = mul i32 %64, %66
%68 = udiv i32 %67, 63
%69 = xor i32 %68, 63
br label %70
70: ; preds = %60, %54
%71 = phi i32 [ %59, %54 ], [ %69, %60 ]
%72 = shl i32 %71, 5
%73 = or i32 %51, %72
%74 = load i32, i32* %9
%75 = icmp ult i32 %74, 16
br i1 %75, label %76, label %82
76: ; preds = %70
%77 = load i32, i32* %9
%78 = shl i32 %77, 1
%79 = load i32, i32* %10
%80 = mul i32 %78, %79
%81 = udiv i32 %80, 31
br label %92
82: ; preds = %70
%83 = load i32, i32* %9
%84 = sub i32 %83, 16
%85 = shl i32 %84, 1
%86 = xor i32 %85, 31
%87 = load i32, i32* %10
%88 = xor i32 %87, 31
%89 = mul i32 %86, %88
%90 = udiv i32 %89, 31
%91 = xor i32 %90, 31
br label %92
92: ; preds = %82, %76
%93 = phi i32 [ %81, %76 ], [ %91, %82 ]
%94 = or i32 %73, %93
%95 = trunc i32 %94 to i16
ret i16 %95
}

LLVM IR (-O2):
define zeroext i16 @blend_hardlight16(i16 zeroext %0, i16 zeroext %1) {
%3 = zext i16 %0 to i32
%4 = zext i16 %1 to i32
%5 = lshr i32 %4, 11
%6 = lshr i32 %3, 5
%7 = and i32 %6, 63
%8 = lshr i32 %4, 5
%9 = and i32 %8, 63
%10 = and i32 %3, 31
%11 = and i32 %4, 31
%12 = icmp sgt i16 %0, -1
%13 = lshr i32 %3, 10
%14 = and i32 %13, 62
br i1 %12, label %15, label %20
15: ; preds = %2
%16 = mul nuw nsw i32 %14, %5
%17 = trunc i32 %16 to i16
%18 = udiv i16 %17, 31
%19 = zext i16 %18 to i32
br label %27
20: ; preds = %2
%21 = add nsw i32 %14, -32
%22 = xor i32 %21, 31
%23 = xor i32 %5, 31
%24 = mul nsw i32 %22, %23
%25 = udiv i32 %24, 31
%26 = xor i32 %25, 31
br label %27
27: ; preds = %20, %15
%28 = phi i32 [ %19, %15 ], [ %26, %20 ]
%29 = icmp ult i32 %7, 32
%30 = shl nuw nsw i32 %7, 1
br i1 %29, label %31, label %36
31: ; preds = %27
%32 = mul nuw nsw i32 %30, %9
%33 = trunc i32 %32 to i16
%34 = udiv i16 %33, 63
%35 = zext i16 %34 to i32
br label %43
36: ; preds = %27
%37 = add nsw i32 %30, -64
%38 = xor i32 %37, 63
%39 = xor i32 %9, 63
%40 = mul nsw i32 %38, %39
%41 = udiv i32 %40, 63
%42 = xor i32 %41, 63
br label %43
43: ; preds = %36, %31
%44 = phi i32 [ %35, %31 ], [ %42, %36 ]
%45 = icmp ult i32 %10, 16
%46 = shl nuw nsw i32 %10, 1
br i1 %45, label %47, label %52
47: ; preds = %43
%48 = mul nuw nsw i32 %46, %11
%49 = trunc i32 %48 to i16
%50 = udiv i16 %49, 31
%51 = zext i16 %50 to i32
br label %59
52: ; preds = %43
%53 = add nsw i32 %46, -32
%54 = xor i32 %53, 31
%55 = xor i32 %11, 31
%56 = mul nsw i32 %54, %55
%57 = udiv i32 %56, 31
%58 = xor i32 %57, 31
br label %59
59: ; preds = %52, %47
%60 = phi i32 [ %51, %47 ], [ %58, %52 ]
%61 = shl i32 %28, 11
%62 = shl i32 %44, 5
%63 = or i32 %62, %61
%64 = or i32 %63, %60
%65 = trunc i32 %64 to i16
ret i16 %65
}