void makecomat1D_2T(int *input1, int *mask1, int xi1, int *input2, int *mask2,
                    int xi2, int *coords, int levels1, int levels2,
                    int *output) {
  int x, xval, i, j;
  for (x = 0; x < xi1; x++) {
    if (mask1[x] == 1) {
      i = input1[x];
      xval = x + coords[0];
      if ((xval >= 0) && (xval < xi1)) {
        if (mask2[xval] == 1) {
          j = input2[xval];
          if (i >= 0 && i < levels1 && j >= 0 && j < levels2) {
            output[i * levels2 + j]++;
          }
        }
      }
    }
  }
}


x86 assembly (clang, -O0):
<makecomat1D_2T>:
push   %rbp
mov    %rsp,%rbp
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1206 <makecomat1D_2T+0x106>
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
cmpl   $0x1,(%rax,%rcx,4)
jne    11f3 <makecomat1D_2T+0xf3>
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x38(%rbp)
mov    -0x30(%rbp),%eax
mov    0x10(%rbp),%rcx
add    (%rcx),%eax
mov    %eax,-0x34(%rbp)
cmpl   $0x0,-0x34(%rbp)
jl     11ee <makecomat1D_2T+0xee>
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11ee <makecomat1D_2T+0xee>
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
cmpl   $0x1,(%rax,%rcx,4)
jne    11e9 <makecomat1D_2T+0xe9>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x3c(%rbp)
cmpl   $0x0,-0x38(%rbp)
jl     11e4 <makecomat1D_2T+0xe4>
mov    -0x38(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    11e4 <makecomat1D_2T+0xe4>
cmpl   $0x0,-0x3c(%rbp)
jl     11e4 <makecomat1D_2T+0xe4>
mov    -0x3c(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    11e4 <makecomat1D_2T+0xe4>
mov    0x28(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   0x20(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%edx
add    $0x1,%edx
mov    %edx,(%rax,%rcx,4)
jmp    11e9 <makecomat1D_2T+0xe9>
jmp    11ee <makecomat1D_2T+0xee>
jmp    11f3 <makecomat1D_2T+0xf3>
jmp    11f8 <makecomat1D_2T+0xf8>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1130 <makecomat1D_2T+0x30>
pop    %rbp
ret

x86 assembly (clang, -O2):
<makecomat1D_2T>:
push   %rbp
push   %r14
push   %rbx
test   %edx,%edx
jle    117d <makecomat1D_2T+0x7d>
mov    0x38(%rsp),%r9
mov    0x30(%rsp),%r11d
mov    0x28(%rsp),%r10d
mov    0x20(%rsp),%r14
mov    %edx,%edx
xor    %ebx,%ebx
jmp    1139 <makecomat1D_2T+0x39>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%rbx
cmp    %rbx,%rdx
je     117d <makecomat1D_2T+0x7d>
cmpl   $0x1,(%rsi,%rbx,4)
jne    1130 <makecomat1D_2T+0x30>
mov    (%r14),%eax
add    %rbx,%rax
test   %eax,%eax
js     1130 <makecomat1D_2T+0x30>
cmp    %edx,%eax
jge    1130 <makecomat1D_2T+0x30>
mov    %eax,%eax
cmpl   $0x1,(%r8,%rax,4)
jne    1130 <makecomat1D_2T+0x30>
mov    (%rdi,%rbx,4),%ebp
test   %ebp,%ebp
js     1130 <makecomat1D_2T+0x30>
cmp    %r10d,%ebp
jge    1130 <makecomat1D_2T+0x30>
mov    (%rcx,%rax,4),%eax
test   %eax,%eax
js     1130 <makecomat1D_2T+0x30>
cmp    %r11d,%eax
jge    1130 <makecomat1D_2T+0x30>
imul   %r11d,%ebp
add    %ebp,%eax
cltq
addl   $0x1,(%r9,%rax,4)
jmp    1130 <makecomat1D_2T+0x30>
pop    %rbx
pop    %r14
pop    %rbp
ret

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

LLVM IR (-O2):
define void @makecomat1D_2T(i32* nocapture %0, i32* nocapture %1, i32 %2, i32* nocapture %3, i32* nocapture %4, i32 %5, i32* nocapture %6, i32 %7, i32 %8, i32* nocapture %9) {
%11 = icmp sgt i32 %2, 0
br i1 %11, label %12, label %53
12: ; preds = %10
%13 = zext i32 %2 to i64
br label %14
14: ; preds = %12, %50
%15 = phi i64 [ 0, %12 ], [ %51, %50 ]
%16 = getelementptr inbounds i32, i32* %1, i64 %15
%17 = load i32, i32* %16
%18 = icmp eq i32 %17, 1
br i1 %18, label %19, label %50
19: ; preds = %14
%20 = getelementptr inbounds i32, i32* %0, i64 %15
%21 = load i32, i32* %20
%22 = load i32, i32* %6
%23 = trunc i64 %15 to i32
%24 = add nsw i32 %22, %23
%25 = icmp sgt i32 %24, -1
%26 = icmp slt i32 %24, %2
%27 = and i1 %25, %26
br i1 %27, label %28, label %50
28: ; preds = %19
%29 = zext i32 %24 to i64
%30 = getelementptr inbounds i32, i32* %4, i64 %29
%31 = load i32, i32* %30
%32 = icmp eq i32 %31, 1
br i1 %32, label %33, label %50
33: ; preds = %28
%34 = getelementptr inbounds i32, i32* %3, i64 %29
%35 = load i32, i32* %34
%36 = icmp sgt i32 %21, -1
br i1 %36, label %37, label %50
37: ; preds = %33
%38 = icmp slt i32 %21, %7
%39 = icmp sgt i32 %35, -1
%40 = select i1 %38, i1 %39, i1 false
%41 = icmp slt i32 %35, %8
%42 = select i1 %40, i1 %41, i1 false
br i1 %42, label %43, label %50
43: ; preds = %37
%44 = mul nsw i32 %21, %8
%45 = add nsw i32 %35, %44
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds i32, i32* %9, i64 %46
%48 = load i32, i32* %47
%49 = add nsw i32 %48, 1
store i32 %49, i32* %47
br label %50
50: ; preds = %14, %28, %43, %37, %33, %19
%51 = add nuw nsw i64 %15, 1
%52 = icmp eq i64 %51, %13
br i1 %52, label %53, label %14
53: ; preds = %50, %10
ret void
}