TC(R, n) unsigned *R;
int n;
{
  register int rowsize;
  register unsigned mask;
  register unsigned *rowj;
  register unsigned *rp;
  register unsigned *rend;
  register unsigned *ccol;
  unsigned *relend;
  unsigned *cword;
  unsigned *rowi;
  rowsize = (((n) + 31) / 32) * sizeof(unsigned);
  relend = (unsigned *)((unsigned)R + n * rowsize);
  cword = R;
  mask = 1;
  rowi = R;
  while (rowi < relend) {
    ccol = cword;
    rowj = R;
    while (rowj < relend) {
      if (*ccol & mask) {
        rp = rowi;
        rend = (unsigned *)((unsigned)rowj + rowsize);
        while (rowj < rend)
          *rowj++ |= *rp++;
      } else {
        rowj = (unsigned *)((unsigned)rowj + rowsize);
      }
      ccol = (unsigned *)((unsigned)ccol + rowsize);
    }
    mask <<= 1;
    if (mask == 0) {
      mask = 1;
      cword++;
    }
    rowi = (unsigned *)((unsigned)rowi + rowsize);
  }
}


x86 assembly (clang, -O0):
<TC>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1f,%eax
mov    $0x20,%ecx
cltd
idiv   %ecx
cltq
shl    $0x2,%rax
mov    %eax,-0x18(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    %ecx,%eax
mov    %eax,%eax
mov    %rax,-0x48(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x50(%rbp)
movl   $0x1,-0x1c(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x58(%rbp)
mov    -0x58(%rbp),%rax
cmp    -0x48(%rbp),%rax
jae    1232 <TC+0x132>
mov    -0x50(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x48(%rbp),%rax
jae    11fa <TC+0xfa>
mov    -0x40(%rbp),%rax
mov    (%rax),%eax
and    -0x1c(%rbp),%eax
cmp    $0x0,%eax
je     11db <TC+0xdb>
mov    -0x58(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x28(%rbp),%rax
add    -0x18(%rbp),%eax
mov    %eax,%eax
mov    %rax,-0x38(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x38(%rbp),%rax
jae    11d6 <TC+0xd6>
mov    -0x30(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x30(%rbp)
mov    (%rax),%ecx
mov    -0x28(%rbp),%rax
mov    %rax,%rdx
add    $0x4,%rdx
mov    %rdx,-0x28(%rbp)
or     (%rax),%ecx
mov    %ecx,(%rax)
jmp    119f <TC+0x9f>
jmp    11e8 <TC+0xe8>
mov    -0x28(%rbp),%rax
add    -0x18(%rbp),%eax
mov    %eax,%eax
mov    %rax,-0x28(%rbp)
mov    -0x40(%rbp),%rax
add    -0x18(%rbp),%eax
mov    %eax,%eax
mov    %rax,-0x40(%rbp)
jmp    116a <TC+0x6a>
mov    -0x1c(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x1c(%rbp)
cmpl   $0x0,-0x1c(%rbp)
jne    1220 <TC+0x120>
movl   $0x1,-0x1c(%rbp)
mov    -0x50(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x50(%rbp)
mov    -0x58(%rbp),%rax
add    -0x18(%rbp),%eax
mov    %eax,%eax
mov    %rax,-0x58(%rbp)
jmp    114c <TC+0x4c>
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<TC>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
lea    0x1f(%rsi),%eax
lea    0x3e(%rsi),%r14d
test   %eax,%eax
cmovns %eax,%r14d
sar    $0x5,%r14d
shl    $0x2,%r14d
imul   %r14d,%esi
add    %edi,%esi
mov    %rdi,-0x8(%rsp)
cmp    %rdi,%rsi
jbe    12d0 <TC+0x1d0>
mov    -0x8(%rsp),%r9
mov    %r9,%rdi
mov    $0x1,%r11d
jmp    1171 <TC+0x71>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
xor    %eax,%eax
add    %r11d,%r11d
mov    $0x1,%ecx
cmove  %ecx,%r11d
sete   %al
lea    (%r9,%rax,4),%r9
add    %r14d,%edi
cmp    %rsi,%rdi
jae    12d0 <TC+0x1d0>
lea    0x4(%rdi),%r10
mov    %r9,%rcx
mov    -0x8(%rsp),%rbx
jmp    118b <TC+0x8b>
nop
mov    %rax,%rbx
add    %r14d,%ecx
cmp    %rsi,%rbx
jae    1150 <TC+0x50>
lea    (%r14,%rbx,1),%eax
test   %r11d,(%rcx)
je     1180 <TC+0x80>
cmp    %rax,%rbx
jae    1183 <TC+0x83>
lea    0x4(%rbx),%rdx
cmp    %rdx,%rax
mov    %rdx,%r15
cmova  %rax,%r15
mov    %rbx,%rbp
not    %rbp
add    %rbp,%r15
cmp    $0x1c,%r15
jae    11be <TC+0xbe>
mov    %rdi,%rdx
jmp    12b0 <TC+0x1b0>
cmp    %rdx,%rax
cmova  %rax,%rdx
add    %rbp,%rdx
and    $0xfffffffffffffffc,%rdx
lea    (%r10,%rdx,1),%rbp
cmp    %rbp,%rbx
jae    11e9 <TC+0xe9>
add    %rbx,%rdx
add    $0x4,%rdx
cmp    %rdx,%rdi
jae    11e9 <TC+0xe9>
mov    %rdi,%rdx
jmp    12b0 <TC+0x1b0>
shr    $0x2,%r15
add    $0x1,%r15
mov    %r15,%r12
and    $0xfffffffffffffff8,%r12
lea    -0x8(%r12),%rdx
mov    %rdx,%r8
shr    $0x3,%r8
add    $0x1,%r8
test   %rdx,%rdx
je     12c6 <TC+0x1c6>
mov    %r8,%r13
and    $0xfffffffffffffffe,%r13
xor    %ebp,%ebp
nopw   0x0(%rax,%rax,1)
movups (%rdi,%rbp,4),%xmm0
movups 0x10(%rdi,%rbp,4),%xmm1
movups (%rbx,%rbp,4),%xmm2
orps   %xmm0,%xmm2
movups 0x10(%rbx,%rbp,4),%xmm0
orps   %xmm1,%xmm0
movups 0x20(%rbx,%rbp,4),%xmm1
movups 0x30(%rbx,%rbp,4),%xmm3
movups %xmm2,(%rbx,%rbp,4)
movups %xmm0,0x10(%rbx,%rbp,4)
movups 0x20(%rdi,%rbp,4),%xmm0
orps   %xmm1,%xmm0
movups 0x30(%rdi,%rbp,4),%xmm1
orps   %xmm3,%xmm1
movups %xmm0,0x20(%rbx,%rbp,4)
movups %xmm1,0x30(%rbx,%rbp,4)
add    $0x10,%rbp
add    $0xfffffffffffffffe,%r13
jne    1220 <TC+0x120>
test   $0x1,%r8b
je     1296 <TC+0x196>
movups (%rdi,%rbp,4),%xmm0
movups 0x10(%rdi,%rbp,4),%xmm1
movups (%rbx,%rbp,4),%xmm2
orps   %xmm0,%xmm2
movups 0x10(%rbx,%rbp,4),%xmm0
orps   %xmm1,%xmm0
movups %xmm2,(%rbx,%rbp,4)
movups %xmm0,0x10(%rbx,%rbp,4)
lea    (%rbx,%r12,4),%rbx
cmp    %r12,%r15
je     1183 <TC+0x83>
lea    (%rdi,%r12,4),%rdx
nopw   0x0(%rax,%rax,1)
mov    (%rdx),%ebp
or     %ebp,(%rbx)
add    $0x4,%rdx
add    $0x4,%rbx
cmp    %rax,%rbx
jb     12b0 <TC+0x1b0>
jmp    1183 <TC+0x83>
xor    %ebp,%ebp
test   $0x1,%r8b
jne    1275 <TC+0x175>
jmp    1296 <TC+0x196>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define i32 @TC(i32* %0, i32 %1) {
%3 = add nsw i32 %1, 31
%4 = sdiv i32 %3, 32
%5 = shl nsw i32 %4, 2
%6 = ptrtoint i32* %0 to i64
%7 = trunc i64 %6 to i32
%8 = mul nsw i32 %5, %1
%9 = add i32 %8, %7
%10 = zext i32 %9 to i64
%11 = inttoptr i64 %10 to i32*
%12 = icmp ugt i32* %11, %0
br i1 %12, label %13, label %155
13: ; preds = %2, %143
%14 = phi i32* [ %153, %143 ], [ %0, %2 ]
%15 = phi i32* [ %148, %143 ], [ %0, %2 ]
%16 = phi i32 [ %146, %143 ], [ 1, %2 ]
%17 = getelementptr i32, i32* %14, i64 1
br label %18
18: ; preds = %13, %135
%19 = phi i32* [ %141, %135 ], [ %15, %13 ]
%20 = phi i32* [ %136, %135 ], [ %0, %13 ]
%21 = load i32, i32* %19
%22 = and i32 %21, %16
%23 = icmp eq i32 %22, 0
%24 = ptrtoint i32* %20 to i64
%25 = trunc i64 %24 to i32
%26 = add i32 %5, %25
%27 = zext i32 %26 to i64
%28 = inttoptr i64 %27 to i32*
br i1 %23, label %135, label %29
29: ; preds = %18
%30 = icmp ult i32* %20, %28
br i1 %30, label %31, label %135
31: ; preds = %29
%32 = add i64 %24, 4
%33 = call i64 @llvm.umax.i64(i64 %27, i64 %32)
%34 = xor i64 %24, -1
%35 = add i64 %33, %34
%36 = lshr i64 %35, 2
%37 = add nuw nsw i64 %36, 1
%38 = icmp ult i64 %35, 28
br i1 %38, label %123, label %39
39: ; preds = %31
%40 = getelementptr i32, i32* %20, i64 1
%41 = add i64 %24, 4
%42 = call i64 @llvm.umax.i64(i64 %27, i64 %41)
%43 = xor i64 %24, -1
%44 = add i64 %42, %43
%45 = lshr i64 %44, 2
%46 = getelementptr i32, i32* %40, i64 %45
%47 = getelementptr i32, i32* %17, i64 %45
%48 = icmp ult i32* %20, %47
%49 = icmp ult i32* %14, %46
%50 = and i1 %48, %49
br i1 %50, label %123, label %51
51: ; preds = %39
%52 = and i64 %37, 9223372036854775800
%53 = getelementptr i32, i32* %20, i64 %52
%54 = getelementptr i32, i32* %14, i64 %52
%55 = add nsw i64 %52, -8
%56 = lshr exact i64 %55, 3
%57 = add nuw nsw i64 %56, 1
%58 = and i64 %57, 1
%59 = icmp eq i64 %55, 0
br i1 %59, label %101, label %60
60: ; preds = %51
%61 = and i64 %57, 4611686018427387902
br label %62
62: ; preds = %62, %60
%63 = phi i64 [ 0, %60 ], [ %98, %62 ]
%64 = phi i64 [ 0, %60 ], [ %99, %62 ]
%65 = getelementptr i32, i32* %20, i64 %63
%66 = getelementptr i32, i32* %14, i64 %63
%67 = bitcast i32* %66 to <4 x i32>*
%68 = load <4 x i32>, <4 x i32>* %67
%69 = getelementptr i32, i32* %66, i64 4
%70 = bitcast i32* %69 to <4 x i32>*
%71 = load <4 x i32>, <4 x i32>* %70
%72 = bitcast i32* %65 to <4 x i32>*
%73 = load <4 x i32>, <4 x i32>* %72
%74 = getelementptr i32, i32* %65, i64 4
%75 = bitcast i32* %74 to <4 x i32>*
%76 = load <4 x i32>, <4 x i32>* %75
%77 = or <4 x i32> %73, %68
%78 = or <4 x i32> %76, %71
%79 = bitcast i32* %65 to <4 x i32>*
store <4 x i32> %77, <4 x i32>* %79
%80 = bitcast i32* %74 to <4 x i32>*
store <4 x i32> %78, <4 x i32>* %80
%81 = or i64 %63, 8
%82 = getelementptr i32, i32* %20, i64 %81
%83 = getelementptr i32, i32* %14, i64 %81
%84 = bitcast i32* %83 to <4 x i32>*
%85 = load <4 x i32>, <4 x i32>* %84
%86 = getelementptr i32, i32* %83, i64 4
%87 = bitcast i32* %86 to <4 x i32>*
%88 = load <4 x i32>, <4 x i32>* %87
%89 = bitcast i32* %82 to <4 x i32>*
%90 = load <4 x i32>, <4 x i32>* %89
%91 = getelementptr i32, i32* %82, i64 4
%92 = bitcast i32* %91 to <4 x i32>*
%93 = load <4 x i32>, <4 x i32>* %92
%94 = or <4 x i32> %90, %85
%95 = or <4 x i32> %93, %88
%96 = bitcast i32* %82 to <4 x i32>*
store <4 x i32> %94, <4 x i32>* %96
%97 = bitcast i32* %91 to <4 x i32>*
store <4 x i32> %95, <4 x i32>* %97
%98 = add nuw i64 %63, 16
%99 = add i64 %64, 2
%100 = icmp eq i64 %99, %61
br i1 %100, label %101, label %62
101: ; preds = %62, %51
%102 = phi i64 [ 0, %51 ], [ %98, %62 ]
%103 = icmp eq i64 %58, 0
br i1 %103, label %121, label %104
104: ; preds = %101
%105 = getelementptr i32, i32* %20, i64 %102
%106 = getelementptr i32, i32* %14, i64 %102
%107 = bitcast i32* %106 to <4 x i32>*
%108 = load <4 x i32>, <4 x i32>* %107
%109 = getelementptr i32, i32* %106, i64 4
%110 = bitcast i32* %109 to <4 x i32>*
%111 = load <4 x i32>, <4 x i32>* %110
%112 = bitcast i32* %105 to <4 x i32>*
%113 = load <4 x i32>, <4 x i32>* %112
%114 = getelementptr i32, i32* %105, i64 4
%115 = bitcast i32* %114 to <4 x i32>*
%116 = load <4 x i32>, <4 x i32>* %115
%117 = or <4 x i32> %113, %108
%118 = or <4 x i32> %116, %111
%119 = bitcast i32* %105 to <4 x i32>*
store <4 x i32> %117, <4 x i32>* %119
%120 = bitcast i32* %114 to <4 x i32>*
store <4 x i32> %118, <4 x i32>* %120
br label %121
121: ; preds = %101, %104
%122 = icmp eq i64 %37, %52
br i1 %122, label %135, label %123
123: ; preds = %39, %31, %121
%124 = phi i32* [ %20, %39 ], [ %20, %31 ], [ %53, %121 ]
%125 = phi i32* [ %14, %39 ], [ %14, %31 ], [ %54, %121 ]
br label %126
126: ; preds = %123, %126
%127 = phi i32* [ %131, %126 ], [ %124, %123 ]
%128 = phi i32* [ %129, %126 ], [ %125, %123 ]
%129 = getelementptr inbounds i32, i32* %128, i64 1
%130 = load i32, i32* %128
%131 = getelementptr inbounds i32, i32* %127, i64 1
%132 = load i32, i32* %127
%133 = or i32 %132, %130
store i32 %133, i32* %127
%134 = icmp ult i32* %131, %28
br i1 %134, label %126, label %135
135: ; preds = %126, %121, %18, %29
%136 = phi i32* [ %20, %29 ], [ %28, %18 ], [ %53, %121 ], [ %131, %126 ]
%137 = ptrtoint i32* %19 to i64
%138 = trunc i64 %137 to i32
%139 = add i32 %5, %138
%140 = zext i32 %139 to i64
%141 = inttoptr i64 %140 to i32*
%142 = icmp ult i32* %136, %11
br i1 %142, label %18, label %143
143: ; preds = %135
%144 = shl i32 %16, 1
%145 = icmp eq i32 %144, 0
%146 = select i1 %145, i32 1, i32 %144
%147 = zext i1 %145 to i64
%148 = getelementptr i32, i32* %15, i64 %147
%149 = ptrtoint i32* %14 to i64
%150 = trunc i64 %149 to i32
%151 = add i32 %5, %150
%152 = zext i32 %151 to i64
%153 = inttoptr i64 %152 to i32*
%154 = icmp ult i32* %153, %11
br i1 %154, label %13, label %155
155: ; preds = %143, %2
ret i32 undef
}