int dp_lookup(double *T, int n, double t) {
  int l, m, r;
  if (t < T[n - 1]) {
    l = 0;
    r = n;
    m = (l + r) / 2;
    while (1) {
      if (t >= T[m + 1])
        l = m;
      else if (t < T[m])
        r = m;
      else
        break;
      m = (r + l) / 2;
    }
    return m;
  } else {
    return n - 2;
  }
}


x86 assembly (clang, -O0):
<dp_lookup>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    11c9 <dp_lookup+0xc9>
movl   $0x0,-0x24(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x24(%rbp),%eax
add    -0x2c(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x28(%rbp)
movsd  -0x20(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
ucomisd (%rax,%rcx,8),%xmm0
jb     1177 <dp_lookup+0x77>
mov    -0x28(%rbp),%eax
mov    %eax,-0x24(%rbp)
jmp    11a8 <dp_lookup+0xa8>
movsd  -0x20(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    119e <dp_lookup+0x9e>
mov    -0x28(%rbp),%eax
mov    %eax,-0x2c(%rbp)
jmp    11a3 <dp_lookup+0xa3>
jmp    11be <dp_lookup+0xbe>
jmp    11a8 <dp_lookup+0xa8>
mov    -0x2c(%rbp),%eax
add    -0x24(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x28(%rbp)
jmp    114f <dp_lookup+0x4f>
mov    -0x28(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    11d2 <dp_lookup+0xd2>
mov    -0x14(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<dp_lookup>:
movslq %esi,%rax
movsd  -0x8(%rdi,%rax,8),%xmm1
ucomisd %xmm0,%xmm1
jbe    1149 <dp_lookup+0x49>
xor    %ecx,%ecx
mov    %esi,%edx
jmp    1129 <dp_lookup+0x29>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %eax,%ecx
mov    %edx,%eax
lea    (%rax,%rcx,1),%esi
mov    %eax,%edx
mov    %esi,%eax
shr    $0x1f,%eax
add    %esi,%eax
sar    %eax
movslq %eax,%rsi
ucomisd 0x8(%rdi,%rsi,8),%xmm0
jae    1120 <dp_lookup+0x20>
movsd  (%rdi,%rsi,8),%xmm1
ucomisd %xmm0,%xmm1
ja     1124 <dp_lookup+0x24>
ret
add    $0xfffffffe,%esi
mov    %esi,%eax
ret

LLVM IR (-O0):
define i32 @dp_lookup(double* %0, i32 %1, double %2) {
%4 = alloca i32
%5 = alloca double*
%6 = alloca i32
%7 = alloca double
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
store double* %0, double** %5
store i32 %1, i32* %6
store double %2, double* %7
%11 = load double, double* %7
%12 = load double*, double** %5
%13 = load i32, i32* %6
%14 = sub nsw i32 %13, 1
%15 = sext i32 %14 to i64
%16 = getelementptr inbounds double, double* %12, i64 %15
%17 = load double, double* %16
%18 = fcmp olt double %11, %17
br i1 %18, label %19, label %55
19: ; preds = %3
store i32 0, i32* %8
%20 = load i32, i32* %6
store i32 %20, i32* %10
%21 = load i32, i32* %8
%22 = load i32, i32* %10
%23 = add nsw i32 %21, %22
%24 = sdiv i32 %23, 2
store i32 %24, i32* %9
br label %25
25: ; preds = %19, %48
%26 = load double, double* %7
%27 = load double*, double** %5
%28 = load i32, i32* %9
%29 = add nsw i32 %28, 1
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds double, double* %27, i64 %30
%32 = load double, double* %31
%33 = fcmp oge double %26, %32
br i1 %33, label %34, label %36
34: ; preds = %25
%35 = load i32, i32* %9
store i32 %35, i32* %8
br label %48
36: ; preds = %25
%37 = load double, double* %7
%38 = load double*, double** %5
%39 = load i32, i32* %9
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds double, double* %38, i64 %40
%42 = load double, double* %41
%43 = fcmp olt double %37, %42
br i1 %43, label %44, label %46
44: ; preds = %36
%45 = load i32, i32* %9
store i32 %45, i32* %10
br label %47
46: ; preds = %36
br label %53
47: ; preds = %44
br label %48
48: ; preds = %47, %34
%49 = load i32, i32* %10
%50 = load i32, i32* %8
%51 = add nsw i32 %49, %50
%52 = sdiv i32 %51, 2
store i32 %52, i32* %9
br label %25
53: ; preds = %46
%54 = load i32, i32* %9
store i32 %54, i32* %4
br label %58
55: ; preds = %3
%56 = load i32, i32* %6
%57 = sub nsw i32 %56, 2
store i32 %57, i32* %4
br label %58
58: ; preds = %55, %53
%59 = load i32, i32* %4
ret i32 %59
}

LLVM IR (-O2):
define i32 @dp_lookup(double* nocapture %0, i32 %1, double %2) {
%4 = add nsw i32 %1, -1
%5 = sext i32 %4 to i64
%6 = getelementptr inbounds double, double* %0, i64 %5
%7 = load double, double* %6
%8 = fcmp ogt double %7, %2
br i1 %8, label %9, label %28
9: ; preds = %3, %24
%10 = phi i32 [ %25, %24 ], [ 0, %3 ]
%11 = phi i32 [ %27, %24 ], [ %1, %3 ]
%12 = phi i32 [ %26, %24 ], [ %1, %3 ]
%13 = sdiv i32 %11, 2
%14 = add nsw i32 %13, 1
%15 = sext i32 %14 to i64
%16 = getelementptr inbounds double, double* %0, i64 %15
%17 = load double, double* %16
%18 = fcmp ugt double %17, %2
br i1 %18, label %19, label %24
19: ; preds = %9
%20 = sext i32 %13 to i64
%21 = getelementptr inbounds double, double* %0, i64 %20
%22 = load double, double* %21
%23 = fcmp ogt double %22, %2
br i1 %23, label %24, label %30
24: ; preds = %19, %9
%25 = phi i32 [ %13, %9 ], [ %10, %19 ]
%26 = phi i32 [ %12, %9 ], [ %13, %19 ]
%27 = add nsw i32 %26, %25
br label %9
28: ; preds = %3
%29 = add nsw i32 %1, -2
br label %30
30: ; preds = %19, %28
%31 = phi i32 [ %29, %28 ], [ %13, %19 ]
ret i32 %31
}