unsigned long find_max_idx_row_number(double max, double *data,
                                      unsigned long *idxs, unsigned long rows) {
  unsigned long l, u, i;
  double val;
  l = 0;
  u = rows;
  while (l < u) {
    i = (l + u) / 2;
    val = data[idxs[i]];
    if (val > max)
      u = i;
    else if (val < max)
      l = i + 1;
    else
      return i;
  }
  return i;
}


x86 assembly (clang, -O0):
<find_max_idx_row_number>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x10(%rbp)
mov    %rdi,-0x18(%rbp)
mov    %rsi,-0x20(%rbp)
mov    %rdx,-0x28(%rbp)
movq   $0x0,-0x30(%rbp)
mov    -0x28(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x38(%rbp),%rax
jae    11b6 <find_max_idx_row_number+0xb6>
mov    -0x30(%rbp),%rax
add    -0x38(%rbp),%rax
shr    $0x1,%rax
mov    %rax,-0x40(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%rcx
mov    -0x40(%rbp),%rdx
mov    (%rcx,%rdx,8),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
ucomisd -0x10(%rbp),%xmm0
jbe    117a <find_max_idx_row_number+0x7a>
mov    -0x40(%rbp),%rax
mov    %rax,-0x38(%rbp)
jmp    11b1 <find_max_idx_row_number+0xb1>
movsd  -0x48(%rbp),%xmm1
movsd  -0x10(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    119f <find_max_idx_row_number+0x9f>
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    11ac <find_max_idx_row_number+0xac>
mov    -0x40(%rbp),%rax
mov    %rax,-0x8(%rbp)
jmp    11be <find_max_idx_row_number+0xbe>
jmp    11b1 <find_max_idx_row_number+0xb1>
jmp    1125 <find_max_idx_row_number+0x25>
mov    -0x40(%rbp),%rax
mov    %rax,-0x8(%rbp)
mov    -0x8(%rbp),%rax
pop    %rbp
ret

x86 assembly (clang, -O2):
<find_max_idx_row_number>:
test   %rdx,%rdx
je     113f <find_max_idx_row_number+0x3f>
xor    %r8d,%r8d
jmp    111f <find_max_idx_row_number+0x1f>
nopw   0x0(%rax,%rax,1)
lea    0x1(%rax),%r8
mov    %rdx,%rcx
mov    %rcx,%rdx
cmp    %rcx,%r8
jae    113e <find_max_idx_row_number+0x3e>
lea    (%rdx,%r8,1),%rax
shr    %rax
mov    (%rsi,%rax,8),%rcx
movsd  (%rdi,%rcx,8),%xmm1
ucomisd %xmm0,%xmm1
mov    %rax,%rcx
ja     1117 <find_max_idx_row_number+0x17>
ucomisd %xmm1,%xmm0
ja     1110 <find_max_idx_row_number+0x10>
ret
ret

LLVM IR (-O0):
define i64 @find_max_idx_row_number(double %0, double* %1, i64* %2, i64 %3) {
%5 = alloca i64
%6 = alloca double
%7 = alloca double*
%8 = alloca i64*
%9 = alloca i64
%10 = alloca i64
%11 = alloca i64
%12 = alloca i64
%13 = alloca double
store double %0, double* %6
store double* %1, double** %7
store i64* %2, i64** %8
store i64 %3, i64* %9
store i64 0, i64* %10
%14 = load i64, i64* %9
store i64 %14, i64* %11
br label %15
15: ; preds = %46, %4
%16 = load i64, i64* %10
%17 = load i64, i64* %11
%18 = icmp ult i64 %16, %17
br i1 %18, label %19, label %47
19: ; preds = %15
%20 = load i64, i64* %10
%21 = load i64, i64* %11
%22 = add i64 %20, %21
%23 = udiv i64 %22, 2
store i64 %23, i64* %12
%24 = load double*, double** %7
%25 = load i64*, i64** %8
%26 = load i64, i64* %12
%27 = getelementptr inbounds i64, i64* %25, i64 %26
%28 = load i64, i64* %27
%29 = getelementptr inbounds double, double* %24, i64 %28
%30 = load double, double* %29
store double %30, double* %13
%31 = load double, double* %13
%32 = load double, double* %6
%33 = fcmp ogt double %31, %32
br i1 %33, label %34, label %36
34: ; preds = %19
%35 = load i64, i64* %12
store i64 %35, i64* %11
br label %46
36: ; preds = %19
%37 = load double, double* %13
%38 = load double, double* %6
%39 = fcmp olt double %37, %38
br i1 %39, label %40, label %43
40: ; preds = %36
%41 = load i64, i64* %12
%42 = add i64 %41, 1
store i64 %42, i64* %10
br label %45
43: ; preds = %36
%44 = load i64, i64* %12
store i64 %44, i64* %5
br label %49
45: ; preds = %40
br label %46
46: ; preds = %45, %34
br label %15
47: ; preds = %15
%48 = load i64, i64* %12
store i64 %48, i64* %5
br label %49
49: ; preds = %47, %43
%50 = load i64, i64* %5
ret i64 %50
}

LLVM IR (-O2):
define i64 @find_max_idx_row_number(double %0, double* nocapture %1, i64* nocapture %2, i64 %3) {
%5 = icmp eq i64 %3, 0
br i1 %5, label %24, label %6
6: ; preds = %4, %20
%7 = phi i64 [ %22, %20 ], [ %3, %4 ]
%8 = phi i64 [ %21, %20 ], [ 0, %4 ]
%9 = add i64 %7, %8
%10 = lshr i64 %9, 1
%11 = getelementptr inbounds i64, i64* %2, i64 %10
%12 = load i64, i64* %11
%13 = getelementptr inbounds double, double* %1, i64 %12
%14 = load double, double* %13
%15 = fcmp ogt double %14, %0
br i1 %15, label %20, label %16
16: ; preds = %6
%17 = fcmp olt double %14, %0
br i1 %17, label %18, label %24
18: ; preds = %16
%19 = add nuw i64 %10, 1
br label %20
20: ; preds = %6, %18
%21 = phi i64 [ %19, %18 ], [ %8, %6 ]
%22 = phi i64 [ %7, %18 ], [ %10, %6 ]
%23 = icmp ult i64 %21, %22
br i1 %23, label %6, label %24
24: ; preds = %16, %20, %4
%25 = phi i64 [ undef, %4 ], [ %10, %20 ], [ %10, %16 ]
ret i64 %25
}