double *linear_interpolate(double *buffer, int x1, int x2) {
  double slope;
  double y1, y2;
  register int i;
  double new_value;
  double *anchor = buffer;
  y1 = buffer[x1];
  y2 = buffer[x2];
  slope = (y2 - y1) / (double)(x2 - x1);
  new_value = y1 + slope;
  buffer += (x1 + 1);
  for (i = x1 + 1; i < x2; i++) {
    *buffer++ = new_value;
    new_value += slope;
  }
  return anchor;
}


x86 assembly (clang, -O0):
<linear_interpolate>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
subsd  -0x20(%rbp),%xmm0
mov    -0x10(%rbp),%eax
sub    -0xc(%rbp),%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x20(%rbp),%xmm0
addsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
mov    -0x8(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x8(%rbp)
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    11c8 <linear_interpolate+0xc8>
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x8(%rbp)
movsd  %xmm0,(%rax)
movsd  -0x18(%rbp),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1187 <linear_interpolate+0x87>
mov    -0x40(%rbp),%rax
pop    %rbp
ret

x86 assembly (clang, -O2):
<linear_interpolate>:
mov    %rdi,%rax
lea    0x1(%rsi),%r8d
cmp    %edx,%r8d
jge    11d7 <linear_interpolate+0xd7>
movslq %esi,%rcx
movsd  (%rax,%rcx,8),%xmm0
movslq %edx,%r9
movsd  (%rax,%r9,8),%xmm1
sub    %ecx,%r9d
cvtsi2sd %r9d,%xmm2
subsd  %xmm0,%xmm1
divsd  %xmm2,%xmm1
movslq %r8d,%rcx
lea    (%rax,%rcx,8),%rcx
not    %esi
add    %edx,%esi
add    $0xfffffffe,%r9d
and    $0x7,%esi
je     1166 <linear_interpolate+0x66>
neg    %esi
xor    %edi,%edi
nopl   0x0(%rax)
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rcx)
add    $0x8,%rcx
add    $0xffffffff,%edi
cmp    %edi,%esi
jne    1150 <linear_interpolate+0x50>
sub    %edi,%r8d
cmp    $0x7,%r9d
jb     11d7 <linear_interpolate+0xd7>
sub    %r8d,%edx
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rcx,%rsi,8)
addsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%rcx,%rsi,8)
addsd  %xmm1,%xmm0
movsd  %xmm0,0x10(%rcx,%rsi,8)
addsd  %xmm1,%xmm0
movsd  %xmm0,0x18(%rcx,%rsi,8)
addsd  %xmm1,%xmm0
movsd  %xmm0,0x20(%rcx,%rsi,8)
addsd  %xmm1,%xmm0
movsd  %xmm0,0x28(%rcx,%rsi,8)
addsd  %xmm1,%xmm0
movsd  %xmm0,0x30(%rcx,%rsi,8)
addsd  %xmm1,%xmm0
movsd  %xmm0,0x38(%rcx,%rsi,8)
add    $0x8,%rsi
cmp    %esi,%edx
jne    1180 <linear_interpolate+0x80>
ret

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

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