#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */
typedef int uint32_t;
typedef int uint16_t;

/* Variables and functions */
void ranged_hist_float(const char *image, uint16_t rows, uint16_t cols,
                       uint32_t r_stride, uint32_t c_stride,
                       uint32_t *histogram, uint16_t n_bins, float hist_min,
                       float hist_max) {
  const char *row_start, *pixel;
  float bin_factor = (float)n_bins / (hist_max - hist_min);
  uint32_t *last_bin = histogram + n_bins - 1;
  for (row_start = image; row_start != image + rows * r_stride;
       row_start += r_stride) {
    for (pixel = row_start; pixel != row_start + cols * c_stride;
         pixel += c_stride) {
      float val = *(float *)pixel;
      if (val >= hist_min && val < hist_max)
        histogram[(uint16_t)(bin_factor * (val - hist_min))]++;
      else if (val == hist_max)
        (*last_bin)++;
    }
  }
}


x86 assembly (clang, -O0):
<ranged_hist_float>:
push   %rbp
mov    %rsp,%rbp
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8d,-0x18(%rbp)
mov    %r9,-0x20(%rbp)
movss  %xmm0,-0x24(%rbp)
movss  %xmm1,-0x28(%rbp)
cvtsi2ssl 0x10(%rbp),%xmm0
movss  -0x28(%rbp),%xmm1
subss  -0x24(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x3c(%rbp)
mov    -0x20(%rbp),%rax
movslq 0x10(%rbp),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
add    $0xfffffffffffffffc,%rax
mov    %rax,-0x48(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x30(%rbp),%rax
mov    -0x8(%rbp),%rcx
mov    -0xc(%rbp),%edx
imul   -0x14(%rbp),%edx
movslq %edx,%rdx
add    %rdx,%rcx
cmp    %rcx,%rax
je     1257 <ranged_hist_float+0x157>
mov    -0x30(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x38(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    -0x10(%rbp),%edx
imul   -0x18(%rbp),%edx
movslq %edx,%rdx
add    %rdx,%rcx
cmp    %rcx,%rax
je     123c <ranged_hist_float+0x13c>
mov    -0x38(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x4c(%rbp),%xmm0
ucomiss -0x24(%rbp),%xmm0
jb     11fc <ranged_hist_float+0xfc>
movss  -0x4c(%rbp),%xmm1
movss  -0x28(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    11fc <ranged_hist_float+0xfc>
mov    -0x20(%rbp),%rax
movss  -0x3c(%rbp),%xmm0
movss  -0x4c(%rbp),%xmm1
subss  -0x24(%rbp),%xmm1
mulss  %xmm1,%xmm0
cvttss2si %xmm0,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%edx
add    $0x1,%edx
mov    %edx,(%rax,%rcx,4)
jmp    1221 <ranged_hist_float+0x121>
movss  -0x4c(%rbp),%xmm0
ucomiss -0x28(%rbp),%xmm0
jne    121c <ranged_hist_float+0x11c>
jp     121c <ranged_hist_float+0x11c>
mov    -0x48(%rbp),%rax
mov    (%rax),%ecx
add    $0x1,%ecx
mov    %ecx,(%rax)
jmp    1221 <ranged_hist_float+0x121>
jmp    1226 <ranged_hist_float+0x126>
mov    -0x18(%rbp),%ecx
mov    -0x38(%rbp),%rax
movslq %ecx,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
jmp    1183 <ranged_hist_float+0x83>
jmp    1241 <ranged_hist_float+0x141>
mov    -0x14(%rbp),%ecx
mov    -0x30(%rbp),%rax
movslq %ecx,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
jmp    115d <ranged_hist_float+0x5d>
pop    %rbp
ret

x86 assembly (clang, -O2):
<ranged_hist_float>:
push   %rbx
imul   %ecx,%esi
test   %esi,%esi
je     1189 <ranged_hist_float+0x89>
movslq 0x10(%rsp),%rax
cvtsi2ss %eax,%xmm2
movaps %xmm1,%xmm3
subss  %xmm0,%xmm3
divss  %xmm3,%xmm2
lea    (%r9,%rax,4),%r11
add    $0xfffffffffffffffc,%r11
movslq %esi,%r10
add    %rdi,%r10
imul   %r8d,%edx
movslq %edx,%rsi
movslq %r8d,%rax
movslq %ecx,%r8
jmp    1148 <ranged_hist_float+0x48>
nopl   (%rax)
add    %r8,%rdi
cmp    %r10,%rdi
je     1189 <ranged_hist_float+0x89>
test   %edx,%edx
je     1140 <ranged_hist_float+0x40>
xor    %ecx,%ecx
jmp    116e <ranged_hist_float+0x6e>
subss  %xmm0,%xmm3
mulss  %xmm2,%xmm3
cvttss2si %xmm3,%ebx
movslq %ebx,%rbx
lea    (%r9,%rbx,4),%rbx
addl   $0x1,(%rbx)
add    %rax,%rcx
cmp    %rcx,%rsi
je     1140 <ranged_hist_float+0x40>
movss  (%rdi,%rcx,1),%xmm3
ucomiss %xmm0,%xmm3
jb     117d <ranged_hist_float+0x7d>
ucomiss %xmm3,%xmm1
ja     1150 <ranged_hist_float+0x50>
ucomiss %xmm1,%xmm3
mov    %r11,%rbx
jne    1166 <ranged_hist_float+0x66>
jnp    1163 <ranged_hist_float+0x63>
jmp    1166 <ranged_hist_float+0x66>
pop    %rbx
ret

LLVM IR (-O0):
define void @ranged_hist_float(i8* %0, i32 %1, i32 %2, i32 %3, i32 %4, i32* %5, i32 %6, float %7, float %8) {
%10 = alloca i8*
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32*
%16 = alloca i32
%17 = alloca float
%18 = alloca float
%19 = alloca i8*
%20 = alloca i8*
%21 = alloca float
%22 = alloca i32*
%23 = alloca float
store i8* %0, i8** %10
store i32 %1, i32* %11
store i32 %2, i32* %12
store i32 %3, i32* %13
store i32 %4, i32* %14
store i32* %5, i32** %15
store i32 %6, i32* %16
store float %7, float* %17
store float %8, float* %18
%24 = load i32, i32* %16
%25 = sitofp i32 %24 to float
%26 = load float, float* %18
%27 = load float, float* %17
%28 = fsub float %26, %27
%29 = fdiv float %25, %28
store float %29, float* %21
%30 = load i32*, i32** %15
%31 = load i32, i32* %16
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds i32, i32* %30, i64 %32
%34 = getelementptr inbounds i32, i32* %33, i64 -1
store i32* %34, i32** %22
%35 = load i8*, i8** %10
store i8* %35, i8** %19
br label %36
36: ; preds = %95, %9
%37 = load i8*, i8** %19
%38 = load i8*, i8** %10
%39 = load i32, i32* %11
%40 = load i32, i32* %13
%41 = mul nsw i32 %39, %40
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds i8, i8* %38, i64 %42
%44 = icmp ne i8* %37, %43
br i1 %44, label %45, label %100
45: ; preds = %36
%46 = load i8*, i8** %19
store i8* %46, i8** %20
br label %47
47: ; preds = %89, %45
%48 = load i8*, i8** %20
%49 = load i8*, i8** %19
%50 = load i32, i32* %12
%51 = load i32, i32* %14
%52 = mul nsw i32 %50, %51
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds i8, i8* %49, i64 %53
%55 = icmp ne i8* %48, %54
br i1 %55, label %56, label %94
56: ; preds = %47
%57 = load i8*, i8** %20
%58 = bitcast i8* %57 to float*
%59 = load float, float* %58
store float %59, float* %23
%60 = load float, float* %23
%61 = load float, float* %17
%62 = fcmp oge float %60, %61
br i1 %62, label %63, label %79
63: ; preds = %56
%64 = load float, float* %23
%65 = load float, float* %18
%66 = fcmp olt float %64, %65
br i1 %66, label %67, label %79
67: ; preds = %63
%68 = load i32*, i32** %15
%69 = load float, float* %21
%70 = load float, float* %23
%71 = load float, float* %17
%72 = fsub float %70, %71
%73 = fmul float %69, %72
%74 = fptosi float %73 to i32
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds i32, i32* %68, i64 %75
%77 = load i32, i32* %76
%78 = add nsw i32 %77, 1
store i32 %78, i32* %76
br label %88
79: ; preds = %63, %56
%80 = load float, float* %23
%81 = load float, float* %18
%82 = fcmp oeq float %80, %81
br i1 %82, label %83, label %87
83: ; preds = %79
%84 = load i32*, i32** %22
%85 = load i32, i32* %84
%86 = add nsw i32 %85, 1
store i32 %86, i32* %84
br label %87
87: ; preds = %83, %79
br label %88
88: ; preds = %87, %67
br label %89
89: ; preds = %88
%90 = load i32, i32* %14
%91 = load i8*, i8** %20
%92 = sext i32 %90 to i64
%93 = getelementptr inbounds i8, i8* %91, i64 %92
store i8* %93, i8** %20
br label %47
94: ; preds = %47
br label %95
95: ; preds = %94
%96 = load i32, i32* %13
%97 = load i8*, i8** %19
%98 = sext i32 %96 to i64
%99 = getelementptr inbounds i8, i8* %97, i64 %98
store i8* %99, i8** %19
br label %36
100: ; preds = %36
ret void
}

LLVM IR (-O2):
define void @ranged_hist_float(i8* %0, i32 %1, i32 %2, i32 %3, i32 %4, i32* nocapture %5, i32 %6, float %7, float %8) {
%10 = sitofp i32 %6 to float
%11 = fsub float %8, %7
%12 = fdiv float %10, %11
%13 = sext i32 %6 to i64
%14 = getelementptr inbounds i32, i32* %5, i64 %13
%15 = getelementptr inbounds i32, i32* %14, i64 -1
%16 = mul nsw i32 %3, %1
%17 = sext i32 %16 to i64
%18 = getelementptr inbounds i8, i8* %0, i64 %17
%19 = icmp eq i32 %16, 0
br i1 %19, label %54, label %20
20: ; preds = %9
%21 = mul nsw i32 %4, %2
%22 = sext i32 %21 to i64
%23 = sext i32 %4 to i64
%24 = sext i32 %3 to i64
%25 = icmp eq i32 %21, 0
br label %26
26: ; preds = %20, %51
%27 = phi i8* [ %0, %20 ], [ %52, %51 ]
%28 = getelementptr inbounds i8, i8* %27, i64 %22
br i1 %25, label %51, label %29
29: ; preds = %26, %48
%30 = phi i8* [ %49, %48 ], [ %27, %26 ]
%31 = bitcast i8* %30 to float*
%32 = load float, float* %31
%33 = fcmp oge float %32, %7
%34 = fcmp olt float %32, %8
%35 = and i1 %33, %34
br i1 %35, label %36, label %42
36: ; preds = %29
%37 = fsub float %32, %7
%38 = fmul float %12, %37
%39 = fptosi float %38 to i32
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds i32, i32* %5, i64 %40
br label %44
42: ; preds = %29
%43 = fcmp oeq float %32, %8
br i1 %43, label %44, label %48
44: ; preds = %42, %36
%45 = phi i32* [ %41, %36 ], [ %15, %42 ]
%46 = load i32, i32* %45
%47 = add nsw i32 %46, 1
store i32 %47, i32* %45
br label %48
48: ; preds = %44, %42
%49 = getelementptr inbounds i8, i8* %30, i64 %23
%50 = icmp eq i8* %49, %28
br i1 %50, label %51, label %29
51: ; preds = %48, %26
%52 = getelementptr inbounds i8, i8* %27, i64 %24
%53 = icmp eq i8* %52, %18
br i1 %53, label %54, label %26
54: ; preds = %51, %9
ret void
}