int findindex(int *indices, double *vals, double target, int nvals) {
  double ratio;
  double vlow, vhigh;
  int low, high;
  int new;
  if (target <= vals[indices[0]]) {
    return (0);
  }
  if (target >= vals[indices[nvals - 1]]) {
    return (nvals - 1);
  }
  low = 0;
  high = nvals - 1;
  while (high - low > 1) {
    vlow = vals[indices[low]];
    vhigh = vals[indices[high]];
    if (vlow == vhigh) {
      return ((vlow + vhigh) / 2);
    }
    ratio = (target - vlow) / (vhigh - vlow);
    new = low + ratio *(high - low);
    if (new == low) {
      ++new;
    } else if (new == high) {
      --new;
    }
    if (vals[indices[new]] < target) {
      low = new;
    } else {
      high = new;
    }
  }
  if (target == vals[indices[high]]) {
    return (high);
  } else {
    return (low);
  }
}


x86 assembly (clang, -O0):
<findindex>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
mov    %edx,-0x24(%rbp)
movsd  -0x20(%rbp),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq (%rcx),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jb     113f <findindex+0x3f>
movl   $0x0,-0x4(%rbp)
jmp    12e2 <findindex+0x1e2>
movsd  -0x20(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x10(%rbp),%rcx
mov    -0x24(%rbp),%edx
sub    $0x1,%edx
movslq %edx,%rdx
movslq (%rcx,%rdx,4),%rcx
ucomisd (%rax,%rcx,8),%xmm0
jb     1172 <findindex+0x72>
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x4(%rbp)
jmp    12e2 <findindex+0x1e2>
movl   $0x0,-0x44(%rbp)
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x48(%rbp)
mov    -0x48(%rbp),%eax
sub    -0x44(%rbp),%eax
cmp    $0x1,%eax
jle    12ab <findindex+0x1ab>
mov    -0x18(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq -0x44(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq -0x48(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm0
ucomisd -0x40(%rbp),%xmm0
jne    11fd <findindex+0xfd>
jp     11fd <findindex+0xfd>
movsd  -0x38(%rbp),%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  0xe13(%rip),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x4(%rbp)
jmp    12e2 <findindex+0x1e2>
movsd  -0x20(%rbp),%xmm0
subsd  -0x38(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm1
subsd  -0x38(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
cvtsi2sdl -0x44(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm0
mov    -0x48(%rbp),%eax
sub    -0x44(%rbp),%eax
cvtsi2sd %eax,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x44(%rbp),%eax
jne    1257 <findindex+0x157>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1271 <findindex+0x171>
mov    -0x4c(%rbp),%eax
cmp    -0x48(%rbp),%eax
jne    126c <findindex+0x16c>
mov    -0x4c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x4c(%rbp)
jmp    1271 <findindex+0x171>
mov    -0x18(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq -0x4c(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x20(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    12a0 <findindex+0x1a0>
mov    -0x4c(%rbp),%eax
mov    %eax,-0x44(%rbp)
jmp    12a6 <findindex+0x1a6>
mov    -0x4c(%rbp),%eax
mov    %eax,-0x48(%rbp)
jmp    1182 <findindex+0x82>
movsd  -0x20(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq -0x48(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
ucomisd (%rax,%rcx,8),%xmm0
jne    12dc <findindex+0x1dc>
jp     12dc <findindex+0x1dc>
mov    -0x48(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    12e2 <findindex+0x1e2>
mov    -0x44(%rbp),%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<findindex>:
movslq (%rdi),%rax
movsd  (%rsi,%rax,8),%xmm2
xor    %eax,%eax
ucomisd %xmm0,%xmm2
jae    1205 <findindex+0x105>
lea    -0x1(%rdx),%eax
movslq %edx,%rcx
movslq -0x4(%rdi,%rcx,4),%rcx
movsd  (%rsi,%rcx,8),%xmm1
ucomisd %xmm1,%xmm0
jae    1205 <findindex+0x105>
xor    %r10d,%r10d
cmp    $0x3,%edx
jl     11f9 <findindex+0xf9>
ucomisd %xmm1,%xmm2
jne    1146 <findindex+0x46>
jnp    11df <findindex+0xdf>
xor    %r10d,%r10d
mov    $0x1,%r8d
mov    %eax,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %eax,%r9d
movapd %xmm0,%xmm3
subsd  %xmm2,%xmm3
subsd  %xmm2,%xmm1
xorps  %xmm2,%xmm2
cvtsi2sd %r10d,%xmm2
divsd  %xmm1,%xmm3
xorps  %xmm1,%xmm1
cvtsi2sd %edx,%xmm1
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
cvttsd2si %xmm1,%eax
xor    %ecx,%ecx
cmp    %eax,%r9d
sete   %cl
neg    %ecx
cmp    %eax,%r10d
cmove  %r8d,%ecx
add    %eax,%ecx
movslq %ecx,%rax
movslq (%rdi,%rax,4),%rcx
ucomisd (%rsi,%rcx,8),%xmm0
cmova  %eax,%r10d
cmova  %r9d,%eax
mov    %eax,%edx
sub    %r10d,%edx
movslq %eax,%r9
cmp    $0x2,%edx
jl     11f0 <findindex+0xf0>
movslq %r10d,%rcx
movslq (%rdi,%rcx,4),%rcx
movsd  (%rsi,%rcx,8),%xmm2
movslq (%rdi,%r9,4),%rcx
movsd  (%rsi,%rcx,8),%xmm1
ucomisd %xmm1,%xmm2
jne    1160 <findindex+0x60>
jp     1160 <findindex+0x60>
addsd  %xmm2,%xmm1
mulsd  0xe15(%rip),%xmm1
cvttsd2si %xmm1,%eax
ret
movslq (%rdi,%r9,4),%rcx
movsd  (%rsi,%rcx,8),%xmm1
ucomisd %xmm0,%xmm1
cmovne %r10d,%eax
cmovp  %r10d,%eax
ret

LLVM IR (-O0):
define i32 @findindex(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 double
%11 = alloca double
%12 = alloca double
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
store i32* %0, i32** %6
store double* %1, double** %7
store double %2, double* %8
store i32 %3, i32* %9
%16 = load double, double* %8
%17 = load double*, double** %7
%18 = load i32*, i32** %6
%19 = getelementptr inbounds i32, i32* %18, i64 0
%20 = load i32, i32* %19
%21 = sext i32 %20 to i64
%22 = getelementptr inbounds double, double* %17, i64 %21
%23 = load double, double* %22
%24 = fcmp ole double %16, %23
br i1 %24, label %25, label %26
25: ; preds = %4
store i32 0, i32* %5
br label %142
26: ; preds = %4
%27 = load double, double* %8
%28 = load double*, double** %7
%29 = load i32*, i32** %6
%30 = load i32, i32* %9
%31 = sub nsw i32 %30, 1
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds i32, i32* %29, i64 %32
%34 = load i32, i32* %33
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds double, double* %28, i64 %35
%37 = load double, double* %36
%38 = fcmp oge double %27, %37
br i1 %38, label %39, label %42
39: ; preds = %26
%40 = load i32, i32* %9
%41 = sub nsw i32 %40, 1
store i32 %41, i32* %5
br label %142
42: ; preds = %26
store i32 0, i32* %13
%43 = load i32, i32* %9
%44 = sub nsw i32 %43, 1
store i32 %44, i32* %14
br label %45
45: ; preds = %125, %42
%46 = load i32, i32* %14
%47 = load i32, i32* %13
%48 = sub nsw i32 %46, %47
%49 = icmp sgt i32 %48, 1
br i1 %49, label %50, label %126
50: ; preds = %45
%51 = load double*, double** %7
%52 = load i32*, i32** %6
%53 = load i32, i32* %13
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds i32, i32* %52, i64 %54
%56 = load i32, i32* %55
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds double, double* %51, i64 %57
%59 = load double, double* %58
store double %59, double* %11
%60 = load double*, double** %7
%61 = load i32*, i32** %6
%62 = load i32, i32* %14
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds i32, i32* %61, i64 %63
%65 = load i32, i32* %64
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds double, double* %60, i64 %66
%68 = load double, double* %67
store double %68, double* %12
%69 = load double, double* %11
%70 = load double, double* %12
%71 = fcmp oeq double %69, %70
br i1 %71, label %72, label %78
72: ; preds = %50
%73 = load double, double* %11
%74 = load double, double* %12
%75 = fadd double %73, %74
%76 = fdiv double %75, 2.000000e+00
%77 = fptosi double %76 to i32
store i32 %77, i32* %5
br label %142
78: ; preds = %50
%79 = load double, double* %8
%80 = load double, double* %11
%81 = fsub double %79, %80
%82 = load double, double* %12
%83 = load double, double* %11
%84 = fsub double %82, %83
%85 = fdiv double %81, %84
store double %85, double* %10
%86 = load i32, i32* %13
%87 = sitofp i32 %86 to double
%88 = load double, double* %10
%89 = load i32, i32* %14
%90 = load i32, i32* %13
%91 = sub nsw i32 %89, %90
%92 = sitofp i32 %91 to double
%93 = call double @llvm.fmuladd.f64(double %88, double %92, double %87)
%94 = fptosi double %93 to i32
store i32 %94, i32* %15
%95 = load i32, i32* %15
%96 = load i32, i32* %13
%97 = icmp eq i32 %95, %96
br i1 %97, label %98, label %101
98: ; preds = %78
%99 = load i32, i32* %15
%100 = add nsw i32 %99, 1
store i32 %100, i32* %15
br label %109
101: ; preds = %78
%102 = load i32, i32* %15
%103 = load i32, i32* %14
%104 = icmp eq i32 %102, %103
br i1 %104, label %105, label %108
105: ; preds = %101
%106 = load i32, i32* %15
%107 = add nsw i32 %106, -1
store i32 %107, i32* %15
br label %108
108: ; preds = %105, %101
br label %109
109: ; preds = %108, %98
%110 = load double*, double** %7
%111 = load i32*, i32** %6
%112 = load i32, i32* %15
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds i32, i32* %111, i64 %113
%115 = load i32, i32* %114
%116 = sext i32 %115 to i64
%117 = getelementptr inbounds double, double* %110, i64 %116
%118 = load double, double* %117
%119 = load double, double* %8
%120 = fcmp olt double %118, %119
br i1 %120, label %121, label %123
121: ; preds = %109
%122 = load i32, i32* %15
store i32 %122, i32* %13
br label %125
123: ; preds = %109
%124 = load i32, i32* %15
store i32 %124, i32* %14
br label %125
125: ; preds = %123, %121
br label %45
126: ; preds = %45
%127 = load double, double* %8
%128 = load double*, double** %7
%129 = load i32*, i32** %6
%130 = load i32, i32* %14
%131 = sext i32 %130 to i64
%132 = getelementptr inbounds i32, i32* %129, i64 %131
%133 = load i32, i32* %132
%134 = sext i32 %133 to i64
%135 = getelementptr inbounds double, double* %128, i64 %134
%136 = load double, double* %135
%137 = fcmp oeq double %127, %136
br i1 %137, label %138, label %140
138: ; preds = %126
%139 = load i32, i32* %14
store i32 %139, i32* %5
br label %142
140: ; preds = %126
%141 = load i32, i32* %13
store i32 %141, i32* %5
br label %142
142: ; preds = %140, %138, %72, %39, %25
%143 = load i32, i32* %5
ret i32 %143
}

LLVM IR (-O2):
define i32 @findindex(i32* nocapture %0, double* nocapture %1, double %2, i32 %3) {
%5 = load i32, i32* %0
%6 = sext i32 %5 to i64
%7 = getelementptr inbounds double, double* %1, i64 %6
%8 = load double, double* %7
%9 = fcmp ult double %8, %2
br i1 %9, label %10, label %85
10: ; preds = %4
%11 = add nsw i32 %3, -1
%12 = sext i32 %11 to i64
%13 = getelementptr inbounds i32, i32* %0, i64 %12
%14 = load i32, i32* %13
%15 = sext i32 %14 to i64
%16 = getelementptr inbounds double, double* %1, i64 %15
%17 = load double, double* %16
%18 = fcmp ugt double %17, %2
br i1 %18, label %19, label %85
19: ; preds = %10
%20 = icmp sgt i32 %3, 2
br i1 %20, label %21, label %79
21: ; preds = %19
%22 = fcmp oeq double %8, %17
br i1 %22, label %23, label %29
23: ; preds = %58, %21
%24 = phi double [ %17, %21 ], [ %70, %58 ]
%25 = phi double [ %8, %21 ], [ %64, %58 ]
%26 = fadd double %25, %24
%27 = fmul double %26, 5.000000e-01
%28 = fptosi double %27 to i32
br label %85
29: ; preds = %21, %58
%30 = phi i32 [ %54, %58 ], [ 0, %21 ]
%31 = phi i32 [ %55, %58 ], [ %11, %21 ]
%32 = phi i32 [ %56, %58 ], [ %11, %21 ]
%33 = phi double [ %64, %58 ], [ %8, %21 ]
%34 = phi double [ %70, %58 ], [ %17, %21 ]
%35 = fsub double %2, %33
%36 = fsub double %34, %33
%37 = fdiv double %35, %36
%38 = sitofp i32 %30 to double
%39 = sitofp i32 %32 to double
%40 = tail call double @llvm.fmuladd.f64(double %37, double %39, double %38)
%41 = fptosi double %40 to i32
%42 = icmp eq i32 %30, %41
%43 = icmp eq i32 %31, %41
%44 = sext i1 %43 to i32
%45 = select i1 %42, i32 1, i32 %44
%46 = add nsw i32 %45, %41
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds i32, i32* %0, i64 %47
%49 = load i32, i32* %48
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds double, double* %1, i64 %50
%52 = load double, double* %51
%53 = fcmp olt double %52, %2
%54 = select i1 %53, i32 %46, i32 %30
%55 = select i1 %53, i32 %31, i32 %46
%56 = sub nsw i32 %55, %54
%57 = icmp sgt i32 %56, 1
br i1 %57, label %58, label %72
58: ; preds = %29
%59 = sext i32 %54 to i64
%60 = getelementptr inbounds i32, i32* %0, i64 %59
%61 = load i32, i32* %60
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds double, double* %1, i64 %62
%64 = load double, double* %63
%65 = sext i32 %55 to i64
%66 = getelementptr inbounds i32, i32* %0, i64 %65
%67 = load i32, i32* %66
%68 = sext i32 %67 to i64
%69 = getelementptr inbounds double, double* %1, i64 %68
%70 = load double, double* %69
%71 = fcmp oeq double %64, %70
br i1 %71, label %23, label %29
72: ; preds = %29
%73 = sext i32 %55 to i64
%74 = getelementptr inbounds i32, i32* %0, i64 %73
%75 = load i32, i32* %74
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds double, double* %1, i64 %76
%78 = load double, double* %77
br label %79
79: ; preds = %72, %19
%80 = phi double [ %78, %72 ], [ %17, %19 ]
%81 = phi i32 [ %54, %72 ], [ 0, %19 ]
%82 = phi i32 [ %55, %72 ], [ %11, %19 ]
%83 = fcmp oeq double %80, %2
%84 = select i1 %83, i32 %82, i32 %81
br label %85
85: ; preds = %79, %10, %4, %23
%86 = phi i32 [ %28, %23 ], [ 0, %4 ], [ %11, %10 ], [ %84, %79 ]
ret i32 %86
}