void interp(float *y, float *x, int n, int r) {
  int i, j;
  float a;
  for (i = 1; i <= n - 1; i++) {
    for (j = 1; j <= r; j++) {
      a = (j - 1) * 1.0 / r;
      y[(i - 1) * r + j] = (1.0 - a) * x[i] + a * x[i + 1];
    }
  }
}


x86 assembly (clang, -O0):
<interp>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     11ec <interp+0xec>
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     11d9 <interp+0xd9>
mov    -0x20(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm0
movsd  0xeb1(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x18(%rbp),%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x24(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe8a(%rip),%xmm0
subsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm2
movss  -0x24(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mulss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x18(%rbp),%ecx
add    -0x20(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1131 <interp+0x31>
jmp    11de <interp+0xde>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1119 <interp+0x19>
pop    %rbp
ret

x86 assembly (clang, -O2):
<interp>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rdi,-0x28(%rsp)
cmp    $0x2,%edx
jl     1324 <interp+0x224>
cvtsi2sd %ecx,%xmm0
lea    0x1(%rcx),%r8d
movslq %ecx,%rdi
mov    %edx,%r12d
mov    -0x28(%rsp),%rax
lea    0x4(%rax),%rdx
mov    %rdx,-0x8(%rsp)
lea    (%rax,%r8,4),%rax
mov    %rax,-0x18(%rsp)
lea    0x8(%rsi),%rax
mov    %rax,-0x10(%rsp)
lea    -0x1(%r8),%r11
mov    %r11,%r14
and    $0xfffffffffffffffc,%r14
lea    0x1(%r14),%rax
mov    %rax,-0x20(%rsp)
movapd %xmm0,%xmm13
unpcklpd %xmm0,%xmm13
mov    $0x1,%edx
movsd  0xebe(%rip),%xmm12
movdqa 0xe85(%rip),%xmm8
pcmpeqd %xmm9,%xmm9
movapd 0xe87(%rip),%xmm11
movdqa 0xe8e(%rip),%xmm10
jmp    11b3 <interp+0xb3>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%rdx
mov    %rdx,%r13
mov    %r13,%rdx
cmp    %r12,%r13
je     1324 <interp+0x224>
test   %ecx,%ecx
jle    11a0 <interp+0xa0>
lea    0x1(%rdx),%r13
lea    -0x1(%rdx),%r9
imul   %rdi,%r9
mov    $0x1,%ebp
cmp    $0x4,%r11
jb     12c0 <interp+0x1c0>
mov    -0x8(%rsp),%rax
lea    (%rax,%r9,4),%r15
mov    -0x10(%rsp),%rax
lea    (%rax,%rdx,4),%rbx
cmp    %rbx,%r15
jae    11ff <interp+0xff>
lea    (%rsi,%rdx,4),%r10
mov    -0x18(%rsp),%rax
lea    (%rax,%r9,4),%rbx
cmp    %rbx,%r10
jb     12c0 <interp+0x1c0>
movss  (%rsi,%rdx,4),%xmm3
movss  0x4(%rsi,%rdx,4),%xmm14
shufps $0x0,%xmm3,%xmm3
cvtps2pd %xmm3,%xmm3
shufps $0x0,%xmm14,%xmm14
xor    %r10d,%r10d
movdqa %xmm8,%xmm6
nop
movdqa %xmm6,%xmm4
paddd  %xmm9,%xmm4
pshufd $0xee,%xmm4,%xmm5
cvtdq2pd %xmm4,%xmm4
cvtdq2pd %xmm5,%xmm5
divpd  %xmm13,%xmm4
divpd  %xmm13,%xmm5
cvtpd2ps %xmm4,%xmm4
cvtpd2ps %xmm5,%xmm5
unpcklpd %xmm5,%xmm4
cvtps2pd %xmm4,%xmm2
cvtps2pd %xmm5,%xmm5
movapd %xmm11,%xmm1
movapd %xmm11,%xmm7
subpd  %xmm2,%xmm1
subpd  %xmm5,%xmm7
mulps  %xmm14,%xmm4
cvtps2pd %xmm4,%xmm2
movhlps %xmm4,%xmm4
mulpd  %xmm3,%xmm1
cvtps2pd %xmm4,%xmm4
addpd  %xmm2,%xmm1
mulpd  %xmm3,%xmm7
cvtpd2ps %xmm1,%xmm1
addpd  %xmm4,%xmm7
cvtpd2ps %xmm7,%xmm2
unpcklpd %xmm2,%xmm1
movupd %xmm1,(%r15,%r10,4)
add    $0x4,%r10
paddd  %xmm10,%xmm6
cmp    %r10,%r14
jne    1220 <interp+0x120>
mov    -0x20(%rsp),%rbp
cmp    %r14,%r11
je     11a7 <interp+0xa7>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    -0x28(%rsp),%rax
lea    (%rax,%r9,4),%rbx
nopl   0x0(%rax)
lea    -0x1(%rbp),%eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
divsd  %xmm0,%xmm1
cvtsd2ss %xmm1,%xmm1
xorps  %xmm2,%xmm2
cvtss2sd %xmm1,%xmm2
movapd %xmm12,%xmm3
subsd  %xmm2,%xmm3
movss  (%rsi,%rdx,4),%xmm2
cvtss2sd %xmm2,%xmm2
mulsd  %xmm3,%xmm2
mulss  0x4(%rsi,%rdx,4),%xmm1
cvtss2sd %xmm1,%xmm1
addsd  %xmm2,%xmm1
cvtsd2ss %xmm1,%xmm1
movss  %xmm1,(%rbx,%rbp,4)
add    $0x1,%rbp
cmp    %rbp,%r8
jne    12d0 <interp+0x1d0>
jmp    11a7 <interp+0xa7>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @interp(float* %0, float* %1, i32 %2, i32 %3) {
%5 = alloca float*
%6 = alloca float*
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca float
store float* %0, float** %5
store float* %1, float** %6
store i32 %2, i32* %7
store i32 %3, i32* %8
store i32 1, i32* %9
br label %12
12: ; preds = %64, %4
%13 = load i32, i32* %9
%14 = load i32, i32* %7
%15 = sub nsw i32 %14, 1
%16 = icmp sle i32 %13, %15
br i1 %16, label %17, label %67
17: ; preds = %12
store i32 1, i32* %10
br label %18
18: ; preds = %60, %17
%19 = load i32, i32* %10
%20 = load i32, i32* %8
%21 = icmp sle i32 %19, %20
br i1 %21, label %22, label %63
22: ; preds = %18
%23 = load i32, i32* %10
%24 = sub nsw i32 %23, 1
%25 = sitofp i32 %24 to double
%26 = fmul double %25, 1.000000e+00
%27 = load i32, i32* %8
%28 = sitofp i32 %27 to double
%29 = fdiv double %26, %28
%30 = fptrunc double %29 to float
store float %30, float* %11
%31 = load float, float* %11
%32 = fpext float %31 to double
%33 = fsub double 1.000000e+00, %32
%34 = load float*, float** %6
%35 = load i32, i32* %9
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds float, float* %34, i64 %36
%38 = load float, float* %37
%39 = fpext float %38 to double
%40 = load float, float* %11
%41 = load float*, float** %6
%42 = load i32, i32* %9
%43 = add nsw i32 %42, 1
%44 = sext i32 %43 to i64
%45 = getelementptr inbounds float, float* %41, i64 %44
%46 = load float, float* %45
%47 = fmul float %40, %46
%48 = fpext float %47 to double
%49 = call double @llvm.fmuladd.f64(double %33, double %39, double %48)
%50 = fptrunc double %49 to float
%51 = load float*, float** %5
%52 = load i32, i32* %9
%53 = sub nsw i32 %52, 1
%54 = load i32, i32* %8
%55 = mul nsw i32 %53, %54
%56 = load i32, i32* %10
%57 = add nsw i32 %55, %56
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds float, float* %51, i64 %58
store float %50, float* %59
br label %60
60: ; preds = %22
%61 = load i32, i32* %10
%62 = add nsw i32 %61, 1
store i32 %62, i32* %10
br label %18
63: ; preds = %18
br label %64
64: ; preds = %63
%65 = load i32, i32* %9
%66 = add nsw i32 %65, 1
store i32 %66, i32* %9
br label %12
67: ; preds = %12
ret void
}

LLVM IR (-O2):
define void @interp(float* nocapture writeonly %0, float* nocapture %1, i32 %2, i32 %3) {
%5 = icmp sgt i32 %2, 1
br i1 %5, label %6, label %94
6: ; preds = %4
%7 = icmp slt i32 %3, 1
%8 = sitofp i32 %3 to double
%9 = add i32 %3, 1
%10 = sext i32 %3 to i64
%11 = zext i32 %2 to i64
%12 = zext i32 %9 to i64
%13 = getelementptr float, float* %0, i64 1
%14 = getelementptr float, float* %0, i64 %12
%15 = getelementptr float, float* %1, i64 2
%16 = add nsw i64 %12, -1
%17 = icmp ult i64 %16, 4
%18 = and i64 %16, -4
%19 = or i64 %18, 1
%20 = insertelement <4 x double> poison, double %8, i64 0
%21 = shufflevector <4 x double> %20, <4 x double> poison, <4 x i32> zeroinitializer
%22 = icmp eq i64 %16, %18
br label %23
23: ; preds = %6, %91
%24 = phi i64 [ 1, %6 ], [ %92, %91 ]
br i1 %7, label %25, label %27
25: ; preds = %23
%26 = add nuw nsw i64 %24, 1
br label %91
27: ; preds = %23
%28 = getelementptr inbounds float, float* %1, i64 %24
%29 = add nuw nsw i64 %24, 1
%30 = getelementptr inbounds float, float* %1, i64 %29
%31 = add nsw i64 %24, -1
%32 = mul nsw i64 %31, %10
br i1 %17, label %69, label %33
33: ; preds = %27
%34 = getelementptr float, float* %13, i64 %32
%35 = getelementptr float, float* %14, i64 %32
%36 = getelementptr float, float* %15, i64 %24
%37 = icmp ult float* %34, %36
%38 = icmp ult float* %28, %35
%39 = and i1 %37, %38
br i1 %39, label %69, label %40
40: ; preds = %33
%41 = load float, float* %28
%42 = insertelement <4 x float> poison, float %41, i64 0
%43 = shufflevector <4 x float> %42, <4 x float> poison, <4 x i32> zeroinitializer
%44 = fpext <4 x float> %43 to <4 x double>
%45 = load float, float* %30
%46 = insertelement <4 x float> poison, float %45, i64 0
%47 = shufflevector <4 x float> %46, <4 x float> poison, <4 x i32> zeroinitializer
br label %48
48: ; preds = %48, %40
%49 = phi i64 [ 0, %40 ], [ %65, %48 ]
%50 = phi <4 x i32> [ <i32 1, i32 2, i32 3, i32 4>, %40 ], [ %66, %48 ]
%51 = or i64 %49, 1
%52 = add <4 x i32> %50, <i32 -1, i32 -1, i32 -1, i32 -1>
%53 = sitofp <4 x i32> %52 to <4 x double>
%54 = fdiv <4 x double> %53, %21
%55 = fptrunc <4 x double> %54 to <4 x float>
%56 = fpext <4 x float> %55 to <4 x double>
%57 = fsub <4 x double> <double 1.000000e+00, double 1.000000e+00, double 1.000000e+00, double 1.000000e+00>, %56
%58 = fmul <4 x float> %47, %55
%59 = fpext <4 x float> %58 to <4 x double>
%60 = call <4 x double> @llvm.fmuladd.v4f64(<4 x double> %57, <4 x double> %44, <4 x double> %59)
%61 = fptrunc <4 x double> %60 to <4 x float>
%62 = add nsw i64 %51, %32
%63 = getelementptr inbounds float, float* %0, i64 %62
%64 = bitcast float* %63 to <4 x float>*
store <4 x float> %61, <4 x float>* %64
%65 = add nuw i64 %49, 4
%66 = add <4 x i32> %50, <i32 4, i32 4, i32 4, i32 4>
%67 = icmp eq i64 %65, %18
br i1 %67, label %68, label %48
68: ; preds = %48
br i1 %22, label %91, label %69
69: ; preds = %33, %27, %68
%70 = phi i64 [ 1, %33 ], [ 1, %27 ], [ %19, %68 ]
br label %71
71: ; preds = %69, %71
%72 = phi i64 [ %89, %71 ], [ %70, %69 ]
%73 = trunc i64 %72 to i32
%74 = add i32 %73, -1
%75 = sitofp i32 %74 to double
%76 = fdiv double %75, %8
%77 = fptrunc double %76 to float
%78 = fpext float %77 to double
%79 = fsub double 1.000000e+00, %78
%80 = load float, float* %28
%81 = fpext float %80 to double
%82 = load float, float* %30
%83 = fmul float %82, %77
%84 = fpext float %83 to double
%85 = tail call double @llvm.fmuladd.f64(double %79, double %81, double %84)
%86 = fptrunc double %85 to float
%87 = add nsw i64 %72, %32
%88 = getelementptr inbounds float, float* %0, i64 %87
store float %86, float* %88
%89 = add nuw nsw i64 %72, 1
%90 = icmp eq i64 %89, %12
br i1 %90, label %91, label %71
91: ; preds = %71, %68, %25
%92 = phi i64 [ %26, %25 ], [ %29, %68 ], [ %29, %71 ]
%93 = icmp eq i64 %92, %11
br i1 %93, label %94, label %23
94: ; preds = %91, %4
ret void
}