void compute_row_thresholds(unsigned char *thresholds_dest,
                            const unsigned char *min_max_a,
                            const unsigned char *min_max_b, int count,
                            int contrast_threshold) {
  int i;
  unsigned char min_a, max_a, min_b, max_b;
  unsigned char test_threshold;
  for (i = 0; i < count; ++i) {
    min_a = min_max_a[0];
    max_a = min_max_a[1];
    if (min_max_a[2] < min_a)
      min_a = min_max_a[2];
    if (min_max_a[3] > max_a)
      max_a = min_max_a[3];
    min_max_a += 2;
    min_b = min_max_b[0];
    max_b = min_max_b[1];
    if (min_max_b[2] < min_b)
      min_b = min_max_b[2];
    if (min_max_b[3] > max_b)
      max_b = min_max_b[3];
    min_max_b += 2;
    if (min_b < min_a)
      min_a = min_b;
    if (max_b > max_a)
      max_a = max_b;
    test_threshold = (unsigned char)((min_a + max_a) / 2);
    if ((max_a - min_a) < contrast_threshold) {
      if (test_threshold < 127)
        *thresholds_dest++ = ((unsigned char)255);
      else
        *thresholds_dest++ = ((unsigned char)0);
    } else {
      *thresholds_dest++ = test_threshold;
    }
  }
}


x86 assembly (clang, -O0):
<compute_row_thresholds>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1296 <compute_row_thresholds+0x196>
mov    -0x10(%rbp),%rax
mov    (%rax),%al
mov    %al,-0x25(%rbp)
mov    -0x10(%rbp),%rax
mov    0x1(%rax),%al
mov    %al,-0x26(%rbp)
mov    -0x10(%rbp),%rax
movzbl 0x2(%rax),%eax
movzbl -0x25(%rbp),%ecx
cmp    %ecx,%eax
jge    115b <compute_row_thresholds+0x5b>
mov    -0x10(%rbp),%rax
mov    0x2(%rax),%al
mov    %al,-0x25(%rbp)
mov    -0x10(%rbp),%rax
movzbl 0x3(%rax),%eax
movzbl -0x26(%rbp),%ecx
cmp    %ecx,%eax
jle    1179 <compute_row_thresholds+0x79>
mov    -0x10(%rbp),%rax
mov    0x3(%rax),%al
mov    %al,-0x26(%rbp)
mov    -0x10(%rbp),%rax
add    $0x2,%rax
mov    %rax,-0x10(%rbp)
mov    -0x18(%rbp),%rax
mov    (%rax),%al
mov    %al,-0x27(%rbp)
mov    -0x18(%rbp),%rax
mov    0x1(%rax),%al
mov    %al,-0x28(%rbp)
mov    -0x18(%rbp),%rax
movzbl 0x2(%rax),%eax
movzbl -0x27(%rbp),%ecx
cmp    %ecx,%eax
jge    11b6 <compute_row_thresholds+0xb6>
mov    -0x18(%rbp),%rax
mov    0x2(%rax),%al
mov    %al,-0x27(%rbp)
mov    -0x18(%rbp),%rax
movzbl 0x3(%rax),%eax
movzbl -0x28(%rbp),%ecx
cmp    %ecx,%eax
jle    11d4 <compute_row_thresholds+0xd4>
mov    -0x18(%rbp),%rax
mov    0x3(%rax),%al
mov    %al,-0x28(%rbp)
mov    -0x18(%rbp),%rax
add    $0x2,%rax
mov    %rax,-0x18(%rbp)
movzbl -0x27(%rbp),%eax
movzbl -0x25(%rbp),%ecx
cmp    %ecx,%eax
jge    11f6 <compute_row_thresholds+0xf6>
mov    -0x27(%rbp),%al
mov    %al,-0x25(%rbp)
movzbl -0x28(%rbp),%eax
movzbl -0x26(%rbp),%ecx
cmp    %ecx,%eax
jle    120c <compute_row_thresholds+0x10c>
mov    -0x28(%rbp),%al
mov    %al,-0x26(%rbp)
movzbl -0x25(%rbp),%eax
movzbl -0x26(%rbp),%ecx
add    %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %al,-0x29(%rbp)
movzbl -0x26(%rbp),%eax
movzbl -0x25(%rbp),%ecx
sub    %ecx,%eax
cmp    -0x20(%rbp),%eax
jge    126f <compute_row_thresholds+0x16f>
movzbl -0x29(%rbp),%eax
cmp    $0x7f,%eax
jge    1258 <compute_row_thresholds+0x158>
mov    -0x8(%rbp),%rax
mov    %rax,%rcx
add    $0x1,%rcx
mov    %rcx,-0x8(%rbp)
movb   $0xff,(%rax)
jmp    126a <compute_row_thresholds+0x16a>
mov    -0x8(%rbp),%rax
mov    %rax,%rcx
add    $0x1,%rcx
mov    %rcx,-0x8(%rbp)
movb   $0x0,(%rax)
jmp    1283 <compute_row_thresholds+0x183>
mov    -0x29(%rbp),%cl
mov    -0x8(%rbp),%rax
mov    %rax,%rdx
add    $0x1,%rdx
mov    %rdx,-0x8(%rbp)
mov    %cl,(%rax)
jmp    1288 <compute_row_thresholds+0x188>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    111e <compute_row_thresholds+0x1e>
pop    %rbp
ret

x86 assembly (clang, -O2):
<compute_row_thresholds>:
push   %rbp
push   %rbx
test   %ecx,%ecx
jle    1189 <compute_row_thresholds+0x89>
mov    %ecx,%r9d
xor    %ecx,%ecx
nop
movzbl (%rsi,%rcx,2),%r11d
movzbl 0x1(%rsi,%rcx,2),%r10d
movzbl 0x2(%rsi,%rcx,2),%eax
cmp    %r11b,%al
cmovb  %eax,%r11d
movzbl 0x3(%rsi,%rcx,2),%eax
cmp    %r10b,%al
cmova  %eax,%r10d
movzbl (%rdx,%rcx,2),%eax
movzbl 0x1(%rdx,%rcx,2),%ebx
movzbl 0x2(%rdx,%rcx,2),%ebp
cmp    %al,%bpl
cmovb  %ebp,%eax
movzbl 0x3(%rdx,%rcx,2),%ebp
cmp    %bl,%bpl
cmova  %ebp,%ebx
cmp    %r11b,%al
cmovae %r11d,%eax
cmp    %r10b,%bl
cmovbe %r10d,%ebx
movzbl %al,%eax
movzbl %bl,%ebp
lea    (%rax,%rbp,1),%ebx
sub    %eax,%ebp
mov    %ebx,%eax
shr    %eax
cmp    $0xfe,%ebx
sbb    %ebx,%ebx
cmp    %r8d,%ebp
cmovge %eax,%ebx
mov    %bl,(%rdi,%rcx,1)
add    $0x1,%rcx
cmp    %ecx,%r9d
jne    1110 <compute_row_thresholds+0x10>
pop    %rbx
pop    %rbp
ret

LLVM IR (-O0):
define void @compute_row_thresholds(i8* %0, i8* %1, i8* %2, i32 %3, i32 %4) {
%6 = alloca i8*
%7 = alloca i8*
%8 = alloca i8*
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i8
%13 = alloca i8
%14 = alloca i8
%15 = alloca i8
%16 = alloca i8
store i8* %0, i8** %6
store i8* %1, i8** %7
store i8* %2, i8** %8
store i32 %3, i32* %9
store i32 %4, i32* %10
store i32 0, i32* %11
br label %17
17: ; preds = %132, %5
%18 = load i32, i32* %11
%19 = load i32, i32* %9
%20 = icmp slt i32 %18, %19
br i1 %20, label %21, label %135
21: ; preds = %17
%22 = load i8*, i8** %7
%23 = getelementptr inbounds i8, i8* %22, i64 0
%24 = load i8, i8* %23
store i8 %24, i8* %12
%25 = load i8*, i8** %7
%26 = getelementptr inbounds i8, i8* %25, i64 1
%27 = load i8, i8* %26
store i8 %27, i8* %13
%28 = load i8*, i8** %7
%29 = getelementptr inbounds i8, i8* %28, i64 2
%30 = load i8, i8* %29
%31 = zext i8 %30 to i32
%32 = load i8, i8* %12
%33 = zext i8 %32 to i32
%34 = icmp slt i32 %31, %33
br i1 %34, label %35, label %39
35: ; preds = %21
%36 = load i8*, i8** %7
%37 = getelementptr inbounds i8, i8* %36, i64 2
%38 = load i8, i8* %37
store i8 %38, i8* %12
br label %39
39: ; preds = %35, %21
%40 = load i8*, i8** %7
%41 = getelementptr inbounds i8, i8* %40, i64 3
%42 = load i8, i8* %41
%43 = zext i8 %42 to i32
%44 = load i8, i8* %13
%45 = zext i8 %44 to i32
%46 = icmp sgt i32 %43, %45
br i1 %46, label %47, label %51
47: ; preds = %39
%48 = load i8*, i8** %7
%49 = getelementptr inbounds i8, i8* %48, i64 3
%50 = load i8, i8* %49
store i8 %50, i8* %13
br label %51
51: ; preds = %47, %39
%52 = load i8*, i8** %7
%53 = getelementptr inbounds i8, i8* %52, i64 2
store i8* %53, i8** %7
%54 = load i8*, i8** %8
%55 = getelementptr inbounds i8, i8* %54, i64 0
%56 = load i8, i8* %55
store i8 %56, i8* %14
%57 = load i8*, i8** %8
%58 = getelementptr inbounds i8, i8* %57, i64 1
%59 = load i8, i8* %58
store i8 %59, i8* %15
%60 = load i8*, i8** %8
%61 = getelementptr inbounds i8, i8* %60, i64 2
%62 = load i8, i8* %61
%63 = zext i8 %62 to i32
%64 = load i8, i8* %14
%65 = zext i8 %64 to i32
%66 = icmp slt i32 %63, %65
br i1 %66, label %67, label %71
67: ; preds = %51
%68 = load i8*, i8** %8
%69 = getelementptr inbounds i8, i8* %68, i64 2
%70 = load i8, i8* %69
store i8 %70, i8* %14
br label %71
71: ; preds = %67, %51
%72 = load i8*, i8** %8
%73 = getelementptr inbounds i8, i8* %72, i64 3
%74 = load i8, i8* %73
%75 = zext i8 %74 to i32
%76 = load i8, i8* %15
%77 = zext i8 %76 to i32
%78 = icmp sgt i32 %75, %77
br i1 %78, label %79, label %83
79: ; preds = %71
%80 = load i8*, i8** %8
%81 = getelementptr inbounds i8, i8* %80, i64 3
%82 = load i8, i8* %81
store i8 %82, i8* %15
br label %83
83: ; preds = %79, %71
%84 = load i8*, i8** %8
%85 = getelementptr inbounds i8, i8* %84, i64 2
store i8* %85, i8** %8
%86 = load i8, i8* %14
%87 = zext i8 %86 to i32
%88 = load i8, i8* %12
%89 = zext i8 %88 to i32
%90 = icmp slt i32 %87, %89
br i1 %90, label %91, label %93
91: ; preds = %83
%92 = load i8, i8* %14
store i8 %92, i8* %12
br label %93
93: ; preds = %91, %83
%94 = load i8, i8* %15
%95 = zext i8 %94 to i32
%96 = load i8, i8* %13
%97 = zext i8 %96 to i32
%98 = icmp sgt i32 %95, %97
br i1 %98, label %99, label %101
99: ; preds = %93
%100 = load i8, i8* %15
store i8 %100, i8* %13
br label %101
101: ; preds = %99, %93
%102 = load i8, i8* %12
%103 = zext i8 %102 to i32
%104 = load i8, i8* %13
%105 = zext i8 %104 to i32
%106 = add nsw i32 %103, %105
%107 = sdiv i32 %106, 2
%108 = trunc i32 %107 to i8
store i8 %108, i8* %16
%109 = load i8, i8* %13
%110 = zext i8 %109 to i32
%111 = load i8, i8* %12
%112 = zext i8 %111 to i32
%113 = sub nsw i32 %110, %112
%114 = load i32, i32* %10
%115 = icmp slt i32 %113, %114
br i1 %115, label %116, label %127
116: ; preds = %101
%117 = load i8, i8* %16
%118 = zext i8 %117 to i32
%119 = icmp slt i32 %118, 127
br i1 %119, label %120, label %123
120: ; preds = %116
%121 = load i8*, i8** %6
%122 = getelementptr inbounds i8, i8* %121, i32 1
store i8* %122, i8** %6
store i8 -1, i8* %121
br label %126
123: ; preds = %116
%124 = load i8*, i8** %6
%125 = getelementptr inbounds i8, i8* %124, i32 1
store i8* %125, i8** %6
store i8 0, i8* %124
br label %126
126: ; preds = %123, %120
br label %131
127: ; preds = %101
%128 = load i8, i8* %16
%129 = load i8*, i8** %6
%130 = getelementptr inbounds i8, i8* %129, i32 1
store i8* %130, i8** %6
store i8 %128, i8* %129
br label %131
131: ; preds = %127, %126
br label %132
132: ; preds = %131
%133 = load i32, i32* %11
%134 = add nsw i32 %133, 1
store i32 %134, i32* %11
br label %17
135: ; preds = %17
ret void
}

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