int threshold_and_find_means(in_image, out_image, hi, low, value, object_mean,
                             background_mean, rows, cols)
short *background_mean, hi, low, **in_image, *object_mean, **out_image, value;
int rows, cols;
{
  int counter = 0, i, j;
  unsigned long object = 0, background = 0;
  for (i = 0; i < rows; i++) {
    for (j = 0; j < cols; j++) {
      if (in_image[i][j] >= low && in_image[i][j] <= hi) {
        out_image[i][j] = value;
        counter++;
        object = object + in_image[i][j];
      } else {
        out_image[i][j] = 0;
        background = background + in_image[i][j];
      }
    }
  }
  object = object / counter;
  background = background / ((rows * cols) - counter);
  *object_mean = (short)(object);
  *background_mean = (short)(background);
  printf("\n\tTAFM> set %d points", counter);
  printf("\n\tTAFM> object=%d background=%d", *object_mean, *background_mean);
  return (1);
}


x86 assembly (clang, -O0):
<threshold_and_find_means>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%rax
mov    %r8w,%ax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %dx,-0x12(%rbp)
mov    %cx,-0x14(%rbp)
mov    %ax,-0x16(%rbp)
mov    %r9,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
movq   $0x0,-0x38(%rbp)
movq   $0x0,-0x40(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    125d <threshold_and_find_means+0x14d>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    124a <threshold_and_find_means+0x13a>
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movswl (%rax,%rcx,2),%eax
movswl -0x14(%rbp),%ecx
cmp    %ecx,%eax
jl     1201 <threshold_and_find_means+0xf1>
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movswl (%rax,%rcx,2),%eax
movswl -0x12(%rbp),%ecx
cmp    %ecx,%eax
jg     1201 <threshold_and_find_means+0xf1>
mov    -0x16(%rbp),%dx
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x38(%rbp),%rax
mov    -0x8(%rbp),%rcx
movslq -0x28(%rbp),%rdx
mov    (%rcx,%rdx,8),%rcx
movslq -0x2c(%rbp),%rdx
movswq (%rcx,%rdx,2),%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
jmp    1237 <threshold_and_find_means+0x127>
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movw   $0x0,(%rax,%rcx,2)
mov    -0x40(%rbp),%rax
mov    -0x8(%rbp),%rcx
movslq -0x28(%rbp),%rdx
mov    (%rcx,%rdx,8),%rcx
movslq -0x2c(%rbp),%rdx
movswq (%rcx,%rdx,2),%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
jmp    123c <threshold_and_find_means+0x12c>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    116f <threshold_and_find_means+0x5f>
jmp    124f <threshold_and_find_means+0x13f>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    115c <threshold_and_find_means+0x4c>
mov    -0x38(%rbp),%rax
movslq -0x24(%rbp),%rcx
xor    %edx,%edx
div    %rcx
mov    %rax,-0x38(%rbp)
mov    -0x40(%rbp),%rax
mov    0x18(%rbp),%ecx
imul   0x20(%rbp),%ecx
sub    -0x24(%rbp),%ecx
movslq %ecx,%rcx
xor    %edx,%edx
div    %rcx
mov    %rax,-0x40(%rbp)
mov    -0x38(%rbp),%rax
mov    %ax,%cx
mov    -0x20(%rbp),%rax
mov    %cx,(%rax)
mov    -0x40(%rbp),%rax
mov    %ax,%cx
mov    0x10(%rbp),%rax
mov    %cx,(%rax)
mov    -0x24(%rbp),%esi
lea    0xd52(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
mov    -0x20(%rbp),%rax
movswl (%rax),%esi
mov    0x10(%rbp),%rax
movswl (%rax),%edx
lea    0xd4c(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
mov    $0x1,%eax
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<threshold_and_find_means>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %r9,(%rsp)
mov    %rsi,0x10(%rsp)
mov    %rdi,0x8(%rsp)
mov    0x60(%rsp),%edi
mov    0x58(%rsp),%eax
test   %eax,%eax
jle    11cc <threshold_and_find_means+0xbc>
movswl %cx,%ebx
movswl %dx,%r11d
mov    %eax,%r15d
mov    %edi,%ecx
xor    %r13d,%r13d
xor    %r12d,%r12d
xor    %eax,%eax
xor    %esi,%esi
jmp    1169 <threshold_and_find_means+0x59>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%r13
cmp    %r15,%r13
je     11d3 <threshold_and_find_means+0xc3>
test   %edi,%edi
jle    1160 <threshold_and_find_means+0x50>
mov    0x8(%rsp),%rdx
mov    (%rdx,%r13,8),%rbp
mov    0x10(%rsp),%rdx
mov    (%rdx,%r13,8),%r10
xor    %r14d,%r14d
jmp    11a9 <threshold_and_find_means+0x99>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movw   $0x0,(%r10,%r14,2)
movswq 0x0(%rbp,%r14,2),%rdx
add    %rdx,%r12
add    $0x1,%r14
cmp    %r14,%rcx
je     1160 <threshold_and_find_means+0x50>
movswl 0x0(%rbp,%r14,2),%r9d
cmp    %r9d,%ebx
jg     1190 <threshold_and_find_means+0x80>
cmp    %r9d,%r11d
jl     1190 <threshold_and_find_means+0x80>
mov    %r8w,(%r10,%r14,2)
add    $0x1,%esi
movswq 0x0(%rbp,%r14,2),%rdx
add    %rdx,%rax
jmp    11a0 <threshold_and_find_means+0x90>
xor    %esi,%esi
xor    %eax,%eax
xor    %r12d,%r12d
movslq %esi,%rcx
mov    %rax,%rdx
or     %rcx,%rdx
shr    $0x20,%rdx
je     11ec <threshold_and_find_means+0xdc>
xor    %edx,%edx
div    %rcx
mov    %rax,%rbp
jmp    11f2 <threshold_and_find_means+0xe2>
xor    %edx,%edx
div    %esi
mov    %eax,%ebp
mov    0x50(%rsp),%r14
imul   0x58(%rsp),%edi
sub    %esi,%edi
movslq %edi,%rcx
mov    %r12,%rax
or     %rcx,%rax
shr    $0x20,%rax
je     1217 <threshold_and_find_means+0x107>
mov    %r12,%rax
xor    %edx,%edx
div    %rcx
jmp    121e <threshold_and_find_means+0x10e>
mov    %r12d,%eax
xor    %edx,%edx
div    %edi
mov    (%rsp),%rbx
mov    %bp,(%rbx)
mov    %ax,(%r14)
lea    0xdd0(%rip),%rdi
xor    %eax,%eax
call   1030 <printf@plt>
movswl (%rbx),%esi
movswl (%r14),%edx
lea    0xdd1(%rip),%rdi
xor    %eax,%eax
call   1030 <printf@plt>
mov    $0x1,%eax
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
@.str = [22 x i8] c"\0A\09TAFM> set %d points\00"
@.str.1 = [32 x i8] c"\0A\09TAFM> object=%d background=%d\00"
define i32 @threshold_and_find_means(i16** nocapture %0, i16** nocapture %1, i32 %2, i32 %3, i32 %4, i16* nocapture %5, i16* nocapture %6, i32 %7, i32 %8) {
%10 = trunc i32 %4 to i16
%11 = icmp sgt i32 %7, 0
br i1 %11, label %12, label %63
12: ; preds = %9
%13 = icmp sgt i32 %8, 0
%14 = shl i32 %3, 16
%15 = ashr exact i32 %14, 16
%16 = shl i32 %2, 16
%17 = ashr exact i32 %16, 16
%18 = zext i32 %7 to i64
%19 = zext i32 %8 to i64
br label %20
20: ; preds = %12, %57
%21 = phi i64 [ 0, %12 ], [ %61, %57 ]
%22 = phi i64 [ 0, %12 ], [ %60, %57 ]
%23 = phi i64 [ 0, %12 ], [ %59, %57 ]
%24 = phi i32 [ 0, %12 ], [ %58, %57 ]
br i1 %13, label %25, label %57
25: ; preds = %20
%26 = getelementptr inbounds i16*, i16** %0, i64 %21
%27 = load i16*, i16** %26
%28 = getelementptr inbounds i16*, i16** %1, i64 %21
%29 = load i16*, i16** %28
br label %30
30: ; preds = %25, %51
%31 = phi i64 [ 0, %25 ], [ %55, %51 ]
%32 = phi i64 [ %22, %25 ], [ %54, %51 ]
%33 = phi i64 [ %23, %25 ], [ %53, %51 ]
%34 = phi i32 [ %24, %25 ], [ %52, %51 ]
%35 = getelementptr inbounds i16, i16* %27, i64 %31
%36 = load i16, i16* %35
%37 = sext i16 %36 to i32
%38 = icmp sgt i32 %15, %37
%39 = icmp slt i32 %17, %37
%40 = select i1 %38, i1 true, i1 %39
%41 = getelementptr inbounds i16, i16* %29, i64 %31
br i1 %40, label %47, label %42
42: ; preds = %30
store i16 %10, i16* %41
%43 = add nsw i32 %34, 1
%44 = load i16, i16* %35
%45 = sext i16 %44 to i64
%46 = add i64 %33, %45
br label %51
47: ; preds = %30
store i16 0, i16* %41
%48 = load i16, i16* %35
%49 = sext i16 %48 to i64
%50 = add i64 %32, %49
br label %51
51: ; preds = %42, %47
%52 = phi i32 [ %43, %42 ], [ %34, %47 ]
%53 = phi i64 [ %46, %42 ], [ %33, %47 ]
%54 = phi i64 [ %32, %42 ], [ %50, %47 ]
%55 = add nuw nsw i64 %31, 1
%56 = icmp eq i64 %55, %19
br i1 %56, label %57, label %30
57: ; preds = %51, %20
%58 = phi i32 [ %24, %20 ], [ %52, %51 ]
%59 = phi i64 [ %23, %20 ], [ %53, %51 ]
%60 = phi i64 [ %22, %20 ], [ %54, %51 ]
%61 = add nuw nsw i64 %21, 1
%62 = icmp eq i64 %61, %18
br i1 %62, label %63, label %20
63: ; preds = %57, %9
%64 = phi i32 [ 0, %9 ], [ %58, %57 ]
%65 = phi i64 [ 0, %9 ], [ %59, %57 ]
%66 = phi i64 [ 0, %9 ], [ %60, %57 ]
%67 = sext i32 %64 to i64
%68 = udiv i64 %65, %67
%69 = mul nsw i32 %8, %7
%70 = sub nsw i32 %69, %64
%71 = sext i32 %70 to i64
%72 = udiv i64 %66, %71
%73 = trunc i64 %68 to i16
store i16 %73, i16* %5
%74 = trunc i64 %72 to i16
store i16 %74, i16* %6
%75 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([22 x i8], [22 x i8]* @.str, i64 0, i64 0), i32 %64)
%76 = load i16, i16* %5
%77 = sext i16 %76 to i32
%78 = load i16, i16* %6
%79 = sext i16 %78 to i32
%80 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([32 x i8], [32 x i8]* @.str.1, i64 0, i64 0), i32 %77, i32 %79)
ret i32 1
}