double interp1Q_cas(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] + ((y[i] - y[i - 1]) / (x[i] - x[i - 1])) * (xi - x[i - 1]);
}


x86 assembly (clang, -O0):
<interp1Q_cas>:
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_cas+0x3d>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    121f <interp1Q_cas+0x11f>
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_cas+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    121f <interp1Q_cas+0x11f>
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_cas+0xa7>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    117d <interp1Q_cas+0x7d>
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%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
movsd  -0x28(%rbp),%xmm2
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%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_cas>:
movsd  (%rdi),%xmm1
ucomisd %xmm0,%xmm1
jae    1171 <interp1Q_cas+0x71>
movslq %edx,%rax
ucomisd -0x8(%rdi,%rax,8),%xmm0
jae    1176 <interp1Q_cas+0x76>
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_cas+0x30>
mov    %ecx,%ecx
sar    $0x1d,%rax
movsd  (%rsi,%rax,1),%xmm2
movsd  (%rsi,%rcx,8),%xmm3
subsd  %xmm2,%xmm3
movsd  (%rdi,%rax,1),%xmm4
subsd  %xmm4,%xmm1
divsd  %xmm1,%xmm3
subsd  %xmm4,%xmm0
mulsd  %xmm3,%xmm0
addsd  %xmm2,%xmm0
ret
movsd  (%rsi),%xmm0
ret
add    $0xffffffffffffffff,%rax
movsd  (%rsi,%rax,8),%xmm0
ret

LLVM IR (-O0):
define double @interp1Q_cas(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 %89
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 %89
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** %7
%56 = load i32, i32* %10
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds double, double* %55, i64 %57
%59 = load double, double* %58
%60 = load double*, double** %7
%61 = load i32, i32* %10
%62 = sub nsw i32 %61, 1
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds double, double* %60, i64 %63
%65 = load double, double* %64
%66 = fsub double %59, %65
%67 = load double*, double** %6
%68 = load i32, i32* %10
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds double, double* %67, i64 %69
%71 = load double, double* %70
%72 = load double*, double** %6
%73 = load i32, i32* %10
%74 = sub nsw i32 %73, 1
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds double, double* %72, i64 %75
%77 = load double, double* %76
%78 = fsub double %71, %77
%79 = fdiv double %66, %78
%80 = load double, double* %9
%81 = load double*, double** %6
%82 = load i32, i32* %10
%83 = sub nsw i32 %82, 1
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds double, double* %81, i64 %84
%86 = load double, double* %85
%87 = fsub double %80, %86
%88 = call double @llvm.fmuladd.f64(double %79, double %87, double %54)
store double %88, double* %5
br label %89
89: ; preds = %48, %29, %16
%90 = load double, double* %5
ret double %90
}

LLVM IR (-O2):
define double @interp1Q_cas(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 %40
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 %40
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* %1, i64 %25
%32 = load double, double* %31
%33 = fsub double %32, %30
%34 = getelementptr inbounds double, double* %0, i64 %28
%35 = load double, double* %34
%36 = fsub double %21, %35
%37 = fdiv double %33, %36
%38 = fsub double %3, %35
%39 = tail call double @llvm.fmuladd.f64(double %37, double %38, double %30)
br label %40
40: ; preds = %24, %15, %7
%41 = phi double [ %8, %7 ], [ %17, %15 ], [ %39, %24 ]
ret double %41
}