int locate(double *xx, int n, double x) {
  int ju = 0, jm = 0, jl = 0, j = 0;
  int monotonicity = 0;
  jl = 0;
  ju = n + 1;
  if (n == 1)
    monotonicity = 1;
  else if (xx[1] > xx[0])
    monotonicity = 1;
  xx -= 1;
  while (ju - jl > 1) {
    jm = (ju + jl) >> 1;
    if ((x >= xx[jm] && monotonicity == 1) ||
        (x <= xx[jm] && monotonicity == 0))
      jl = jm;
    else
      ju = jm;
  }
  if (x == xx[1]) {
    j = 1;
  } else {
    if (x == xx[n])
      j = n - 1;
    else
      j = jl;
  }
  j -= 1;
  return j;
}


x86 assembly (clang, -O0):
<locate>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movsd  %xmm0,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
movl   $0x0,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
cmpl   $0x1,-0xc(%rbp)
jne    1159 <locate+0x59>
movl   $0x1,-0x2c(%rbp)
jmp    117c <locate+0x7c>
mov    -0x8(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x8(%rbp),%rax
ucomisd (%rax),%xmm0
jbe    1177 <locate+0x77>
movl   $0x1,-0x2c(%rbp)
jmp    117c <locate+0x7c>
mov    -0x8(%rbp),%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,-0x8(%rbp)
mov    -0x1c(%rbp),%eax
sub    -0x24(%rbp),%eax
cmp    $0x1,%eax
jle    1201 <locate+0x101>
mov    -0x1c(%rbp),%eax
add    -0x24(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x20(%rbp)
movsd  -0x18(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
jb     11c5 <locate+0xc5>
cmpl   $0x1,-0x2c(%rbp)
je     11eb <locate+0xeb>
movsd  -0x18(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jb     11f6 <locate+0xf6>
cmpl   $0x0,-0x2c(%rbp)
jne    11f6 <locate+0xf6>
mov    -0x20(%rbp),%eax
mov    %eax,-0x24(%rbp)
jmp    11fc <locate+0xfc>
mov    -0x20(%rbp),%eax
mov    %eax,-0x1c(%rbp)
jmp    1188 <locate+0x88>
movsd  -0x18(%rbp),%xmm0
mov    -0x8(%rbp),%rax
ucomisd 0x8(%rax),%xmm0
jne    1227 <locate+0x127>
jp     1227 <locate+0x127>
movl   $0x1,-0x28(%rbp)
jmp    125e <locate+0x15e>
movsd  -0x18(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
jne    1253 <locate+0x153>
jp     1253 <locate+0x153>
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1259 <locate+0x159>
mov    -0x24(%rbp),%eax
mov    %eax,-0x28(%rbp)
jmp    125e <locate+0x15e>
mov    -0x28(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<locate>:
push   %rbx
cmp    $0x1,%esi
jne    1112 <locate+0x12>
lea    -0x8(%rdi),%r8
mov    $0x1,%r9b
xor    %r10d,%r10d
jmp    112b <locate+0x2b>
movsd  0x8(%rdi),%xmm1
ucomisd (%rdi),%xmm1
seta   %r9b
setbe  %r10b
lea    -0x8(%rdi),%r8
test   %esi,%esi
jle    1178 <locate+0x78>
lea    0x1(%rsi),%eax
xor    %r11d,%r11d
jmp    114c <locate+0x4c>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %edx,%eax
mov    %eax,%ecx
sub    %r11d,%ecx
cmp    $0x1,%ecx
jle    117b <locate+0x7b>
lea    (%r11,%rax,1),%edx
sar    %edx
movslq %edx,%rcx
movsd  -0x8(%rdi,%rcx,8),%xmm1
ucomisd %xmm1,%xmm0
setae  %cl
and    %r9b,%cl
ucomisd %xmm0,%xmm1
setae  %bl
and    %r10b,%bl
or     %cl,%bl
je     1140 <locate+0x40>
mov    %edx,%r11d
jmp    1142 <locate+0x42>
xor    %r11d,%r11d
movsd  (%rdi),%xmm1
mov    $0x1,%eax
ucomisd %xmm0,%xmm1
jne    118c <locate+0x8c>
jnp    11a5 <locate+0xa5>
movslq %esi,%rax
movsd  (%r8,%rax,8),%xmm1
ucomisd %xmm0,%xmm1
mov    %r11d,%eax
jne    11a5 <locate+0xa5>
jp     11a5 <locate+0xa5>
add    $0xffffffff,%esi
mov    %esi,%eax
add    $0xffffffff,%eax
pop    %rbx
ret

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

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