divdif(vec, k, diffn) double vec[];
double *diffn;
int k;
{
  double diftbl[19];
  double *p, *q;
  double y;
  int i, o;
  p = diftbl;
  q = vec;
  for (i = 0; i <= k; i++)
    *p++ = *q++;
  o = k;
  for (o = k; o > 0; o--) {
    p = diftbl;
    q = p;
    for (i = 0; i < o; i++) {
      y = *p++;
      *q++ = *p - y;
    }
    *diffn++ = *p;
  }
  *diffn++ = *(q - 1);
}


x86 assembly (clang, -O0):
<divdif>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %rdx,-0x20(%rbp)
lea    -0xc0(%rbp),%rax
mov    %rax,-0xc8(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0xd0(%rbp)
movl   $0x0,-0xdc(%rbp)
mov    -0xdc(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     118b <divdif+0x8b>
mov    -0xd0(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0xd0(%rbp)
movsd  (%rax),%xmm0
mov    -0xc8(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0xc8(%rbp)
movsd  %xmm0,(%rax)
mov    -0xdc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xdc(%rbp)
jmp    1136 <divdif+0x36>
mov    -0x14(%rbp),%eax
mov    %eax,-0xe0(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0xe0(%rbp)
cmpl   $0x0,-0xe0(%rbp)
jle    1275 <divdif+0x175>
lea    -0xc0(%rbp),%rax
mov    %rax,-0xc8(%rbp)
mov    -0xc8(%rbp),%rax
mov    %rax,-0xd0(%rbp)
movl   $0x0,-0xdc(%rbp)
mov    -0xdc(%rbp),%eax
cmp    -0xe0(%rbp),%eax
jge    1243 <divdif+0x143>
mov    -0xc8(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0xc8(%rbp)
movsd  (%rax),%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
mov    -0xc8(%rbp),%rax
movsd  (%rax),%xmm0
subsd  -0xd8(%rbp),%xmm0
ff
mov    -0xd0(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0xd0(%rbp)
movsd  %xmm0,(%rax)
mov    -0xdc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xdc(%rbp)
jmp    11d0 <divdif+0xd0>
mov    -0xc8(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x20(%rbp)
movsd  %xmm0,(%rax)
mov    -0xe0(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0xe0(%rbp)
jmp    119d <divdif+0x9d>
mov    -0xd0(%rbp),%rax
movsd  -0x8(%rax),%xmm0
mov    -0x20(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x20(%rbp)
movsd  %xmm0,(%rax)
mov    -0x4(%rbp),%eax
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<divdif>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xa0,%rsp
mov    %rdx,%r14
mov    %rdi,%rbx
test   %esi,%esi
js     12f1 <divdif+0x1e1>
mov    %esi,%r12d
mov    %esi,%r13d
lea    0x8(,%r13,8),%rdx
mov    %rsp,%r15
mov    %r15,%rdi
mov    %rbx,%rsi
call   1030 <memcpy@plt>
test   %r12d,%r12d
je     12e9 <divdif+0x1d9>
lea    -0x1(%r12),%eax
lea    0x1(%rax),%r8
add    $0xfffffffffffffffd,%rax
xor    %edx,%edx
jmp    1193 <divdif+0x83>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  %xmm0,(%r14)
add    $0x8,%r14
lea    -0x1(%r12),%ecx
add    $0x1,%rdx
add    $0xffffffffffffffff,%rax
cmp    $0x1,%r12d
mov    %ecx,%r12d
jle    12f1 <divdif+0x1e1>
mov    %r8,%rcx
sub    %rdx,%rcx
movsd  (%rsp),%xmm1
cmp    $0x4,%rcx
jae    11b0 <divdif+0xa0>
xor    %esi,%esi
mov    %r15,%rbx
jmp    12b4 <divdif+0x1a4>
xchg   %ax,%ax
lea    -0x4(%rcx),%rsi
mov    %rsi,%rdi
shr    $0x2,%rdi
add    $0x1,%rdi
unpcklpd %xmm1,%xmm0
cmp    $0x4,%rsi
jae    11d0 <divdif+0xc0>
xor    %ebx,%ebx
jmp    1261 <divdif+0x151>
mov    %rax,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
and    $0xfffffffffffffffe,%rsi
xor    %ebx,%ebx
movapd %xmm0,%xmm1
cs nopw 0x0(%rax,%rax,1)
nop
movupd 0x8(%rsp,%rbx,8),%xmm2
movupd 0x18(%rsp,%rbx,8),%xmm3
movupd 0x28(%rsp,%rbx,8),%xmm4
movupd 0x38(%rsp,%rbx,8),%xmm0
shufpd $0x1,%xmm2,%xmm1
movapd %xmm2,%xmm5
shufpd $0x1,%xmm3,%xmm5
subpd  %xmm1,%xmm2
movapd %xmm3,%xmm1
subpd  %xmm5,%xmm1
movapd %xmm2,(%rsp,%rbx,8)
movapd %xmm1,0x10(%rsp,%rbx,8)
shufpd $0x1,%xmm4,%xmm3
movapd %xmm4,%xmm1
shufpd $0x1,%xmm0,%xmm1
subpd  %xmm3,%xmm4
movapd %xmm0,%xmm2
subpd  %xmm1,%xmm2
movapd %xmm4,0x20(%rsp,%rbx,8)
movapd %xmm2,0x30(%rsp,%rbx,8)
add    $0x8,%rbx
movapd %xmm0,%xmm1
add    $0xfffffffffffffffe,%rsi
jne    11f0 <divdif+0xe0>
mov    %rcx,%rsi
and    $0xfffffffffffffffc,%rsi
test   $0x1,%dil
je     129f <divdif+0x18f>
movupd 0x8(%rsp,%rbx,8),%xmm1
movupd 0x18(%rsp,%rbx,8),%xmm2
shufpd $0x1,%xmm1,%xmm0
movapd %xmm1,%xmm3
shufpd $0x1,%xmm2,%xmm3
subpd  %xmm0,%xmm1
movapd %xmm2,%xmm0
subpd  %xmm3,%xmm2
movapd %xmm1,(%rsp,%rbx,8)
movapd %xmm2,0x10(%rsp,%rbx,8)
lea    (%rsp,%rsi,8),%rbx
unpckhpd %xmm0,%xmm0
cmp    %rsi,%rcx
je     1170 <divdif+0x60>
movapd %xmm0,%xmm1
mov    %rbx,%rcx
nopw   0x0(%rax,%rax,1)
lea    0x8(%rcx),%rbx
movsd  0x8(%rcx),%xmm0
movapd %xmm0,%xmm2
subsd  %xmm1,%xmm2
movsd  %xmm2,(%rcx)
add    $0x1,%esi
movapd %xmm0,%xmm1
mov    %rbx,%rcx
cmp    %esi,%r12d
jne    12c0 <divdif+0x1b0>
jmp    1170 <divdif+0x60>
lea    (%rbx,%r13,8),%rbx
add    $0x8,%rbx
movsd  -0x8(%rbx),%xmm0
movsd  %xmm0,(%r14)
add    $0xa0,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define i32 @divdif(double* %0, i32 %1, double* %2) {
%4 = alloca i32
%5 = alloca double*
%6 = alloca i32
%7 = alloca double*
%8 = alloca [19 x double]
%9 = alloca double*
%10 = alloca double*
%11 = alloca double
%12 = alloca i32
%13 = alloca i32
store double* %0, double** %5
store i32 %1, i32* %6
store double* %2, double** %7
%14 = getelementptr inbounds [19 x double], [19 x double]* %8, i64 0, i64 0
store double* %14, double** %9
%15 = load double*, double** %5
store double* %15, double** %10
store i32 0, i32* %12
br label %16
16: ; preds = %26, %3
%17 = load i32, i32* %12
%18 = load i32, i32* %6
%19 = icmp sle i32 %17, %18
br i1 %19, label %20, label %29
20: ; preds = %16
%21 = load double*, double** %10
%22 = getelementptr inbounds double, double* %21, i32 1
store double* %22, double** %10
%23 = load double, double* %21
%24 = load double*, double** %9
%25 = getelementptr inbounds double, double* %24, i32 1
store double* %25, double** %9
store double %23, double* %24
br label %26
26: ; preds = %20
%27 = load i32, i32* %12
%28 = add nsw i32 %27, 1
store i32 %28, i32* %12
br label %16
29: ; preds = %16
%30 = load i32, i32* %6
store i32 %30, i32* %13
%31 = load i32, i32* %6
store i32 %31, i32* %13
br label %32
32: ; preds = %60, %29
%33 = load i32, i32* %13
%34 = icmp sgt i32 %33, 0
br i1 %34, label %35, label %63
35: ; preds = %32
%36 = getelementptr inbounds [19 x double], [19 x double]* %8, i64 0, i64 0
store double* %36, double** %9
%37 = load double*, double** %9
store double* %37, double** %10
store i32 0, i32* %12
br label %38
38: ; preds = %52, %35
%39 = load i32, i32* %12
%40 = load i32, i32* %13
%41 = icmp slt i32 %39, %40
br i1 %41, label %42, label %55
42: ; preds = %38
%43 = load double*, double** %9
%44 = getelementptr inbounds double, double* %43, i32 1
store double* %44, double** %9
%45 = load double, double* %43
store double %45, double* %11
%46 = load double*, double** %9
%47 = load double, double* %46
%48 = load double, double* %11
%49 = fsub double %47, %48
%50 = load double*, double** %10
%51 = getelementptr inbounds double, double* %50, i32 1
store double* %51, double** %10
store double %49, double* %50
br label %52
52: ; preds = %42
%53 = load i32, i32* %12
%54 = add nsw i32 %53, 1
store i32 %54, i32* %12
br label %38
55: ; preds = %38
%56 = load double*, double** %9
%57 = load double, double* %56
%58 = load double*, double** %7
%59 = getelementptr inbounds double, double* %58, i32 1
store double* %59, double** %7
store double %57, double* %58
br label %60
60: ; preds = %55
%61 = load i32, i32* %13
%62 = add nsw i32 %61, -1
store i32 %62, i32* %13
br label %32
63: ; preds = %32
%64 = load double*, double** %10
%65 = getelementptr inbounds double, double* %64, i64 -1
%66 = load double, double* %65
%67 = load double*, double** %7
%68 = getelementptr inbounds double, double* %67, i32 1
store double* %68, double** %7
store double %66, double* %67
%69 = load i32, i32* %4
ret i32 %69
}

LLVM IR (-O2):
define i32 @divdif(double* nocapture %0, i32 %1, double* nocapture writeonly %2) {
%4 = alloca [19 x double]
%5 = bitcast [19 x double]* %4 to i8*
%6 = getelementptr inbounds [19 x double], [19 x double]* %4, i64 0, i64 0
%7 = icmp slt i32 %1, 0
br i1 %7, label %121, label %8
8: ; preds = %3
%9 = bitcast double* %0 to i8*
%10 = zext i32 %1 to i64
%11 = shl nuw nsw i64 %10, 3
%12 = add nuw nsw i64 %11, 8
%13 = add nuw nsw i64 %10, 1
%14 = getelementptr double, double* %0, i64 %13
%15 = icmp eq i32 %1, 0
br i1 %15, label %121, label %16
16: ; preds = %8
%17 = add i32 %1, -1
%18 = zext i32 %17 to i64
%19 = add nuw nsw i64 %18, 1
br label %20
20: ; preds = %16, %114
%21 = phi i64 [ 0, %16 ], [ %120, %114 ]
%22 = phi i32 [ %1, %16 ], [ %118, %114 ]
%23 = phi double* [ %2, %16 ], [ %117, %114 ]
%24 = sub i64 %19, %21
%25 = add i64 %24, -4
%26 = lshr i64 %25, 2
%27 = add nuw nsw i64 %26, 1
%28 = sub i64 %19, %21
%29 = load double, double* %6
%30 = icmp ult i64 %28, 4
br i1 %30, label %101, label %31
31: ; preds = %20
%32 = and i64 %28, -4
%33 = trunc i64 %32 to i32
%34 = getelementptr [19 x double], [19 x double]* %4, i64 0, i64 %32
%35 = insertelement <2 x double> poison, double %29, i64 1
%36 = and i64 %27, 1
%37 = icmp ult i64 %25, 4
br i1 %37, label %76, label %38
38: ; preds = %31
%39 = and i64 %27, 9223372036854775806
br label %40
40: ; preds = %40, %38
%41 = phi i64 [ 0, %38 ], [ %73, %40 ]
%42 = phi <2 x double> [ %35, %38 ], [ %65, %40 ]
%43 = phi i64 [ 0, %38 ], [ %74, %40 ]
%44 = getelementptr [19 x double], [19 x double]* %4, i64 0, i64 %41
%45 = getelementptr inbounds double, double* %44, i64 1
%46 = bitcast double* %45 to <2 x double>*
%47 = load <2 x double>, <2 x double>* %46
%48 = getelementptr inbounds double, double* %44, i64 3
%49 = bitcast double* %48 to <2 x double>*
%50 = load <2 x double>, <2 x double>* %49
%51 = shufflevector <2 x double> %42, <2 x double> %47, <2 x i32> <i32 1, i32 2>
%52 = shufflevector <2 x double> %47, <2 x double> %50, <2 x i32> <i32 1, i32 2>
%53 = fsub <2 x double> %47, %51
%54 = fsub <2 x double> %50, %52
%55 = bitcast double* %44 to <2 x double>*
store <2 x double> %53, <2 x double>* %55
%56 = getelementptr double, double* %44, i64 2
%57 = bitcast double* %56 to <2 x double>*
store <2 x double> %54, <2 x double>* %57
%58 = or i64 %41, 4
%59 = getelementptr [19 x double], [19 x double]* %4, i64 0, i64 %58
%60 = getelementptr inbounds double, double* %59, i64 1
%61 = bitcast double* %60 to <2 x double>*
%62 = load <2 x double>, <2 x double>* %61
%63 = getelementptr inbounds double, double* %59, i64 3
%64 = bitcast double* %63 to <2 x double>*
%65 = load <2 x double>, <2 x double>* %64
%66 = shufflevector <2 x double> %50, <2 x double> %62, <2 x i32> <i32 1, i32 2>
%67 = shufflevector <2 x double> %62, <2 x double> %65, <2 x i32> <i32 1, i32 2>
%68 = fsub <2 x double> %62, %66
%69 = fsub <2 x double> %65, %67
%70 = bitcast double* %59 to <2 x double>*
store <2 x double> %68, <2 x double>* %70
%71 = getelementptr double, double* %59, i64 2
%72 = bitcast double* %71 to <2 x double>*
store <2 x double> %69, <2 x double>* %72
%73 = add nuw i64 %41, 8
%74 = add i64 %43, 2
%75 = icmp eq i64 %74, %39
br i1 %75, label %76, label %40
76: ; preds = %40, %31
%77 = phi <2 x double> [ undef, %31 ], [ %65, %40 ]
%78 = phi i64 [ 0, %31 ], [ %73, %40 ]
%79 = phi <2 x double> [ %35, %31 ], [ %65, %40 ]
%80 = icmp eq i64 %36, 0
br i1 %80, label %96, label %81
81: ; preds = %76
%82 = getelementptr [19 x double], [19 x double]* %4, i64 0, i64 %78
%83 = getelementptr inbounds double, double* %82, i64 1
%84 = bitcast double* %83 to <2 x double>*
%85 = load <2 x double>, <2 x double>* %84
%86 = getelementptr inbounds double, double* %82, i64 3
%87 = bitcast double* %86 to <2 x double>*
%88 = load <2 x double>, <2 x double>* %87
%89 = shufflevector <2 x double> %79, <2 x double> %85, <2 x i32> <i32 1, i32 2>
%90 = shufflevector <2 x double> %85, <2 x double> %88, <2 x i32> <i32 1, i32 2>
%91 = fsub <2 x double> %85, %89
%92 = fsub <2 x double> %88, %90
%93 = bitcast double* %82 to <2 x double>*
store <2 x double> %91, <2 x double>* %93
%94 = getelementptr double, double* %82, i64 2
%95 = bitcast double* %94 to <2 x double>*
store <2 x double> %92, <2 x double>* %95
br label %96
96: ; preds = %76, %81
%97 = phi <2 x double> [ %77, %76 ], [ %88, %81 ]
%98 = icmp eq i64 %28, %32
%99 = extractelement <2 x double> %97, i64 1
%100 = extractelement <2 x double> %97, i64 1
br i1 %98, label %114, label %101
101: ; preds = %20, %96
%102 = phi double [ %99, %96 ], [ %29, %20 ]
%103 = phi i32 [ %33, %96 ], [ 0, %20 ]
%104 = phi double* [ %34, %96 ], [ %6, %20 ]
br label %105
105: ; preds = %101, %105
%106 = phi double [ %110, %105 ], [ %102, %101 ]
%107 = phi i32 [ %112, %105 ], [ %103, %101 ]
%108 = phi double* [ %109, %105 ], [ %104, %101 ]
%109 = getelementptr inbounds double, double* %108, i64 1
%110 = load double, double* %109
%111 = fsub double %110, %106
store double %111, double* %108
%112 = add nuw nsw i32 %107, 1
%113 = icmp eq i32 %112, %22
br i1 %113, label %114, label %105
114: ; preds = %105, %96
%115 = phi double* [ %34, %96 ], [ %109, %105 ]
%116 = phi double [ %100, %96 ], [ %110, %105 ]
%117 = getelementptr inbounds double, double* %23, i64 1
store double %116, double* %23
%118 = add nsw i32 %22, -1
%119 = icmp sgt i32 %22, 1
%120 = add i64 %21, 1
br i1 %119, label %20, label %121
121: ; preds = %114, %3, %8
%122 = phi double* [ %14, %8 ], [ %0, %3 ], [ %115, %114 ]
%123 = phi double* [ %2, %8 ], [ %2, %3 ], [ %117, %114 ]
%124 = getelementptr inbounds double, double* %122, i64 -1
%125 = load double, double* %124
store double %125, double* %123
ret i32 undef
}