int findspan(float xi, const float knots[], const int lenght_knots) {
  int index_up = lenght_knots;
  int index_low = 0;
  int index_middle;
  while (index_up - index_low > 1) {
    index_middle = ((index_up + index_low) / 2);
    if (xi < knots[index_middle]) {
      index_up = index_middle;
    } else {
      index_low = index_middle;
    }
  }
  return index_low;
}


x86 assembly (clang, -O0):
<findspan>:
push   %rbp
mov    %rsp,%rbp
movss  %xmm0,-0x4(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x18(%rbp),%eax
sub    -0x1c(%rbp),%eax
cmp    $0x1,%eax
jle    116e <findspan+0x6e>
mov    -0x18(%rbp),%eax
add    -0x1c(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x20(%rbp)
movss  -0x4(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss %xmm1,%xmm0
jbe    1163 <findspan+0x63>
mov    -0x20(%rbp),%eax
mov    %eax,-0x18(%rbp)
jmp    1169 <findspan+0x69>
mov    -0x20(%rbp),%eax
mov    %eax,-0x1c(%rbp)
jmp    111d <findspan+0x1d>
mov    -0x1c(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<findspan>:
xor    %eax,%eax
cmp    $0x2,%esi
jl     1138 <findspan+0x38>
xor    %eax,%eax
jmp    111b <findspan+0x1b>
nopl   0x0(%rax,%rax,1)
mov    %ecx,%eax
mov    %esi,%ecx
sub    %eax,%ecx
cmp    $0x1,%ecx
jle    1138 <findspan+0x38>
lea    (%rsi,%rax,1),%edx
mov    %edx,%ecx
shr    $0x1f,%ecx
add    %edx,%ecx
sar    %ecx
movslq %ecx,%rdx
movss  (%rdi,%rdx,4),%xmm1
ucomiss %xmm0,%xmm1
jbe    1110 <findspan+0x10>
mov    %ecx,%esi
jmp    1112 <findspan+0x12>
ret

LLVM IR (-O0):
define i32 @findspan(float %0, float* %1, i32 %2) {
%4 = alloca float
%5 = alloca float*
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
store float %0, float* %4
store float* %1, float** %5
store i32 %2, i32* %6
%10 = load i32, i32* %6
store i32 %10, i32* %7
store i32 0, i32* %8
br label %11
11: ; preds = %32, %3
%12 = load i32, i32* %7
%13 = load i32, i32* %8
%14 = sub nsw i32 %12, %13
%15 = icmp sgt i32 %14, 1
br i1 %15, label %16, label %33
16: ; preds = %11
%17 = load i32, i32* %7
%18 = load i32, i32* %8
%19 = add nsw i32 %17, %18
%20 = sdiv i32 %19, 2
store i32 %20, i32* %9
%21 = load float, float* %4
%22 = load float*, float** %5
%23 = load i32, i32* %9
%24 = sext i32 %23 to i64
%25 = getelementptr inbounds float, float* %22, i64 %24
%26 = load float, float* %25
%27 = fcmp olt float %21, %26
br i1 %27, label %28, label %30
28: ; preds = %16
%29 = load i32, i32* %9
store i32 %29, i32* %7
br label %32
30: ; preds = %16
%31 = load i32, i32* %9
store i32 %31, i32* %8
br label %32
32: ; preds = %30, %28
br label %11
33: ; preds = %11
%34 = load i32, i32* %8
ret i32 %34
}

LLVM IR (-O2):
define i32 @findspan(float %0, float* nocapture %1, i32 %2) {
%4 = icmp sgt i32 %2, 1
br i1 %4, label %5, label %18
5: ; preds = %3, %5
%6 = phi i32 [ %15, %5 ], [ %2, %3 ]
%7 = phi i32 [ %14, %5 ], [ 0, %3 ]
%8 = add nsw i32 %6, %7
%9 = sdiv i32 %8, 2
%10 = sext i32 %9 to i64
%11 = getelementptr inbounds float, float* %1, i64 %10
%12 = load float, float* %11
%13 = fcmp ogt float %12, %0
%14 = select i1 %13, i32 %7, i32 %9
%15 = select i1 %13, i32 %9, i32 %6
%16 = sub nsw i32 %15, %14
%17 = icmp sgt i32 %16, 1
br i1 %17, label %5, label %18
18: ; preds = %5, %3
%19 = phi i32 [ 0, %3 ], [ %14, %5 ]
ret i32 %19
}