void difdiv(int n, double *x, double *f) {
  int k = 0;
  int i = 0;
  for (k = 1; k < n; k++) {
    for (i = n - 1; i >= k; i--) {
      f[i] = (f[i] - f[i - 1]) / (x[i] - x[i - k]);
    }
  }
}


x86 assembly (clang, -O0):
<difdiv>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11b5 <difdiv+0xb5>
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jl     11a2 <difdiv+0xa2>
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm1
divsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x20(%rbp)
jmp    1139 <difdiv+0x39>
jmp    11a7 <difdiv+0xa7>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1124 <difdiv+0x24>
pop    %rbp
ret

x86 assembly (clang, -O2):
<difdiv>:
cmp    $0x2,%edi
jl     116b <difdiv+0x6b>
mov    %edi,%r8d
lea    -0x10(%rsi),%r9
mov    $0x1,%edi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %r8,%rax
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  -0x8(%rdx,%rax,8),%xmm0
subsd  -0x10(%rdx,%rax,8),%xmm0
movsd  -0x8(%rsi,%rax,8),%xmm1
subsd  (%r9,%rax,8),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rdx,%rax,8)
lea    -0x1(%rax),%rcx
mov    %rcx,%rax
cmp    %rdi,%rcx
jg     1130 <difdiv+0x30>
add    $0x1,%rdi
add    $0xfffffffffffffff8,%r9
cmp    %r8,%rdi
jne    1120 <difdiv+0x20>
ret

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

LLVM IR (-O2):
define void @difdiv(i32 %0, double* nocapture %1, double* nocapture %2) {
%4 = icmp sgt i32 %0, 1
br i1 %4, label %5, label %30
5: ; preds = %3
%6 = zext i32 %0 to i64
%7 = zext i32 %0 to i64
br label %8
8: ; preds = %27, %5
%9 = phi i64 [ 1, %5 ], [ %28, %27 ]
br label %10
10: ; preds = %8, %10
%11 = phi i64 [ %6, %8 ], [ %12, %10 ]
%12 = add nsw i64 %11, -1
%13 = getelementptr inbounds double, double* %2, i64 %12
%14 = load double, double* %13
%15 = add nsw i64 %11, -2
%16 = getelementptr inbounds double, double* %2, i64 %15
%17 = load double, double* %16
%18 = fsub double %14, %17
%19 = getelementptr inbounds double, double* %1, i64 %12
%20 = load double, double* %19
%21 = sub nsw i64 %12, %9
%22 = getelementptr inbounds double, double* %1, i64 %21
%23 = load double, double* %22
%24 = fsub double %20, %23
%25 = fdiv double %18, %24
store double %25, double* %13
%26 = icmp sgt i64 %12, %9
br i1 %26, label %10, label %27
27: ; preds = %10
%28 = add nuw nsw i64 %9, 1
%29 = icmp eq i64 %28, %7
br i1 %29, label %30, label %8
30: ; preds = %27, %3
ret void
}