double interpolate_linear(int tCount, double *tValues, double *yValues,
                          double t) {
  int i, prevI = 0, nextI = 0;
  double prevT = tValues[0], nextT = 0;
  if (tCount == 1) {
    return yValues[0];
  }
  if (t <= tValues[0]) {
    return yValues[0];
  } else if (t >= tValues[tCount - 1]) {
    return yValues[tCount - 1];
  }
  for (i = 1; i < tCount; i++) {
    nextI = i;
    nextT = tValues[i];
    if (t >= prevT && t <= nextT)
      break;
    prevT = nextT;
    prevI = nextI;
  }
  return yValues[prevI] +
         (t - prevT) * (yValues[nextI] - yValues[prevI]) / (nextT - prevT);
}


x86 assembly (clang, -O0):
<interpolate_linear>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0xc(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
movsd  %xmm0,-0x28(%rbp)
movl   $0x0,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
cmpl   $0x1,-0xc(%rbp)
jne    1153 <interpolate_linear+0x53>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1278 <interpolate_linear+0x178>
movsd  -0x28(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
ucomisd %xmm1,%xmm0
jb     117c <interpolate_linear+0x7c>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1278 <interpolate_linear+0x178>
movsd  -0x28(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
ucomisd (%rax,%rcx,8),%xmm0
jb     11b5 <interpolate_linear+0xb5>
mov    -0x20(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1278 <interpolate_linear+0x178>
jmp    11ba <interpolate_linear+0xba>
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    122c <interpolate_linear+0x12c>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x28(%rbp),%xmm0
ucomisd -0x40(%rbp),%xmm0
jb     120e <interpolate_linear+0x10e>
movsd  -0x28(%rbp),%xmm1
movsd  -0x48(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jb     120e <interpolate_linear+0x10e>
jmp    122c <interpolate_linear+0x12c>
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11c1 <interpolate_linear+0xc1>
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x28(%rbp),%xmm1
subsd  -0x40(%rbp),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm1
movsd  -0x48(%rbp),%xmm2
subsd  -0x40(%rbp),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<interpolate_linear>:
cmp    $0x1,%edi
je     113f <interpolate_linear+0x3f>
movsd  (%rsi),%xmm2
ucomisd %xmm0,%xmm2
jae    113f <interpolate_linear+0x3f>
lea    -0x1(%rdi),%eax
movslq %eax,%rcx
ucomisd (%rsi,%rcx,8),%xmm0
jae    1144 <interpolate_linear+0x44>
xorpd  %xmm3,%xmm3
xor    %ecx,%ecx
cmp    $0x2,%edi
jl     114a <interpolate_linear+0x4a>
movsd  0x8(%rsi),%xmm3
ucomisd %xmm2,%xmm0
jb     114e <interpolate_linear+0x4e>
ucomisd %xmm0,%xmm3
jb     114e <interpolate_linear+0x4e>
mov    $0x1,%eax
jmp    118a <interpolate_linear+0x8a>
movsd  (%rdx),%xmm0
ret
movsd  (%rdx,%rcx,8),%xmm0
ret
xor    %eax,%eax
jmp    118a <interpolate_linear+0x8a>
mov    %edi,%edi
add    $0xfffffffffffffffe,%rdi
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
movapd %xmm3,%xmm2
cmp    %rcx,%rdi
je     1184 <interpolate_linear+0x84>
movsd  0x10(%rsi,%rcx,8),%xmm3
add    $0x1,%rcx
ucomisd %xmm2,%xmm0
jb     1160 <interpolate_linear+0x60>
ucomisd %xmm0,%xmm3
jb     1160 <interpolate_linear+0x60>
lea    0x1(%rcx),%eax
jmp    118a <interpolate_linear+0x8a>
mov    %eax,%ecx
movapd %xmm2,%xmm3
mov    %ecx,%ecx
movsd  (%rdx,%rcx,8),%xmm4
subsd  %xmm2,%xmm0
mov    %eax,%eax
movsd  (%rdx,%rax,8),%xmm1
subsd  %xmm4,%xmm1
mulsd  %xmm0,%xmm1
subsd  %xmm2,%xmm3
divsd  %xmm3,%xmm1
addsd  %xmm4,%xmm1
movapd %xmm1,%xmm0
ret

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

LLVM IR (-O2):
define double @interpolate_linear(i32 %0, double* nocapture %1, double* nocapture %2, double %3) {
%5 = load double, double* %1
%6 = icmp eq i32 %0, 1
br i1 %6, label %7, label %9
7: ; preds = %4
%8 = load double, double* %2
br label %62
9: ; preds = %4
%10 = fcmp ult double %5, %3
br i1 %10, label %13, label %11
11: ; preds = %9
%12 = load double, double* %2
br label %62
13: ; preds = %9
%14 = add i32 %0, -1
%15 = sext i32 %14 to i64
%16 = getelementptr inbounds double, double* %1, i64 %15
%17 = load double, double* %16
%18 = fcmp ugt double %17, %3
br i1 %18, label %19, label %28
19: ; preds = %13
%20 = icmp sgt i32 %0, 1
br i1 %20, label %21, label %45
21: ; preds = %19
%22 = zext i32 %0 to i64
%23 = getelementptr inbounds double, double* %1, i64 1
%24 = load double, double* %23
%25 = fcmp ugt double %5, %3
%26 = fcmp ult double %24, %3
%27 = select i1 %25, i1 true, i1 %26
br i1 %27, label %37, label %45
28: ; preds = %13
%29 = getelementptr inbounds double, double* %2, i64 %15
%30 = load double, double* %29
br label %62
31: ; preds = %37
%32 = getelementptr inbounds double, double* %1, i64 %40
%33 = load double, double* %32
%34 = fcmp ugt double %38, %3
%35 = fcmp ult double %33, %3
%36 = select i1 %34, i1 true, i1 %35
br i1 %36, label %37, label %42
37: ; preds = %21, %31
%38 = phi double [ %33, %31 ], [ %24, %21 ]
%39 = phi i64 [ %40, %31 ], [ 1, %21 ]
%40 = add nuw nsw i64 %39, 1
%41 = icmp eq i64 %40, %22
br i1 %41, label %45, label %31
42: ; preds = %31
%43 = trunc i64 %39 to i32
%44 = trunc i64 %40 to i32
br label %45
45: ; preds = %37, %21, %42, %19
%46 = phi i32 [ 0, %19 ], [ %43, %42 ], [ 0, %21 ], [ %14, %37 ]
%47 = phi double [ %5, %19 ], [ %38, %42 ], [ %5, %21 ], [ %38, %37 ]
%48 = phi i32 [ 0, %19 ], [ %44, %42 ], [ 1, %21 ], [ %14, %37 ]
%49 = phi double [ 0.000000e+00, %19 ], [ %33, %42 ], [ %24, %21 ], [ %38, %37 ]
%50 = zext i32 %46 to i64
%51 = getelementptr inbounds double, double* %2, i64 %50
%52 = load double, double* %51
%53 = fsub double %3, %47
%54 = zext i32 %48 to i64
%55 = getelementptr inbounds double, double* %2, i64 %54
%56 = load double, double* %55
%57 = fsub double %56, %52
%58 = fmul double %53, %57
%59 = fsub double %49, %47
%60 = fdiv double %58, %59
%61 = fadd double %52, %60
br label %62
62: ; preds = %45, %28, %11, %7
%63 = phi double [ %8, %7 ], [ %12, %11 ], [ %30, %28 ], [ %61, %45 ]
ret double %63
}