void BWTIncPutPackedTextToRank(const unsigned int *packedText,
                               unsigned int *__restrict rank,
                               unsigned int *__restrict cumulativeCount,
                               const unsigned int numChar) {
  unsigned int i, j;
  unsigned int c, t;
  unsigned int packedMask;
  unsigned int rankIndex;
  unsigned int lastWord, numCharInLastWord;
  lastWord = (numChar - 1) / 16;
  numCharInLastWord = numChar - lastWord * 16;
  packedMask = 0xFFFFFFFF >> (32 - 2);
  rankIndex = numChar - 1;
  t = packedText[lastWord] >> (32 - numCharInLastWord * 2);
  for (i = 0; i < numCharInLastWord; i++) {
    c = t & packedMask;
    cumulativeCount[c + 1]++;
    rank[rankIndex] = c;
    rankIndex--;
    t >>= 2;
  }
  for (i = lastWord; i--;) {
    t = packedText[i];
    for (j = 0; j < 16; j++) {
      c = t & packedMask;
      cumulativeCount[c + 1]++;
      rank[rankIndex] = c;
      rankIndex--;
      t >>= 2;
    }
  }
  cumulativeCount[2] += cumulativeCount[1];
  cumulativeCount[3] += cumulativeCount[2];
  cumulativeCount[4] += cumulativeCount[3];
}


x86 assembly (clang, -O0):
<BWTIncPutPackedTextToRank>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
sub    $0x1,%eax
shr    $0x4,%eax
mov    %eax,-0x38(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x38(%rbp),%ecx
shl    $0x4,%ecx
sub    %ecx,%eax
mov    %eax,-0x3c(%rbp)
movl   $0x3,-0x30(%rbp)
mov    -0x1c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    (%rax,%rcx,4),%eax
mov    -0x3c(%rbp),%edx
shl    $0x1,%edx
mov    $0x20,%ecx
sub    %edx,%ecx
shr    %cl,%eax
mov    %eax,-0x2c(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jae    11b7 <BWTIncPutPackedTextToRank+0xb7>
mov    -0x2c(%rbp),%eax
and    -0x30(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,%ecx
mov    (%rax,%rcx,4),%edx
add    $0x1,%edx
mov    %edx,(%rax,%rcx,4)
mov    -0x28(%rbp),%edx
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    %edx,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x34(%rbp)
mov    -0x2c(%rbp),%eax
shr    $0x2,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1160 <BWTIncPutPackedTextToRank+0x60>
mov    -0x38(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    %eax,%ecx
add    $0xffffffff,%ecx
mov    %ecx,-0x20(%rbp)
cmp    $0x0,%eax
je     123f <BWTIncPutPackedTextToRank+0x13f>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x2c(%rbp)
movl   $0x0,-0x24(%rbp)
cmpl   $0x10,-0x24(%rbp)
jae    123a <BWTIncPutPackedTextToRank+0x13a>
mov    -0x2c(%rbp),%eax
and    -0x30(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,%ecx
mov    (%rax,%rcx,4),%edx
add    $0x1,%edx
mov    %edx,(%rax,%rcx,4)
mov    -0x28(%rbp),%edx
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    %edx,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x34(%rbp)
mov    -0x2c(%rbp),%eax
shr    $0x2,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11e5 <BWTIncPutPackedTextToRank+0xe5>
jmp    11bd <BWTIncPutPackedTextToRank+0xbd>
mov    -0x18(%rbp),%rax
mov    0x4(%rax),%ecx
mov    -0x18(%rbp),%rax
add    0x8(%rax),%ecx
mov    %ecx,0x8(%rax)
mov    -0x18(%rbp),%rax
mov    0x8(%rax),%ecx
mov    -0x18(%rbp),%rax
add    0xc(%rax),%ecx
mov    %ecx,0xc(%rax)
mov    -0x18(%rbp),%rax
mov    0xc(%rax),%ecx
mov    -0x18(%rbp),%rax
add    0x10(%rax),%ecx
mov    %ecx,0x10(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<BWTIncPutPackedTextToRank>:
push   %r15
push   %r14
push   %rbx
mov    %ecx,%r8d
lea    -0x1(%r8),%r9d
mov    %r9d,%r10d
shr    $0x4,%r10d
mov    %r9d,%eax
and    $0xfffffff0,%eax
mov    %ecx,%r11d
mov    %r9d,%r15d
sub    %eax,%r11d
je     119f <BWTIncPutPackedTextToRank+0x9f>
mov    %r10d,%ecx
mov    (%rdi,%rcx,4),%r14d
lea    (%r11,%r11,1),%ecx
neg    %cl
shr    %cl,%r14d
not    %eax
add    %r8d,%eax
jne    1140 <BWTIncPutPackedTextToRank+0x40>
mov    %r9d,%r15d
jmp    1184 <BWTIncPutPackedTextToRank+0x84>
mov    %r11d,%r8d
and    $0xfffffffe,%r8d
mov    %r9d,%r15d
nopw   0x0(%rax,%rax,1)
mov    %r14d,%ebx
and    $0x3,%ebx
addl   $0x1,0x4(%rdx,%rbx,4)
mov    %r15d,%eax
mov    %ebx,(%rsi,%rax,4)
lea    -0x1(%r15),%eax
mov    %r14d,%ebx
shr    $0x2,%ebx
and    $0x3,%ebx
addl   $0x1,0x4(%rdx,%rbx,4)
mov    %ebx,(%rsi,%rax,4)
add    $0xfffffffe,%r15d
shr    $0x4,%r14d
add    $0xfffffffe,%r8d
jne    1150 <BWTIncPutPackedTextToRank+0x50>
test   $0x1,%r11b
je     119f <BWTIncPutPackedTextToRank+0x9f>
and    $0x3,%r14d
addl   $0x1,0x4(%rdx,%r14,4)
mov    %r15d,%eax
mov    %r14d,(%rsi,%rax,4)
add    $0xffffffff,%r15d
cmp    $0x10,%r9d
jb     11f6 <BWTIncPutPackedTextToRank+0xf6>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0xffffffff,%r10d
mov    (%rdi,%r10,4),%eax
mov    $0xfffffff0,%r8d
xchg   %ax,%ax
mov    %eax,%ecx
and    $0x3,%ecx
addl   $0x1,0x4(%rdx,%rcx,4)
mov    %r15d,%ebx
mov    %ecx,(%rsi,%rbx,4)
lea    -0x1(%r15),%ecx
mov    %eax,%ebx
shr    $0x2,%ebx
and    $0x3,%ebx
addl   $0x1,0x4(%rdx,%rbx,4)
mov    %ebx,(%rsi,%rcx,4)
add    $0xfffffffe,%r15d
shr    $0x4,%eax
add    $0x2,%r8d
jne    11c0 <BWTIncPutPackedTextToRank+0xc0>
test   %r10d,%r10d
jne    11b0 <BWTIncPutPackedTextToRank+0xb0>
mov    0x8(%rdx),%eax
add    0x4(%rdx),%eax
mov    %eax,0x8(%rdx)
add    0xc(%rdx),%eax
mov    %eax,0xc(%rdx)
add    %eax,0x10(%rdx)
pop    %rbx
pop    %r14
pop    %r15
ret

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

LLVM IR (-O2):
define void @BWTIncPutPackedTextToRank(i32* nocapture %0, i32* noalias nocapture writeonly %1, i32* noalias nocapture %2, i32 %3) {
%5 = add i32 %3, -1
%6 = lshr i32 %5, 4
%7 = and i32 %5, -16
%8 = sub i32 %3, %7
%9 = icmp eq i32 %8, 0
br i1 %9, label %38, label %10
10: ; preds = %4
%11 = zext i32 %6 to i64
%12 = getelementptr inbounds i32, i32* %0, i64 %11
%13 = load i32, i32* %12
%14 = mul i32 %8, -2
%15 = add i32 %14, 32
%16 = lshr i32 %13, %15
%17 = xor i32 %7, -1
%18 = and i32 %8, 1
%19 = sub i32 0, %3
%20 = icmp eq i32 %17, %19
br i1 %20, label %23, label %21
21: ; preds = %10
%22 = and i32 %8, -2
br label %41
23: ; preds = %41, %10
%24 = phi i32 [ undef, %10 ], [ %63, %41 ]
%25 = phi i32 [ %5, %10 ], [ %63, %41 ]
%26 = phi i32 [ %16, %10 ], [ %64, %41 ]
%27 = icmp eq i32 %18, 0
br i1 %27, label %38, label %28
28: ; preds = %23
%29 = and i32 %26, 3
%30 = add nuw nsw i32 %29, 1
%31 = zext i32 %30 to i64
%32 = getelementptr inbounds i32, i32* %2, i64 %31
%33 = load i32, i32* %32
%34 = add i32 %33, 1
store i32 %34, i32* %32
%35 = zext i32 %25 to i64
%36 = getelementptr inbounds i32, i32* %1, i64 %35
store i32 %29, i32* %36
%37 = add i32 %25, -1
br label %38
38: ; preds = %28, %23, %4
%39 = phi i32 [ %5, %4 ], [ %24, %23 ], [ %37, %28 ]
%40 = icmp ult i32 %5, 16
br i1 %40, label %102, label %69
41: ; preds = %41, %21
%42 = phi i32 [ %5, %21 ], [ %63, %41 ]
%43 = phi i32 [ %16, %21 ], [ %64, %41 ]
%44 = phi i32 [ 0, %21 ], [ %65, %41 ]
%45 = and i32 %43, 3
%46 = add nuw nsw i32 %45, 1
%47 = zext i32 %46 to i64
%48 = getelementptr inbounds i32, i32* %2, i64 %47
%49 = load i32, i32* %48
%50 = add i32 %49, 1
store i32 %50, i32* %48
%51 = zext i32 %42 to i64
%52 = getelementptr inbounds i32, i32* %1, i64 %51
store i32 %45, i32* %52
%53 = add i32 %42, -1
%54 = lshr i32 %43, 2
%55 = and i32 %54, 3
%56 = add nuw nsw i32 %55, 1
%57 = zext i32 %56 to i64
%58 = getelementptr inbounds i32, i32* %2, i64 %57
%59 = load i32, i32* %58
%60 = add i32 %59, 1
store i32 %60, i32* %58
%61 = zext i32 %53 to i64
%62 = getelementptr inbounds i32, i32* %1, i64 %61
store i32 %55, i32* %62
%63 = add i32 %42, -2
%64 = lshr i32 %43, 4
%65 = add i32 %44, 2
%66 = icmp eq i32 %65, %22
br i1 %66, label %23, label %41
67: ; preds = %76
%68 = icmp eq i32 %72, 0
br i1 %68, label %102, label %69
69: ; preds = %38, %67
%70 = phi i32 [ %72, %67 ], [ %6, %38 ]
%71 = phi i32 [ %98, %67 ], [ %39, %38 ]
%72 = add i32 %70, -1
%73 = zext i32 %72 to i64
%74 = getelementptr inbounds i32, i32* %0, i64 %73
%75 = load i32, i32* %74
br label %76
76: ; preds = %76, %69
%77 = phi i32 [ 0, %69 ], [ %100, %76 ]
%78 = phi i32 [ %71, %69 ], [ %98, %76 ]
%79 = phi i32 [ %75, %69 ], [ %99, %76 ]
%80 = and i32 %79, 3
%81 = add nuw nsw i32 %80, 1
%82 = zext i32 %81 to i64
%83 = getelementptr inbounds i32, i32* %2, i64 %82
%84 = load i32, i32* %83
%85 = add i32 %84, 1
store i32 %85, i32* %83
%86 = zext i32 %78 to i64
%87 = getelementptr inbounds i32, i32* %1, i64 %86
store i32 %80, i32* %87
%88 = add i32 %78, -1
%89 = lshr i32 %79, 2
%90 = and i32 %89, 3
%91 = add nuw nsw i32 %90, 1
%92 = zext i32 %91 to i64
%93 = getelementptr inbounds i32, i32* %2, i64 %92
%94 = load i32, i32* %93
%95 = add i32 %94, 1
store i32 %95, i32* %93
%96 = zext i32 %88 to i64
%97 = getelementptr inbounds i32, i32* %1, i64 %96
store i32 %90, i32* %97
%98 = add i32 %78, -2
%99 = lshr i32 %79, 4
%100 = add nuw nsw i32 %77, 2
%101 = icmp eq i32 %100, 16
br i1 %101, label %67, label %76
102: ; preds = %67, %38
%103 = getelementptr inbounds i32, i32* %2, i64 1
%104 = load i32, i32* %103
%105 = getelementptr inbounds i32, i32* %2, i64 2
%106 = load i32, i32* %105
%107 = add i32 %106, %104
store i32 %107, i32* %105
%108 = getelementptr inbounds i32, i32* %2, i64 3
%109 = load i32, i32* %108
%110 = add i32 %109, %107
store i32 %110, i32* %108
%111 = getelementptr inbounds i32, i32* %2, i64 4
%112 = load i32, i32* %111
%113 = add i32 %112, %110
store i32 %113, i32* %111
ret void
}