int findX(double *x, double xval, int i1, int N) {
  int u;
  if (x[i1] >= xval)
    return (i1);
  if (x[N - 1] < xval)
    return (N);
  while (1) {
    u = (N + i1) / 2;
    if (u == N || u == i1)
      return (u);
    if (x[u] >= xval) {
      if (x[u - 1] < xval)
        return (u);
      N = u;
    }
    if (x[u] <= xval) {
      if (u == N - 1)
        return (N);
      if (x[u + 1] >= xval)
        return (u + 1);
      i1 = u;
    }
  }
}


x86 assembly (clang, -O0):
<findX>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
mov    %esi,-0x1c(%rbp)
mov    %edx,-0x20(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x18(%rbp),%xmm0
jb     1136 <findX+0x36>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    1253 <findX+0x153>
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x18(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1162 <findX+0x62>
mov    -0x20(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    1253 <findX+0x153>
jmp    1167 <findX+0x67>
mov    -0x20(%rbp),%eax
add    -0x1c(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x20(%rbp),%eax
je     1190 <findX+0x90>
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jne    119b <findX+0x9b>
mov    -0x24(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    1253 <findX+0x153>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x18(%rbp),%xmm0
jb     11e5 <findX+0xe5>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x18(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    11df <findX+0xdf>
mov    -0x24(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    1253 <findX+0x153>
mov    -0x24(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x18(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jb     124e <findX+0x14e>
mov    -0x24(%rbp),%eax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jne    121d <findX+0x11d>
mov    -0x20(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    1253 <findX+0x153>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x18(%rbp),%xmm0
jb     1248 <findX+0x148>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4(%rbp)
jmp    1253 <findX+0x153>
mov    -0x24(%rbp),%eax
mov    %eax,-0x1c(%rbp)
jmp    1167 <findX+0x67>
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<findX>:
movslq %esi,%rax
movsd  (%rdi,%rax,8),%xmm1
ucomisd %xmm0,%xmm1
jae    1169 <findX+0x69>
movslq %edx,%rax
ucomisd -0x8(%rdi,%rax,8),%xmm0
jbe    1120 <findX+0x20>
mov    %edx,%eax
ret
nopl   0x0(%rax)
lea    (%rdx,%rsi,1),%ecx
mov    %ecx,%eax
shr    $0x1f,%eax
add    %ecx,%eax
sar    %eax
cmp    %edx,%eax
je     116f <findX+0x6f>
cmp    %esi,%eax
je     116f <findX+0x6f>
movslq %eax,%rcx
movsd  (%rdi,%rcx,8),%xmm1
ucomisd %xmm0,%xmm1
jb     114c <findX+0x4c>
ucomisd -0x8(%rdi,%rcx,8),%xmm0
mov    %eax,%edx
ja     116f <findX+0x6f>
ucomisd %xmm1,%xmm0
jb     1120 <findX+0x20>
lea    -0x1(%rdx),%esi
cmp    %esi,%eax
je     1119 <findX+0x19>
movsd  0x8(%rdi,%rcx,8),%xmm1
ucomisd %xmm0,%xmm1
jae    116c <findX+0x6c>
mov    %eax,%esi
jmp    1120 <findX+0x20>
mov    %esi,%eax
ret
add    $0x1,%eax
ret

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

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