double interp1Q(double x[], double y[], int nx, double xi) {
  int i;
  if (xi <= x[0]) {
    return y[0];
  }
  if (xi >= x[nx - 1]) {
    return y[nx - 1];
  }
  i = 1;
  while (x[i] < xi) {
    i = i + 1;
  }
  return y[i - 1] + (xi - x[i - 1]) / (x[i] - x[i - 1]) * y[i];
}


x86 assembly (clang, -O0):
<interp1Q>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
ucomisd %xmm1,%xmm0
jb     113d <interp1Q+0x3d>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    120d <interp1Q+0x10d>
movsd  -0x28(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
ucomisd (%rax,%rcx,8),%xmm0
jb     1176 <interp1Q+0x76>
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    120d <interp1Q+0x10d>
movl   $0x1,-0x2c(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x28(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    11a7 <interp1Q+0xa7>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    117d <interp1Q+0x7d>
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x28(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm2
divsd  %xmm2,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<interp1Q>:
movsd  (%rdi),%xmm1
ucomisd %xmm0,%xmm1
jae    1165 <interp1Q+0x65>
movslq %edx,%rax
ucomisd -0x8(%rdi,%rax,8),%xmm0
jae    116a <interp1Q+0x6a>
movabs $0xffffffff00000000,%rax
ff ff ff
xor    %ecx,%ecx
movabs $0x100000000,%rdx
nopl   0x0(%rax,%rax,1)
add    %rdx,%rax
movsd  0x8(%rdi,%rcx,8),%xmm1
add    $0x1,%rcx
ucomisd %xmm1,%xmm0
ja     1130 <interp1Q+0x30>
mov    %ecx,%ecx
sar    $0x1d,%rax
movsd  (%rdi,%rax,1),%xmm2
subsd  %xmm2,%xmm0
subsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mulsd  (%rsi,%rcx,8),%xmm0
addsd  (%rsi,%rax,1),%xmm0
ret
movsd  (%rsi),%xmm0
ret
add    $0xffffffffffffffff,%rax
movsd  (%rsi,%rax,8),%xmm0
ret

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

LLVM IR (-O2):
define double @interp1Q(double* nocapture %0, double* nocapture %1, i32 %2, double %3) {
%5 = load double, double* %0
%6 = fcmp ult double %5, %3
br i1 %6, label %9, label %7
7: ; preds = %4
%8 = load double, double* %1
br label %39
9: ; preds = %4
%10 = add nsw i32 %2, -1
%11 = sext i32 %10 to i64
%12 = getelementptr inbounds double, double* %0, i64 %11
%13 = load double, double* %12
%14 = fcmp ugt double %13, %3
br i1 %14, label %18, label %15
15: ; preds = %9
%16 = getelementptr inbounds double, double* %1, i64 %11
%17 = load double, double* %16
br label %39
18: ; preds = %9, %18
%19 = phi i64 [ %23, %18 ], [ 1, %9 ]
%20 = getelementptr inbounds double, double* %0, i64 %19
%21 = load double, double* %20
%22 = fcmp olt double %21, %3
%23 = add nuw i64 %19, 1
br i1 %22, label %18, label %24
24: ; preds = %18
%25 = and i64 %19, 4294967295
%26 = shl i64 %19, 32
%27 = add i64 %26, -4294967296
%28 = ashr exact i64 %27, 32
%29 = getelementptr inbounds double, double* %1, i64 %28
%30 = load double, double* %29
%31 = getelementptr inbounds double, double* %0, i64 %28
%32 = load double, double* %31
%33 = fsub double %3, %32
%34 = fsub double %21, %32
%35 = fdiv double %33, %34
%36 = getelementptr inbounds double, double* %1, i64 %25
%37 = load double, double* %36
%38 = tail call double @llvm.fmuladd.f64(double %35, double %37, double %30)
br label %39
39: ; preds = %24, %15, %7
%40 = phi double [ %8, %7 ], [ %17, %15 ], [ %38, %24 ]
ret double %40
}