#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 masked_ranged_hist_float(const char *image, uint16_t rows, uint16_t cols,
                              uint32_t r_stride, uint32_t c_stride,
                              const uint16_t *starts, const uint16_t *ends,
                              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, starts++, ends++) {
    for (pixel = row_start + (*starts) * c_stride;
         pixel != row_start + (*ends) * 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):
<masked_ranged_hist_float>:
push   %rbp
mov    %rsp,%rbp
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
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 0x20(%rbp),%xmm0
movss  -0x28(%rbp),%xmm1
subss  -0x24(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x3c(%rbp)
mov    0x18(%rbp),%rax
movslq 0x20(%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     128a <masked_ranged_hist_float+0x18a>
mov    -0x30(%rbp),%rax
mov    -0x20(%rbp),%rcx
mov    (%rcx),%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
mov    -0x38(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    0x10(%rbp),%rdx
mov    (%rdx),%edx
imul   -0x18(%rbp),%edx
movslq %edx,%rdx
add    %rdx,%rcx
cmp    %rcx,%rax
je     1257 <masked_ranged_hist_float+0x157>
mov    -0x38(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x4c(%rbp),%xmm0
ucomiss -0x24(%rbp),%xmm0
jb     1217 <masked_ranged_hist_float+0x117>
movss  -0x4c(%rbp),%xmm1
movss  -0x28(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    1217 <masked_ranged_hist_float+0x117>
mov    0x18(%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    123c <masked_ranged_hist_float+0x13c>
movss  -0x4c(%rbp),%xmm0
ucomiss -0x28(%rbp),%xmm0
jne    1237 <masked_ranged_hist_float+0x137>
jp     1237 <masked_ranged_hist_float+0x137>
mov    -0x48(%rbp),%rax
mov    (%rax),%ecx
add    $0x1,%ecx
mov    %ecx,(%rax)
jmp    123c <masked_ranged_hist_float+0x13c>
jmp    1241 <masked_ranged_hist_float+0x141>
mov    -0x18(%rbp),%ecx
mov    -0x38(%rbp),%rax
movslq %ecx,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
jmp    119b <masked_ranged_hist_float+0x9b>
jmp    125c <masked_ranged_hist_float+0x15c>
mov    -0x14(%rbp),%ecx
mov    -0x30(%rbp),%rax
movslq %ecx,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
mov    -0x20(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x20(%rbp)
mov    0x10(%rbp),%rax
add    $0x4,%rax
mov    %rax,0x10(%rbp)
jmp    1165 <masked_ranged_hist_float+0x65>
pop    %rbp
ret

x86 assembly (clang, -O2):
<masked_ranged_hist_float>:
push   %rbx
imul   %ecx,%esi
test   %esi,%esi
je     11b2 <masked_ranged_hist_float+0xb2>
movslq 0x20(%rsp),%rax
mov    0x18(%rsp),%r10
cvtsi2ss %eax,%xmm2
mov    0x10(%rsp),%rdx
movaps %xmm1,%xmm3
subss  %xmm0,%xmm3
divss  %xmm3,%xmm2
lea    (%r10,%rax,4),%rax
add    $0xfffffffffffffffc,%rax
movslq %esi,%r11
add    %rdi,%r11
movslq %r8d,%rsi
movslq %ecx,%r8
jmp    1150 <masked_ranged_hist_float+0x50>
add    %r8,%rdi
add    $0x4,%r9
add    $0x4,%rdx
cmp    %r11,%rdi
je     11b2 <masked_ranged_hist_float+0xb2>
movslq (%r9),%rcx
imul   %rsi,%rcx
mov    (%rdx),%ebx
imul   %esi,%ebx
cmp    %ebx,%ecx
je     1140 <masked_ranged_hist_float+0x40>
add    %rdi,%rcx
jmp    1198 <masked_ranged_hist_float+0x98>
cs nopw 0x0(%rax,%rax,1)
nop
subss  %xmm0,%xmm3
mulss  %xmm2,%xmm3
cvttss2si %xmm3,%ebx
movslq %ebx,%rbx
lea    (%r10,%rbx,4),%rbx
addl   $0x1,(%rbx)
add    %rsi,%rcx
movslq (%rdx),%rbx
imul   %rsi,%rbx
add    %rdi,%rbx
cmp    %rbx,%rcx
je     1140 <masked_ranged_hist_float+0x40>
movss  (%rcx),%xmm3
ucomiss %xmm0,%xmm3
jb     11a6 <masked_ranged_hist_float+0xa6>
ucomiss %xmm3,%xmm1
ja     1170 <masked_ranged_hist_float+0x70>
ucomiss %xmm1,%xmm3
mov    %rax,%rbx
jne    1186 <masked_ranged_hist_float+0x86>
jnp    1183 <masked_ranged_hist_float+0x83>
jmp    1186 <masked_ranged_hist_float+0x86>
pop    %rbx
ret

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

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