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 = count; i > 0; --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)
mov    -0x1c(%rbp),%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jle    1293 <compute_row_thresholds+0x193>
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    1158 <compute_row_thresholds+0x58>
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    1176 <compute_row_thresholds+0x76>
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    11b3 <compute_row_thresholds+0xb3>
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    11d1 <compute_row_thresholds+0xd1>
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    11f3 <compute_row_thresholds+0xf3>
mov    -0x27(%rbp),%al
mov    %al,-0x25(%rbp)
movzbl -0x28(%rbp),%eax
movzbl -0x26(%rbp),%ecx
cmp    %ecx,%eax
jle    1209 <compute_row_thresholds+0x109>
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    126c <compute_row_thresholds+0x16c>
movzbl -0x29(%rbp),%eax
cmp    $0x7f,%eax
jge    1255 <compute_row_thresholds+0x155>
mov    -0x8(%rbp),%rax
mov    %rax,%rcx
add    $0x1,%rcx
mov    %rcx,-0x8(%rbp)
movb   $0xff,(%rax)
jmp    1267 <compute_row_thresholds+0x167>
mov    -0x8(%rbp),%rax
mov    %rax,%rcx
add    $0x1,%rcx
mov    %rcx,-0x8(%rbp)
movb   $0x0,(%rax)
jmp    1280 <compute_row_thresholds+0x180>
mov    -0x29(%rbp),%cl
mov    -0x8(%rbp),%rax
mov    %rax,%rdx
add    $0x1,%rdx
mov    %rdx,-0x8(%rbp)
mov    %cl,(%rax)
jmp    1285 <compute_row_thresholds+0x185>
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    111d <compute_row_thresholds+0x1d>
pop    %rbp
ret

x86 assembly (clang, -O2):
<compute_row_thresholds>:
push   %r14
push   %rbx
test   %ecx,%ecx
jle    11b0 <compute_row_thresholds+0xb0>
add    $0x1,%ecx
xor    %r14d,%r14d
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movzbl (%rsi,%r14,2),%r10d
movzbl 0x1(%rsi,%r14,2),%r9d
movzbl 0x2(%rsi,%r14,2),%r11d
cmp    %r10b,%r11b
cmovb  %r11d,%r10d
movzbl 0x3(%rsi,%r14,2),%r11d
cmp    %r9b,%r11b
cmova  %r11d,%r9d
movzbl (%rdx,%r14,2),%eax
movzbl 0x1(%rdx,%r14,2),%ebx
movzbl 0x2(%rdx,%r14,2),%r11d
cmp    %al,%r11b
cmovb  %r11d,%eax
movzbl 0x3(%rdx,%r14,2),%r11d
cmp    %bl,%r11b
cmova  %r11d,%ebx
cmp    %r10b,%al
cmovae %r10d,%eax
cmp    %r9b,%bl
cmovbe %r9d,%ebx
movzbl %al,%r9d
movzbl %bl,%ebx
lea    (%rbx,%r9,1),%r10d
sub    %r9d,%ebx
mov    %r10d,%r9d
shr    %r9d
cmp    $0xfe,%r10d
sbb    %eax,%eax
cmp    %r8d,%ebx
cmovge %r9d,%eax
mov    %al,(%rdi,%r14,1)
add    $0x1,%r14
add    $0xffffffff,%ecx
cmp    $0x1,%ecx
jg     1120 <compute_row_thresholds+0x20>
pop    %rbx
pop    %r14
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
%17 = load i32, i32* %9
store i32 %17, i32* %11
br label %18
18: ; preds = %132, %5
%19 = load i32, i32* %11
%20 = icmp sgt i32 %19, 0
br i1 %20, label %21, label %135
21: ; preds = %18
%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 %18
135: ; preds = %18
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 ], [ %3, %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 nsw i32 %9, -1
%50 = icmp sgt i32 %9, 1
br i1 %50, label %7, label %51
51: ; preds = %7, %5
ret void
}