void init_splinen(float *a, int n) {
  a[0] = 1;
  for (int k = 2; k <= n; k++)
    a[0] /= k;
  for (int k = 1; k <= n + 1; k++)
    a[k] = -a[k - 1] * (n + 2 - k) / k;
}


x86 assembly (clang, -O0):
<init_splinen>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    -0x8(%rbp),%rax
movss  0xee9(%rip),%xmm0
movss  %xmm0,(%rax)
movl   $0x2,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     1151 <init_splinen+0x51>
cvtsi2ssl -0x10(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,(%rax)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1122 <init_splinen+0x22>
movl   $0x1,-0x14(%rbp)
mov    -0x14(%rbp),%eax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     11bd <init_splinen+0xbd>
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mov    -0xc(%rbp),%eax
add    $0x2,%eax
sub    -0x14(%rbp),%eax
cvtsi2ss %eax,%xmm1
mulss  %xmm1,%xmm0
cvtsi2ssl -0x14(%rbp),%xmm1
divss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1158 <init_splinen+0x58>
pop    %rbp
ret

x86 assembly (clang, -O2):
<init_splinen>:
movl   $0x3f800000,(%rdi)
cmp    $0x2,%esi
jge    111c <init_splinen+0x1c>
movss  0xeed(%rip),%xmm0
test   %esi,%esi
jns    11c0 <init_splinen+0xc0>
ret
lea    -0x1(%rsi),%r8d
lea    -0x2(%rsi),%ecx
mov    %r8d,%r9d
and    $0x3,%r9d
cmp    $0x3,%ecx
jae    113e <init_splinen+0x3e>
mov    $0x2,%ecx
movss  0xec4(%rip),%xmm0
jmp    118d <init_splinen+0x8d>
and    $0xfffffffc,%r8d
mov    $0x2,%ecx
movss  0xeb1(%rip),%xmm0
nop
xorps  %xmm1,%xmm1
cvtsi2ss %ecx,%xmm1
divss  %xmm1,%xmm0
lea    0x1(%rcx),%edx
xorps  %xmm1,%xmm1
cvtsi2ss %edx,%xmm1
divss  %xmm1,%xmm0
lea    0x2(%rcx),%edx
xorps  %xmm1,%xmm1
cvtsi2ss %edx,%xmm1
lea    0x3(%rcx),%eax
xorps  %xmm2,%xmm2
cvtsi2ss %eax,%xmm2
divss  %xmm1,%xmm0
divss  %xmm2,%xmm0
add    $0x4,%ecx
cmp    %r8d,%edx
jne    1150 <init_splinen+0x50>
test   %r9d,%r9d
je     11b4 <init_splinen+0xb4>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
xorps  %xmm1,%xmm1
cvtsi2ss %ecx,%xmm1
divss  %xmm1,%xmm0
add    $0x1,%ecx
add    $0xffffffff,%r9d
jne    11a0 <init_splinen+0xa0>
movss  %xmm0,(%rdi)
test   %esi,%esi
js     111b <init_splinen+0x1b>
lea    0x1(%rsi),%r8d
test   %esi,%esi
je     123c <init_splinen+0x13c>
mov    %r8d,%r9d
and    $0xfffffffe,%r9d
mov    $0x1,%eax
movaps 0xe35(%rip),%xmm1
mov    %esi,%ecx
nopl   (%rax)
xorps  %xmm1,%xmm0
lea    0x1(%rcx),%edx
xorps  %xmm2,%xmm2
cvtsi2ss %edx,%xmm2
mulss  %xmm0,%xmm2
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
divss  %xmm0,%xmm2
xorps  %xmm0,%xmm0
cvtsi2ss %ecx,%xmm0
movss  %xmm2,(%rdi,%rax,4)
xorps  %xmm1,%xmm2
mulss  %xmm2,%xmm0
lea    0x1(%rax),%rdx
xorps  %xmm2,%xmm2
cvtsi2ss %edx,%xmm2
divss  %xmm2,%xmm0
movss  %xmm0,0x4(%rdi,%rax,4)
add    $0x2,%rax
add    $0xfffffffe,%ecx
cmp    %r9,%rdx
jne    11e0 <init_splinen+0xe0>
test   $0x1,%r8b
je     111b <init_splinen+0x1b>
jmp    124b <init_splinen+0x14b>
mov    $0x1,%eax
test   $0x1,%r8b
je     111b <init_splinen+0x1b>
add    $0x2,%esi
xorps  0xdbb(%rip),%xmm0
sub    %eax,%esi
xorps  %xmm1,%xmm1
cvtsi2ss %esi,%xmm1
xorps  %xmm2,%xmm2
cvtsi2ss %eax,%xmm2
mulss  %xmm0,%xmm1
divss  %xmm2,%xmm1
movss  %xmm1,(%rdi,%rax,4)
ret

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

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