int SearchGrid(double x, const double xval[], int Mxval) {
  int below;
  double xlow, xhigh;
  xlow = xval[0];
  xhigh = xval[Mxval - 1];
  below = (int)((double)(Mxval - 2) * (x - xlow) / (xhigh - xlow));
  if (below < 0)
    below = 0;
  else if (below > Mxval - 2)
    below = Mxval - 2;
  while (x > xval[below + 1])
    below = below + 1;
  while (x < xval[below])
    below = below - 1;
  return below;
}


x86 assembly (clang, -O0):
<SearchGrid>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x14(%rbp),%eax
sub    $0x2,%eax
cvtsi2sd %eax,%xmm0
movsd  -0x8(%rbp),%xmm1
subsd  -0x20(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x28(%rbp),%xmm1
subsd  -0x20(%rbp),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x18(%rbp)
cmpl   $0x0,-0x18(%rbp)
jge    1177 <SearchGrid+0x77>
movl   $0x0,-0x18(%rbp)
jmp    1196 <SearchGrid+0x96>
mov    -0x18(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x2,%ecx
cmp    %ecx,%eax
jle    1191 <SearchGrid+0x91>
mov    -0x14(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x18(%rbp)
jmp    1196 <SearchGrid+0x96>
jmp    119b <SearchGrid+0x9b>
movsd  -0x8(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x18(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    11c6 <SearchGrid+0xc6>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    119b <SearchGrid+0x9b>
jmp    11cb <SearchGrid+0xcb>
movsd  -0x8(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    11f5 <SearchGrid+0xf5>
mov    -0x18(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    11cb <SearchGrid+0xcb>
mov    -0x18(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<SearchGrid>:
movslq %esi,%rcx
movsd  -0x8(%rdi,%rcx,8),%xmm1
add    $0xfffffffe,%ecx
cvtsi2sd %ecx,%xmm2
movsd  (%rdi),%xmm3
movapd %xmm0,%xmm4
subsd  %xmm3,%xmm4
mulsd  %xmm2,%xmm4
subsd  %xmm3,%xmm1
divsd  %xmm1,%xmm4
cvttsd2si %xmm4,%edx
cmp    %edx,%ecx
cmovge %edx,%ecx
xor    %eax,%eax
test   %edx,%edx
cmovns %ecx,%eax
movslq %eax,%rcx
lea    (%rdi,%rcx,8),%rcx
add    $0x8,%rcx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%eax
ucomisd (%rcx),%xmm0
lea    0x8(%rcx),%rcx
ja     1150 <SearchGrid+0x50>
lea    -0x1(%rax),%ecx
movslq %ecx,%rcx
lea    (%rdi,%rcx,8),%rcx
nopw   0x0(%rax,%rax,1)
movsd  (%rcx),%xmm1
add    $0xffffffff,%eax
add    $0xfffffffffffffff8,%rcx
ucomisd %xmm0,%xmm1
ja     1170 <SearchGrid+0x70>
ret

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

LLVM IR (-O2):
define i32 @SearchGrid(double %0, double* nocapture %1, i32 %2) {
%4 = load double, double* %1
%5 = add nsw i32 %2, -1
%6 = sext i32 %5 to i64
%7 = getelementptr inbounds double, double* %1, i64 %6
%8 = load double, double* %7
%9 = add nsw i32 %2, -2
%10 = sitofp i32 %9 to double
%11 = fsub double %0, %4
%12 = fmul double %11, %10
%13 = fsub double %8, %4
%14 = fdiv double %12, %13
%15 = fptosi double %14 to i32
%16 = icmp slt i32 %15, 0
%17 = icmp slt i32 %9, %15
%18 = select i1 %17, i32 %9, i32 %15
%19 = select i1 %16, i32 0, i32 %18
%20 = sext i32 %19 to i64
br label %21
21: ; preds = %21, %3
%22 = phi i64 [ %23, %21 ], [ %20, %3 ]
%23 = add nsw i64 %22, 1
%24 = getelementptr inbounds double, double* %1, i64 %23
%25 = load double, double* %24
%26 = fcmp olt double %25, %0
br i1 %26, label %21, label %27
27: ; preds = %21
%28 = shl i64 %22, 32
%29 = ashr exact i64 %28, 32
br label %30
30: ; preds = %27, %30
%31 = phi i64 [ %29, %27 ], [ %35, %30 ]
%32 = getelementptr inbounds double, double* %1, i64 %31
%33 = load double, double* %32
%34 = fcmp ogt double %33, %0
%35 = add i64 %31, -1
br i1 %34, label %30, label %36
36: ; preds = %30
%37 = trunc i64 %31 to i32
ret i32 %37
}