int color_feature_histogram(unsigned char type, int **obj_id, int width,
                            int height, int n_object, int *area,
                            unsigned char *image, double **color_histogram) {
  unsigned int i = 0, j = 0;
  for (i = 0; i < height; i++) {
    for (j = 0; j < width; j++) {
      if (!obj_id[i][j])
        continue;
      register unsigned int idx = (j + i * width) * 3;
      unsigned char val = ((unsigned char)(0.2126 * (double)(image[idx])) +
                           (unsigned char)(0.7152 * (double)(image[idx + 1])) +
                           (unsigned char)(0.0722 * (double)(image[idx + 2])));
      color_histogram[(obj_id[i][j] - 1)][val]++;
    }
  }
  return 1;
}


x86 assembly (clang, -O0):
<color_feature_histogram>:
push   %rbp
mov    %rsp,%rbp
mov    %dil,%al
mov    0x18(%rbp),%rdi
mov    0x10(%rbp),%rdi
mov    %al,-0x1(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
mov    %r9,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
movl   $0x0,-0x30(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jae    124b <color_feature_histogram+0x14b>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jae    1238 <color_feature_histogram+0x138>
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    (%rax,%rcx,8),%rax
mov    -0x30(%rbp),%ecx
cmpl   $0x0,(%rax,%rcx,4)
jne    1175 <color_feature_histogram+0x75>
jmp    122a <color_feature_histogram+0x12a>
mov    -0x30(%rbp),%ecx
mov    -0x2c(%rbp),%eax
mov    -0x14(%rbp),%edx
imul   %edx,%eax
add    %eax,%ecx
mov    %ecx,%eax
lea    (%rax,%rax,2),%eax
mov    %eax,-0x34(%rbp)
mov    0x10(%rbp),%rcx
mov    -0x34(%rbp),%eax
mov    %eax,%edx
movzbl (%rcx,%rax,1),%eax
cvtsi2sd %eax,%xmm0
movsd  0xe64(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
movzbl %al,%eax
mov    %edx,%esi
add    $0x1,%esi
movzbl (%rcx,%rsi,1),%esi
cvtsi2sd %esi,%xmm0
movsd  0xe4c(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%esi
movzbl %sil,%esi
add    %esi,%eax
add    $0x2,%edx
movzbl (%rcx,%rdx,1),%ecx
cvtsi2sd %ecx,%xmm0
movsd  0xe33(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%ecx
movzbl %cl,%ecx
add    %ecx,%eax
mov    %al,-0x35(%rbp)
mov    0x18(%rbp),%rax
mov    -0x10(%rbp),%rcx
mov    -0x2c(%rbp),%edx
mov    (%rcx,%rdx,8),%rcx
mov    -0x30(%rbp),%edx
mov    (%rcx,%rdx,4),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movzbl -0x35(%rbp),%ecx
movsd  0xde0(%rip),%xmm0
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    114c <color_feature_histogram+0x4c>
jmp    123d <color_feature_histogram+0x13d>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1139 <color_feature_histogram+0x39>
mov    $0x1,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<color_feature_histogram>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %ecx,%ecx
je     11e8 <color_feature_histogram+0xe8>
mov    0x40(%rsp),%r10
mov    0x38(%rsp),%r15
mov    %ecx,%r8d
mov    %edx,%r12d
lea    (%rdx,%rdx,2),%r9d
xor    %r11d,%r11d
movsd  0xecf(%rip),%xmm0
movsd  0xecf(%rip),%xmm1
movsd  0xecf(%rip),%xmm2
movsd  0xecf(%rip),%xmm3
xor    %r14d,%r14d
jmp    1160 <color_feature_histogram+0x60>
xchg   %ax,%ax
add    $0x1,%r14
add    %r9d,%r11d
cmp    %r8,%r14
je     11e8 <color_feature_histogram+0xe8>
test   %edx,%edx
je     1150 <color_feature_histogram+0x50>
mov    (%rsi,%r14,8),%r13
mov    %r11d,%eax
xor    %edi,%edi
jmp    117c <color_feature_histogram+0x7c>
nop
add    $0x1,%rdi
add    $0x3,%eax
cmp    %rdi,%r12
je     1150 <color_feature_histogram+0x50>
movslq 0x0(%r13,%rdi,4),%rcx
test   %rcx,%rcx
je     1170 <color_feature_histogram+0x70>
mov    %eax,%ebx
movzbl (%r15,%rbx,1),%ebx
xorps  %xmm4,%xmm4
cvtsi2sd %ebx,%xmm4
mulsd  %xmm0,%xmm4
cvttsd2si %xmm4,%ebx
lea    0x1(%rax),%ebp
movzbl (%r15,%rbp,1),%ebp
xorps  %xmm4,%xmm4
cvtsi2sd %ebp,%xmm4
mulsd  %xmm1,%xmm4
cvttsd2si %xmm4,%ebp
add    %bl,%bpl
lea    0x2(%rax),%ebx
movzbl (%r15,%rbx,1),%ebx
xorps  %xmm4,%xmm4
cvtsi2sd %ebx,%xmm4
mulsd  %xmm2,%xmm4
cvttsd2si %xmm4,%ebx
add    %bpl,%bl
mov    -0x8(%r10,%rcx,8),%rcx
movzbl %bl,%ebx
movsd  (%rcx,%rbx,8),%xmm4
addsd  %xmm3,%xmm4
movsd  %xmm4,(%rcx,%rbx,8)
jmp    1170 <color_feature_histogram+0x70>
mov    $0x1,%eax
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @color_feature_histogram(i8 zeroext %0, i32** %1, i32 %2, i32 %3, i32 %4, i32* %5, i8* %6, double** %7) {
%9 = alloca i8
%10 = alloca i32**
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32*
%15 = alloca i8*
%16 = alloca double**
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca i8
store i8 %0, i8* %9
store i32** %1, i32*** %10
store i32 %2, i32* %11
store i32 %3, i32* %12
store i32 %4, i32* %13
store i32* %5, i32** %14
store i8* %6, i8** %15
store double** %7, double*** %16
store i32 0, i32* %17
store i32 0, i32* %18
store i32 0, i32* %17
br label %21
21: ; preds = %104, %8
%22 = load i32, i32* %17
%23 = load i32, i32* %12
%24 = icmp ult i32 %22, %23
br i1 %24, label %25, label %107
25: ; preds = %21
store i32 0, i32* %18
br label %26
26: ; preds = %100, %25
%27 = load i32, i32* %18
%28 = load i32, i32* %11
%29 = icmp ult i32 %27, %28
br i1 %29, label %30, label %103
30: ; preds = %26
%31 = load i32**, i32*** %10
%32 = load i32, i32* %17
%33 = zext i32 %32 to i64
%34 = getelementptr inbounds i32*, i32** %31, i64 %33
%35 = load i32*, i32** %34
%36 = load i32, i32* %18
%37 = zext i32 %36 to i64
%38 = getelementptr inbounds i32, i32* %35, i64 %37
%39 = load i32, i32* %38
%40 = icmp ne i32 %39, 0
br i1 %40, label %42, label %41
41: ; preds = %30
br label %100
42: ; preds = %30
%43 = load i32, i32* %18
%44 = load i32, i32* %17
%45 = load i32, i32* %11
%46 = mul i32 %44, %45
%47 = add i32 %43, %46
%48 = mul i32 %47, 3
store i32 %48, i32* %19
%49 = load i8*, i8** %15
%50 = load i32, i32* %19
%51 = zext i32 %50 to i64
%52 = getelementptr inbounds i8, i8* %49, i64 %51
%53 = load i8, i8* %52
%54 = uitofp i8 %53 to double
%55 = fmul double 2.126000e-01, %54
%56 = fptoui double %55 to i8
%57 = zext i8 %56 to i32
%58 = load i8*, i8** %15
%59 = load i32, i32* %19
%60 = add i32 %59, 1
%61 = zext i32 %60 to i64
%62 = getelementptr inbounds i8, i8* %58, i64 %61
%63 = load i8, i8* %62
%64 = uitofp i8 %63 to double
%65 = fmul double 7.152000e-01, %64
%66 = fptoui double %65 to i8
%67 = zext i8 %66 to i32
%68 = add nsw i32 %57, %67
%69 = load i8*, i8** %15
%70 = load i32, i32* %19
%71 = add i32 %70, 2
%72 = zext i32 %71 to i64
%73 = getelementptr inbounds i8, i8* %69, i64 %72
%74 = load i8, i8* %73
%75 = uitofp i8 %74 to double
%76 = fmul double 7.220000e-02, %75
%77 = fptoui double %76 to i8
%78 = zext i8 %77 to i32
%79 = add nsw i32 %68, %78
%80 = trunc i32 %79 to i8
store i8 %80, i8* %20
%81 = load double**, double*** %16
%82 = load i32**, i32*** %10
%83 = load i32, i32* %17
%84 = zext i32 %83 to i64
%85 = getelementptr inbounds i32*, i32** %82, i64 %84
%86 = load i32*, i32** %85
%87 = load i32, i32* %18
%88 = zext i32 %87 to i64
%89 = getelementptr inbounds i32, i32* %86, i64 %88
%90 = load i32, i32* %89
%91 = sub nsw i32 %90, 1
%92 = sext i32 %91 to i64
%93 = getelementptr inbounds double*, double** %81, i64 %92
%94 = load double*, double** %93
%95 = load i8, i8* %20
%96 = zext i8 %95 to i64
%97 = getelementptr inbounds double, double* %94, i64 %96
%98 = load double, double* %97
%99 = fadd double %98, 1.000000e+00
store double %99, double* %97
br label %100
100: ; preds = %42, %41
%101 = load i32, i32* %18
%102 = add i32 %101, 1
store i32 %102, i32* %18
br label %26
103: ; preds = %26
br label %104
104: ; preds = %103
%105 = load i32, i32* %17
%106 = add i32 %105, 1
store i32 %106, i32* %17
br label %21
107: ; preds = %21
ret i32 1
}

LLVM IR (-O2):
define i32 @color_feature_histogram(i8 zeroext %0, i32** nocapture %1, i32 %2, i32 %3, i32 %4, i32* nocapture readnone %5, i8* nocapture %6, double** nocapture %7) {
%9 = icmp eq i32 %3, 0
br i1 %9, label %66, label %10
10: ; preds = %8
%11 = icmp eq i32 %2, 0
%12 = zext i32 %3 to i64
%13 = zext i32 %2 to i64
br label %14
14: ; preds = %10, %63
%15 = phi i64 [ 0, %10 ], [ %64, %63 ]
br i1 %11, label %63, label %16
16: ; preds = %14
%17 = getelementptr inbounds i32*, i32** %1, i64 %15
%18 = load i32*, i32** %17
%19 = trunc i64 %15 to i32
%20 = mul i32 %19, %2
br label %21
21: ; preds = %16, %60
%22 = phi i64 [ 0, %16 ], [ %61, %60 ]
%23 = getelementptr inbounds i32, i32* %18, i64 %22
%24 = load i32, i32* %23
%25 = icmp eq i32 %24, 0
br i1 %25, label %60, label %26
26: ; preds = %21
%27 = trunc i64 %22 to i32
%28 = add i32 %20, %27
%29 = mul i32 %28, 3
%30 = zext i32 %29 to i64
%31 = getelementptr inbounds i8, i8* %6, i64 %30
%32 = load i8, i8* %31
%33 = uitofp i8 %32 to double
%34 = fmul double %33, 2.126000e-01
%35 = fptoui double %34 to i8
%36 = add i32 %29, 1
%37 = zext i32 %36 to i64
%38 = getelementptr inbounds i8, i8* %6, i64 %37
%39 = load i8, i8* %38
%40 = uitofp i8 %39 to double
%41 = fmul double %40, 7.152000e-01
%42 = fptoui double %41 to i8
%43 = add i8 %42, %35
%44 = add i32 %29, 2
%45 = zext i32 %44 to i64
%46 = getelementptr inbounds i8, i8* %6, i64 %45
%47 = load i8, i8* %46
%48 = uitofp i8 %47 to double
%49 = fmul double %48, 7.220000e-02
%50 = fptoui double %49 to i8
%51 = add i8 %43, %50
%52 = add nsw i32 %24, -1
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds double*, double** %7, i64 %53
%55 = load double*, double** %54
%56 = zext i8 %51 to i64
%57 = getelementptr inbounds double, double* %55, i64 %56
%58 = load double, double* %57
%59 = fadd double %58, 1.000000e+00
store double %59, double* %57
br label %60
60: ; preds = %21, %26
%61 = add nuw nsw i64 %22, 1
%62 = icmp eq i64 %61, %13
br i1 %62, label %63, label %21
63: ; preds = %60, %14
%64 = add nuw nsw i64 %15, 1
%65 = icmp eq i64 %64, %12
br i1 %65, label %66, label %14
66: ; preds = %63, %8
ret i32 1
}