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


x86 assembly (clang, -O0):
<k_locate>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movss  %xmm0,-0x14(%rbp)
mov    %rdx,-0x20(%rbp)
movq   $0xffffffffffffffff,-0x38(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%rcx
sub    $0x1,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
ucomiss (%rax),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x28(%rbp),%rax
sub    -0x38(%rbp),%rax
cmp    $0x1,%rax
jbe    11a6 <k_locate+0xa6>
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
seta   %al
and    $0x1,%al
movzbl %al,%eax
cmp    -0x3c(%rbp),%eax
jne    1199 <k_locate+0x99>
mov    -0x30(%rbp),%rax
mov    %rax,-0x38(%rbp)
jmp    11a1 <k_locate+0xa1>
mov    -0x30(%rbp),%rax
mov    %rax,-0x28(%rbp)
jmp    1148 <k_locate+0x48>
mov    -0x38(%rbp),%rcx
mov    -0x20(%rbp),%rax
mov    %rcx,(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<k_locate>:
lea    -0x1(%rsi),%rcx
mov    $0xffffffffffffffff,%r8
cmp    $0xfffffffffffffffd,%rcx
ja     115e <k_locate+0x5e>
movss  -0x4(%rdi,%rsi,4),%xmm1
movss  (%rdi),%xmm2
mov    $0xffffffffffffffff,%r8
jmp    113c <k_locate+0x3c>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %rcx,%rsi
sub    %r8,%rcx
cmp    $0x1,%rcx
jbe    115e <k_locate+0x5e>
ucomiss %xmm2,%xmm1
seta   %r9b
lea    (%rsi,%r8,1),%rcx
shr    %rcx
ucomiss (%rdi,%rcx,4),%xmm0
setbe  %al
xor    %r9b,%al
je     1130 <k_locate+0x30>
mov    %rcx,%r8
mov    %rsi,%rcx
jmp    1133 <k_locate+0x33>
mov    %r8,(%rdx)
ret

LLVM IR (-O0):
define void @k_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 -1, i64* %11
%13 = load i64, i64* %6
store i64 %13, i64* %9
%14 = load float*, float** %5
%15 = load i64, i64* %6
%16 = sub i64 %15, 1
%17 = getelementptr inbounds float, float* %14, i64 %16
%18 = load float, float* %17
%19 = load float*, float** %5
%20 = getelementptr inbounds float, float* %19, i64 0
%21 = load float, float* %20
%22 = fcmp ogt 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 ogt 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 i64, i64* %11
%50 = load i64*, i64** %8
store i64 %49, i64* %50
ret void
}

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