unsigned short crc16poly(char *data_p, unsigned short length,
                         unsigned short poly) {
  unsigned char i;
  unsigned int data;
  unsigned int crc = 0xffff;
  if (length == 0)
    return (~crc);
  do {
    for (i = 0, data = (unsigned int)0xff & *data_p++; i < 8; i++, data >>= 1) {
      if ((crc & 0x0001) ^ (data & 0x0001))
        crc = (crc >> 1) ^ poly;
      else
        crc >>= 1;
    }
  } while (--length);
  crc = ~crc;
  data = crc;
  crc = (crc << 8) | ((data >> 8) & 0xff);
  return (crc);
}


x86 assembly (clang, -O0):
<crc16poly>:
push   %rbp
mov    %rsp,%rbp
mov    %dx,%ax
mov    %si,%cx
mov    %rdi,-0x10(%rbp)
mov    %cx,-0x12(%rbp)
mov    %ax,-0x14(%rbp)
movl   $0xffff,-0x20(%rbp)
movzwl -0x12(%rbp),%eax
cmp    $0x0,%eax
jne    1139 <crc16poly+0x39>
mov    -0x20(%rbp),%eax
xor    $0xffffffff,%eax
mov    %ax,-0x2(%rbp)
jmp    1200 <crc16poly+0x100>
jmp    113e <crc16poly+0x3e>
movb   $0x0,-0x15(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,%rcx
add    $0x1,%rcx
mov    %rcx,-0x10(%rbp)
movsbl (%rax),%eax
and    $0xff,%eax
mov    %eax,-0x1c(%rbp)
movzbl -0x15(%rbp),%eax
cmp    $0x8,%eax
jge    11b8 <crc16poly+0xb8>
mov    -0x20(%rbp),%eax
and    $0x1,%eax
mov    -0x1c(%rbp),%ecx
and    $0x1,%ecx
xor    %ecx,%eax
cmp    $0x0,%eax
je     1194 <crc16poly+0x94>
mov    -0x20(%rbp),%eax
shr    $0x1,%eax
movzwl -0x14(%rbp),%ecx
xor    %ecx,%eax
mov    %eax,-0x20(%rbp)
jmp    119d <crc16poly+0x9d>
mov    -0x20(%rbp),%eax
shr    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    11a2 <crc16poly+0xa2>
mov    -0x15(%rbp),%al
add    $0x1,%al
mov    %al,-0x15(%rbp)
mov    -0x1c(%rbp),%eax
shr    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    115c <crc16poly+0x5c>
jmp    11bd <crc16poly+0xbd>
mov    -0x12(%rbp),%ax
add    $0xffff,%ax
mov    %ax,-0x12(%rbp)
cmp    $0x0,%ax
jne    113e <crc16poly+0x3e>
mov    -0x20(%rbp),%eax
xor    $0xffffffff,%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x20(%rbp),%eax
shl    $0x8,%eax
mov    -0x1c(%rbp),%ecx
shr    $0x8,%ecx
and    $0xff,%ecx
or     %ecx,%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    %ax,-0x2(%rbp)
movzwl -0x2(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<crc16poly>:
xor    %eax,%eax
test   %esi,%esi
je     11d5 <crc16poly+0xd5>
push   %rbx
movzwl %dx,%r9d
movzwl %si,%r8d
mov    $0xffff,%ebx
xor    %r10d,%r10d
nopl   0x0(%rax,%rax,1)
movzbl (%rdi,%r10,1),%esi
mov    %ebx,%ecx
xor    %esi,%ecx
and    $0x1,%ecx
neg    %ecx
shr    %ebx
and    %r9d,%ecx
xor    %ebx,%ecx
mov    %esi,%edx
shr    %edx
xor    %ecx,%edx
and    $0x1,%edx
neg    %edx
shr    %ecx
and    %r9d,%edx
xor    %ecx,%edx
mov    %esi,%ecx
shr    $0x2,%ecx
xor    %edx,%ecx
and    $0x1,%ecx
neg    %ecx
shr    %edx
and    %r9d,%ecx
xor    %edx,%ecx
mov    %esi,%edx
shr    $0x3,%edx
xor    %ecx,%edx
and    $0x1,%edx
neg    %edx
shr    %ecx
and    %r9d,%edx
xor    %ecx,%edx
mov    %esi,%ecx
shr    $0x4,%ecx
xor    %edx,%ecx
and    $0x1,%ecx
neg    %ecx
shr    %edx
and    %r9d,%ecx
xor    %edx,%ecx
mov    %esi,%edx
shr    $0x5,%edx
xor    %ecx,%edx
and    $0x1,%edx
neg    %edx
shr    %ecx
and    %r9d,%edx
xor    %ecx,%edx
mov    %esi,%ecx
shr    $0x6,%ecx
xor    %edx,%ecx
and    $0x1,%ecx
neg    %ecx
shr    %edx
and    %r9d,%ecx
xor    %edx,%ecx
shr    $0x7,%esi
mov    %ecx,%edx
and    $0x1,%edx
shr    %ecx
cmp    %esi,%edx
mov    %r9d,%ebx
cmove  %eax,%ebx
xor    %ecx,%ebx
add    $0x1,%r10
cmp    %r10w,%r8w
jne    1120 <crc16poly+0x20>
not    %ebx
mov    %ebx,%ecx
shl    $0x8,%ecx
movzbl %bh,%eax
or     %ecx,%eax
pop    %rbx
ret

LLVM IR (-O0):
define zeroext i16 @crc16poly(i8* %0, i16 zeroext %1, i16 zeroext %2) {
%4 = alloca i16
%5 = alloca i8*
%6 = alloca i16
%7 = alloca i16
%8 = alloca i8
%9 = alloca i32
%10 = alloca i32
store i8* %0, i8** %5
store i16 %1, i16* %6
store i16 %2, i16* %7
store i32 65535, i32* %10
%11 = load i16, i16* %6
%12 = zext i16 %11 to i32
%13 = icmp eq i32 %12, 0
br i1 %13, label %14, label %18
14: ; preds = %3
%15 = load i32, i32* %10
%16 = xor i32 %15, -1
%17 = trunc i32 %16 to i16
store i16 %17, i16* %4
br label %68
18: ; preds = %3
br label %19
19: ; preds = %52, %18
store i8 0, i8* %8
%20 = load i8*, i8** %5
%21 = getelementptr inbounds i8, i8* %20, i32 1
store i8* %21, i8** %5
%22 = load i8, i8* %20
%23 = sext i8 %22 to i32
%24 = and i32 255, %23
store i32 %24, i32* %9
br label %25
25: ; preds = %46, %19
%26 = load i8, i8* %8
%27 = zext i8 %26 to i32
%28 = icmp slt i32 %27, 8
br i1 %28, label %29, label %51
29: ; preds = %25
%30 = load i32, i32* %10
%31 = and i32 %30, 1
%32 = load i32, i32* %9
%33 = and i32 %32, 1
%34 = xor i32 %31, %33
%35 = icmp ne i32 %34, 0
br i1 %35, label %36, label %42
36: ; preds = %29
%37 = load i32, i32* %10
%38 = lshr i32 %37, 1
%39 = load i16, i16* %7
%40 = zext i16 %39 to i32
%41 = xor i32 %38, %40
store i32 %41, i32* %10
br label %45
42: ; preds = %29
%43 = load i32, i32* %10
%44 = lshr i32 %43, 1
store i32 %44, i32* %10
br label %45
45: ; preds = %42, %36
br label %46
46: ; preds = %45
%47 = load i8, i8* %8
%48 = add i8 %47, 1
store i8 %48, i8* %8
%49 = load i32, i32* %9
%50 = lshr i32 %49, 1
store i32 %50, i32* %9
br label %25
51: ; preds = %25
br label %52
52: ; preds = %51
%53 = load i16, i16* %6
%54 = add i16 %53, -1
store i16 %54, i16* %6
%55 = icmp ne i16 %54, 0
br i1 %55, label %19, label %56
56: ; preds = %52
%57 = load i32, i32* %10
%58 = xor i32 %57, -1
store i32 %58, i32* %10
%59 = load i32, i32* %10
store i32 %59, i32* %9
%60 = load i32, i32* %10
%61 = shl i32 %60, 8
%62 = load i32, i32* %9
%63 = lshr i32 %62, 8
%64 = and i32 %63, 255
%65 = or i32 %61, %64
store i32 %65, i32* %10
%66 = load i32, i32* %10
%67 = trunc i32 %66 to i16
store i16 %67, i16* %4
br label %68
68: ; preds = %56, %14
%69 = load i16, i16* %4
ret i16 %69
}

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