void aan_ecualizar_histograma(float *h, float *e, int *f) {
  int k;
  float sum_e = 0.0;
  for (k = 0; k < 256; k++) {
    int i;
    float sum_h = 0.0;
    sum_e = sum_e + e[k];
    for (i = 0; i < 256; i++) {
      if ((sum_h + h[i]) > sum_e)
        break;
      sum_h = sum_h + h[i];
    }
    f[k] = i - 1;
  }
}


x86 assembly (clang, -O0):
<aan_ecualizar_histograma>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
movl   $0x0,-0x1c(%rbp)
cmpl   $0x100,-0x1c(%rbp)
jge    11c4 <aan_ecualizar_histograma+0xc4>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x20(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
cmpl   $0x100,-0x24(%rbp)
jge    11a5 <aan_ecualizar_histograma+0xa5>
movss  -0x28(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
ucomiss -0x20(%rbp),%xmm0
jbe    1180 <aan_ecualizar_histograma+0x80>
jmp    11a5 <aan_ecualizar_histograma+0xa5>
movss  -0x28(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x28(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1152 <aan_ecualizar_histograma+0x52>
mov    -0x24(%rbp),%edx
sub    $0x1,%edx
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    111f <aan_ecualizar_histograma+0x1f>
pop    %rbp
ret

x86 assembly (clang, -O2):
<aan_ecualizar_histograma>:
xorps  %xmm0,%xmm0
xor    %eax,%eax
jmp    1126 <aan_ecualizar_histograma+0x26>
nopw   0x0(%rax,%rax,1)
add    $0x1,%rcx
add    $0xffffffff,%ecx
mov    %ecx,(%rdx,%rax,4)
add    $0x1,%rax
cmp    $0x100,%rax
je     117c <aan_ecualizar_histograma+0x7c>
addss  (%rsi,%rax,4),%xmm0
xor    %ecx,%ecx
xorps  %xmm1,%xmm1
addss  (%rdi,%rcx,4),%xmm1
ucomiss %xmm0,%xmm1
ja     1114 <aan_ecualizar_histograma+0x14>
addss  0x4(%rdi,%rcx,4),%xmm1
ucomiss %xmm0,%xmm1
ja     1110 <aan_ecualizar_histograma+0x10>
addss  0x8(%rdi,%rcx,4),%xmm1
ucomiss %xmm0,%xmm1
ja     1170 <aan_ecualizar_histograma+0x70>
addss  0xc(%rdi,%rcx,4),%xmm1
ucomiss %xmm0,%xmm1
ja     1176 <aan_ecualizar_histograma+0x76>
add    $0x4,%rcx
cmp    $0x100,%rcx
jne    1130 <aan_ecualizar_histograma+0x30>
mov    $0x100,%ecx
jmp    1114 <aan_ecualizar_histograma+0x14>
nop
add    $0x2,%rcx
jmp    1114 <aan_ecualizar_histograma+0x14>
add    $0x3,%rcx
jmp    1114 <aan_ecualizar_histograma+0x14>
ret

LLVM IR (-O0):
define void @aan_ecualizar_histograma(float* %0, float* %1, i32* %2) {
%4 = alloca float*
%5 = alloca float*
%6 = alloca i32*
%7 = alloca i32
%8 = alloca float
%9 = alloca i32
%10 = alloca float
store float* %0, float** %4
store float* %1, float** %5
store i32* %2, i32** %6
store float 0.000000e+00, float* %8
store i32 0, i32* %7
br label %11
11: ; preds = %54, %3
%12 = load i32, i32* %7
%13 = icmp slt i32 %12, 256
br i1 %13, label %14, label %57
14: ; preds = %11
store float 0.000000e+00, float* %10
%15 = load float, float* %8
%16 = load float*, float** %5
%17 = load i32, i32* %7
%18 = sext i32 %17 to i64
%19 = getelementptr inbounds float, float* %16, i64 %18
%20 = load float, float* %19
%21 = fadd float %15, %20
store float %21, float* %8
store i32 0, i32* %9
br label %22
22: ; preds = %44, %14
%23 = load i32, i32* %9
%24 = icmp slt i32 %23, 256
br i1 %24, label %25, label %47
25: ; preds = %22
%26 = load float, float* %10
%27 = load float*, float** %4
%28 = load i32, i32* %9
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds float, float* %27, i64 %29
%31 = load float, float* %30
%32 = fadd float %26, %31
%33 = load float, float* %8
%34 = fcmp ogt float %32, %33
br i1 %34, label %35, label %36
35: ; preds = %25
br label %47
36: ; preds = %25
%37 = load float, float* %10
%38 = load float*, float** %4
%39 = load i32, i32* %9
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds float, float* %38, i64 %40
%42 = load float, float* %41
%43 = fadd float %37, %42
store float %43, float* %10
br label %44
44: ; preds = %36
%45 = load i32, i32* %9
%46 = add nsw i32 %45, 1
store i32 %46, i32* %9
br label %22
47: ; preds = %35, %22
%48 = load i32, i32* %9
%49 = sub nsw i32 %48, 1
%50 = load i32*, i32** %6
%51 = load i32, i32* %7
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds i32, i32* %50, i64 %52
store i32 %49, i32* %53
br label %54
54: ; preds = %47
%55 = load i32, i32* %7
%56 = add nsw i32 %55, 1
store i32 %56, i32* %7
br label %11
57: ; preds = %11
ret void
}

LLVM IR (-O2):
define void @aan_ecualizar_histograma(float* nocapture %0, float* nocapture %1, i32* nocapture writeonly %2) {
br label %4
4: ; preds = %3, %41
%5 = phi i64 [ 0, %3 ], [ %45, %41 ]
%6 = phi float [ 0.000000e+00, %3 ], [ %9, %41 ]
%7 = getelementptr inbounds float, float* %1, i64 %5
%8 = load float, float* %7
%9 = fadd float %6, %8
br label %10
10: ; preds = %35, %4
%11 = phi i64 [ 0, %4 ], [ %36, %35 ]
%12 = phi float [ 0.000000e+00, %4 ], [ %33, %35 ]
%13 = getelementptr inbounds float, float* %0, i64 %11
%14 = load float, float* %13
%15 = fadd float %12, %14
%16 = fcmp ogt float %15, %9
br i1 %16, label %38, label %17
17: ; preds = %10
%18 = or i64 %11, 1
%19 = getelementptr inbounds float, float* %0, i64 %18
%20 = load float, float* %19
%21 = fadd float %15, %20
%22 = fcmp ogt float %21, %9
br i1 %22, label %38, label %23
23: ; preds = %17
%24 = or i64 %11, 2
%25 = getelementptr inbounds float, float* %0, i64 %24
%26 = load float, float* %25
%27 = fadd float %21, %26
%28 = fcmp ogt float %27, %9
br i1 %28, label %38, label %29
29: ; preds = %23
%30 = or i64 %11, 3
%31 = getelementptr inbounds float, float* %0, i64 %30
%32 = load float, float* %31
%33 = fadd float %27, %32
%34 = fcmp ogt float %33, %9
br i1 %34, label %38, label %35
35: ; preds = %29
%36 = add nuw nsw i64 %11, 4
%37 = icmp eq i64 %36, 256
br i1 %37, label %41, label %10
38: ; preds = %29, %23, %17, %10
%39 = phi i64 [ %11, %10 ], [ %18, %17 ], [ %24, %23 ], [ %30, %29 ]
%40 = trunc i64 %39 to i32
br label %41
41: ; preds = %35, %38
%42 = phi i32 [ %40, %38 ], [ 256, %35 ]
%43 = add nsw i32 %42, -1
%44 = getelementptr inbounds i32, i32* %2, i64 %5
store i32 %43, i32* %44
%45 = add nuw nsw i64 %5, 1
%46 = icmp eq i64 %45, 256
br i1 %46, label %47, label %4
47: ; preds = %41
ret void
}