void transitive_closure(unsigned *R, int n) {
  int rowsize;
  unsigned i;
  unsigned *rowj;
  unsigned *rp;
  unsigned *rend;
  unsigned *ccol;
  unsigned *relend;
  unsigned *cword;
  unsigned *rowi;
  rowsize = (((n) + (32 - 1) / 32));
  relend = R + n * rowsize;
  cword = R;
  i = 0;
  rowi = R;
  while (rowi < relend) {
    ccol = cword;
    rowj = R;
    while (rowj < relend) {
      if (*ccol & (1 << i)) {
        rp = rowi;
        rend = rowj + rowsize;
        while (rowj < rend)
          *rowj++ |= *rp++;
      } else {
        rowj += rowsize;
      }
      ccol += rowsize;
    }
    if (++i >= 32) {
      i = 0;
      cword++;
    }
    rowi += rowsize;
  }
}


x86 assembly (clang, -O0):
<transitive_closure>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    -0xc(%rbp),%eax
add    $0x0,%eax
mov    %eax,-0x10(%rbp)
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
imul   -0x10(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x48(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x50(%rbp),%rax
cmp    -0x40(%rbp),%rax
jae    1252 <transitive_closure+0x152>
mov    -0x48(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x20(%rbp)
mov    -0x20(%rbp),%rax
cmp    -0x40(%rbp),%rax
jae    1213 <transitive_closure+0x113>
mov    -0x38(%rbp),%rax
mov    (%rax),%eax
mov    -0x14(%rbp),%ecx
mov    $0x1,%edx
shl    %cl,%edx
mov    %edx,%ecx
and    %ecx,%eax
cmp    $0x0,%eax
je     11e4 <transitive_closure+0xe4>
mov    -0x50(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x10(%rbp),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
mov    -0x20(%rbp),%rax
cmp    -0x30(%rbp),%rax
jae    11df <transitive_closure+0xdf>
mov    -0x28(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x28(%rbp)
mov    (%rax),%ecx
mov    -0x20(%rbp),%rax
mov    %rax,%rdx
add    $0x4,%rdx
mov    %rdx,-0x20(%rbp)
or     (%rax),%ecx
mov    %ecx,(%rax)
jmp    11a8 <transitive_closure+0xa8>
jmp    11f9 <transitive_closure+0xf9>
mov    -0x10(%rbp),%ecx
mov    -0x20(%rbp),%rax
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x20(%rbp)
mov    -0x10(%rbp),%ecx
mov    -0x38(%rbp),%rax
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
jmp    1162 <transitive_closure+0x62>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
cmp    $0x20,%eax
jb     1238 <transitive_closure+0x138>
movl   $0x0,-0x14(%rbp)
mov    -0x48(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x48(%rbp)
mov    -0x10(%rbp),%ecx
mov    -0x50(%rbp),%rax
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x50(%rbp)
jmp    1144 <transitive_closure+0x44>
pop    %rbp
ret

x86 assembly (clang, -O2):
<transitive_closure>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rdi,-0x28(%rsp)
mov    %esi,%eax
imul   %eax,%eax
test   %eax,%eax
je     133a <transitive_closure+0x23a>
mov    -0x28(%rsp),%r10
lea    (%r10,%rax,4),%r11
movslq %esi,%r13
lea    0x0(,%r13,4),%r14
mov    %r10,%r8
xor    %eax,%eax
xor    %ecx,%ecx
mov    %esi,-0x2c(%rsp)
mov    %r11,-0x18(%rsp)
mov    %r14,-0x20(%rsp)
jmp    1181 <transitive_closure+0x81>
nopw   0x0(%rax,%rax,1)
add    $0x1,%ecx
xor    %eax,%eax
cmp    $0x20,%ecx
setae  %al
mov    $0x0,%edx
cmovae %edx,%ecx
mov    -0x10(%rsp),%r10
lea    (%r10,%rax,4),%r10
add    %r14,%r8
mov    -0x8(%rsp),%rax
add    $0x1,%rax
cmp    %r11,%r8
jae    133a <transitive_closure+0x23a>
mov    %rax,-0x8(%rsp)
imul   %r13,%rax
mov    -0x28(%rsp),%rdx
lea    (%rdx,%rax,4),%r15
add    $0x4,%r15
mov    $0x1,%r12d
shl    %cl,%r12d
mov    %r10,-0x10(%rsp)
mov    %rdx,%rdi
jmp    11cc <transitive_closure+0xcc>
mov    %rdx,%rdi
mov    -0x2c(%rsp),%esi
mov    -0x18(%rsp),%r11
mov    -0x20(%rsp),%r14
nopl   0x0(%rax,%rax,1)
lea    (%r10,%r13,4),%r10
mov    %rdi,%rdx
cmp    %r11,%rdi
jae    1150 <transitive_closure+0x50>
lea    (%rdx,%r13,4),%rbp
test   %r12d,(%r10)
cmove  %rbp,%rdi
je     11c0 <transitive_closure+0xc0>
test   %esi,%esi
jle    11c0 <transitive_closure+0xc0>
lea    (%r14,%rdx,1),%rdi
lea    0x4(%rdx),%rax
cmp    %rax,%rdi
mov    %rax,%r9
cmova  %rdi,%r9
mov    %rdx,%rbx
not    %rbx
add    %rbx,%r9
cmp    $0x1c,%r9
jb     1221 <transitive_closure+0x121>
cmp    %rax,%rdi
cmova  %rdi,%rax
add    %rbx,%rax
and    $0xfffffffffffffffc,%rax
lea    (%r15,%rax,1),%rdi
cmp    %rdi,%rdx
jae    1246 <transitive_closure+0x146>
add    %rdx,%rax
add    $0x4,%rax
cmp    %rax,%r8
jae    1246 <transitive_closure+0x146>
mov    %r8,%rax
mov    %rdx,%rdi
nopw   0x0(%rax,%rax,1)
mov    (%rax),%edx
or     %edx,(%rdi)
add    $0x4,%rax
add    $0x4,%rdi
cmp    %rbp,%rdi
jb     1230 <transitive_closure+0x130>
jmp    11c0 <transitive_closure+0xc0>
mov    %r15,%rsi
shr    $0x2,%r9
add    $0x1,%r9
mov    %r9,%r11
and    $0xfffffffffffffff8,%r11
lea    -0x8(%r11),%rax
mov    %rax,%r15
shr    $0x3,%r15
add    $0x1,%r15
test   %rax,%rax
je     132f <transitive_closure+0x22f>
mov    %r15,%rdi
and    $0xfffffffffffffffe,%rdi
xor    %r14d,%r14d
nopw   0x0(%rax,%rax,1)
movups (%r8,%r14,4),%xmm0
movups 0x10(%r8,%r14,4),%xmm1
movups (%rdx,%r14,4),%xmm2
orps   %xmm0,%xmm2
movups 0x10(%rdx,%r14,4),%xmm0
orps   %xmm1,%xmm0
movups 0x20(%rdx,%r14,4),%xmm1
movups 0x30(%rdx,%r14,4),%xmm3
movups %xmm2,(%rdx,%r14,4)
movups %xmm0,0x10(%rdx,%r14,4)
movups 0x20(%r8,%r14,4),%xmm0
orps   %xmm1,%xmm0
movups 0x30(%r8,%r14,4),%xmm1
orps   %xmm3,%xmm1
movups %xmm0,0x20(%rdx,%r14,4)
movups %xmm1,0x30(%rdx,%r14,4)
add    $0x10,%r14
add    $0xfffffffffffffffe,%rdi
jne    1280 <transitive_closure+0x180>
test   $0x1,%r15b
je     1308 <transitive_closure+0x208>
movups (%r8,%r14,4),%xmm0
movups 0x10(%r8,%r14,4),%xmm1
movups (%rdx,%r14,4),%xmm2
orps   %xmm0,%xmm2
movups 0x10(%rdx,%r14,4),%xmm0
orps   %xmm1,%xmm0
movups %xmm2,(%rdx,%r14,4)
movups %xmm0,0x10(%rdx,%r14,4)
lea    (%rdx,%r11,4),%rdx
cmp    %r11,%r9
mov    %rsi,%r15
je     11aa <transitive_closure+0xaa>
lea    (%r8,%r11,4),%rax
mov    -0x2c(%rsp),%esi
mov    -0x18(%rsp),%r11
mov    -0x20(%rsp),%r14
jmp    1224 <transitive_closure+0x124>
xor    %r14d,%r14d
test   $0x1,%r15b
jne    12e1 <transitive_closure+0x1e1>
jmp    1308 <transitive_closure+0x208>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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