void kcof2acof(int order, float *kcof, float *acof) {
  float atmp[20 + 2];
  int i, p;
  acof[0] = 1.0;
  atmp[0] = 1.0;
  for (p = 1; p <= order; p++) {
    acof[p] = 0.0;
    for (i = 1; i < p; i++) {
      atmp[i] = acof[i] + kcof[p - 1] * acof[p - i];
    }
    for (i = 1; i < p; i++) {
      acof[i] = atmp[i];
    }
    acof[p] = kcof[p - 1];
  }
}


x86 assembly (clang, -O0):
<kcof2acof>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x18(%rbp),%rax
movss  0xee5(%rip),%xmm0
movss  %xmm0,(%rax)
movss  0xed9(%rip),%xmm0
movss  %xmm0,-0x70(%rbp)
movl   $0x1,-0x78(%rbp)
mov    -0x78(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     1218 <kcof2acof+0x118>
mov    -0x18(%rbp),%rax
movslq -0x78(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movl   $0x1,-0x74(%rbp)
mov    -0x74(%rbp),%eax
cmp    -0x78(%rbp),%eax
jge    11b3 <kcof2acof+0xb3>
mov    -0x18(%rbp),%rax
movslq -0x74(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x78(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x78(%rbp),%ecx
sub    -0x74(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movslq -0x74(%rbp),%rax
movss  %xmm0,-0x70(%rbp,%rax,4)
mov    -0x74(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x74(%rbp)
jmp    1156 <kcof2acof+0x56>
movl   $0x1,-0x74(%rbp)
mov    -0x74(%rbp),%eax
cmp    -0x78(%rbp),%eax
jge    11eb <kcof2acof+0xeb>
movslq -0x74(%rbp),%rax
movss  -0x70(%rbp,%rax,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x74(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x74(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x74(%rbp)
jmp    11ba <kcof2acof+0xba>
mov    -0x10(%rbp),%rax
mov    -0x78(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x78(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x78(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x78(%rbp)
jmp    1133 <kcof2acof+0x33>
pop    %rbp
ret

x86 assembly (clang, -O2):
<kcof2acof>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
movl   $0x3f800000,(%rdx)
test   %edi,%edi
jle    1277 <kcof2acof+0x167>
mov    %rdx,%rbx
mov    %rsi,%r14
lea    0x4(%rdx),%rbp
mov    %edi,%r8d
lea    -0xc(%rdx),%r12
mov    $0x1,%r13d
xor    %r15d,%r15d
mov    %rbp,0x8(%rsp)
mov    %r8,(%rsp)
jmp    1186 <kcof2acof+0x76>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movss  -0x4(%r14,%r13,4),%xmm0
movss  %xmm0,(%rbx,%r13,4)
add    $0x1,%r13
add    $0x1,%r15
add    $0x4,%r12
add    $0x4,%rbp
cmp    %r8,%r15
je     1277 <kcof2acof+0x167>
movl   $0x0,(%rbx,%r13,4)
cmp    $0x2,%r13
jb     1160 <kcof2acof+0x50>
movss  -0x4(%r14,%r13,4),%xmm0
mov    $0x1,%eax
cmp    $0x8,%r15
jb     1211 <kcof2acof+0x101>
mov    %r15,%rcx
and    $0xfffffffffffffff8,%rcx
mov    %r15,%rdx
movabs $0x7ffffffffffffff8,%rax
ff ff 7f
and    %rax,%rdx
lea    0x1(%rdx),%rax
movaps %xmm0,%xmm1
shufps $0x0,%xmm0,%xmm1
mov    %r12,%rsi
xor    %edi,%edi
nopl   (%rax)
movups 0x4(%rbx,%rdi,4),%xmm2
movups 0x14(%rbx,%rdi,4),%xmm3
movups -0x10(%rsi),%xmm4
movups (%rsi),%xmm5
shufps $0x1b,%xmm5,%xmm5
shufps $0x1b,%xmm4,%xmm4
mulps  %xmm1,%xmm5
addps  %xmm2,%xmm5
mulps  %xmm1,%xmm4
addps  %xmm3,%xmm4
movups %xmm5,0x14(%rsp,%rdi,4)
movups %xmm4,0x24(%rsp,%rdi,4)
add    $0x8,%rdi
add    $0xffffffffffffffe0,%rsi
cmp    %rdi,%rcx
jne    11d0 <kcof2acof+0xc0>
cmp    %rdx,%r15
je     1240 <kcof2acof+0x130>
lea    0x0(,%rax,4),%rdx
mov    %rbp,%rcx
sub    %rdx,%rcx
nop
movss  (%rcx),%xmm1
mulss  %xmm0,%xmm1
addss  (%rbx,%rax,4),%xmm1
movss  %xmm1,0x10(%rsp,%rax,4)
add    $0x1,%rax
add    $0xfffffffffffffffc,%rcx
cmp    %rax,%r13
jne    1220 <kcof2acof+0x110>
cmp    $0x2,%r13
jb     1160 <kcof2acof+0x50>
movabs $0x3fffffffc,%rax
lea    (%rax,%r15,4),%rdx
and    %rax,%rdx
add    $0x4,%rdx
mov    0x8(%rsp),%rdi
lea    0x14(%rsp),%rsi
call   1030 <memcpy@plt>
mov    (%rsp),%r8
jmp    1160 <kcof2acof+0x50>
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @kcof2acof(i32 %0, float* %1, float* %2) {
%4 = alloca i32
%5 = alloca float*
%6 = alloca float*
%7 = alloca [22 x float]
%8 = alloca i32
%9 = alloca i32
store i32 %0, i32* %4
store float* %1, float** %5
store float* %2, float** %6
%10 = load float*, float** %6
%11 = getelementptr inbounds float, float* %10, i64 0
store float 1.000000e+00, float* %11
%12 = getelementptr inbounds [22 x float], [22 x float]* %7, i64 0, i64 0
store float 1.000000e+00, float* %12
store i32 1, i32* %9
br label %13
13: ; preds = %80, %3
%14 = load i32, i32* %9
%15 = load i32, i32* %4
%16 = icmp sle i32 %14, %15
br i1 %16, label %17, label %83
17: ; preds = %13
%18 = load float*, float** %6
%19 = load i32, i32* %9
%20 = sext i32 %19 to i64
%21 = getelementptr inbounds float, float* %18, i64 %20
store float 0.000000e+00, float* %21
store i32 1, i32* %8
br label %22
22: ; preds = %49, %17
%23 = load i32, i32* %8
%24 = load i32, i32* %9
%25 = icmp slt i32 %23, %24
br i1 %25, label %26, label %52
26: ; preds = %22
%27 = load float*, float** %6
%28 = load i32, i32* %8
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds float, float* %27, i64 %29
%31 = load float, float* %30
%32 = load float*, float** %5
%33 = load i32, i32* %9
%34 = sub nsw i32 %33, 1
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds float, float* %32, i64 %35
%37 = load float, float* %36
%38 = load float*, float** %6
%39 = load i32, i32* %9
%40 = load i32, i32* %8
%41 = sub nsw i32 %39, %40
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds float, float* %38, i64 %42
%44 = load float, float* %43
%45 = call float @llvm.fmuladd.f32(float %37, float %44, float %31)
%46 = load i32, i32* %8
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds [22 x float], [22 x float]* %7, i64 0, i64 %47
store float %45, float* %48
br label %49
49: ; preds = %26
%50 = load i32, i32* %8
%51 = add nsw i32 %50, 1
store i32 %51, i32* %8
br label %22
52: ; preds = %22
store i32 1, i32* %8
br label %53
53: ; preds = %66, %52
%54 = load i32, i32* %8
%55 = load i32, i32* %9
%56 = icmp slt i32 %54, %55
br i1 %56, label %57, label %69
57: ; preds = %53
%58 = load i32, i32* %8
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds [22 x float], [22 x float]* %7, i64 0, i64 %59
%61 = load float, float* %60
%62 = load float*, float** %6
%63 = load i32, i32* %8
%64 = sext i32 %63 to i64
%65 = getelementptr inbounds float, float* %62, i64 %64
store float %61, float* %65
br label %66
66: ; preds = %57
%67 = load i32, i32* %8
%68 = add nsw i32 %67, 1
store i32 %68, i32* %8
br label %53
69: ; preds = %53
%70 = load float*, float** %5
%71 = load i32, i32* %9
%72 = sub nsw i32 %71, 1
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds float, float* %70, i64 %73
%75 = load float, float* %74
%76 = load float*, float** %6
%77 = load i32, i32* %9
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds float, float* %76, i64 %78
store float %75, float* %79
br label %80
80: ; preds = %69
%81 = load i32, i32* %9
%82 = add nsw i32 %81, 1
store i32 %82, i32* %9
br label %13
83: ; preds = %13
ret void
}

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