void locate(float xx[], unsigned long n, float x, unsigned long *j) {
  unsigned long ju, jm, jl;
  int ascnd;
  jl = 0;
  ju = n + 1;
  ascnd = (xx[n] >= xx[1]);
  while (ju - jl > 1) {
    jm = (ju + jl) >> 1;
    if (x >= xx[jm] == ascnd)
      jl = jm;
    else
      ju = jm;
  }
  if (x == xx[1])
    *j = 1;
  else if (x == xx[n])
    *j = n - 1;
  else
    *j = jl;
}


x86 assembly (clang, -O0):
<locate>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movss  %xmm0,-0x14(%rbp)
mov    %rdx,-0x20(%rbp)
movq   $0x0,-0x38(%rbp)
mov    -0x10(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
ucomiss 0x4(%rax),%xmm0
setae  %al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x28(%rbp),%rax
sub    -0x38(%rbp),%rax
cmp    $0x1,%rax
jbe    11a7 <locate+0xa7>
mov    -0x28(%rbp),%rax
add    -0x38(%rbp),%rax
shr    $0x1,%rax
mov    %rax,-0x30(%rbp)
movss  -0x14(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
ucomiss (%rax,%rcx,4),%xmm0
setae  %al
and    $0x1,%al
movzbl %al,%eax
cmp    -0x3c(%rbp),%eax
jne    119a <locate+0x9a>
mov    -0x30(%rbp),%rax
mov    %rax,-0x38(%rbp)
jmp    11a2 <locate+0xa2>
mov    -0x30(%rbp),%rax
mov    %rax,-0x28(%rbp)
jmp    1149 <locate+0x49>
movss  -0x14(%rbp),%xmm0
mov    -0x8(%rbp),%rax
ucomiss 0x4(%rax),%xmm0
jne    11d0 <locate+0xd0>
jp     11d0 <locate+0xd0>
mov    -0x20(%rbp),%rax
movq   $0x1,(%rax)
jmp    1211 <locate+0x111>
movss  -0x14(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%rcx
ucomiss (%rax,%rcx,4),%xmm0
jne    1201 <locate+0x101>
jp     1201 <locate+0x101>
mov    -0x10(%rbp),%rcx
sub    $0x1,%rcx
mov    -0x20(%rbp),%rax
mov    %rcx,(%rax)
jmp    120c <locate+0x10c>
mov    -0x38(%rbp),%rcx
mov    -0x20(%rbp),%rax
mov    %rcx,(%rax)
jmp    1211 <locate+0x111>
pop    %rbp
ret

x86 assembly (clang, -O2):
<locate>:
lea    0x1(%rsi),%r8
movss  (%rdi,%rsi,4),%xmm1
movss  0x4(%rdi),%xmm2
xor    %r9d,%r9d
cmp    $0x2,%r8
jae    113c <locate+0x3c>
mov    $0x1,%ecx
ucomiss %xmm0,%xmm2
jne    115e <locate+0x5e>
jp     115e <locate+0x5e>
mov    %rcx,(%rdx)
ret
nopw   0x0(%rax,%rax,1)
mov    %rcx,%r8
sub    %r9,%rcx
cmp    $0x1,%rcx
jbe    1117 <locate+0x17>
ucomiss %xmm2,%xmm1
setae  %r10b
lea    (%r8,%r9,1),%rcx
shr    %rcx
ucomiss (%rdi,%rcx,4),%xmm0
setb   %al
xor    %r10b,%al
je     1130 <locate+0x30>
mov    %rcx,%r9
mov    %r8,%rcx
jmp    1133 <locate+0x33>
ucomiss %xmm0,%xmm1
mov    %r9,%rcx
jne    1123 <locate+0x23>
jp     1123 <locate+0x23>
add    $0xffffffffffffffff,%rsi
mov    %rsi,%rcx
mov    %rcx,(%rdx)
ret

LLVM IR (-O0):
define void @locate(float* %0, i64 %1, float %2, i64* %3) {
%5 = alloca float*
%6 = alloca i64
%7 = alloca float
%8 = alloca i64*
%9 = alloca i64
%10 = alloca i64
%11 = alloca i64
%12 = alloca i32
store float* %0, float** %5
store i64 %1, i64* %6
store float %2, float* %7
store i64* %3, i64** %8
store i64 0, i64* %11
%13 = load i64, i64* %6
%14 = add i64 %13, 1
store i64 %14, i64* %9
%15 = load float*, float** %5
%16 = load i64, i64* %6
%17 = getelementptr inbounds float, float* %15, i64 %16
%18 = load float, float* %17
%19 = load float*, float** %5
%20 = getelementptr inbounds float, float* %19, i64 1
%21 = load float, float* %20
%22 = fcmp oge float %18, %21
%23 = zext i1 %22 to i32
store i32 %23, i32* %12
br label %24
24: ; preds = %47, %4
%25 = load i64, i64* %9
%26 = load i64, i64* %11
%27 = sub i64 %25, %26
%28 = icmp ugt i64 %27, 1
br i1 %28, label %29, label %48
29: ; preds = %24
%30 = load i64, i64* %9
%31 = load i64, i64* %11
%32 = add i64 %30, %31
%33 = lshr i64 %32, 1
store i64 %33, i64* %10
%34 = load float, float* %7
%35 = load float*, float** %5
%36 = load i64, i64* %10
%37 = getelementptr inbounds float, float* %35, i64 %36
%38 = load float, float* %37
%39 = fcmp oge float %34, %38
%40 = zext i1 %39 to i32
%41 = load i32, i32* %12
%42 = icmp eq i32 %40, %41
br i1 %42, label %43, label %45
43: ; preds = %29
%44 = load i64, i64* %10
store i64 %44, i64* %11
br label %47
45: ; preds = %29
%46 = load i64, i64* %10
store i64 %46, i64* %9
br label %47
47: ; preds = %45, %43
br label %24
48: ; preds = %24
%49 = load float, float* %7
%50 = load float*, float** %5
%51 = getelementptr inbounds float, float* %50, i64 1
%52 = load float, float* %51
%53 = fcmp oeq float %49, %52
br i1 %53, label %54, label %56
54: ; preds = %48
%55 = load i64*, i64** %8
store i64 1, i64* %55
br label %71
56: ; preds = %48
%57 = load float, float* %7
%58 = load float*, float** %5
%59 = load i64, i64* %6
%60 = getelementptr inbounds float, float* %58, i64 %59
%61 = load float, float* %60
%62 = fcmp oeq float %57, %61
br i1 %62, label %63, label %67
63: ; preds = %56
%64 = load i64, i64* %6
%65 = sub i64 %64, 1
%66 = load i64*, i64** %8
store i64 %65, i64* %66
br label %70
67: ; preds = %56
%68 = load i64, i64* %11
%69 = load i64*, i64** %8
store i64 %68, i64* %69
br label %70
70: ; preds = %67, %63
br label %71
71: ; preds = %70, %54
ret void
}

LLVM IR (-O2):
define void @locate(float* nocapture %0, i64 %1, float %2, i64* nocapture writeonly %3) {
%5 = add i64 %1, 1
%6 = getelementptr inbounds float, float* %0, i64 %1
%7 = load float, float* %6
%8 = getelementptr inbounds float, float* %0, i64 1
%9 = load float, float* %8
%10 = fcmp oge float %7, %9
%11 = icmp ugt i64 %5, 1
br i1 %11, label %12, label %25
12: ; preds = %4, %12
%13 = phi i64 [ %22, %12 ], [ %5, %4 ]
%14 = phi i64 [ %21, %12 ], [ 0, %4 ]
%15 = add i64 %13, %14
%16 = lshr i64 %15, 1
%17 = getelementptr inbounds float, float* %0, i64 %16
%18 = load float, float* %17
%19 = fcmp ugt float %18, %2
%20 = xor i1 %10, %19
%21 = select i1 %20, i64 %16, i64 %14
%22 = select i1 %20, i64 %13, i64 %16
%23 = sub i64 %22, %21
%24 = icmp ugt i64 %23, 1
br i1 %24, label %12, label %25
25: ; preds = %12, %4
%26 = phi i64 [ 0, %4 ], [ %21, %12 ]
%27 = fcmp oeq float %9, %2
br i1 %27, label %32, label %28
28: ; preds = %25
%29 = fcmp oeq float %7, %2
br i1 %29, label %30, label %32
30: ; preds = %28
%31 = add i64 %1, -1
br label %32
32: ; preds = %28, %25, %30
%33 = phi i64 [ %31, %30 ], [ 1, %25 ], [ %26, %28 ]
store i64 %33, i64* %3
ret void
}