int bch(int totalBits, int dataBits, int data, int gpb) {
  int ecBits = totalBits - dataBits;
  int top = 1 << (totalBits - 1);
  data <<= ecBits;
  gpb <<= (dataBits - 1);
  while (top > 1 << (ecBits - 1)) {
    if ((top & data)) {
      data ^= gpb;
    }
    top >>= 1;
    gpb >>= 1;
  }
  return data;
}


x86 assembly (clang, -O0):
<bch>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %ecx,-0x10(%rbp)
mov    -0x4(%rbp),%eax
sub    -0x8(%rbp),%eax
mov    %eax,-0x14(%rbp)
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
mov    $0x1,%eax
shl    %cl,%eax
mov    %eax,-0x18(%rbp)
mov    -0x14(%rbp),%ecx
mov    -0xc(%rbp),%eax
shl    %cl,%eax
mov    %eax,-0xc(%rbp)
mov    -0x8(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x10(%rbp),%eax
shl    %cl,%eax
mov    %eax,-0x10(%rbp)
mov    -0x18(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
mov    $0x1,%edx
shl    %cl,%edx
mov    %edx,%ecx
cmp    %ecx,%eax
jle    118b <bch+0x8b>
mov    -0x18(%rbp),%eax
and    -0xc(%rbp),%eax
cmp    $0x0,%eax
je     1174 <bch+0x74>
mov    -0x10(%rbp),%eax
xor    -0xc(%rbp),%eax
mov    %eax,-0xc(%rbp)
mov    -0x18(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x18(%rbp)
mov    -0x10(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1142 <bch+0x42>
mov    -0xc(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<bch>:
mov    %ecx,%r8d
mov    %edx,%eax
lea    -0x1(%rdi),%ecx
mov    %edi,%r9d
sub    %esi,%r9d
mov    $0x1,%edi
mov    $0x1,%edx
shl    %cl,%edx
mov    %r9d,%ecx
shl    %cl,%eax
add    $0xff,%r9b
mov    %r9d,%ecx
shl    %cl,%edi
cmp    %edi,%edx
jle    1153 <bch+0x53>
add    $0xff,%sil
mov    %esi,%ecx
shl    %cl,%r8d
cs nopw 0x0(%rax,%rax,1)
nop
mov    %edx,%ecx
and    %eax,%ecx
cmovne %r8d,%ecx
xor    %ecx,%eax
sar    %edx
sar    %r8d
cmp    %edi,%edx
jg     1140 <bch+0x40>
ret

LLVM IR (-O0):
define i32 @bch(i32 %0, i32 %1, i32 %2, i32 %3) {
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
store i32 %0, i32* %5
store i32 %1, i32* %6
store i32 %2, i32* %7
store i32 %3, i32* %8
%11 = load i32, i32* %5
%12 = load i32, i32* %6
%13 = sub nsw i32 %11, %12
store i32 %13, i32* %9
%14 = load i32, i32* %5
%15 = sub nsw i32 %14, 1
%16 = shl i32 1, %15
store i32 %16, i32* %10
%17 = load i32, i32* %9
%18 = load i32, i32* %7
%19 = shl i32 %18, %17
store i32 %19, i32* %7
%20 = load i32, i32* %6
%21 = sub nsw i32 %20, 1
%22 = load i32, i32* %8
%23 = shl i32 %22, %21
store i32 %23, i32* %8
br label %24
24: ; preds = %39, %4
%25 = load i32, i32* %10
%26 = load i32, i32* %9
%27 = sub nsw i32 %26, 1
%28 = shl i32 1, %27
%29 = icmp sgt i32 %25, %28
br i1 %29, label %30, label %44
30: ; preds = %24
%31 = load i32, i32* %10
%32 = load i32, i32* %7
%33 = and i32 %31, %32
%34 = icmp ne i32 %33, 0
br i1 %34, label %35, label %39
35: ; preds = %30
%36 = load i32, i32* %8
%37 = load i32, i32* %7
%38 = xor i32 %37, %36
store i32 %38, i32* %7
br label %39
39: ; preds = %35, %30
%40 = load i32, i32* %10
%41 = ashr i32 %40, 1
store i32 %41, i32* %10
%42 = load i32, i32* %8
%43 = ashr i32 %42, 1
store i32 %43, i32* %8
br label %24
44: ; preds = %24
%45 = load i32, i32* %7
ret i32 %45
}

LLVM IR (-O2):
define i32 @bch(i32 %0, i32 %1, i32 %2, i32 %3) {
%5 = sub nsw i32 %0, %1
%6 = add nsw i32 %0, -1
%7 = shl nuw i32 1, %6
%8 = shl i32 %2, %5
%9 = add nsw i32 %5, -1
%10 = shl nuw i32 1, %9
%11 = icmp sgt i32 %7, %10
br i1 %11, label %12, label %26
12: ; preds = %4
%13 = add nsw i32 %1, -1
%14 = shl i32 %3, %13
br label %15
15: ; preds = %12, %15
%16 = phi i32 [ %23, %15 ], [ %7, %12 ]
%17 = phi i32 [ %22, %15 ], [ %8, %12 ]
%18 = phi i32 [ %24, %15 ], [ %14, %12 ]
%19 = and i32 %16, %17
%20 = icmp eq i32 %19, 0
%21 = select i1 %20, i32 0, i32 %18
%22 = xor i32 %21, %17
%23 = ashr i32 %16, 1
%24 = ashr i32 %18, 1
%25 = icmp sgt i32 %23, %10
br i1 %25, label %15, label %26
26: ; preds = %15, %4
%27 = phi i32 [ %8, %4 ], [ %22, %15 ]
ret i32 %27
}