unsigned int FindSpanIncSeq(int n, double u, double *U, int PrevSpan) {
  unsigned int low, high, mid;
  if (u < U[PrevSpan + 1]) {
    return (PrevSpan);
  } else if (u < U[PrevSpan + 2]) {
    return (PrevSpan + 1);
  } else if (u >= U[n - 1]) {
    return (n - 1);
  } else {
    low = PrevSpan + 2;
    high = n - 1;
    mid = (low + high) / 2;
    while (u < U[mid] || u >= U[mid + 1]) {
      if (u < U[mid]) {
        high = mid;
      } else {
        low = mid + 1;
      }
      mid = (low + high) / 2;
    }
    return (mid);
  }
}


x86 assembly (clang, -O0):
<FindSpanIncSeq>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
movsd  -0x10(%rbp),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    113f <FindSpanIncSeq+0x3f>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    1249 <FindSpanIncSeq+0x149>
movsd  -0x10(%rbp),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    116e <FindSpanIncSeq+0x6e>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4(%rbp)
jmp    1249 <FindSpanIncSeq+0x149>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x8(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
ucomisd (%rax,%rcx,8),%xmm0
jb     1199 <FindSpanIncSeq+0x99>
mov    -0x8(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x4(%rbp)
jmp    1249 <FindSpanIncSeq+0x149>
mov    -0x1c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x20(%rbp)
mov    -0x8(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x20(%rbp),%eax
add    -0x24(%rbp),%eax
shr    $0x1,%eax
mov    %eax,-0x28(%rbp)
movsd  -0x10(%rbp),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
mov    $0x1,%al
ucomisd %xmm1,%xmm0
mov    %al,-0x29(%rbp)
ja     11f3 <FindSpanIncSeq+0xf3>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,%ecx
ucomisd (%rax,%rcx,8),%xmm0
setae  %al
mov    %al,-0x29(%rbp)
mov    -0x29(%rbp),%al
test   $0x1,%al
jne    1203 <FindSpanIncSeq+0x103>
jmp    1243 <FindSpanIncSeq+0x143>
movsd  -0x10(%rbp),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    1229 <FindSpanIncSeq+0x129>
mov    -0x28(%rbp),%eax
mov    %eax,-0x24(%rbp)
jmp    1232 <FindSpanIncSeq+0x132>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
add    -0x24(%rbp),%eax
shr    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    11b7 <FindSpanIncSeq+0xb7>
mov    -0x28(%rbp),%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<FindSpanIncSeq>:
mov    %edx,%eax
movslq %edx,%rcx
movsd  0x8(%rsi,%rcx,8),%xmm1
ucomisd %xmm0,%xmm1
ja     1120 <FindSpanIncSeq+0x20>
movsd  0x10(%rsi,%rcx,8),%xmm1
ucomisd %xmm0,%xmm1
jbe    1121 <FindSpanIncSeq+0x21>
add    $0x1,%eax
ret
movslq %edi,%rcx
add    $0xffffffff,%edi
ucomisd -0x8(%rsi,%rcx,8),%xmm0
jae    1166 <FindSpanIncSeq+0x66>
add    $0x2,%eax
mov    %eax,%ecx
add    %edi,%eax
jmp    1145 <FindSpanIncSeq+0x45>
nopl   0x0(%rax,%rax,1)
mov    %eax,%edi
lea    (%rdi,%rcx,1),%eax
shr    %eax
movsd  (%rsi,%rax,8),%xmm1
ucomisd %xmm0,%xmm1
ja     1140 <FindSpanIncSeq+0x40>
mov    %eax,%ecx
ucomisd 0x8(%rsi,%rcx,8),%xmm0
jb     1120 <FindSpanIncSeq+0x20>
add    $0x1,%eax
mov    %eax,%ecx
lea    (%rdi,%rcx,1),%eax
jmp    1145 <FindSpanIncSeq+0x45>
mov    %edi,%eax
ret

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

LLVM IR (-O2):
define i32 @FindSpanIncSeq(i32 %0, double %1, double* nocapture %2, i32 %3) {
%5 = add nsw i32 %3, 1
%6 = sext i32 %5 to i64
%7 = getelementptr inbounds double, double* %2, i64 %6
%8 = load double, double* %7
%9 = fcmp ogt double %8, %1
br i1 %9, label %43, label %10
10: ; preds = %4
%11 = add nsw i32 %3, 2
%12 = sext i32 %11 to i64
%13 = getelementptr inbounds double, double* %2, i64 %12
%14 = load double, double* %13
%15 = fcmp ogt double %14, %1
br i1 %15, label %43, label %16
16: ; preds = %10
%17 = add nsw i32 %0, -1
%18 = sext i32 %17 to i64
%19 = getelementptr inbounds double, double* %2, i64 %18
%20 = load double, double* %19
%21 = fcmp ugt double %20, %1
br i1 %21, label %22, label %43
22: ; preds = %16
%23 = add i32 %11, %17
br label %24
24: ; preds = %39, %22
%25 = phi i32 [ %11, %22 ], [ %40, %39 ]
%26 = phi i32 [ %17, %22 ], [ %41, %39 ]
%27 = phi i32 [ %23, %22 ], [ %42, %39 ]
%28 = lshr i32 %27, 1
%29 = zext i32 %28 to i64
%30 = getelementptr inbounds double, double* %2, i64 %29
%31 = load double, double* %30
%32 = fcmp ogt double %31, %1
br i1 %32, label %39, label %33
33: ; preds = %24
%34 = add nuw i32 %28, 1
%35 = zext i32 %34 to i64
%36 = getelementptr inbounds double, double* %2, i64 %35
%37 = load double, double* %36
%38 = fcmp ugt double %37, %1
br i1 %38, label %43, label %39
39: ; preds = %24, %33
%40 = phi i32 [ %34, %33 ], [ %25, %24 ]
%41 = phi i32 [ %26, %33 ], [ %28, %24 ]
%42 = add i32 %41, %40
br label %24
43: ; preds = %33, %16, %10, %4
%44 = phi i32 [ %3, %4 ], [ %5, %10 ], [ %17, %16 ], [ %28, %33 ]
ret i32 %44
}