void makecomat1D(int *input, int *mask, int xi, int *coords, int levels,
                 int *output) {
  int x, xval, i, j;
  for (x = 0; x < xi; x++) {
    if (mask[x] == 1) {
      i = input[x];
      xval = x + coords[0];
      if ((xval >= 0) && (xval < xi)) {
        if (mask[xval] == 1) {
          j = input[xval];
          if (i >= 0 && i < levels && j >= 0 && j < levels) {
            output[i * levels + j]++;
          }
        }
      }
    }
  }
}


x86 assembly (clang, -O0):
<makecomat1D>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
mov    %r9,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11f8 <makecomat1D+0xf8>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
cmpl   $0x1,(%rax,%rcx,4)
jne    11e5 <makecomat1D+0xe5>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x20(%rbp),%rcx
add    (%rcx),%eax
mov    %eax,-0x38(%rbp)
cmpl   $0x0,-0x38(%rbp)
jl     11e0 <makecomat1D+0xe0>
mov    -0x38(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11e0 <makecomat1D+0xe0>
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
cmpl   $0x1,(%rax,%rcx,4)
jne    11db <makecomat1D+0xdb>
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x40(%rbp)
cmpl   $0x0,-0x3c(%rbp)
jl     11d6 <makecomat1D+0xd6>
mov    -0x3c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    11d6 <makecomat1D+0xd6>
cmpl   $0x0,-0x40(%rbp)
jl     11d6 <makecomat1D+0xd6>
mov    -0x40(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    11d6 <makecomat1D+0xd6>
mov    -0x30(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x24(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%edx
add    $0x1,%edx
mov    %edx,(%rax,%rcx,4)
jmp    11db <makecomat1D+0xdb>
jmp    11e0 <makecomat1D+0xe0>
jmp    11e5 <makecomat1D+0xe5>
jmp    11ea <makecomat1D+0xea>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1122 <makecomat1D+0x22>
pop    %rbp
ret

x86 assembly (clang, -O2):
<makecomat1D>:
push   %rbx
test   %edx,%edx
jle    1161 <makecomat1D+0x61>
mov    %edx,%ebx
xor    %edx,%edx
jmp    1119 <makecomat1D+0x19>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rdx
cmp    %rdx,%rbx
je     1161 <makecomat1D+0x61>
cmpl   $0x1,(%rsi,%rdx,4)
jne    1110 <makecomat1D+0x10>
mov    (%rcx),%eax
add    %rdx,%rax
test   %eax,%eax
js     1110 <makecomat1D+0x10>
cmp    %ebx,%eax
jge    1110 <makecomat1D+0x10>
mov    %eax,%r11d
cmpl   $0x1,(%rsi,%r11,4)
jne    1110 <makecomat1D+0x10>
mov    (%rdi,%rdx,4),%r10d
test   %r10d,%r10d
js     1110 <makecomat1D+0x10>
cmp    %r8d,%r10d
jge    1110 <makecomat1D+0x10>
mov    (%rdi,%r11,4),%eax
test   %eax,%eax
js     1110 <makecomat1D+0x10>
cmp    %r8d,%eax
jge    1110 <makecomat1D+0x10>
imul   %r8d,%r10d
add    %r10d,%eax
cltq
addl   $0x1,(%r9,%rax,4)
jmp    1110 <makecomat1D+0x10>
pop    %rbx
ret

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

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