double interp(double x, double *xs, double *ys, int count) {
  int i;
  double dx, dy;
  if (x < xs[0]) {
    return 1.0;
  }
  if (x > xs[count - 1]) {
    return 0.0;
  }
  for (i = 0; i < count - 1; i++) {
    if (xs[i + 1] > x) {
      break;
    }
  }
  dx = xs[i + 1] - xs[i];
  dy = ys[i + 1] - ys[i];
  return ys[i] + (x - xs[i]) * dy / dx;
}


x86 assembly (clang, -O0):
<interp>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x10(%rbp)
mov    %rdi,-0x18(%rbp)
mov    %rsi,-0x20(%rbp)
mov    %edx,-0x24(%rbp)
movsd  -0x10(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
ucomisd %xmm1,%xmm0
jbe    113d <interp+0x3d>
movsd  0xecd(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    122e <interp+0x12e>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    1167 <interp+0x67>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    122e <interp+0x12e>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    11b4 <interp+0xb4>
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x10(%rbp),%xmm0
jbe    11a1 <interp+0xa1>
jmp    11b4 <interp+0xb4>
jmp    11a6 <interp+0xa6>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    116e <interp+0x6e>
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x10(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
mulsd  -0x38(%rbp),%xmm1
divsd  -0x30(%rbp),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<interp>:
movsd  (%rdi),%xmm1
ucomisd %xmm0,%xmm1
jbe    1113 <interp+0x13>
movsd  0xeee(%rip),%xmm0
ret
movslq %edx,%rax
xorpd  %xmm1,%xmm1
ucomisd -0x8(%rdi,%rax,8),%xmm0
ja     1197 <interp+0x97>
add    $0xffffffff,%edx
xor    %r8d,%r8d
test   %edx,%edx
cmovg  %edx,%r8d
mov    %r8,%rdx
neg    %rdx
mov    $0xffffffffffffffff,%rax
nopl   0x0(%rax,%rax,1)
lea    (%rdx,%rax,1),%rcx
cmp    $0xffffffffffffffff,%rcx
je     115f <interp+0x5f>
movsd  0x10(%rdi,%rax,8),%xmm1
add    $0x1,%rax
ucomisd %xmm0,%xmm1
jbe    1140 <interp+0x40>
mov    %eax,%r8d
jmp    1162 <interp+0x62>
mov    %r8,%rax
mov    %r8d,%ecx
movsd  0x8(%rdi,%rcx,8),%xmm1
movsd  (%rdi,%rax,8),%xmm2
subsd  %xmm2,%xmm1
movsd  0x8(%rsi,%rcx,8),%xmm3
movsd  (%rsi,%rax,8),%xmm4
subsd  %xmm4,%xmm3
subsd  %xmm2,%xmm0
mulsd  %xmm3,%xmm0
divsd  %xmm1,%xmm0
addsd  %xmm4,%xmm0
movapd %xmm0,%xmm1
movapd %xmm1,%xmm0
ret

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

LLVM IR (-O2):
define double @interp(double %0, double* nocapture %1, double* nocapture %2, i32 %3) {
%5 = load double, double* %1
%6 = fcmp ogt double %5, %0
br i1 %6, label %46, label %7
7: ; preds = %4
%8 = add nsw i32 %3, -1
%9 = sext i32 %8 to i64
%10 = getelementptr inbounds double, double* %1, i64 %9
%11 = load double, double* %10
%12 = fcmp olt double %11, %0
br i1 %12, label %46, label %13
13: ; preds = %7
%14 = call i32 @llvm.smax.i32(i32 %8, i32 0)
%15 = zext i32 %14 to i64
br label %16
16: ; preds = %13, %19
%17 = phi i64 [ 0, %13 ], [ %20, %19 ]
%18 = icmp eq i64 %17, %15
br i1 %18, label %27, label %19
19: ; preds = %16
%20 = add nuw nsw i64 %17, 1
%21 = getelementptr inbounds double, double* %1, i64 %20
%22 = load double, double* %21
%23 = fcmp ogt double %22, %0
br i1 %23, label %24, label %16
24: ; preds = %19
%25 = trunc i64 %17 to i32
%26 = and i64 %17, 4294967295
br label %27
27: ; preds = %16, %24
%28 = phi i64 [ %26, %24 ], [ %15, %16 ]
%29 = phi i32 [ %25, %24 ], [ %14, %16 ]
%30 = add nuw nsw i32 %29, 1
%31 = zext i32 %30 to i64
%32 = getelementptr inbounds double, double* %1, i64 %31
%33 = load double, double* %32
%34 = getelementptr inbounds double, double* %1, i64 %28
%35 = load double, double* %34
%36 = fsub double %33, %35
%37 = getelementptr inbounds double, double* %2, i64 %31
%38 = load double, double* %37
%39 = getelementptr inbounds double, double* %2, i64 %28
%40 = load double, double* %39
%41 = fsub double %38, %40
%42 = fsub double %0, %35
%43 = fmul double %42, %41
%44 = fdiv double %43, %36
%45 = fadd double %40, %44
br label %46
46: ; preds = %7, %4, %27
%47 = phi double [ %45, %27 ], [ 1.000000e+00, %4 ], [ 0.000000e+00, %7 ]
ret double %47
}