void rfctolpc(float *buf) {
  float a, b;
  register float *cptr;
  register int n, k;
  cptr = buf + 12;
  for (n = -1; n >= -12; n--) {
    *(cptr + n) = -(*(cptr + n));
    for (k = -1; 2 * k >= n; k--) {
      a = *(cptr + k);
      b = *(cptr + n - k);
      *(cptr + k) = a - b * *(cptr + n);
      *(cptr + n - k) = b - a * *(cptr + n);
    }
  }
}


x86 assembly (clang, -O0):
<rfctolpc>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    -0x8(%rbp),%rax
add    $0x30,%rax
mov    %rax,-0x18(%rbp)
movl   $0xffffffff,-0x1c(%rbp)
cmpl   $0xfffffff4,-0x1c(%rbp)
jl     1239 <rfctolpc+0x139>
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movl   $0xffffffff,-0x20(%rbp)
mov    -0x20(%rbp),%eax
shl    $0x1,%eax
cmp    -0x1c(%rbp),%eax
jl     1226 <rfctolpc+0x126>
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0xc(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rdx
xor    %ecx,%ecx
sub    %rdx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x10(%rbp)
movss  -0xc(%rbp),%xmm1
movss  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x10(%rbp),%xmm1
movss  -0xc(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rdx
xor    %ecx,%ecx
sub    %rdx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x20(%rbp)
jmp    1153 <rfctolpc+0x53>
jmp    122b <rfctolpc+0x12b>
mov    -0x1c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x1c(%rbp)
jmp    111b <rfctolpc+0x1b>
pop    %rbp
ret

x86 assembly (clang, -O2):
<rfctolpc>:
xor    %eax,%eax
movaps 0xef7(%rip),%xmm0
jmp    1120 <rfctolpc+0x20>
nopl   0x0(%rax,%rax,1)
add    $0xffffffffffffffff,%rax
lea    -0x1(%rax),%ecx
cmp    $0xfffffff3,%ecx
jbe    1286 <rfctolpc+0x186>
movss  0x2c(%rdi,%rax,4),%xmm1
movaps %xmm1,%xmm2
xorps  %xmm0,%xmm2
movss  %xmm2,0x2c(%rdi,%rax,4)
test   %rax,%rax
je     1110 <rfctolpc+0x10>
movss  0x2c(%rdi),%xmm2
movss  0x30(%rdi,%rax,4),%xmm3
mulss  %xmm3,%xmm1
addss  %xmm2,%xmm1
movss  %xmm1,0x2c(%rdi)
mulss  0x2c(%rdi,%rax,4),%xmm2
lea    -0x1(%rax),%rcx
subss  %xmm2,%xmm3
movss  %xmm3,0x30(%rdi,%rax,4)
cmp    $0xfffffffc,%ecx
ja     1110 <rfctolpc+0x10>
movss  0x28(%rdi),%xmm1
movss  0x34(%rdi,%rax,4),%xmm2
movss  0x2c(%rdi,%rax,4),%xmm3
mulss  %xmm2,%xmm3
movaps %xmm1,%xmm4
subss  %xmm3,%xmm4
movss  %xmm4,0x28(%rdi)
mulss  0x2c(%rdi,%rax,4),%xmm1
subss  %xmm1,%xmm2
movss  %xmm2,0x34(%rdi,%rax,4)
cmp    $0xfffffffa,%ecx
ja     1110 <rfctolpc+0x10>
movss  0x24(%rdi),%xmm1
movss  0x38(%rdi,%rax,4),%xmm2
movss  0x2c(%rdi,%rax,4),%xmm3
mulss  %xmm2,%xmm3
movaps %xmm1,%xmm4
subss  %xmm3,%xmm4
movss  %xmm4,0x24(%rdi)
mulss  0x2c(%rdi,%rax,4),%xmm1
subss  %xmm1,%xmm2
movss  %xmm2,0x38(%rdi,%rax,4)
cmp    $0xfffffff8,%ecx
ja     1110 <rfctolpc+0x10>
movss  0x20(%rdi),%xmm1
movss  0x3c(%rdi,%rax,4),%xmm2
movss  0x2c(%rdi,%rax,4),%xmm3
mulss  %xmm2,%xmm3
movaps %xmm1,%xmm4
subss  %xmm3,%xmm4
movss  %xmm4,0x20(%rdi)
mulss  0x2c(%rdi,%rax,4),%xmm1
subss  %xmm1,%xmm2
movss  %xmm2,0x3c(%rdi,%rax,4)
cmp    $0xfffffff6,%ecx
ja     1110 <rfctolpc+0x10>
movss  0x1c(%rdi),%xmm1
movss  0x40(%rdi,%rax,4),%xmm2
movss  0x2c(%rdi,%rax,4),%xmm3
mulss  %xmm2,%xmm3
movaps %xmm1,%xmm4
subss  %xmm3,%xmm4
movss  %xmm4,0x1c(%rdi)
mulss  0x2c(%rdi,%rax,4),%xmm1
subss  %xmm1,%xmm2
movss  %xmm2,0x40(%rdi,%rax,4)
cmp    $0xfffffff4,%ecx
ja     1110 <rfctolpc+0x10>
movss  0x18(%rdi),%xmm1
movss  0x44(%rdi,%rax,4),%xmm2
movss  0x2c(%rdi,%rax,4),%xmm3
mulss  %xmm2,%xmm3
movaps %xmm1,%xmm4
subss  %xmm3,%xmm4
movss  %xmm4,0x18(%rdi)
mulss  0x2c(%rdi,%rax,4),%xmm1
subss  %xmm1,%xmm2
movss  %xmm2,0x44(%rdi,%rax,4)
jmp    1110 <rfctolpc+0x10>
ret

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

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