unsigned extract_collapse_mask(int *iy, int N, int B) {
  unsigned collapse_mask;
  int N0;
  int i;
  if (B <= 1)
    return 1;
  N0 = N / B;
  collapse_mask = 0;
  i = 0;
  do {
    int j;
    j = 0;
    do {
      collapse_mask |= (iy[i * N0 + j] != 0) << i;
    } while (++j < N0);
  } while (++i < B);
  return collapse_mask;
}


x86 assembly (clang, -O0):
<extract_collapse_mask>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %edx,-0x18(%rbp)
cmpl   $0x1,-0x18(%rbp)
jg     1124 <extract_collapse_mask+0x24>
movl   $0x1,-0x4(%rbp)
jmp    119a <extract_collapse_mask+0x9a>
mov    -0x14(%rbp),%eax
cltd
idivl  -0x18(%rbp)
mov    %eax,-0x20(%rbp)
movl   $0x0,-0x1c(%rbp)
movl   $0x0,-0x24(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x20(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
cmpl   $0x0,(%rax,%rcx,4)
setne  %al
and    $0x1,%al
movzbl %al,%eax
mov    -0x24(%rbp),%ecx
shl    %cl,%eax
or     -0x1c(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
cmp    -0x20(%rbp),%eax
jl     1143 <extract_collapse_mask+0x43>
jmp    1182 <extract_collapse_mask+0x82>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
cmp    -0x18(%rbp),%eax
jl     113c <extract_collapse_mask+0x3c>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<extract_collapse_mask>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    $0x1,%eax
cmp    $0x2,%edx
jge    111f <extract_collapse_mask+0x1f>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %edx,%ecx
mov    %esi,%eax
cltd
idiv   %ecx
cmp    $0x2,%eax
mov    $0x1,%r13d
cmovge %eax,%r13d
movslq %eax,%r8
mov    %ecx,%r9d
mov    %r13d,%r15d
and    $0x7ffffff8,%r15d
lea    -0x8(%r15),%r10
mov    %r10,%r11
shr    $0x3,%r11
add    $0x1,%r11
mov    %r11,%r14
and    $0xfffffffffffffffe,%r14
lea    0x30(%rdi),%rbx
lea    0x0(,%r8,4),%r12
xor    %ecx,%ecx
pxor   %xmm8,%xmm8
movdqa 0xe8c(%rip),%xmm1
mov    %rdi,%rsi
xor    %eax,%eax
jmp    118f <extract_collapse_mask+0x8f>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rcx
add    %r12,%rbx
add    %r12,%rsi
cmp    %r9,%rcx
je     1114 <extract_collapse_mask+0x14>
cmp    $0x8,%r13d
jae    11a0 <extract_collapse_mask+0xa0>
xor    %edx,%edx
jmp    12a0 <extract_collapse_mask+0x1a0>
nopl   0x0(%rax)
movd   %eax,%xmm2
movd   %ecx,%xmm0
xorps  %xmm3,%xmm3
movss  %xmm0,%xmm3
test   %r10,%r10
je     12bc <extract_collapse_mask+0x1bc>
pxor   %xmm5,%xmm5
mov    %r14,%rdx
xor    %eax,%eax
movdqa %xmm2,%xmm0
cs nopw 0x0(%rax,%rax,1)
nop
movdqu -0x30(%rbx,%rax,4),%xmm7
movdqu -0x20(%rbx,%rax,4),%xmm6
movdqu -0x10(%rbx,%rax,4),%xmm2
movdqu (%rbx,%rax,4),%xmm4
pcmpeqd %xmm8,%xmm7
pandn  %xmm1,%xmm7
pcmpeqd %xmm8,%xmm6
pandn  %xmm1,%xmm6
pslld  %xmm3,%xmm7
por    %xmm0,%xmm7
pslld  %xmm3,%xmm6
por    %xmm5,%xmm6
pcmpeqd %xmm8,%xmm2
pandn  %xmm1,%xmm2
pcmpeqd %xmm8,%xmm4
pandn  %xmm1,%xmm4
pslld  %xmm3,%xmm2
pslld  %xmm3,%xmm4
por    %xmm7,%xmm2
por    %xmm6,%xmm4
add    $0x10,%rax
movdqa %xmm2,%xmm0
movdqa %xmm4,%xmm5
add    $0xfffffffffffffffe,%rdx
jne    11d0 <extract_collapse_mask+0xd0>
test   $0x1,%r11b
je     127a <extract_collapse_mask+0x17a>
mov    %rcx,%rdx
imul   %r8,%rdx
add    %rdx,%rax
movdqu (%rdi,%rax,4),%xmm0
movdqu 0x10(%rdi,%rax,4),%xmm5
pcmpeqd %xmm8,%xmm0
pandn  %xmm1,%xmm0
pcmpeqd %xmm8,%xmm5
pandn  %xmm1,%xmm5
pslld  %xmm3,%xmm0
por    %xmm0,%xmm2
pslld  %xmm3,%xmm5
por    %xmm5,%xmm4
por    %xmm4,%xmm2
pshufd $0xee,%xmm2,%xmm0
por    %xmm2,%xmm0
pshufd $0x55,%xmm0,%xmm2
por    %xmm0,%xmm2
movd   %xmm2,%eax
mov    %r15,%rdx
cmp    %r13,%r15
je     1180 <extract_collapse_mask+0x80>
xor    %ebp,%ebp
cmpl   $0x0,(%rsi,%rdx,4)
setne  %bpl
shl    %cl,%ebp
or     %ebp,%eax
add    $0x1,%rdx
cmp    %rdx,%r13
jne    12a0 <extract_collapse_mask+0x1a0>
jmp    1180 <extract_collapse_mask+0x80>
pxor   %xmm4,%xmm4
xor    %eax,%eax
test   $0x1,%r11b
jne    1243 <extract_collapse_mask+0x143>
jmp    127a <extract_collapse_mask+0x17a>

LLVM IR (-O0):
define i32 @extract_collapse_mask(i32* %0, i32 %1, i32 %2) {
%4 = alloca i32
%5 = alloca i32*
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
store i32* %0, i32** %5
store i32 %1, i32* %6
store i32 %2, i32* %7
%12 = load i32, i32* %7
%13 = icmp sle i32 %12, 1
br i1 %13, label %14, label %15
14: ; preds = %3
store i32 1, i32* %4
br label %49
15: ; preds = %3
%16 = load i32, i32* %6
%17 = load i32, i32* %7
%18 = sdiv i32 %16, %17
store i32 %18, i32* %9
store i32 0, i32* %8
store i32 0, i32* %10
br label %19
19: ; preds = %42, %15
store i32 0, i32* %11
br label %20
20: ; preds = %36, %19
%21 = load i32*, i32** %5
%22 = load i32, i32* %10
%23 = load i32, i32* %9
%24 = mul nsw i32 %22, %23
%25 = load i32, i32* %11
%26 = add nsw i32 %24, %25
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds i32, i32* %21, i64 %27
%29 = load i32, i32* %28
%30 = icmp ne i32 %29, 0
%31 = zext i1 %30 to i32
%32 = load i32, i32* %10
%33 = shl i32 %31, %32
%34 = load i32, i32* %8
%35 = or i32 %34, %33
store i32 %35, i32* %8
br label %36
36: ; preds = %20
%37 = load i32, i32* %11
%38 = add nsw i32 %37, 1
store i32 %38, i32* %11
%39 = load i32, i32* %9
%40 = icmp slt i32 %38, %39
br i1 %40, label %20, label %41
41: ; preds = %36
br label %42
42: ; preds = %41
%43 = load i32, i32* %10
%44 = add nsw i32 %43, 1
store i32 %44, i32* %10
%45 = load i32, i32* %7
%46 = icmp slt i32 %44, %45
br i1 %46, label %19, label %47
47: ; preds = %42
%48 = load i32, i32* %8
store i32 %48, i32* %4
br label %49
49: ; preds = %47, %14
%50 = load i32, i32* %4
ret i32 %50
}

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