int findIndexBin(double *CDF, int beginIndex, int endIndex, double value) {
  if (endIndex < beginIndex)
    return -1;
  int middleIndex = beginIndex + ((endIndex - beginIndex) / 2);
  if (CDF[middleIndex] >= value) {
    if (middleIndex == 0)
      return middleIndex;
    else if (CDF[middleIndex - 1] < value)
      return middleIndex;
    else if (CDF[middleIndex - 1] == value) {
      while (middleIndex > 0 && CDF[middleIndex - 1] == value)
        middleIndex--;
      return middleIndex;
    }
  }
  if (CDF[middleIndex] > value)
    return findIndexBin(CDF, beginIndex, middleIndex + 1, value);
  return findIndexBin(CDF, middleIndex - 1, endIndex, value);
}


x86 assembly (clang, -O0):
<findIndexBin>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %edx,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    113f <findIndexBin+0x2f>
movl   $0xffffffff,-0x4(%rbp)
jmp    1298 <findIndexBin+0x188>
mov    -0x14(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x18(%rbp),%eax
sub    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x28(%rbp),%eax
add    %ecx,%eax
mov    %eax,-0x24(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x20(%rbp),%xmm0
jb     1247 <findIndexBin+0x137>
cmpl   $0x0,-0x24(%rbp)
jne    118a <findIndexBin+0x7a>
mov    -0x24(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    1298 <findIndexBin+0x188>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x20(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    11b6 <findIndexBin+0xa6>
mov    -0x24(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    1298 <findIndexBin+0x188>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x20(%rbp),%xmm0
jne    1238 <findIndexBin+0x128>
jp     1238 <findIndexBin+0x128>
jmp    11de <findIndexBin+0xce>
xor    %eax,%eax
cmpl   $0x0,-0x24(%rbp)
mov    %al,-0x29(%rbp)
jle    120f <findIndexBin+0xff>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x20(%rbp),%xmm0
sete   %al
setnp  %cl
and    %cl,%al
mov    %al,-0x29(%rbp)
mov    -0x29(%rbp),%al
test   $0x1,%al
jne    121f <findIndexBin+0x10f>
jmp    122d <findIndexBin+0x11d>
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    11de <findIndexBin+0xce>
mov    -0x24(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    1298 <findIndexBin+0x188>
jmp    123d <findIndexBin+0x12d>
jmp    1242 <findIndexBin+0x132>
jmp    1247 <findIndexBin+0x137>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x20(%rbp),%xmm0
jbe    127e <findIndexBin+0x16e>
mov    -0x10(%rbp),%rdi
mov    -0x14(%rbp),%esi
mov    -0x24(%rbp),%edx
add    $0x1,%edx
movsd  -0x20(%rbp),%xmm0
call   1030 <findIndexBin@plt>
mov    %eax,-0x4(%rbp)
jmp    1298 <findIndexBin+0x188>
mov    -0x10(%rbp),%rdi
mov    -0x24(%rbp),%esi
sub    $0x1,%esi
mov    -0x18(%rbp),%edx
movsd  -0x20(%rbp),%xmm0
call   1030 <findIndexBin@plt>
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<findIndexBin>:
mov    $0xffffffff,%eax
cmp    %esi,%edx
jge    1110 <findIndexBin+0x10>
ret
nopw   0x0(%rax,%rax,1)
mov    %edx,%r8d
sub    %esi,%r8d
mov    %r8d,%ecx
shr    $0x1f,%ecx
add    %r8d,%ecx
sar    %ecx
add    %esi,%ecx
movslq %ecx,%r8
movsd  (%rdi,%r8,8),%xmm1
ucomisd %xmm0,%xmm1
jb     114b <findIndexBin+0x4b>
test   %ecx,%ecx
je     1167 <findIndexBin+0x67>
movsd  -0x8(%rdi,%r8,8),%xmm2
ucomisd %xmm2,%xmm0
ja     116a <findIndexBin+0x6a>
ucomisd %xmm0,%xmm2
jne    114b <findIndexBin+0x4b>
jnp    116d <findIndexBin+0x6d>
ucomisd %xmm0,%xmm1
jbe    115c <findIndexBin+0x5c>
add    $0x1,%ecx
mov    %ecx,%edx
cmp    %esi,%ecx
jge    1110 <findIndexBin+0x10>
jmp    1109 <findIndexBin+0x9>
add    $0xffffffff,%ecx
mov    %ecx,%esi
cmp    %ecx,%edx
jge    1110 <findIndexBin+0x10>
jmp    1109 <findIndexBin+0x9>
xor    %eax,%eax
ret
mov    %ecx,%eax
ret
mov    %ecx,%edx
xor    %eax,%eax
test   %ecx,%ecx
cmovs  %ecx,%eax
cs nopw 0x0(%rax,%rax,1)
test   %edx,%edx
jle    1109 <findIndexBin+0x9>
movsd  -0x8(%rdi,%rdx,8),%xmm1
add    $0xffffffffffffffff,%rdx
ucomisd %xmm0,%xmm1
jne    1196 <findIndexBin+0x96>
jnp    1180 <findIndexBin+0x80>
add    $0x1,%edx
mov    %edx,%eax
ret

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

LLVM IR (-O2):
define i32 @findIndexBin(double* %0, i32 %1, i32 %2, double %3) {
%5 = icmp slt i32 %2, %1
br i1 %5, label %50, label %6
6: ; preds = %4, %45
%7 = phi i32 [ %10, %45 ], [ %2, %4 ]
%8 = phi i32 [ %46, %45 ], [ %1, %4 ]
br label %9
9: ; preds = %6, %42
%10 = phi i32 [ %7, %6 ], [ %43, %42 ]
%11 = sub nsw i32 %10, %8
%12 = sdiv i32 %11, 2
%13 = add nsw i32 %12, %8
%14 = sext i32 %13 to i64
%15 = getelementptr inbounds double, double* %0, i64 %14
%16 = load double, double* %15
%17 = fcmp ult double %16, %3
br i1 %17, label %40, label %18
18: ; preds = %9
%19 = icmp eq i32 %13, 0
br i1 %19, label %50, label %20
20: ; preds = %18
%21 = add nsw i32 %13, -1
%22 = sext i32 %21 to i64
%23 = getelementptr inbounds double, double* %0, i64 %22
%24 = load double, double* %23
%25 = fcmp olt double %24, %3
br i1 %25, label %50, label %26
26: ; preds = %20
%27 = fcmp oeq double %24, %3
br i1 %27, label %28, label %40
28: ; preds = %26
%29 = zext i32 %13 to i64
%30 = call i32 @llvm.smin.i32(i32 %13, i32 0)
br label %31
31: ; preds = %28, %35
%32 = phi i64 [ %29, %28 ], [ %36, %35 ]
%33 = trunc i64 %32 to i32
%34 = icmp sgt i32 %33, 0
br i1 %34, label %35, label %50
35: ; preds = %31
%36 = add nsw i64 %32, -1
%37 = getelementptr inbounds double, double* %0, i64 %36
%38 = load double, double* %37
%39 = fcmp oeq double %38, %3
br i1 %39, label %31, label %48
40: ; preds = %26, %9
%41 = fcmp ogt double %16, %3
br i1 %41, label %42, label %45
42: ; preds = %40
%43 = add nsw i32 %13, 1
%44 = icmp slt i32 %43, %8
br i1 %44, label %50, label %9
45: ; preds = %40
%46 = add nsw i32 %13, -1
%47 = icmp slt i32 %10, %46
br i1 %47, label %50, label %6
48: ; preds = %35
%49 = trunc i64 %32 to i32
br label %50
50: ; preds = %45, %42, %20, %18, %31, %48, %4
%51 = phi i32 [ -1, %4 ], [ %49, %48 ], [ %30, %31 ], [ -1, %42 ], [ %13, %20 ], [ 0, %18 ], [ -1, %45 ]
ret i32 %51
}