int find_starting_index(int window_start_coordinate, int *snp_locations,
                        int start_index, int end_index) {
  int current_index = 0;
  if (start_index == end_index || start_index + 1 == end_index) {
    return start_index;
  } else if (end_index < start_index) {
    return end_index;
  }
  if (snp_locations[start_index] >= window_start_coordinate) {
    return start_index;
  }
  current_index = (int)((end_index - start_index) / 2) + start_index;
  if (current_index >= end_index) {
    current_index = end_index - 1;
  }
  if (current_index < start_index) {
    current_index = start_index + 1;
  }
  if (snp_locations[current_index] < window_start_coordinate) {
    start_index = current_index;
  } else if (snp_locations[current_index] > window_start_coordinate) {
    end_index = current_index;
  } else {
    return current_index;
  }
  return find_starting_index(window_start_coordinate, snp_locations,
                             start_index, end_index);
}


x86 assembly (clang, -O0):
<find_starting_index>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %edi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x18(%rbp),%eax
je     1147 <find_starting_index+0x37>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
cmp    -0x18(%rbp),%eax
jne    1152 <find_starting_index+0x42>
mov    -0x14(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    122e <find_starting_index+0x11e>
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1169 <find_starting_index+0x59>
mov    -0x18(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    122e <find_starting_index+0x11e>
jmp    116e <find_starting_index+0x5e>
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
cmp    -0x8(%rbp),%eax
jl     118d <find_starting_index+0x7d>
mov    -0x14(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    122e <find_starting_index+0x11e>
mov    -0x18(%rbp),%eax
sub    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
add    -0x14(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jl     11b6 <find_starting_index+0xa6>
mov    -0x18(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11cb <find_starting_index+0xbb>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
cmp    -0x8(%rbp),%eax
jge    11ea <find_starting_index+0xda>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x14(%rbp)
jmp    1219 <find_starting_index+0x109>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
cmp    -0x8(%rbp),%eax
jle    1209 <find_starting_index+0xf9>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x18(%rbp)
jmp    1214 <find_starting_index+0x104>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    122e <find_starting_index+0x11e>
jmp    1219 <find_starting_index+0x109>
mov    -0x8(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x14(%rbp),%edx
mov    -0x18(%rbp),%ecx
call   1030 <find_starting_index@plt>
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<find_starting_index>:
cmp    %edx,%ecx
jne    111c <find_starting_index+0x1c>
mov    %ecx,%eax
ret
nopw   0x0(%rax,%rax,1)
mov    %eax,%edx
mov    %ecx,%eax
mov    %eax,%ecx
cmp    %edx,%eax
mov    %edx,%eax
je     1160 <find_starting_index+0x60>
lea    0x1(%rdx),%r8d
cmp    %ecx,%r8d
je     1161 <find_starting_index+0x61>
mov    %ecx,%r9d
sub    %edx,%r9d
jl     1104 <find_starting_index+0x4>
movslq %edx,%rax
cmp    %edi,(%rsi,%rax,4)
jge    1161 <find_starting_index+0x61>
mov    %r9d,%r10d
shr    $0x1f,%r10d
add    %r9d,%r10d
sar    %r10d
add    %edx,%r10d
lea    -0x1(%rcx),%eax
cmp    %ecx,%r10d
cmovl  %r10d,%eax
cmp    %edx,%eax
cmovl  %r8d,%eax
movslq %eax,%r8
cmp    %edi,(%rsi,%r8,4)
jl     1110 <find_starting_index+0x10>
jg     1114 <find_starting_index+0x14>
ret
mov    %edx,%eax
ret

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

LLVM IR (-O2):
define i32 @find_starting_index(i32 %0, i32* nocapture %1, i32 %2, i32 %3) {
%5 = icmp eq i32 %3, %2
br i1 %5, label %37, label %6
6: ; preds = %4, %33
%7 = phi i32 [ %35, %33 ], [ %3, %4 ]
%8 = phi i32 [ %34, %33 ], [ %2, %4 ]
%9 = add nsw i32 %8, 1
%10 = icmp eq i32 %9, %7
br i1 %10, label %37, label %11
11: ; preds = %6
%12 = icmp slt i32 %7, %8
br i1 %12, label %37, label %13
13: ; preds = %11
%14 = sext i32 %8 to i64
%15 = getelementptr inbounds i32, i32* %1, i64 %14
%16 = load i32, i32* %15
%17 = icmp slt i32 %16, %0
br i1 %17, label %18, label %37
18: ; preds = %13
%19 = sub nsw i32 %7, %8
%20 = sdiv i32 %19, 2
%21 = add nsw i32 %20, %8
%22 = icmp slt i32 %21, %7
%23 = add nsw i32 %7, -1
%24 = select i1 %22, i32 %21, i32 %23
%25 = icmp slt i32 %24, %8
%26 = select i1 %25, i32 %9, i32 %24
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds i32, i32* %1, i64 %27
%29 = load i32, i32* %28
%30 = icmp slt i32 %29, %0
br i1 %30, label %33, label %31
31: ; preds = %18
%32 = icmp sgt i32 %29, %0
br i1 %32, label %33, label %37
33: ; preds = %31, %18
%34 = phi i32 [ %26, %18 ], [ %8, %31 ]
%35 = phi i32 [ %7, %18 ], [ %26, %31 ]
%36 = icmp eq i32 %35, %34
br i1 %36, label %37, label %6
37: ; preds = %6, %33, %11, %13, %31, %4
%38 = phi i32 [ %3, %4 ], [ %26, %31 ], [ %8, %13 ], [ %7, %11 ], [ %34, %33 ], [ %8, %6 ]
ret i32 %38
}