dupdate(diffn, k, f) register double *diffn;
int k;
double f;
{
  double new, old;
  int i;
  new = f;
  for (i = 0; i < k; i++) {
    old = *diffn;
    *diffn++ = new;
    new = new - old;
  }
  *diffn++ = new;
}


x86 assembly (clang, -O0):
<dupdate>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    116f <dupdate+0x6f>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x28(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x10(%rbp)
movsd  %xmm0,(%rax)
movsd  -0x28(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1121 <dupdate+0x21>
movsd  -0x28(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x10(%rbp)
movsd  %xmm0,(%rax)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<dupdate>:
test   %esi,%esi
jle    1113 <dupdate+0x13>
lea    -0x1(%rsi),%ecx
mov    %esi,%eax
and    $0x3,%eax
cmp    $0x3,%ecx
jae    1118 <dupdate+0x18>
jmp    1162 <dupdate+0x62>
movsd  %xmm0,(%rdi)
ret
and    $0xfffffffc,%esi
nopl   0x0(%rax,%rax,1)
movapd %xmm0,%xmm1
subsd  (%rdi),%xmm1
movsd  %xmm0,(%rdi)
movapd %xmm1,%xmm0
subsd  0x8(%rdi),%xmm0
movsd  %xmm1,0x8(%rdi)
movapd %xmm0,%xmm1
subsd  0x10(%rdi),%xmm1
movsd  %xmm0,0x10(%rdi)
movapd %xmm1,%xmm0
subsd  0x18(%rdi),%xmm0
movsd  %xmm1,0x18(%rdi)
add    $0x20,%rdi
add    $0xfffffffc,%esi
jne    1120 <dupdate+0x20>
mov    %rdi,%rcx
movapd %xmm0,%xmm1
test   %eax,%eax
je     1188 <dupdate+0x88>
movapd %xmm0,%xmm1
xchg   %ax,%ax
lea    0x8(%rdi),%rcx
subsd  (%rdi),%xmm1
movsd  %xmm0,(%rdi)
mov    %rcx,%rdi
movapd %xmm1,%xmm0
add    $0xffffffff,%eax
jne    1170 <dupdate+0x70>
movsd  %xmm1,(%rcx)
ret

LLVM IR (-O0):
define i32 @dupdate(double* %0, i32 %1, double %2) {
%4 = alloca i32
%5 = alloca double*
%6 = alloca i32
%7 = alloca double
%8 = alloca double
%9 = alloca double
%10 = alloca i32
store double* %0, double** %5
store i32 %1, i32* %6
store double %2, double* %7
%11 = load double, double* %7
store double %11, double* %8
store i32 0, i32* %10
br label %12
12: ; preds = %25, %3
%13 = load i32, i32* %10
%14 = load i32, i32* %6
%15 = icmp slt i32 %13, %14
br i1 %15, label %16, label %28
16: ; preds = %12
%17 = load double*, double** %5
%18 = load double, double* %17
store double %18, double* %9
%19 = load double, double* %8
%20 = load double*, double** %5
%21 = getelementptr inbounds double, double* %20, i32 1
store double* %21, double** %5
store double %19, double* %20
%22 = load double, double* %8
%23 = load double, double* %9
%24 = fsub double %22, %23
store double %24, double* %8
br label %25
25: ; preds = %16
%26 = load i32, i32* %10
%27 = add nsw i32 %26, 1
store i32 %27, i32* %10
br label %12
28: ; preds = %12
%29 = load double, double* %8
%30 = load double*, double** %5
%31 = getelementptr inbounds double, double* %30, i32 1
store double* %31, double** %5
store double %29, double* %30
%32 = load i32, i32* %4
ret i32 %32
}

LLVM IR (-O2):
define i32 @dupdate(double* nocapture %0, i32 %1, double %2) {
%4 = icmp sgt i32 %1, 0
br i1 %4, label %5, label %44
5: ; preds = %3
%6 = add i32 %1, -1
%7 = and i32 %1, 3
%8 = icmp ult i32 %6, 3
br i1 %8, label %29, label %9
9: ; preds = %5
%10 = and i32 %1, -4
br label %11
11: ; preds = %11, %9
%12 = phi double* [ %0, %9 ], [ %25, %11 ]
%13 = phi double [ %2, %9 ], [ %26, %11 ]
%14 = phi i32 [ 0, %9 ], [ %27, %11 ]
%15 = load double, double* %12
%16 = getelementptr inbounds double, double* %12, i64 1
store double %13, double* %12
%17 = fsub double %13, %15
%18 = load double, double* %16
%19 = getelementptr inbounds double, double* %12, i64 2
store double %17, double* %16
%20 = fsub double %17, %18
%21 = load double, double* %19
%22 = getelementptr inbounds double, double* %12, i64 3
store double %20, double* %19
%23 = fsub double %20, %21
%24 = load double, double* %22
%25 = getelementptr inbounds double, double* %12, i64 4
store double %23, double* %22
%26 = fsub double %23, %24
%27 = add i32 %14, 4
%28 = icmp eq i32 %27, %10
br i1 %28, label %29, label %11
29: ; preds = %11, %5
%30 = phi double* [ undef, %5 ], [ %25, %11 ]
%31 = phi double [ undef, %5 ], [ %26, %11 ]
%32 = phi double* [ %0, %5 ], [ %25, %11 ]
%33 = phi double [ %2, %5 ], [ %26, %11 ]
%34 = icmp eq i32 %7, 0
br i1 %34, label %44, label %35
35: ; preds = %29, %35
%36 = phi double* [ %40, %35 ], [ %32, %29 ]
%37 = phi double [ %41, %35 ], [ %33, %29 ]
%38 = phi i32 [ %42, %35 ], [ 0, %29 ]
%39 = load double, double* %36
%40 = getelementptr inbounds double, double* %36, i64 1
store double %37, double* %36
%41 = fsub double %37, %39
%42 = add i32 %38, 1
%43 = icmp eq i32 %42, %7
br i1 %43, label %44, label %35
44: ; preds = %29, %35, %3
%45 = phi double [ %2, %3 ], [ %31, %29 ], [ %41, %35 ]
%46 = phi double* [ %0, %3 ], [ %30, %29 ], [ %40, %35 ]
store double %45, double* %46
ret i32 undef
}